A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite...A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.展开更多
CO_(2) pre-injection during hydraulic fracturing is an important method for the development of medium to deep heavy oil reservoirs.It reduces the interfacial tension and viscosity of crude oil,enhances its flowability...CO_(2) pre-injection during hydraulic fracturing is an important method for the development of medium to deep heavy oil reservoirs.It reduces the interfacial tension and viscosity of crude oil,enhances its flowability,maintains reservoir pressure,and increases reservoir drainage capacity.Taking the Badaowan Formation as an example,in this study a detailed three-dimensional geomechanical model based on static data from well logging interpretations is elaborated,which can take into account both vertical and horizontal geological variations and mechanical characteristics.A comprehensive analysis of the impact of key construction parameters on Pre-CO_(2) based fracturing(such as cluster spacing and injection volume),is therefore conducted.Thereafter,using optimized construction parameters,a non-structured grid for dynamic development prediction is introduced,and the capacity variations of different production scenarios are assessed.On the basis of the simulation results,reasonable fracturing parameters are finally determined,including cluster spacing,fracturing fluid volume,proppant concentration,and well spacing.展开更多
Krylov subspace projection methods are known to be highly efficient for solving large linear systems. Many different versions arise from different choices to the left and right subspaces. These methods were classified...Krylov subspace projection methods are known to be highly efficient for solving large linear systems. Many different versions arise from different choices to the left and right subspaces. These methods were classified into two groups in terms of the different forms of matrix H-m, the main properties in applications and the new versions of these two types of methods were briefly reviewed, then one of the most efficient versions, GMRES method was applied to oil reservoir simulation. The block Pseudo-Elimination method was used to generate the preconditioned matrix. Numerical results show much better performance of this preconditioned techniques and the GMRES method than that of preconditioned ORTHMIN method, which is now in use in oil reservoir simulation. Finally, some limitations of Krylov subspace methods and some potential improvements to this type of methods are further presented.展开更多
It is known that the pore media characteristics of glutenite reservoirs are different from those of conventional sandstone reservoirs.Low reservoir permeability and naturally developed microfractures make water inject...It is known that the pore media characteristics of glutenite reservoirs are different from those of conventional sandstone reservoirs.Low reservoir permeability and naturally developed microfractures make water injection in this kind of reservoir very difficult.In this study,new exploitation methods are explored.Using a real glutenite reservoir as a basis,a three-dimensional fine geological model is elaborated.Then,combining the model with reservoir performance information,and through a historical fitting analysis,the saturation abundance distribution of remaining oil in the reservoir is determined.It is shown that,using this information,predictions can be made about whether the considered reservoir is suitable for horizontal well fracturing or not.The direction,well length,well spacing and productivity of horizontal well are also obtained.展开更多
Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data ...Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data (such as the distribution of fractures, low permeability zones, faults and major stress) and simulation experiences to design the grid direction of simulation model qualitatively. The research of the paper indicates that the key to determine the grid direction is to determine the principal permeability direction. Under the circumstances of few static materials, a new grid direction determination method has been developed by using field data (well location map and inter-well permeability) on the bases of Darcy’s law and tensor analysis theory. The grid direction of WZ11-7 Oilfield simulation model has been determined using four production wells and two production zones (L1 and L3) in WZ11-7-2 well group, the results are in conformity with the geological studied major stress. Therefore, this method can give insights into the numerical simulation study.展开更多
The exploitation of an oil field is a complex and multidisciplinary task, which demands a lot of prior knowledge, time, and money. A good reservoir characterization is deemed essential in the accomplishment of Enhance...The exploitation of an oil field is a complex and multidisciplinary task, which demands a lot of prior knowledge, time, and money. A good reservoir characterization is deemed essential in the accomplishment of Enhanced Oil Recovery (EOR) processes in order to estimate accurately the properties of the porous medium affecting the flow properties. Several techniques at a field scale are currently being used to determine these properties, which are time and money consuming. But these alone do not guarantee the success of the project. Reservoir simulation and numerical techniques were then included in the pre-development and follow-up studies as an effective tool to determine the productivity and future behavior of the oil field. As the computational power increased, more advanced and detailed models were developed, including different chemical and physical phenomena. But alongside this process, there was an active research in the area of reservoir simulation, improving the accuracy and efficiency of the numerical schemes used for the flow, transport, and energy equations. The aim of this review is to address the topics described. Firstly, the origin of an oil recovery process, the economic factors and field tests involved are introduced. Secondly, the oil and porous medium origin and characterization as well as an introduction to the fundamental concepts and equations are associated to reservoir simulation. Finally, a brief description and analysis of the techniques are used in reservoir simulation employing finite difference methods, their downsides and possible ways to overcome these problems.展开更多
In the last decades, the numerical analysis has abridged the time and cost of oil extractive operations, and thus their implementation has been promoted so much that today is used to solve a large number of phenomena ...In the last decades, the numerical analysis has abridged the time and cost of oil extractive operations, and thus their implementation has been promoted so much that today is used to solve a large number of phenomena that occur during the extraction of hydrocarbons, not only economizing the operation, but also shedding light to the phenomenon, and helping to have a better understanding. Nevertheless in spite of the success of the applications of numerical simulations, it has been until today, looking for their optimization, because the discretization methods of the domain and the model are through meshes, which require modifying the physics of the problem so that this matches the numerical method, and not in backwards.展开更多
Three dimensional geophysical models were abstracted and established according to characteristics of oil and gas reservoir.Then direct current fields for different models were simulated with finite element software(fi...Three dimensional geophysical models were abstracted and established according to characteristics of oil and gas reservoir.Then direct current fields for different models were simulated with finite element software(finite element program generator) by hole-to-surface resistivity method.Numerical solution was compared with analytical one for the homogeneity earth model.And a new parameter of deviation rate was proposed by analyzing different plot curves.The results show that the relative error of solution for homogeneity earth model may attain to 0.043%.And deviation rate decreases from 18% to 1% and its anomaly range becomes wide gradually when the depth of oil and gas reservoir increases from 200 to 1 500 m.If resistivity ratio of oil and gas reservoir to sur-rounding rock decreases from 100 to 10 for the resistive oil and gas reservoir,the amplitude attenuation of deviation rate nearly reaches 8%.When there exists stratum above oil and gas reservoir,and influence of resistive stratum may be eliminated or weakened and anomaly of oil and gas reservoir can be strengthened.展开更多
After long-term waterflooding in unconsolidated sandstone reservoir, the high-permeability channels are easy to evolve, which leads to a significant reduction in water flooding efficiency and a poor oilfield developme...After long-term waterflooding in unconsolidated sandstone reservoir, the high-permeability channels are easy to evolve, which leads to a significant reduction in water flooding efficiency and a poor oilfield development effect. The current researches on the formation parameters variation are mainly based on the experiment analysis or field statistics, while lacking quantitative research of combining microcosmic and macroscopic mechanism. A network model was built after taking the detachment and entrapment mechanisms of particles in unconsolidated sandstone reservoir into consideration. Then a coupled mathematical model for the formation parameters variation was established based on the network modeling and the model of fluids flowing in porous media. The model was solved by a finite-difference method and the Gauss-Seidel iterative technique. A novel field-scale reservoir numerical simulator was written in Fortran 90 and it can be used to predict 1) the evolvement of high-permeability channels caused by particles release and migration in the long-term water flooding process, and 2) well production performances and remaining oil distribution. In addition, a series of oil field examples with inverted nine-spot pattern was made on the new numerical simulator. The results show that the high-permeability channels are more likely to develop along the main streamlines between the injection and production wells, and the formation parameters variation has an obvious influence on the remaining oil distribution.展开更多
The oil production of the multi-fractured horizontal wells(MFHWs) declines quickly in unconventional oil reservoirs due to the fast depletion of natural energy. Gas injection has been acknowledged as an effective meth...The oil production of the multi-fractured horizontal wells(MFHWs) declines quickly in unconventional oil reservoirs due to the fast depletion of natural energy. Gas injection has been acknowledged as an effective method to improve oil recovery factor from unconventional oil reservoirs. Hydrocarbon gas huff-n-puff becomes preferable when the CO_(2) source is limited. However, the impact of complex fracture networks and well interference on the EOR performance of multiple MFHWs is still unclear. The optimal gas huff-n-puff parameters are significant for enhancing oil recovery. This work aims to optimize the hydrocarbon gas injection and production parameters for multiple MFHWs with complex fracture networks in unconventional oil reservoirs. Firstly, the numerical model based on unstructured grids is developed to characterize the complex fracture networks and capture the dynamic fracture features.Secondly, the PVT phase behavior simulation was carried out to provide the fluid model for numerical simulation. Thirdly, the optimal parameters for hydrocarbon gas huff-n-puff were obtained. Finally, the dominant factors of hydrocarbon gas huff-n-puff under complex fracture networks are obtained by fuzzy mathematical method. Results reveal that the current pressure of hydrocarbon gas injection can achieve miscible displacement. The optimal injection and production parameters are obtained by single-factor analysis to analyze the effect of individual parameter. Gas injection time is the dominant factor of hydrocarbon gas huff-n-puff in unconventional oil reservoirs with complex fracture networks. This work can offer engineers guidance for hydrocarbon gas huff-n-puff of multiple MFHWs considering the complex fracture networks.展开更多
Based on the typical dissection of various onshore tight oil fields in China,the tight oil migration and accumulation mechanism and enrichment-controlling factors in continental lake basins are analyzed through nuclea...Based on the typical dissection of various onshore tight oil fields in China,the tight oil migration and accumulation mechanism and enrichment-controlling factors in continental lake basins are analyzed through nuclear magnetic resonance(NMR)displacement physical simulation and Lattice Boltzmann numerical simulation by using the samples of source rock,reservoir rock and crude oil.In continental lake basins,the dynamic forces driving hydrocarbon generation and expulsion of high-quality source rocks are the foundational power that determines the charging efficiency and accumulation effect of tight oil,the oil migration resistance is a key element that influences the charging efficiency and accumulation effect of tight oil,and the coupling of charging force with pore-throat resistance in tight reservoir controls the tight oil accumulation and sweet spot enrichment.The degree of tight oil enrichment in continental lake basins is controlled by four factors:source rock,reservoir pore-throat size,anisotropy of reservoir structure,and fractures.The high-quality source rocks control the near-source distribution of tight oil,reservoir physical properties and pore-throat size are positively correlated with the degree of tight oil enrichment,the anisotropy of reservoir structure reveals that the parallel migration rate is the highest,and intralayer fractures can improve the migration and accumulation efficiency and the oil saturation.展开更多
In east China, almost all the old fields have stepped into high water cut stage, with main development blocks into phase of super high water\|cut and high ratio of total oil produced to OOIP. Because of the complicate...In east China, almost all the old fields have stepped into high water cut stage, with main development blocks into phase of super high water\|cut and high ratio of total oil produced to OOIP. Because of the complicated distribution of remaining oil and great difficulties of potential excavation, effective succeed development methods are in urgent. Through physical modeling in lab, field experiments and industrial deployment, a set of necessary industrialized technology is posted, which is suitable for terrestrial fluvial delta reservoirs during super high water bearing stage, including the detailed description of major reservoir, remaining oil monitoring, fine potential exploitation by water flooding and IOR technology by chemical flooding. This technology has been applied to Shengli, Zhongyuan and Henan oilfields with remarkable benefits and social effects.展开更多
基金Supported by the National Science and Technology Major Project (2017ZX05013-005)。
文摘A mathematical model, fully coupling multiple porous media deformation and fluid flow, was established based on the elastic theory of porous media and fluid-solid coupling mechanism in tight oil reservoirs. The finite element method was used to determine the numerical solution and the accuracy of the model was verified. On this basis, the model was used to simulate productivity of multistage fractured horizontal wells in tight oil reservoirs. The results show that during the production of tight oil wells, the reservoir region close to artificial fractures deteriorated in physical properties significantly, e.g. the aperture and conductivity of artificial fractures dropped by 52.12% and 89.02% respectively. The simulations of 3000-day production of a horizontal well in tight oil reservoir showed that the predicted productivity by the uncoupled model had an error of 38.30% from that by the fully-coupled model. Apparently, ignoring the influence of fluid-solid interaction effect led to serious deviations of the productivity prediction results. The productivity of horizontal well in tight oil reservoir was most sensitive to the start-up pressure gradient, and second most sensitive to the opening of artificial fractures. Enhancing the initial conductivity of artificial fractures was helpful to improve the productivity of tight oil wells. The influence of conductivity, spacing, number and length of artificial fractures should be considered comprehensively in fracturing design. Increasing the number of artificial fractures unilaterally could not achieve the expected increase in production.
基金supported by the Cutting-Edge Project Foundation of Petro-China(Cold-Based Method to Enhance Heavy Oil Recovery)(Grant No.2021DJ1406)Open Fund(PLN201802)of National Key Laboratory of Oil and Gas Reservoir Geology and Exploitation(Southwest Petroleum University).
文摘CO_(2) pre-injection during hydraulic fracturing is an important method for the development of medium to deep heavy oil reservoirs.It reduces the interfacial tension and viscosity of crude oil,enhances its flowability,maintains reservoir pressure,and increases reservoir drainage capacity.Taking the Badaowan Formation as an example,in this study a detailed three-dimensional geomechanical model based on static data from well logging interpretations is elaborated,which can take into account both vertical and horizontal geological variations and mechanical characteristics.A comprehensive analysis of the impact of key construction parameters on Pre-CO_(2) based fracturing(such as cluster spacing and injection volume),is therefore conducted.Thereafter,using optimized construction parameters,a non-structured grid for dynamic development prediction is introduced,and the capacity variations of different production scenarios are assessed.On the basis of the simulation results,reasonable fracturing parameters are finally determined,including cluster spacing,fracturing fluid volume,proppant concentration,and well spacing.
文摘Krylov subspace projection methods are known to be highly efficient for solving large linear systems. Many different versions arise from different choices to the left and right subspaces. These methods were classified into two groups in terms of the different forms of matrix H-m, the main properties in applications and the new versions of these two types of methods were briefly reviewed, then one of the most efficient versions, GMRES method was applied to oil reservoir simulation. The block Pseudo-Elimination method was used to generate the preconditioned matrix. Numerical results show much better performance of this preconditioned techniques and the GMRES method than that of preconditioned ORTHMIN method, which is now in use in oil reservoir simulation. Finally, some limitations of Krylov subspace methods and some potential improvements to this type of methods are further presented.
文摘It is known that the pore media characteristics of glutenite reservoirs are different from those of conventional sandstone reservoirs.Low reservoir permeability and naturally developed microfractures make water injection in this kind of reservoir very difficult.In this study,new exploitation methods are explored.Using a real glutenite reservoir as a basis,a three-dimensional fine geological model is elaborated.Then,combining the model with reservoir performance information,and through a historical fitting analysis,the saturation abundance distribution of remaining oil in the reservoir is determined.It is shown that,using this information,predictions can be made about whether the considered reservoir is suitable for horizontal well fracturing or not.The direction,well length,well spacing and productivity of horizontal well are also obtained.
文摘Grid direction selection and grid size design are two important elements that need to be considered in the grid direction design in reservoir numerical simulation. Reservoir engineers normally utilize geological data (such as the distribution of fractures, low permeability zones, faults and major stress) and simulation experiences to design the grid direction of simulation model qualitatively. The research of the paper indicates that the key to determine the grid direction is to determine the principal permeability direction. Under the circumstances of few static materials, a new grid direction determination method has been developed by using field data (well location map and inter-well permeability) on the bases of Darcy’s law and tensor analysis theory. The grid direction of WZ11-7 Oilfield simulation model has been determined using four production wells and two production zones (L1 and L3) in WZ11-7-2 well group, the results are in conformity with the geological studied major stress. Therefore, this method can give insights into the numerical simulation study.
文摘The exploitation of an oil field is a complex and multidisciplinary task, which demands a lot of prior knowledge, time, and money. A good reservoir characterization is deemed essential in the accomplishment of Enhanced Oil Recovery (EOR) processes in order to estimate accurately the properties of the porous medium affecting the flow properties. Several techniques at a field scale are currently being used to determine these properties, which are time and money consuming. But these alone do not guarantee the success of the project. Reservoir simulation and numerical techniques were then included in the pre-development and follow-up studies as an effective tool to determine the productivity and future behavior of the oil field. As the computational power increased, more advanced and detailed models were developed, including different chemical and physical phenomena. But alongside this process, there was an active research in the area of reservoir simulation, improving the accuracy and efficiency of the numerical schemes used for the flow, transport, and energy equations. The aim of this review is to address the topics described. Firstly, the origin of an oil recovery process, the economic factors and field tests involved are introduced. Secondly, the oil and porous medium origin and characterization as well as an introduction to the fundamental concepts and equations are associated to reservoir simulation. Finally, a brief description and analysis of the techniques are used in reservoir simulation employing finite difference methods, their downsides and possible ways to overcome these problems.
基金financial support from the Escuela Superior de Ingenieria Mecanica y Electrica Unidad Ticoman,Instituto Politecnico Nacional,through Project no.20150465.
文摘In the last decades, the numerical analysis has abridged the time and cost of oil extractive operations, and thus their implementation has been promoted so much that today is used to solve a large number of phenomena that occur during the extraction of hydrocarbons, not only economizing the operation, but also shedding light to the phenomenon, and helping to have a better understanding. Nevertheless in spite of the success of the applications of numerical simulations, it has been until today, looking for their optimization, because the discretization methods of the domain and the model are through meshes, which require modifying the physics of the problem so that this matches the numerical method, and not in backwards.
基金Projects(2006AA06Z105,2007AA06Z134) supported by the National High-Tech Research and Development Program of China
文摘Three dimensional geophysical models were abstracted and established according to characteristics of oil and gas reservoir.Then direct current fields for different models were simulated with finite element software(finite element program generator) by hole-to-surface resistivity method.Numerical solution was compared with analytical one for the homogeneity earth model.And a new parameter of deviation rate was proposed by analyzing different plot curves.The results show that the relative error of solution for homogeneity earth model may attain to 0.043%.And deviation rate decreases from 18% to 1% and its anomaly range becomes wide gradually when the depth of oil and gas reservoir increases from 200 to 1 500 m.If resistivity ratio of oil and gas reservoir to sur-rounding rock decreases from 100 to 10 for the resistive oil and gas reservoir,the amplitude attenuation of deviation rate nearly reaches 8%.When there exists stratum above oil and gas reservoir,and influence of resistive stratum may be eliminated or weakened and anomaly of oil and gas reservoir can be strengthened.
文摘After long-term waterflooding in unconsolidated sandstone reservoir, the high-permeability channels are easy to evolve, which leads to a significant reduction in water flooding efficiency and a poor oilfield development effect. The current researches on the formation parameters variation are mainly based on the experiment analysis or field statistics, while lacking quantitative research of combining microcosmic and macroscopic mechanism. A network model was built after taking the detachment and entrapment mechanisms of particles in unconsolidated sandstone reservoir into consideration. Then a coupled mathematical model for the formation parameters variation was established based on the network modeling and the model of fluids flowing in porous media. The model was solved by a finite-difference method and the Gauss-Seidel iterative technique. A novel field-scale reservoir numerical simulator was written in Fortran 90 and it can be used to predict 1) the evolvement of high-permeability channels caused by particles release and migration in the long-term water flooding process, and 2) well production performances and remaining oil distribution. In addition, a series of oil field examples with inverted nine-spot pattern was made on the new numerical simulator. The results show that the high-permeability channels are more likely to develop along the main streamlines between the injection and production wells, and the formation parameters variation has an obvious influence on the remaining oil distribution.
基金funded by the National Natural Science Foundation of China(No.51974268)Open Fund of Key Laboratory of Ministry of Education for Improving Oil and Gas Recovery(NEPUEOR-2022-03)Research and Innovation Fund for Graduate Students of Southwest Petroleum University(No.2022KYCX005)。
文摘The oil production of the multi-fractured horizontal wells(MFHWs) declines quickly in unconventional oil reservoirs due to the fast depletion of natural energy. Gas injection has been acknowledged as an effective method to improve oil recovery factor from unconventional oil reservoirs. Hydrocarbon gas huff-n-puff becomes preferable when the CO_(2) source is limited. However, the impact of complex fracture networks and well interference on the EOR performance of multiple MFHWs is still unclear. The optimal gas huff-n-puff parameters are significant for enhancing oil recovery. This work aims to optimize the hydrocarbon gas injection and production parameters for multiple MFHWs with complex fracture networks in unconventional oil reservoirs. Firstly, the numerical model based on unstructured grids is developed to characterize the complex fracture networks and capture the dynamic fracture features.Secondly, the PVT phase behavior simulation was carried out to provide the fluid model for numerical simulation. Thirdly, the optimal parameters for hydrocarbon gas huff-n-puff were obtained. Finally, the dominant factors of hydrocarbon gas huff-n-puff under complex fracture networks are obtained by fuzzy mathematical method. Results reveal that the current pressure of hydrocarbon gas injection can achieve miscible displacement. The optimal injection and production parameters are obtained by single-factor analysis to analyze the effect of individual parameter. Gas injection time is the dominant factor of hydrocarbon gas huff-n-puff in unconventional oil reservoirs with complex fracture networks. This work can offer engineers guidance for hydrocarbon gas huff-n-puff of multiple MFHWs considering the complex fracture networks.
基金This work is supported in partial by Major State Basic Research Project （No. G19990328, Parallel Computations of the Large-Scale Reservoir Simulation （2003-2004） （Cooperated with China National 0ffshore 0il Corporation）, and National Natural Science Foundation Project （No. 60303020, 2004.1-2006.12）.
基金Supported by the National Science and Technology Major Project of China(2016ZX05046-001).
文摘Based on the typical dissection of various onshore tight oil fields in China,the tight oil migration and accumulation mechanism and enrichment-controlling factors in continental lake basins are analyzed through nuclear magnetic resonance(NMR)displacement physical simulation and Lattice Boltzmann numerical simulation by using the samples of source rock,reservoir rock and crude oil.In continental lake basins,the dynamic forces driving hydrocarbon generation and expulsion of high-quality source rocks are the foundational power that determines the charging efficiency and accumulation effect of tight oil,the oil migration resistance is a key element that influences the charging efficiency and accumulation effect of tight oil,and the coupling of charging force with pore-throat resistance in tight reservoir controls the tight oil accumulation and sweet spot enrichment.The degree of tight oil enrichment in continental lake basins is controlled by four factors:source rock,reservoir pore-throat size,anisotropy of reservoir structure,and fractures.The high-quality source rocks control the near-source distribution of tight oil,reservoir physical properties and pore-throat size are positively correlated with the degree of tight oil enrichment,the anisotropy of reservoir structure reveals that the parallel migration rate is the highest,and intralayer fractures can improve the migration and accumulation efficiency and the oil saturation.
文摘In east China, almost all the old fields have stepped into high water cut stage, with main development blocks into phase of super high water\|cut and high ratio of total oil produced to OOIP. Because of the complicated distribution of remaining oil and great difficulties of potential excavation, effective succeed development methods are in urgent. Through physical modeling in lab, field experiments and industrial deployment, a set of necessary industrialized technology is posted, which is suitable for terrestrial fluvial delta reservoirs during super high water bearing stage, including the detailed description of major reservoir, remaining oil monitoring, fine potential exploitation by water flooding and IOR technology by chemical flooding. This technology has been applied to Shengli, Zhongyuan and Henan oilfields with remarkable benefits and social effects.
