Determining fluid flow characteristics of porous mediums

1. A computer implemented method for simulating aging and a wettability alteration of a physical rock sample in the presence of at least two fluids including water and oil, with the method comprising, comprises:

retrieving by a computing system a representation of a physical rock sample, the representation including chemical composition, pore space, and grain space data corresponding to the physical rock sample;

iteratively simulating by the computing system aging and wettability alteration fora surface location in the pore space of the representation of the physical rock sample;

calculating by the computing system, a change in a local curvature for the surface location in the pore space of the representation of the physical rock sample from a prior iteration of the wettability alteration, with calculating by the computing system, the change in local curvature comprising:

calculating by the computing system, local capillary pressure and local disjoining pressure of the surface location in the pore space; and the method further comprises:

determining from the calculated change in local curvature whether water-oil film breakage will occur;

for some of the surface locations, determine whether water-oil film breakage will occur by comparing the capillary pressure to the sum of local curvature-based pressure and local disjoining pressure; and

for at least some of the surface locations with water-oil film breakage, determine the degree of wettability alteration based on local surface properties, local fluid chemical composition, and aging;

determining by the computer system a modified surface wettability of the surface location corresponding to the representation of the physical rock based on the change in the calculated local curvature; and

verifying by the computer system a degree of aging in relation to an experimental value of aging of the physical rock sample, and terminating the iterative simulating, upon verification of the aging.

2. The method of claim 1 further comprises:

assigning input parameters needed for wettability simulation to the computer simulation environment.

3. The method of claim 1 wherein the representation of the rock sample is modeled as a triangulated mesh and fluid is represented by voxels with voxels assigned attributes.

4. The method of claim 3 wherein voxels are cut to conform the fluid to the triangulated mesh.

5. The method of claim 1 further comprises:

performing a drainage simulation to establish a fluid phase distribution according to user-specified conditions.

6. The method of claim 5 wherein the user-specified conditions are selected from the group consisting of representative laboratory and reservoir conditions.

7. The method of claim 6 wherein performing drainage simulation further comprises:

simulating drainage using a fluid computational technique to establish fluid phase distribution for the user-specified conditions through the rock sample representation.

8. The method of claim 5 wherein performing drainage simulation further comprises:

calculating the change in local curvature at each surface in the rock sample representation; and

applying results of the change in the local curvature calculation to calculate the local capillary pressure at each surface location.

9. The method of claim 1 wherein after performing an iteration of determining the numerical representation of wettability alteration, the method further comprises:

allowing the at least two fluids to redistribute due to the modified surface wettability.

10. The method of claim 1 wherein the method provides a numerical prediction of surface wettability alteration and resulting spatial distribution of wettability in the form of the calculated contact angle at each surface location.

11. The method of claim 1 wherein a discrete surface element and one or more fluid elements are in communication, and local wettability determinations are made that depend on the local properties of mineral type, fluid properties, fluid composition, capillary pressure, and disjoining pressure; and wherein wettability alteration is iteratively performed followed by fluid phase redistribution, with the wettability alteration repeated until the computing system detects a convergence measured by a subsequent repetition not significantly changing values of determined local curvature.

12. A computer system comprises:

one or more processor devices;

memory coupled to the one or more processor devices;

storage storing executable computer instructions to simulate aging and a wettability alteration of a physical rock sample in the presence of at least two fluids that are oil and water, the instructions to configure the one or more processors to:

retrieve a representation of a physical rock sample, the representation including chemical composition, pore space, and grain space data corresponding to the physical rock sample;

iteratively simulate by the computing system aging and wettability alteration for a surface location in the pore space of the representation of the physical rock sample;

calculate a change in a local curvature for the surface location in the pore space of the representation of the physical rock sample from a prior iteration of the wettability alteration, with instructions to calculate the change in local curvature comprising:

calculate local capillary pressure and local disjoining pressure of the surface location in the pore space; and the method further comprises:

determine from the calculated change in local curvature whether water-oil film breakage will occur;

for some of the surface locations, determine whether water-oil film breakage will occur by comparing the capillary pressure to the sum of local curvature-based pressure and local disjoining pressure; and

for at least some of the surface locations with water-oil film breakage, determine the degree of wettability alteration based on local surface properties, local fluid chemical composition, and aging;

determine a modified surface wettability of the surface location corresponding to the representation of the physical rock based on the change in the calculated local curvature; and

verify a degree of aging in relation to an experimental value of aging of the physical rock sample, and terminate the iterative simulation, upon verification of the aging.

13. The computer system of claim 12 wherein instructions further configure the one or more processors to:

assign input parameters needed for wettability simulation to the computer simulation environment.

14. The computer system of claim 12 wherein the representation of the rock sample is modeled as a triangulated mesh and fluid is represented by voxels with voxels assigned attributes.

