The pore structure characteristic parameters of vuggy carbonate reservoirs were extracted,and matrix resistivity of vuggy reservoir was calculated by the percolation network simulation.A segmented cross-scale resistiv...The pore structure characteristic parameters of vuggy carbonate reservoirs were extracted,and matrix resistivity of vuggy reservoir was calculated by the percolation network simulation.A segmented cross-scale resistivity calculation method was established,in which the finite element method was used to simulate the resistivity of vuggy reservoirs.The mathematical models of vug porosity and water saturation with rock resistivity in vuggy carbonate reservoir were established,and the relationships between them were obtained.Experimental results verified the reliability of the simulation results.The method presented provides new technical means and research method for the resistivity log interpretation of vuggy carbonate reservoirs.The matrix porosity,vug porosity and matrix pore water saturation are the key factors determining the resistivity of reservoir rocks.展开更多
We present a path morphology method to separate total rock pore space into matrix, fractures and vugs and derive their pore structure spectrum. Thus, we can achieve fine pore evaluation in fracture–vug reservoirs bas...We present a path morphology method to separate total rock pore space into matrix, fractures and vugs and derive their pore structure spectrum. Thus, we can achieve fine pore evaluation in fracture–vug reservoirs based on electric imaging logging data. We automatically identify and extract fracture–vug information from the electric imaging images by adopting a path morphological operator that remains flexible enough to fit rectilinear and slightly curved structures because they are independent of the structuring element shape. The Otsu method was used to extract fracture–vug information from the background noise caused by the matrix. To accommodate the differences in scale and form of the different target regions,including fracture and vug path, operators with different lengths were selected for their recognition and extraction at the corresponding scale. Polynomial and elliptic functions are used to fit the extracted fractures and vugs, respectively, and the fracture–vug parameters are deduced from the fitted edge. Finally, test examples of numerical simulation data and several measured well data have been provided for the verification of the effectiveness and adaptability of the path morphology method in the application of electric imaging logging data processing. This also provides algorithm support for the fine evaluation of fracture–vug reservoirs.展开更多
A review of the pressure transient analysis of flow in reservoirs having natural fractures,vugs and/or caves is presented to provide an insight into how much knowledge has been acquired about this phenomenon and to hi...A review of the pressure transient analysis of flow in reservoirs having natural fractures,vugs and/or caves is presented to provide an insight into how much knowledge has been acquired about this phenomenon and to highlight the gaps still open for further research.A comparison-based approach is adopted which involved the review of works by several authors and identifying the limiting assumptions,model restrictions and applicability.Pressure transient analysis provides information to aid the identification of important features of reservoirs.It also provides an explanation to complex reservoir pressuredependent variations which have led to improved understanding and optimization of the reservoir dynamics.Pressure transient analysis techniques,however,have limitations as not all its models find application in naturally fractured and vuggy reservoirs as the flow dynamics differ considerably.Pollard’s model presented in 1953 provided the foundation for existing pressure transient analysis in these types of reservoirs,and since then,several authors have modified this basic model and come up with more accurate models to characterize the dynamic pressure behavior in reservoirs with natural fractures,vugs and/or caves,with most having inherent limitations.This paper summarizes what has been done,what knowledge is considered established and the gaps left to be researched on.展开更多
Considering the fluid flow non-darcy characteristics in fracture-vug carbonate reservoirs, a new multi-scale conduit flow model production prediction method for fracture-vug carbonate reservoirs was presented using im...Considering the fluid flow non-darcy characteristics in fracture-vug carbonate reservoirs, a new multi-scale conduit flow model production prediction method for fracture-vug carbonate reservoirs was presented using image segmentation technique of electric imaging logging data. Firstly, based on Hagen-Poiseuille's law of incompressible fluid flow and the different cross-sectional areas in single fractures and vugs in carbonate reservoirs, a multi-scale conduit flow model for fracture-vug carbonate reservoir was established, and a multi-scale conduit radial fluid flow equation was deduced. Then, conduit flow production index was introduced. The conduit flow production index was calculated using fracture-vug area extracted from the result of electrical image segmentation. Finally, production prediction of fracture-vug carbonate reservoir was realized by using electric imaging logging data. The method has been applied to Ordovician fracture-vug carbonate reservoirs in the Tabei area, and the predicted results are in good agreement with the oil testing data.展开更多
Based on the similarity criterion, volcanic rock samples were taken from outcrops to make experimental models. Water flooding experiments of five-spot well pattern, nine-spot well pattern, five-spot to nine-spot well ...Based on the similarity criterion, volcanic rock samples were taken from outcrops to make experimental models. Water flooding experiments of five-spot well pattern, nine-spot well pattern, five-spot to nine-spot well pattern, the relationship between relative well and fracture positions, and injection rate in dissolution vug-cave reservoirs with/without fractures were carried out to analyze variation regularities of development indexes, find out development characteristics of water flooding with different well patterns and sort out the optimal water flooding development mode. For dissolution vug-cave reservoirs without fractures, five-spot well pattern waterflooding has very small sweeping area, serious water channeling and low oil recovery. When the well pattern was adjusted from five-spot to nine-spot well pattern, oil recovery could be largely improved, but the corner well far from the injector is little affected. In dissolution vug-cave reservoirs with fractures, when the injector and producer are not connected by fractures, the fractures could effectively connect the poorly linked vugs to improve the development effect of water flooding. Whether there are fractures or not in dissolution vug-cave reservoirs, the development effect of nine-spot well-pattern is much better than that of five-spot well pattern and five-spot to nine-spot well pattern, this is more evident when there are fractures, and the edge well has better development indexes than corner well. At the high-water cut stage of water flooding with nine-spot well pattern, the oil recovery can be further improved with staggered line-drive pattern by converting the corner well into injection well. It is helpful to increase the oil production of corner well of nine-spot well pattern by increasing the injection rate, and improve ultimate oil recovery, but the water-free production period would be greatly shortened and water-free recovery would decrease.展开更多
Based on the characteristics of fractured vuggy porous media,a novel mathematical model was proposed to model fluid flow in such media on fine scale,i.e.,the discrete fracture-vug network model.The new model consists ...Based on the characteristics of fractured vuggy porous media,a novel mathematical model was proposed to model fluid flow in such media on fine scale,i.e.,the discrete fracture-vug network model.The new model consists of three systems:porous rock system,fracture system,and vug system.The fractures and vugs are embedded in porous rock,and the isolated vugs could be connected via the discrete fracture network.The flow in porous rock and fractures follows Darcy's law,and the vugs system is free fluid region.Using a two-scale homogenization limit theory,we obtained a macroscopic Darcy's law governing the media on coarse scale.The theoretical formula of the equivalent permeability of the fractured vuggy porous media was derived.The model and method of this paper were verified by some numerical examples.At the end the permeability of some fractured vuggy porous media with typical fracture-vug structures was analyzed.展开更多
A key aspect in extracting quantitative information from FMI logs is to segment the FMI image to get images of pores, vugs and fractures. A segmentation method based on the dyadic wavelet transform in 2-D is introduce...A key aspect in extracting quantitative information from FMI logs is to segment the FMI image to get images of pores, vugs and fractures. A segmentation method based on the dyadic wavelet transform in 2-D is introduced in this paper. The first step is to find all the edge pixels of the FMI image using the 2-D wavelet transform. The second step is to calculate a segmentation threshold based on the average value of the edge pixels. Field data processing examples show that sub-images of vugs and fractures can be correctly separated from original FMI data continuously and automatically along the depth axis. The image segmentation lays the foundation for in-situ parameter calculation.展开更多
Fractured karst geothermal reservoir is a kind of typical geothermal reservoirs with the advantages of abundant storage water and easy reinjection of tail water during the period of geothermal utilization.Such geother...Fractured karst geothermal reservoir is a kind of typical geothermal reservoirs with the advantages of abundant storage water and easy reinjection of tail water during the period of geothermal utilization.Such geothermal system is also one of the geothermal reservoirs with the greatest potential for the development and utilization of geothermal energy in China.However,its geological structures are diverse(e.g.pore,fracture and vug),exhibiting complex characteristics of multiple scales,strong heterogeneity and various flow regimes.Therefore,the fluid-heat transfer processes and geothermal production performance of fractured karst geothermal reservoirs are not clarified.In this paper,a numerical model considering thermoehydraulic coupling processes based on the discrete fractureevug network approach is put forward,according to the characteristics of fracturedevuggy geothermal reservoirs.In addition,the accuracy of the numerical model is verified.The results obtained from this research are as follows.First,the numerical model considering the thermoehydraulic coupling process is put forward,in which the Darcy's law is used to describe the flow zone of porous medium,the NaviereStokes equation is used to illustrate the free flow zone of vugs,and the BeaverseJosepheSaffman boundary condition is used to couple the fluid flow between these two zones.Second,the connectivity of fracture network is the key parameter to control and evaluate the flow and heat transfer effects in fractured vuggy geothermal reservoirs.The existence of vugs plays an important role in the fluid flow and heat transfer in geothermal reservoirs.Third,the thermoehydraulic coupling model based on the discrete fractureevug network can effectively describe the fluid flow and heat transfer processes in fractured vuggy geothermal reservoirs.The connectivity of fracture networks controls the thermoehydraulic coupling processes in fractured vuggy geothermal reservoirs.Fourth,the existence of vugs seriously impacts the thermoehydraulic coupling processes in geothermal reservoirs.For instance,on the one hand,it increases the number of high-speed flow channels spanning across the system and even makes the system get connected.On the other hand,it increases the speed of local flow channels inside the system.In conclusion,this proposed method is of great significance for studying the development characteristics and optimizing their geothermal production performance of fractured vuggy geothermal reservoirs.展开更多
基金Supported by the National Natural Science Foundation of China(U1562217)
文摘The pore structure characteristic parameters of vuggy carbonate reservoirs were extracted,and matrix resistivity of vuggy reservoir was calculated by the percolation network simulation.A segmented cross-scale resistivity calculation method was established,in which the finite element method was used to simulate the resistivity of vuggy reservoirs.The mathematical models of vug porosity and water saturation with rock resistivity in vuggy carbonate reservoir were established,and the relationships between them were obtained.Experimental results verified the reliability of the simulation results.The method presented provides new technical means and research method for the resistivity log interpretation of vuggy carbonate reservoirs.The matrix porosity,vug porosity and matrix pore water saturation are the key factors determining the resistivity of reservoir rocks.
基金granted access to projects supported by the National Major Fundamental Research Program of China ‘‘On basic research problems in applied geophysics for deep oil and gas fields’’(Grant Number 2013CB228605)CNPC Science and Technology Project(Grant Number 2016A-3303)and CNPC Logging Project(Grant Number 2017E-15)
文摘We present a path morphology method to separate total rock pore space into matrix, fractures and vugs and derive their pore structure spectrum. Thus, we can achieve fine pore evaluation in fracture–vug reservoirs based on electric imaging logging data. We automatically identify and extract fracture–vug information from the electric imaging images by adopting a path morphological operator that remains flexible enough to fit rectilinear and slightly curved structures because they are independent of the structuring element shape. The Otsu method was used to extract fracture–vug information from the background noise caused by the matrix. To accommodate the differences in scale and form of the different target regions,including fracture and vug path, operators with different lengths were selected for their recognition and extraction at the corresponding scale. Polynomial and elliptic functions are used to fit the extracted fractures and vugs, respectively, and the fracture–vug parameters are deduced from the fitted edge. Finally, test examples of numerical simulation data and several measured well data have been provided for the verification of the effectiveness and adaptability of the path morphology method in the application of electric imaging logging data processing. This also provides algorithm support for the fine evaluation of fracture–vug reservoirs.
基金the financial support received from the College of Petroleum Engineering and Geosciences at KFUPM through the project SF20006 toward the completion of this work。
文摘A review of the pressure transient analysis of flow in reservoirs having natural fractures,vugs and/or caves is presented to provide an insight into how much knowledge has been acquired about this phenomenon and to highlight the gaps still open for further research.A comparison-based approach is adopted which involved the review of works by several authors and identifying the limiting assumptions,model restrictions and applicability.Pressure transient analysis provides information to aid the identification of important features of reservoirs.It also provides an explanation to complex reservoir pressuredependent variations which have led to improved understanding and optimization of the reservoir dynamics.Pressure transient analysis techniques,however,have limitations as not all its models find application in naturally fractured and vuggy reservoirs as the flow dynamics differ considerably.Pollard’s model presented in 1953 provided the foundation for existing pressure transient analysis in these types of reservoirs,and since then,several authors have modified this basic model and come up with more accurate models to characterize the dynamic pressure behavior in reservoirs with natural fractures,vugs and/or caves,with most having inherent limitations.This paper summarizes what has been done,what knowledge is considered established and the gaps left to be researched on.
基金Supported by the China National Science and Technology Major Project(2011ZX05020-008)
文摘Considering the fluid flow non-darcy characteristics in fracture-vug carbonate reservoirs, a new multi-scale conduit flow model production prediction method for fracture-vug carbonate reservoirs was presented using image segmentation technique of electric imaging logging data. Firstly, based on Hagen-Poiseuille's law of incompressible fluid flow and the different cross-sectional areas in single fractures and vugs in carbonate reservoirs, a multi-scale conduit flow model for fracture-vug carbonate reservoir was established, and a multi-scale conduit radial fluid flow equation was deduced. Then, conduit flow production index was introduced. The conduit flow production index was calculated using fracture-vug area extracted from the result of electrical image segmentation. Finally, production prediction of fracture-vug carbonate reservoir was realized by using electric imaging logging data. The method has been applied to Ordovician fracture-vug carbonate reservoirs in the Tabei area, and the predicted results are in good agreement with the oil testing data.
基金Supported by the China National Science and Technology Major Project(2016ZX05014-003-004)
文摘Based on the similarity criterion, volcanic rock samples were taken from outcrops to make experimental models. Water flooding experiments of five-spot well pattern, nine-spot well pattern, five-spot to nine-spot well pattern, the relationship between relative well and fracture positions, and injection rate in dissolution vug-cave reservoirs with/without fractures were carried out to analyze variation regularities of development indexes, find out development characteristics of water flooding with different well patterns and sort out the optimal water flooding development mode. For dissolution vug-cave reservoirs without fractures, five-spot well pattern waterflooding has very small sweeping area, serious water channeling and low oil recovery. When the well pattern was adjusted from five-spot to nine-spot well pattern, oil recovery could be largely improved, but the corner well far from the injector is little affected. In dissolution vug-cave reservoirs with fractures, when the injector and producer are not connected by fractures, the fractures could effectively connect the poorly linked vugs to improve the development effect of water flooding. Whether there are fractures or not in dissolution vug-cave reservoirs, the development effect of nine-spot well-pattern is much better than that of five-spot well pattern and five-spot to nine-spot well pattern, this is more evident when there are fractures, and the edge well has better development indexes than corner well. At the high-water cut stage of water flooding with nine-spot well pattern, the oil recovery can be further improved with staggered line-drive pattern by converting the corner well into injection well. It is helpful to increase the oil production of corner well of nine-spot well pattern by increasing the injection rate, and improve ultimate oil recovery, but the water-free production period would be greatly shortened and water-free recovery would decrease.
基金supported by the National Basic Research Program of China("973"Program)(Grant No.2006CB202404)
文摘Based on the characteristics of fractured vuggy porous media,a novel mathematical model was proposed to model fluid flow in such media on fine scale,i.e.,the discrete fracture-vug network model.The new model consists of three systems:porous rock system,fracture system,and vug system.The fractures and vugs are embedded in porous rock,and the isolated vugs could be connected via the discrete fracture network.The flow in porous rock and fractures follows Darcy's law,and the vugs system is free fluid region.Using a two-scale homogenization limit theory,we obtained a macroscopic Darcy's law governing the media on coarse scale.The theoretical formula of the equivalent permeability of the fractured vuggy porous media was derived.The model and method of this paper were verified by some numerical examples.At the end the permeability of some fractured vuggy porous media with typical fracture-vug structures was analyzed.
基金The research was supported by the FifteenthNational Scientific and Technological Project (2001-BA605A-03-02)
文摘A key aspect in extracting quantitative information from FMI logs is to segment the FMI image to get images of pores, vugs and fractures. A segmentation method based on the dyadic wavelet transform in 2-D is introduced in this paper. The first step is to find all the edge pixels of the FMI image using the 2-D wavelet transform. The second step is to calculate a segmentation threshold based on the average value of the edge pixels. Field data processing examples show that sub-images of vugs and fractures can be correctly separated from original FMI data continuously and automatically along the depth axis. The image segmentation lays the foundation for in-situ parameter calculation.
基金Major science and technology projects of PetroChina“Study on production mechanism and reservoir numerical simulation technology of deep carbonate reservoir”(Serial number:ZD2019-183-008)Project sponsored by the National Natural Science Foundation of China“Study on fluid-solid coupling seepage theory and numerical simulation of condensate gas reservoir in fractured-vuggy carbonate rocks”(Serial number:52074336)Key project of the National Natural Science Foundation of China“Study on heat and mass transfer mechanism in coupled heat,fluid,dynamic and chemical fields of enhanced geothermal system”(Serial number:51936001).
文摘Fractured karst geothermal reservoir is a kind of typical geothermal reservoirs with the advantages of abundant storage water and easy reinjection of tail water during the period of geothermal utilization.Such geothermal system is also one of the geothermal reservoirs with the greatest potential for the development and utilization of geothermal energy in China.However,its geological structures are diverse(e.g.pore,fracture and vug),exhibiting complex characteristics of multiple scales,strong heterogeneity and various flow regimes.Therefore,the fluid-heat transfer processes and geothermal production performance of fractured karst geothermal reservoirs are not clarified.In this paper,a numerical model considering thermoehydraulic coupling processes based on the discrete fractureevug network approach is put forward,according to the characteristics of fracturedevuggy geothermal reservoirs.In addition,the accuracy of the numerical model is verified.The results obtained from this research are as follows.First,the numerical model considering the thermoehydraulic coupling process is put forward,in which the Darcy's law is used to describe the flow zone of porous medium,the NaviereStokes equation is used to illustrate the free flow zone of vugs,and the BeaverseJosepheSaffman boundary condition is used to couple the fluid flow between these two zones.Second,the connectivity of fracture network is the key parameter to control and evaluate the flow and heat transfer effects in fractured vuggy geothermal reservoirs.The existence of vugs plays an important role in the fluid flow and heat transfer in geothermal reservoirs.Third,the thermoehydraulic coupling model based on the discrete fractureevug network can effectively describe the fluid flow and heat transfer processes in fractured vuggy geothermal reservoirs.The connectivity of fracture networks controls the thermoehydraulic coupling processes in fractured vuggy geothermal reservoirs.Fourth,the existence of vugs seriously impacts the thermoehydraulic coupling processes in geothermal reservoirs.For instance,on the one hand,it increases the number of high-speed flow channels spanning across the system and even makes the system get connected.On the other hand,it increases the speed of local flow channels inside the system.In conclusion,this proposed method is of great significance for studying the development characteristics and optimizing their geothermal production performance of fractured vuggy geothermal reservoirs.
