An appropriate potential function, which satisfies the boundary condition is found by using arrays of infinite sources and sinks.

The fluid flow among the fluids was not considered in fracture reservoir in the past, and only the source and sink was considered in the fluid exchange between the fracture and matrix.

However, as a source and a sink as well, soil is being changed because of human activities, such as changing land use and applying fertilizers.