Abstract:
This invention relates to methods for determination of thermophysical properties of porous media filled with fluid, gas, or other mineral medium, and can be used, in particular, in the oil and gas industry. When implementing the method, it is necessary to initially determine a composition of a saturated porous medium and thermal conductivity coefficients of its components. A three-dimensional image of a sample of the porous medium is obtained by X-ray scanning. A thermal conductivity coefficient of the medium is determined by solving thermal conductivity problems, based on decomposition of a computational domain followed by composition.

Description:
CROSS-REFERENCE TO RELATED APPLICATION 
       [0001]    This application is a United States National Stage Application under 35 U.S.C. §371 and claims priority to Patent Cooperation Treaty Application Number PCT/RU2010/000804 filed Dec. 30, 2010, which is incorporated herein by reference in its entirety. 
     
    
     FIELD OF THE DISCLOSURE 
       [0002]    This invention relates to methods for determination of thermophysical properties of porous media filled with fluid, gas, or other mineral medium. The invention can be used for such porous media as composite materials, geological formations, bone tissues. 
         [0003]    The materials are considered to be “porous” if a medium (skeleton) they are composed of, is porous (has pores or discontinuities). These pores can be void, but most often they are filled with some fluids, such as water, oil, gas, etc. 
         [0004]    The areas which can be filled with fluid include, in particular, pores, layer fractures, cavities, etc. Water, oil and other liquids can serve as the fluid here. Natural gas, air, or any other gaseous substance can be regarded as the gas. Any solid formation in fixed or bulk condition, having a mineralogical composition, can serve as the mineral medium. 
         [0005]    The method can be used in oil and gas industry, in particular, for determination of thermal conductivity of an oil-saturated porous geological medium, in the development of oil and gas fields (injection wells, production wells). 
       BACKGROUND OF THE DISCLOSURE 
       [0006]    Thermal methods for enhanced oil recovery provide for a preliminary modeling of the heat-and-mass exchange processes in a reservoir (an oil reservoir) and in wells, as well as determination of thermal conditions for downhole equipment. So information on thermophysical properties of rock (such as the thermal conductivity, heat capacity and thermal diffusivity) is of primary importance. 
         [0007]    There are known methods for determination of the thermal conductivity of porous rock, such as described in S. L. Lee, J. H. Yang, Modeling of effective thermal conductivity for a non-homogeneous anisotropic porous medium, Int. J. Heat Mass Transfer, Vol. 41, Issue 6-7, pp. 931-937, 1998, or described in H. Ye, Effective thermal conductivity of two-scale porous media, APPLIED PHYSICS LETTERS 89, 081908 (2006). 
         [0008]    However, the above-mentioned methods allow for determining thermal conductivity of only small-sized samples and are characterized by increased labor requirements and increased time consumption. 
       SUMMARY OF THE DISCLOSURE 
       [0009]    The claimed method for determination of thermal conductivity of a saturated porous medium allows for improving the accuracy of determination of thermal conductivity as compared to known methods, and consists in determining a microstructure of a sample of a porous medium and in solving a thermal conductivity problem, based on decomposition of a computational domain with following composition. 
         [0010]    The claimed method comprises the following steps:
       determining a composition of a saturated porous medium and thermal conductivity coefficients of determined components,   performing an X-ray scanning of a sample of the saturated porous medium and obtaining a three-dimensional image of the sample which is a first single computational domain and consists of N voxels, with N=N x ·N y ·N z , where N x , N y , N z  are numbers of voxels along x, y, z axes, respectively,   combining the N voxels into n coarse cells, with n=n x ·n y ·n z , where n x , n y , n z  are numbers of coarse cells along the x, y, z axes, respectively, with 2&lt;n x &lt;N x , 2&lt;n y &lt;N y , 2&lt;n z &lt;N z ,   solving a first thermal conductivity problem for each of the coarse cells in each of the x, y, z directions, and determining an effective thermal conductivity coefficient in each direction,   combining the coarse cells into a second single computational domain,   solving a second thermal conductivity problem for the second single computational domain which is composed of the coarse cells, and determining an effective thermal conductivity coefficient for the whole computational domain in three directions,   determining a thermal conductivity coefficient of the saturated porous medium as an average of the determined thermal conductivity coefficients along the axes.       
 
     
    
     
       BRIEF DESCRIPTION OF THE FIGURES 
         [0018]      FIG. 1  shows a microstructure of a parallelepiped-shaped sample; 
           [0019]      FIG. 2  shows decomposition of a computational domain; 
           [0020]      FIG. 3  shows a coarse cell; 
           [0021]      FIG. 4  shows boundaries for x, y, z directions; and 
           [0022]      FIG. 5  shows a rock sample 240×240×240 mm in size. 
       
    
    
     DETAILED DESCRIPTION 
       [0023]    At first it is necessary to determine a composition of a porous medium and, consequently, thermal conductivity coefficients of determined components of the porous medium. The components can be determined by known methods, such as methods described in the paper by V. N. Shvanov, Petrography of Sandy Rock (Composition, Systematics, and Description of Mineral Species)—Leningrad: Nedra Publishing House, 1987, p. 269. The thermal conductivity coefficients of the determined components are determined from reference literature, e.g., from Clauser, C and E. Huenges, 1995, Thermal Conductivity of Rocks and Minerals, Rock Physics and Phase Relations. 
         [0024]    Then, it is necessary to obtain information on microstructure of the porous medium. A rectangular-parallelepiped-shaped sample of the saturated porous medium is taken for tests. Results of this test are usually represented as data on filling of each unit of a computational domain (the sample of the medium represents the computational domain) with one or another material. The unit of the computational domain is selected depending on the test type. 
         [0025]    A microstructure of the sample is determined by X-ray scanning of the sample. The X-ray computer tomographs which are used for obtaining the data on the microstructure are described, for example, in http://www.ngi.no/no/Innholdsbokser/Referansjeprosjekter-LISTER-/Referanser/Rock-structure-visualisation-with-CT/A or in the paper by Sasov, Microtomography, Journal of Microscopy, Vol. 147 (2): 169-192, 1987. 
         [0026]    The scanning gives data on the microstructure of the sample having a height H, a length L and a width B in the form of a first three-dimensional computational domain consisting of N x ·N y ·N z  three-dimensional cells (voxels), where N x , N y , N z  are numbers of voxels along x, y, z axes, respectively, with each of these voxels having an edge length equal to a scanning resolution value ( FIG. 1 ). In the microstructure determination, each of these voxels is considered to be filled; that is to say, it is considered to correspond to a certain material from the two-component or multi-component medium being studied, e.g., to a skeleton or a fluid. For example, white color in  FIG. 1  corresponds to the skeleton of the porous medium, while black color corresponds to the fluid in pores. Consequently, the thermal conductivity coefficients of each of the voxels are considered to be known. 
         [0027]      FIG. 2  shows decomposition of the first computational domain having a length δx, a width δy and a height δz, in the course of which N x ·N y ·N z  voxels of the first computational domain are combined into n x ·n y ·n z  of coarse cells having dimensions of a cell equal to D x ·D y ·D z  in x-, y-, z-axis directions, respectively, so that a number of the coarse cells along a respective axis meets the following inequalities: 2&lt;n x &lt;N x , 2&lt;n y &lt;N y , 2&lt;n z &lt;N z . 
         [0028]    For each of these coarse cells (with a number ijk) in each of the x, y, z directions, a first thermal conductivity problem is solved, based on the established assumptions ( FIG. 3 ). A mathematical model of this first thermal conductivity problem consists of a thermal conductivity equation (1) and boundary conditions (2): 
         [0000]      div( l   a grad T )=0,   (1)
 
         [0000]    where l a  is a thermal conductivity coefficient which is defined according to the X-ray scanning results for each voxel; a is a voxel number in a coarse cell (for example, in case of a two-component medium, we assume that a=1 corresponds to a skeleton, while a=0 corresponds to a pore-filling material); T is a temperature; div, grad are mathematical operators. 
         [0029]    The boundary conditions are established as follows ( FIG. 4 ). In each of the x, y, z directions, a first cell surface perpendicular to an axis under consideration is considered to be B 1 , and an opposite surface is considered to be B 2 . Temperatures on the above-mentioned surfaces are taken equal to T 1           T 2 , respectively. Other surfaces of the coarse cell are considered to be thermally insulated. So, the boundary conditions are expressed as: 
         [0000]      On B 1  boundary: T=T 1 , 
         [0000]      On B 2  boundary: T=T 2 , 
         [0000]      On all other boundaries: ( l   a  grad T,n )=0,   (2)
 
         [0000]    where n is a normal to the boundary surface. 
         [0030]    It should be noted at the same time that the mathematical model of the first thermal conductivity problem for the first computational domain is built, based on the following important assumptions: a heat flux propagates in a direction under consideration only, i.e., we assume that the heat flux propagates along the coordinate axis (x, y, z) for which the current equation is set up; principal axes of anisotropy coincide with the coordinate axes (i.e. a length of irregularity in properties coincides with the x, y, z coordinate axes). 
         [0031]    Temperature fields, from which fields of normal temperature derivatives with respect to the boundary surfaces of a coarse cell are calculated, are the solution to the mathematical model consisting of the thermal conductivity equations (1) and the boundary conditions (2) based on the adopted assumptions. 
         [0032]    A heat flux that passes across a boundary of a coarse cell in the x, y, z directions is determined through an area integral of a temperature derivative in the corresponding directions (axes) and is equal to: 
         [0000]    
       
         
           
             
               
                 Q 
                 x 
                 ijk 
               
               = 
               
                 
                   ∫ 
                   
                     S 
                     yz 
                   
                 
                  
                 
                   
                     l 
                     α 
                   
                    
                   
                     
                       ∂ 
                       T 
                     
                     
                       ∂ 
                       x 
                     
                   
                    
                   
                       
                   
                    
                   
                      
                     S 
                   
                 
               
             
             , 
             
               
                 Q 
                 y 
                 ijk 
               
               = 
               
                 
                   ∫ 
                   
                     S 
                     xz 
                   
                 
                  
                 
                   
                     l 
                     α 
                   
                    
                   
                     
                       ∂ 
                       T 
                     
                     
                       ∂ 
                       y 
                     
                   
                    
                   
                       
                   
                    
                   
                      
                     S 
                   
                 
               
             
             , 
             
               
                 Q 
                 z 
                 ijk 
               
               = 
               
                 
                   ∫ 
                   
                     S 
                     xy 
                   
                 
                  
                 
                   
                     l 
                     α 
                   
                    
                   
                     
                       ∂ 
                       T 
                     
                     
                       ∂ 
                       z 
                     
                   
                    
                   
                       
                   
                    
                   
                      
                     S 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    where S yz  is an area of the sample surface perpendicular to the x axis, S xz  is an area of the sample surface perpendicular to the y axis, S xy , is an area of the sample surface perpendicular to the z axis. 
         [0033]    In case of a macroscopic description of the porous medium, the heat propagation is described by Fourier&#39;s law, wherein a thermal conductivity tensor is replaced with an effective thermal conductivity tensor. 
         [0034]    For each of the coarse cells, a heat propagation described by Fourier&#39;s law is given by: 
         [0000]    
       
         
           
             
               = 
               
                 
                   
                     Q 
                     x 
                     ijk 
                   
                    
                   Dx 
                 
                 
                   
                     ( 
                     
                       
                         T 
                         2 
                       
                       - 
                       
                         T 
                         1 
                       
                     
                     ) 
                   
                    
                   DyDz 
                 
               
             
             , 
             
               
 
             
              
             
               = 
               
                 
                   
                     Q 
                     y 
                     ijk 
                   
                    
                   Dy 
                 
                 
                   
                     ( 
                     
                       
                         T 
                         2 
                       
                       - 
                       
                         T 
                         1 
                       
                     
                     ) 
                   
                    
                   DxDz 
                 
               
             
             , 
             
               
 
             
              
             
               = 
               
                 
                   
                     Q 
                     z 
                     ijk 
                   
                    
                   Dz 
                 
                 
                   
                     ( 
                     
                       
                         T 
                         2 
                       
                       - 
                       
                         T 
                         1 
                       
                     
                     ) 
                   
                    
                   DxDy 
                 
               
             
             , 
           
         
       
     
         [0000]    where          ,          ,           are effective thermal conductivity coefficients in the x, y, z directions, respectively. 
         [0035]    A next operation of the method is to combine the coarse cells into a second single computational domain, while considering that each of the coarse cells is homogeneous and has the thermal conductivity coefficient equal to the effective thermal conductivity coefficient in each of such cells. 
         [0036]    Then, a second thermal conductivity problem is solved for the second combined computational domain which is composed of the coarse cells. A mathematical model of this second thermal conductivity problem consists of a thermal conductivity equation for the second combined computational domain and is given by: 
         [0000]      div([l           ijk ]grad T )=0, 
         [0000]    and of boundary conditions which represent given T 1  and T 2  temperatures on the first surface perpendicular to the axis under consideration (x, y, z) and on the opposite surface, respectively, while all other surfaces are considered to be thermally insulated. 
         [0037]    The assumptions that the heat flux propagates in the direction under consideration only (along the x, y, z axis) and that the main axes of anisotropy match the coordinate axes, remain in force. 
         [0038]    A heat flux that passes across the boundaries of the second combined computation domain in the x, y, z directions is equal to: 
         [0000]    
       
         
           
             
               
                 Q 
                 x 
               
               = 
               
                 
                   ∫ 
                   
                     S 
                     yz 
                     ′ 
                   
                 
                  
                 
                    
                   
                     
                       ∂ 
                       T 
                     
                     
                       ∂ 
                       x 
                     
                   
                    
                   
                       
                   
                    
                   
                      
                     S 
                   
                 
               
             
             , 
             
               
                 Q 
                 y 
               
               = 
               
                 
                   ∫ 
                   
                     S 
                     xz 
                     ′ 
                   
                 
                  
                 
                    
                   
                     
                       ∂ 
                       T 
                     
                     
                       ∂ 
                       y 
                     
                   
                    
                   
                       
                   
                    
                   
                      
                     S 
                   
                 
               
             
             , 
             
               
                 Q 
                 z 
               
               = 
               
                 
                   ∫ 
                   
                     S 
                     xy 
                     ′ 
                   
                 
                  
                 
                    
                   
                     
                       ∂ 
                       T 
                     
                     
                       ∂ 
                       z 
                     
                   
                    
                   
                       
                   
                    
                   
                      
                     S 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    where S′ yz  is an area of the sample surface perpendicular to the x axis, S′ xz  is an area of the sample surface perpendicular to the y axis, S′ xy  is an area of the sample surface perpendicular to the z axis. 
         [0039]    When solving the second thermal conductivity problem in combination with the boundary conditions based on the above-mentioned assumptions, we determine effective thermal conductivity coefficients for the whole computational domain in three (x, y, z) directions: 
         [0000]    
       
         
           
             
               
                 l 
                 
                   x 
                   ∋ 
                   φ 
                 
               
               = 
               
                 
                   
                     Q 
                     x 
                   
                    
                   
                      
                     x 
                   
                 
                 
                   
                     ( 
                     
                       
                         T 
                         2 
                       
                       - 
                       
                         T 
                         1 
                       
                     
                     ) 
                   
                    
                   
                      
                     y 
                   
                    
                   
                      
                     z 
                   
                 
               
             
             , 
             
               
 
             
              
             
               
                 l 
                 
                   y 
                   ∋ 
                   φ 
                 
               
               = 
               
                 
                   
                     Q 
                     y 
                   
                    
                   
                      
                     y 
                   
                 
                 
                   
                     ( 
                     
                       
                         T 
                         2 
                       
                       - 
                       
                         T 
                         1 
                       
                     
                     ) 
                   
                    
                   
                      
                     x 
                   
                    
                   
                      
                     z 
                   
                 
               
             
             , 
             
               
 
             
              
             
               
                 l 
                 
                   z 
                   ∋ 
                   φ 
                 
               
               = 
               
                 
                   
                     Q 
                     z 
                   
                    
                   
                      
                     z 
                   
                 
                 
                   
                     ( 
                     
                       
                         T 
                         2 
                       
                       - 
                       
                         T 
                         1 
                       
                     
                     ) 
                   
                    
                   
                      
                     x 
                   
                    
                   
                      
                     y 
                   
                 
               
             
             , 
           
         
       
     
         [0000]    where δx, δy, δz are sample dimensions. 
         [0040]    A thermal conductivity coefficient of the medium is equal to the average of the thermal conductivity coefficients determined along each of the coordinate axes, namely: 
         [0000]        l =⅓*(         +         +         ).
 
         [0041]    The suggested method was tested on samples 240×240×240 mm in size and was applied to samples 1,800×1,800×1,800 mm in size. 
         [0042]    For example, a thermal conductivity coefficient was determined by the suggested method for a sample 240×240×240 mm in size. 
         [0043]    The number of coarse cells was equal to n=4 along each of coordinate axes, and a total number of coarse cells was equal to 4×4×4=64. 
         [0044]    Effective thermal conductivity coefficients were determined for each of the coarse cells in x, y, z directions. The coarse cells were then combined into a single computational domain where each of the coarse cells was represented as a homogeneous material having the thermal conductivity coefficients determined in the previous step. Then, an effective thermal conductivity coefficient of the whole sample was determined in the x, y, z directions. Table 1 shows the values obtained for different filling fluids. 
         [0045]    Comparison of the thermal conductivity coefficient determination results with the exact solution shows that the error of determination in the x, y, z directions does not exceed 0.56% for the sample under consideration (Table 1). 
         [0000]    
       
         
               
             
               
               
               
               
               
               
               
               
             
               
             
               
               
               
               
               
               
               
               
             
               
             
               
               
               
               
               
               
               
               
             
           
               
                 TABLE 1 
               
               
                   
               
             
             
               
                 X direction 
               
             
          
           
               
                 Ratio of the thermal conductivity of the 
                 0.01 
                 0.1 
                 0.247 
                 0.5 
                 0.75 
                 0.99 
                 1 
               
               
                 saturating fluid to the thermal 
               
               
                 conductivity of rock, λ f /λ s   
               
               
                 Suggested method, λ x eff /λ s   
                 0.555 
                 0.650 
                 0.747 
                 0.856 
                 0.936 
                 0.998 
                 1 
               
               
                 Exact solution, λ x eff /λ s   
                 0.556 
                 0.649 
                 0.747 
                 0.857 
                 0.936 
                 0.998 
                 1 
               
               
                 Error, % 
                 0.135 
                 0.081 
                 0.043 
                 0.006 
                 0.001 
                 0.008 
                 0 
               
             
          
           
               
                 Y direction 
               
             
          
           
               
                 Ratio of the thermal conductivity of the 
                 0.01 
                 0.1 
                 0.247 
                 0.5 
                 0.75 
                 0.99 
                 1 
               
               
                 saturating fluid to the thermal 
               
               
                 conductivity of rock, λ f /λ s   
               
               
                 Suggested method, λ y eff /λ s   
                 0.549 
                 0.649 
                 0.747 
                 0.857 
                 0.936 
                 0.998 
                 1 
               
               
                 Exact solution, λ y eff /λ s   
                 0.553 
                 0.649 
                 0.747 
                 0.857 
                 0.936 
                 0.998 
                 1 
               
               
                 Error, % 
                 0.559 
                 0.094 
                 0.018 
                 0.008 
                 0.004 
                 0.008 
                 0 
               
             
          
           
               
                 Z direction 
               
             
          
           
               
                 Ratio of the thermal conductivity of the 
                 0.01 
                 0.1 
                 0.247 
                 0.5 
                 0.75 
                 0.99 
                 1 
               
               
                 saturating fluid to the thermal 
               
               
                 conductivity of rock, λ f /λ s   
               
               
                 Suggested method, λ z eff /λ s   
                 0.607 
                 0.684 
                 0.765 
                 0.863 
                 0.937 
                 0.998 
                 1 
               
               
                 Exact solution, λ z eff /λ s   
                 0.608 
                 0.684 
                 0.766 
                 0.863 
                 0.937 
                 0.998 
                 1 
               
               
                 Error, % 
                 0.171 
                 0.074 
                 0.052 
                 0.023 
                 0.007 
                 0.009 
                 0 
               
               
                   
               
             
          
         
       
     
         [0046]    Table 2 shows the results of the thermal conductivity coefficient determination for a sample 1,800×1,800×1,800 mm in size, by the suggested method. 
         [0000]    
       
         
               
               
               
             
               
               
               
             
               
               
               
             
               
               
               
               
             
           
               
                   
                 TABLE 2 
               
             
             
               
                   
                   
               
               
                   
                 Ratio of the thermal 
                   
               
               
                   
                 conductivity of the 
               
               
                   
                 saturating fluid to the 
               
               
                   
                 thermal conductivity 
               
               
                   
                 of rock, λ f /λ s   
               
             
          
           
               
                   
                 0.01 
                 0.247 
               
             
          
           
               
                   
                 λ eff , (W/m · K) 
                   
               
               
                   
                   
               
             
          
           
               
                   
                 X direction 
                 0.562 
                 0.750 
               
               
                   
                 Y direction 
                 0.548 
                 0.747 
               
               
                   
                 Z direction 
                 0.601 
                 0.764