Due to the tremendous amount of high-resolution measurement information,array laterolog is widely used in evaluations of deviated anisotropic reservoirs.However,the precision of a complementary numerical simulation sh...Due to the tremendous amount of high-resolution measurement information,array laterolog is widely used in evaluations of deviated anisotropic reservoirs.However,the precision of a complementary numerical simulation should be improved as high as the core of fine-scale reservoir evaluation.Therefore,the 3D finite element method(3D-FEM)is presented to simulate the array laterolog responses.Notably,a downscaled physical simulation system is introduced to validate and calibrate the precision of the 3D-FEM.First,the size of the downscaled system is determined by COMSOL.Then,the surrounding and investigated beds are represented by a sodium chloride solution and planks soaked in solution,respectively.Finally,a half-space measurement scheme is presented to improve the experimental efficiency.Moreover,the corresponding sensitivity function and separation factor are established to analyze the effects of the formation anisotro py and dipping angle on the array laterolog responses.The numerical and experimental results indicate that the half-space method is practical,and the mean relative error between the numerical and experimental results is less than 5%,which indicates that the numerical simulation is accurate.With the proposed approach,the reversal angle of array laterolog response curves in anisotropic formations can be observed,and this range is determined to be 50°-62°.展开更多
A scaling-down experiment system of array laterolog resistivity was developed, and a corresponding formation model was built by 3 D finite element numerical method to study the effect of different factors on the loggi...A scaling-down experiment system of array laterolog resistivity was developed, and a corresponding formation model was built by 3 D finite element numerical method to study the effect of different factors on the logging response quantitatively. The error between the experimental and numerical results was less than 5%, validating the reliability of the numerical simulation method. The single factor analysis of the formation relative dip, resistivity anisotropy and drilling fluid invasion was carried out by numerical simulation method, and the results show that:(1) The increase of relative dip can lead to the increase of formation resistivity, but the increasing value is relatively small, and the values of five array resistivity curves will reverse when the relative dip angle reaches a certain degree.(2) The increase of anisotropy coefficient λ can also cause the formation resistivity to rise, and the resistivity will increase by about 10% when λ increases from 1.0 to 1.5 in vertical wells.(3) Drilling fluid invasion has a more significant effect on the logging response than the former two factors. The order of the five curves will change due to drilling fluid invasion in anisotropic formation and the change rule is contrary to resistivity anisotropy. Taking the logging data of the Yingxi oilfield in the Qaidam Basin as an example, an anisotropic formation model considering drilling fluid invasion was built, and the numerical simulation results from the above methods were basically consistent with the logging data, which verified the accuracy of the method again. The results of this study lay a theoretical foundation for multiple-parameter inversion in anisotropic formation under complex well conditions.展开更多
The estimation of fractures is key to evaluating fractured carbonate reservoirs. It is difficult to evaluate this kind of reservoir because of its heterogeneously distributed fractures and anisotropy. A three-dimensio...The estimation of fractures is key to evaluating fractured carbonate reservoirs. It is difficult to evaluate this kind of reservoir because of its heterogeneously distributed fractures and anisotropy. A three-dimensional numerical model was used to simulate the responses of the dual laterolog (DLL) in a fractured formation based on a macro-isotropic anisotropic model. Accordingly, a fast fracture- computing method was developed. First, the apparent conductivity of the DLL is linearly related to the porosity of the fracture and the conductivity of pore fluid. Second, the amplitude difference of the deep and shallow apparent resistivity logs is mainly dependent on the dip angle of the fracture. Then the response of the DLL to a formation with dip angle fractures is approximately depicted as a function of the bulk resistivity of the rock, the porosity of the fractures and the conductivity of fracture fluid. This function can be used to compute the porosity of fracture quickly. The actual data show that the fracture parameters determined by the DLL closely coincide with the formation micro imager log.展开更多
In order to investigate the influences of caliper, formation thickness and invaded zone on the form of dual laterologs, forward modeling technique were applied to calculate the dual laterologs for different cases. The...In order to investigate the influences of caliper, formation thickness and invaded zone on the form of dual laterologs, forward modeling technique were applied to calculate the dual laterologs for different cases. The result shows that the resistivity logs become smoother and lower as the borehole diameter increases, the increase of the contrast between mud resistivity and formation resistivity induce the logs to be more pointed. When the formation thickness is less than 1m, the two-peak on the logs for resistive invasion vanished, and for thickness between 1 m and 4 m, the form of logs does not vary significantly. If the formation thickness is greater than 4 m, a platform appears on the logs at the middle of the formation. The thinner the invaded zone is, the more obvious the invasion feature on the laterologs is. For thick invaded zone the form of logs tend to be that of an uninvaded resistive formation. The form and amplitude of logs depend on the resistivity contrast between invaded zone, uninvaded formation and adjacentlayers.展开更多
Fracture porosity is one of the key parameters for characterizing fractured reservoirs.However,fracture porosity calculation is difficult with conventional logging data due to severe anisotropy of the reservoirs.To de...Fracture porosity is one of the key parameters for characterizing fractured reservoirs.However,fracture porosity calculation is difficult with conventional logging data due to severe anisotropy of the reservoirs.To deal with the problem,the equivalent macroscopic anisotropic formation model based on dual laterolog(DLL)data is adopted to cyclically assign such parameters as bedrock resistivity(RB),fluid resistivity in fractures(RFL),fracture dip angle(FDA)and fracture thickness as well as fracture spacing,and to produce massive data for formation modeling.A large number of training data obtained through three dimensional finite element forward modeling and the functional relationship between DLL responses and fracture parameters that are trained and summarized by deep neural network,are combined to establish a new fast forward model for calculating DLL responses in fractured formations.A new fracture porosity inversion model for fractured reservoirs based on gradient optimization inversion algorithm combined with multi-initial inversion strategy is then proposed.While running the model,formation is divided into eight intervals according to bedrock resistivity and fracture dip angle from 0°to 90°is divided every 0.5°to improve the operation speed and efficiency.The results of numerical verification show that when bedrock resistivity is greater than 1000Ωm,the mean absolute error(MAE)of fracture porosity inversion is 0.001658%for horizontal fractures,0.00413%for intermediate fractures and 0.0027%for quasi-vertical fractures.When bedrock resistivity is between 100Ωm and 1000Ωm,MAE of fracture porosity inversion is 0.003%for horizontal fractures,0.0034%for intermediate fractures and 0.00348%for quasi-vertical fractures.Fracture parameters determined by the fracture porosity inversion model with actual data are in good agreement with the results of micro resistivity imaging logging.展开更多
基金funded by the National Natural Science Foundation of China(41974146,42074134)the Graduate Innovation Project of China University of Petroleum(East China)(YCX2021005)。
文摘Due to the tremendous amount of high-resolution measurement information,array laterolog is widely used in evaluations of deviated anisotropic reservoirs.However,the precision of a complementary numerical simulation should be improved as high as the core of fine-scale reservoir evaluation.Therefore,the 3D finite element method(3D-FEM)is presented to simulate the array laterolog responses.Notably,a downscaled physical simulation system is introduced to validate and calibrate the precision of the 3D-FEM.First,the size of the downscaled system is determined by COMSOL.Then,the surrounding and investigated beds are represented by a sodium chloride solution and planks soaked in solution,respectively.Finally,a half-space measurement scheme is presented to improve the experimental efficiency.Moreover,the corresponding sensitivity function and separation factor are established to analyze the effects of the formation anisotro py and dipping angle on the array laterolog responses.The numerical and experimental results indicate that the half-space method is practical,and the mean relative error between the numerical and experimental results is less than 5%,which indicates that the numerical simulation is accurate.With the proposed approach,the reversal angle of array laterolog response curves in anisotropic formations can be observed,and this range is determined to be 50°-62°.
基金Supported by the Scientific Research and Technological Development Project of CNPC(2019A-3608).
文摘A scaling-down experiment system of array laterolog resistivity was developed, and a corresponding formation model was built by 3 D finite element numerical method to study the effect of different factors on the logging response quantitatively. The error between the experimental and numerical results was less than 5%, validating the reliability of the numerical simulation method. The single factor analysis of the formation relative dip, resistivity anisotropy and drilling fluid invasion was carried out by numerical simulation method, and the results show that:(1) The increase of relative dip can lead to the increase of formation resistivity, but the increasing value is relatively small, and the values of five array resistivity curves will reverse when the relative dip angle reaches a certain degree.(2) The increase of anisotropy coefficient λ can also cause the formation resistivity to rise, and the resistivity will increase by about 10% when λ increases from 1.0 to 1.5 in vertical wells.(3) Drilling fluid invasion has a more significant effect on the logging response than the former two factors. The order of the five curves will change due to drilling fluid invasion in anisotropic formation and the change rule is contrary to resistivity anisotropy. Taking the logging data of the Yingxi oilfield in the Qaidam Basin as an example, an anisotropic formation model considering drilling fluid invasion was built, and the numerical simulation results from the above methods were basically consistent with the logging data, which verified the accuracy of the method again. The results of this study lay a theoretical foundation for multiple-parameter inversion in anisotropic formation under complex well conditions.
基金supported by State Key Laboratory of Petroleum Resources and Prospecting (Grant No. 2009009)Basic Scientific Research Special Fund of Central Colleges (Grant No. 2010ZY28)National Major Projects (Grant No.2008ZX05023-005 and 2008ZX05035-002)
基金This paper is supported by the National High Technology Research andDevelopment Programof China (No .863-306-ZT01-02-3) .
文摘The estimation of fractures is key to evaluating fractured carbonate reservoirs. It is difficult to evaluate this kind of reservoir because of its heterogeneously distributed fractures and anisotropy. A three-dimensional numerical model was used to simulate the responses of the dual laterolog (DLL) in a fractured formation based on a macro-isotropic anisotropic model. Accordingly, a fast fracture- computing method was developed. First, the apparent conductivity of the DLL is linearly related to the porosity of the fracture and the conductivity of pore fluid. Second, the amplitude difference of the deep and shallow apparent resistivity logs is mainly dependent on the dip angle of the fracture. Then the response of the DLL to a formation with dip angle fractures is approximately depicted as a function of the bulk resistivity of the rock, the porosity of the fractures and the conductivity of fracture fluid. This function can be used to compute the porosity of fracture quickly. The actual data show that the fracture parameters determined by the DLL closely coincide with the formation micro imager log.
文摘In order to investigate the influences of caliper, formation thickness and invaded zone on the form of dual laterologs, forward modeling technique were applied to calculate the dual laterologs for different cases. The result shows that the resistivity logs become smoother and lower as the borehole diameter increases, the increase of the contrast between mud resistivity and formation resistivity induce the logs to be more pointed. When the formation thickness is less than 1m, the two-peak on the logs for resistive invasion vanished, and for thickness between 1 m and 4 m, the form of logs does not vary significantly. If the formation thickness is greater than 4 m, a platform appears on the logs at the middle of the formation. The thinner the invaded zone is, the more obvious the invasion feature on the laterologs is. For thick invaded zone the form of logs tend to be that of an uninvaded resistive formation. The form and amplitude of logs depend on the resistivity contrast between invaded zone, uninvaded formation and adjacentlayers.
基金This work was financially supported by the National Natural Science Foundation of China(NSFC)Basic Research Program on Deep Petroleum Resource Accumulation and Key Engineering Technologies(Grant No.U19B6003-04-03-03)State Key Laboratory of Shale Oil and Gas Enrichment Mechanisms and Effective Development Projects(No.20-YYGZ-KF-GC-11)+1 种基金the Strategic Priority Research program of the Chinese Academy of Sciences(Grant No.XDA14010101)the National Science and Technology Major Project(Grant No.2017ZX05005005-005 and 2016ZX05014002-001).
文摘Fracture porosity is one of the key parameters for characterizing fractured reservoirs.However,fracture porosity calculation is difficult with conventional logging data due to severe anisotropy of the reservoirs.To deal with the problem,the equivalent macroscopic anisotropic formation model based on dual laterolog(DLL)data is adopted to cyclically assign such parameters as bedrock resistivity(RB),fluid resistivity in fractures(RFL),fracture dip angle(FDA)and fracture thickness as well as fracture spacing,and to produce massive data for formation modeling.A large number of training data obtained through three dimensional finite element forward modeling and the functional relationship between DLL responses and fracture parameters that are trained and summarized by deep neural network,are combined to establish a new fast forward model for calculating DLL responses in fractured formations.A new fracture porosity inversion model for fractured reservoirs based on gradient optimization inversion algorithm combined with multi-initial inversion strategy is then proposed.While running the model,formation is divided into eight intervals according to bedrock resistivity and fracture dip angle from 0°to 90°is divided every 0.5°to improve the operation speed and efficiency.The results of numerical verification show that when bedrock resistivity is greater than 1000Ωm,the mean absolute error(MAE)of fracture porosity inversion is 0.001658%for horizontal fractures,0.00413%for intermediate fractures and 0.0027%for quasi-vertical fractures.When bedrock resistivity is between 100Ωm and 1000Ωm,MAE of fracture porosity inversion is 0.003%for horizontal fractures,0.0034%for intermediate fractures and 0.00348%for quasi-vertical fractures.Fracture parameters determined by the fracture porosity inversion model with actual data are in good agreement with the results of micro resistivity imaging logging.