15. The computer system of claim 14 wherein voxels are cut to conform the fluid to the triangulated mesh.

16. The computer system of claim 12 further configure the one or more processors:

perform a drainage simulation to establish a fluid phase distribution according to user-specified conditions.

17. The computer system of claim 16 wherein the user-specified conditions are selected from the group consisting of representative laboratory and reservoir conditions.

18. The computer system of claim 17 wherein instructions to perform drainage simulation further configure the one or more processors:

simulate drainage using a fluid computational technique to establish fluid phase distribution for the user-specified conditions through the rock sample representation.

19. The computer system of claim 16 wherein instructions to perform drainage simulation further configure the one or more processors:

calculate the change in local curvature calculation at each surface in the rock sample representation; and

apply results of the change in the local curvature calculation to calculate the local capillary pressure at each every surface location.

20. The computer system of claim 12 wherein after performing an iteration of determining the numerical representation of wettability alteration, the instructions further configure the one or more processors to:

allow the at least two fluids to redistribute due to the modified surface wettability.

21. The computer system of claim 12 wherein the computer system provides a numerical prediction of surface wettability alteration and resulting spatial distribution of wettability in the form of the calculated contact angle at each surface location.

22. The computer system of claim 12 wherein a discrete surface element and one or more fluid elements are in communication, and local wettability determinations are made that depend on the local properties of mineral type, fluid properties, fluid composition, capillary pressure, and disjoining pressure.

23. The computer system of claim 12 wherein wettability alteration is iteratively performed followed by fluid phase redistribution, with the wettability alteration repeated until the computer system detects a convergence measured by a subsequent repetition not significantly changing values of determined local curvature.

24. A computer program product tangibly stored on a computer readable non-transitory storage device that stores executable computer instructions to simulate aging and a wettability alteration of a physical rock sample in the presence of at least two fluids, the instructions for causing a computing system to:

retrieve a representation of a physical rock sample, the representation including chemical composition, pore space, and grain space data corresponding to the physical rock sample;

iteratively simulate by the computing system aging and wettability alteration for a surface location in the pore space of the representation of the physical rock sample;

calculate a change in a local curvature for the surface location in the pore space of the representation of the physical rock sample from a prior iteration of the wettability alteration, with instructions to calculate the change in local curvature comprising:

calculate local capillary pressure and local disjoining pressure of the surface location in the pore space; and the method further comprises:

determine from the calculated change in local curvature whether water-oil film breakage will occur;

for some of the surface locations, determine whether water-oil film breakage will occur by comparing the capillary pressure to the sum of local curvature-based pressure and local disjoining pressure; and

for at least some of the surface locations with water-oil film breakage, determine the degree of wettability alteration based on local surface properties, local fluid chemical composition, and aging;

determine a modified surface wettability of the surface location corresponding to the representation of the physical rock based on the change in the calculated local curvature; and

verify a degree of aging in relation to an experimental value of aging of the physical rock sample, and terminate the iterative simulation, upon verification of the aging.

25. The computer program product of claim 24 wherein instructions further comprises instructions to:

assign input parameters needed for wettability simulation to the computer simulation environment.

26. The computer program product of claim 24 wherein the representation of the rock sample is modeled as a triangulated mesh and fluid is represented by voxels with voxels assigned attributes.

27. The computer program product of claim 26 wherein voxels are cut to conform the fluid to the triangulated mesh.

28. The computer program product of claim 24 further comprises instructions to:

receive user specified conditions; and

perform a drainage simulation to establish a fluid phase distribution according to the received user-specified conditions.

29. The computer program product of claim 28 wherein instructions to perform drainage simulation further comprises instructions to

calculate the change in a local curvature calculation at each surface in the rock sample representation; and

apply results of the change in the local curvature calculation to calculate the local capillary pressure at each surface location.

30. The computer program product of claim 24 wherein after performing an iteration of determining the numerical representation of wettability alteration, the computing system is further configured to:

allow the at least two fluids to redistribute due to the modified surface wettability.

31. The computer program product of claim 24 wherein the instructions cause the computing system to provide a numerical prediction of surface wettability alteration and resulting spatial distribution of wettability in the form of the calculated contact angle at each surface location.

32. The computer program product of claim 24 wherein a discrete surface element and one or more fluid elements are in communication, and local wettability determinations are made that depend on the local properties of mineral type, fluid properties, fluid composition, capillary pressure, and disjoining pressure.

33. The computer program product of claim 24 wherein wettability alteration is iteratively performed followed by fluid phase redistribution, with the wettability alteration repeated until the computing system detects a convergence measured by a subsequent repetition not significantly changing values of determined local curvature.

34. The method of claim 1 wherein wettability alteration is iteratively performed followed by fluid phase redistribution, with the wettability alteration repeated until the computing system detects a convergence measured by a subsequent repetition not significantly changing values of determined local curvature.