The positive and negative difference of deep and shallow resisitivity in formation without invasion is caused mainly by the fixed two tool constants from numerical simulations. A dynamic calibration method for tool co...The positive and negative difference of deep and shallow resisitivity in formation without invasion is caused mainly by the fixed two tool constants from numerical simulations. A dynamic calibration method for tool constants is proposed based on the effects of the mud and formation resisitivity ratio on the two constants calculated using the finite element method （FEM）. Finally, four specific examples are given to validate the dynamic calibration method. It is an automatic borehole correction method and can give more accurate formation resistivity. The method is useful for dual laterolog logging.展开更多
Cave carbonate formations are characterized by heterogeneity, which makes electrical log prediction difficult. It is currently important to know how to use the dual laterolog to accurately identify and quantitatively ...Cave carbonate formations are characterized by heterogeneity, which makes electrical log prediction difficult. It is currently important to know how to use the dual laterolog to accurately identify and quantitatively evaluate caves. Using numerical simulation to calculate electrical log responses can provide a theoretical basis for cave identification and evaluation. In this paper, based on the dual laterolog principles, we first study different size spherical cave models using the finite element method （FEM）, determine a relation between resistivity and cave filling after comprehensively studying the log responses of cave models with different filling material, and finally study the dual laterolog responses on caves filled with shale, limestone, conglomerate, and thin laminated formation of sand and shale. The numerical results provide a theoretical basis for identification and evaluation of carbonate cave reservoirs.展开更多
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 lm, 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.展开更多
There exist different response characteristics in the resistivity measurements of dual laterolog （DLL） and logging while drilling （LWD） electromagnetic wave propagation logging in highly deviated and horizontal we...There exist different response characteristics in the resistivity measurements of dual laterolog （DLL） and logging while drilling （LWD） electromagnetic wave propagation logging in highly deviated and horizontal wells due to the difference in their measuring principles. In this study, we first use the integral equation method simulated the response characteristics of LWD resistivity and use the three dimensional finite element method （3D-FEM） simulated the response characteristics of DLL resistivity in horizontal wells, and then analyzed the response differences between the DLL and LWD resistivity. The comparative analysis indicated that the response differences may be caused by different factors such as differences in the angle of instrument inclination, anisotropy, formation interface, and mud intrusion. In the interface, the curves of the LWD resistivity become sharp with increases in the deviation while those of the DLL resistivity gradually become smooth. Both curves are affected by the anisotropy although the effect on DLL resistivity is lower than the LWD resistivity. These differences aid in providing a reasonable explanation in the horizontal well. However, this can also simultaneously lead to false results. At the end of the study, we explain the effects of the differences in the interpretation of the horizontal well based on the results and actual data analysis.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.11134011)
文摘The positive and negative difference of deep and shallow resisitivity in formation without invasion is caused mainly by the fixed two tool constants from numerical simulations. A dynamic calibration method for tool constants is proposed based on the effects of the mud and formation resisitivity ratio on the two constants calculated using the finite element method （FEM）. Finally, four specific examples are given to validate the dynamic calibration method. It is an automatic borehole correction method and can give more accurate formation resistivity. The method is useful for dual laterolog logging.
基金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)
文摘Cave carbonate formations are characterized by heterogeneity, which makes electrical log prediction difficult. It is currently important to know how to use the dual laterolog to accurately identify and quantitatively evaluate caves. Using numerical simulation to calculate electrical log responses can provide a theoretical basis for cave identification and evaluation. In this paper, based on the dual laterolog principles, we first study different size spherical cave models using the finite element method （FEM）, determine a relation between resistivity and cave filling after comprehensively studying the log responses of cave models with different filling material, and finally study the dual laterolog responses on caves filled with shale, limestone, conglomerate, and thin laminated formation of sand and shale. The numerical results provide a theoretical basis for identification and evaluation of carbonate cave reservoirs.
基金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 lm, 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.
基金supported by the National Science and Technology Major Project of China(Nos.2016ZX05014-002-001,2016ZX05002-005-001,and 2017ZX05005-005-005)
文摘There exist different response characteristics in the resistivity measurements of dual laterolog （DLL） and logging while drilling （LWD） electromagnetic wave propagation logging in highly deviated and horizontal wells due to the difference in their measuring principles. In this study, we first use the integral equation method simulated the response characteristics of LWD resistivity and use the three dimensional finite element method （3D-FEM） simulated the response characteristics of DLL resistivity in horizontal wells, and then analyzed the response differences between the DLL and LWD resistivity. The comparative analysis indicated that the response differences may be caused by different factors such as differences in the angle of instrument inclination, anisotropy, formation interface, and mud intrusion. In the interface, the curves of the LWD resistivity become sharp with increases in the deviation while those of the DLL resistivity gradually become smooth. Both curves are affected by the anisotropy although the effect on DLL resistivity is lower than the LWD resistivity. These differences aid in providing a reasonable explanation in the horizontal well. However, this can also simultaneously lead to false results. At the end of the study, we explain the effects of the differences in the interpretation of the horizontal well based on the results and actual data analysis.
