CN110488367B - Resistivity inversion initial value selection method based on array laterolog data - Google Patents
Resistivity inversion initial value selection method based on array laterolog data Download PDFInfo
- Publication number
- CN110488367B CN110488367B CN201910786024.7A CN201910786024A CN110488367B CN 110488367 B CN110488367 B CN 110488367B CN 201910786024 A CN201910786024 A CN 201910786024A CN 110488367 B CN110488367 B CN 110488367B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- initial value
- mlr4
- mlr1
- logging
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000010187 selection method Methods 0.000 title claims abstract description 14
- 238000011010 flushing procedure Methods 0.000 claims abstract description 21
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 18
- 238000000034 method Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims abstract description 4
- 238000012417 linear regression Methods 0.000 claims abstract description 4
- 238000012937 correction Methods 0.000 claims description 7
- 239000011435 rock Substances 0.000 claims description 7
- 101000605076 Homo sapiens Ligand-dependent nuclear receptor corepressor-like protein Proteins 0.000 abstract description 4
- 102100038259 Ligand-dependent nuclear receptor corepressor-like protein Human genes 0.000 abstract description 4
- 230000009286 beneficial effect Effects 0.000 abstract description 4
- 230000001502 supplementing effect Effects 0.000 abstract description 3
- 238000004422 calculation algorithm Methods 0.000 description 14
- 230000009545 invasion Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 230000009977 dual effect Effects 0.000 description 4
- 238000011156 evaluation Methods 0.000 description 4
- 230000006698 induction Effects 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000013528 artificial neural network Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002068 genetic effect Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 239000004005 microsphere Substances 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000002922 simulated annealing Methods 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000005406 washing Methods 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/18—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation specially adapted for well-logging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Mining & Mineral Resources (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Fluid Mechanics (AREA)
- Geochemistry & Mineralogy (AREA)
- Developing Agents For Electrophotography (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a resistivity inversion initial value selection method based on array lateral logging data, which comprises the following steps of: forward calculation is carried out based on the radial step formation model, resistivity amplitude difference coefficients S1, S2, S3, S4, S5 and S6 are defined, and MLR1 in array lateral logging aimed at by the radial step formation model is obtaineda、MLR2a、MLR3a、MLR4aEstablishing a chart of variation of resistivity amplitude difference coefficients along with the radius of a flushing zone according to the apparent resistivities of the four basic logging curves; establishing an initial value r of the radius of the flushing zone through linear regression based on a chart of the variation of the resistivity amplitude difference coefficient along with the radius of the flushing zonexoczA quantitative relationship with a resistivity magnitude difference coefficient; MLR4 for basic well logging curve provided by array lateral well loggingaApparent resistivity of the earth layer is used as an initial value R of the formation resistivitytczThen, the initial value R of the resistivity of the flushing zone is obtained based on the pseudo-geometric factor definitionxoczThe method can effectively reduce the deviation between the initial value and the true value without supplementing other logging information, and is beneficial to improving the inversion precision and speed of resistivity logging.
Description
Technical Field
The invention relates to the technical field of logging evaluation in exploration and development of oil and gas resources, in particular to a resistivity inversion initial value selection method based on array lateral logging information.
Background
The reservoir resistivity data has important application in the aspect of reservoir evaluation, and the resistivity can be used for not only qualitatively identifying oil and gas reservoirs (general reservoirs show high-resistance characteristics) but also quantitatively calculating the oil saturation (by using formulas such as Archie or Indonesia and the like). Various resistivity logging methods are widely applied to oil field sites, such as dual lateral logging, dual induction logging, array lateral logging, array induction logging, electrical imaging logging and the like, and with increasingly refined oil and gas reservoir evaluation, the array type resistivity logging instrument can provide a plurality of apparent resistivity curves with higher longitudinal resolution compared with the traditional dual lateral and dual induction logging instrument, so that the array type resistivity logging instrument is increasingly applied to the oil field sites, and the processing work of corresponding array type resistivity logging information becomes a key problem to be solved urgently.
The response of the lateral logging instrument is influenced by factors such as mud invasion, well bores, surrounding rocks and the like, so that the measured apparent resistivity curve cannot truly reflect the resistivity of a reservoir, and inversion processing needs to be carried out on the apparent resistivity curve to recover the resistivity information of a mud invasion section. At present, the processing of the array lateral logging data in the oilfield field is mainly divided into two types. One is to use the apparent resistivity directly as the true resistivity of the reservoir, which has an evaluation error within an acceptable range when the mud invasion of the reservoir is not deep, but causes a serious error at a horizon where the mud invasion is severe. The other processing method is to approximate estimate the resistivity of the invaded zone by using microspheres or micro lateral logging, then correct the resistivity of the invaded zone by using a mud invasion correction plate to obtain the radius of the invaded zone and the resistivity of the stratum.
With the development of computer technology, people introduce inversion ideas into the field of electrical logging inversion. The idea of inversion of the electrical logging is that a resistivity logging inversion model is constructed firstly, repeated iteration is carried out through a corresponding inversion algorithm, and when the error between a forward operator and an actual value is smaller than a preset threshold value, the corresponding forward model parameter at the moment is taken as an inversion parameter to be output. The influence of mud invasion is eliminated by carrying out quick inversion processing on the logging data measured in real time on site, the original characteristics of the stratum can be preliminarily recovered, basic data are provided for later fine interpretation, and the method is also an important basis for logging personnel to judge the quality of logging data. The marquitet algorithm is the most widely applied resistivity inversion method in the oil field. The inversion algorithm has extremely high convergence speed and is beneficial to the real-time inversion of logging information, but the defects of the algorithm are also exposed in the practical application process, that is, the inversion result is greatly influenced by an initial value, and if the initial value is unreasonable in setting (is far away from the true value), the true resistivity value of the stratum is difficult to accurately invert, and the inversion precision is seriously influenced. In order to solve the problems of the marquit algorithm, a learner introduces global optimization algorithms such as a simulated annealing method, a genetic algorithm, an artificial neural network algorithm, an immune algorithm, a differential evolution algorithm, a particle swarm and the like into inversion processing of electric logging information, the algorithm has strong global search capability and can effectively improve inversion accuracy, but has generally slow convergence speed, is not suitable for fast processing of the electric logging information and is limited in use, and the marquit algorithm is still the first choice for real-time inversion of field array lateral logging information.
Disclosure of Invention
The invention aims to solve the technical problem of providing a resistivity inversion initial value selection method based on array lateral logging data.
The technical scheme adopted by the invention for solving the technical problems is as follows: a resistivity inversion initial value selection method based on array laterolog data is constructed, and the method comprises the following steps:
step S100: forward calculation is carried out based on a radial step stratum model, resistivity amplitude difference coefficients S1, S2, S3, S4, S5 and S6 are defined, and MLR1 in array lateral logging aimed at by the radial step stratum model is obtaineda、MLR2a、MLR3a、MLR4aCorrecting the apparent resistivities of the four basic logging curves, performing borehole and surrounding rock correction, and establishing a chart of variation of resistivity amplitude difference coefficients along with the radius of a flushing zone;
step S200: establishing an initial value r of the radius of the flushing zone through linear regression based on a chart of the variation of the resistivity amplitude difference coefficient along with the radius of the flushing zonexoczA quantitative relationship with a resistivity magnitude difference coefficient;
step S300: MLR4 for basic well logging curve provided by array lateral well loggingaApparent resistivity of the earth layer is used as an initial value R of the formation resistivitytczThen, the initial value R of the resistivity of the flushing zone is obtained based on the pseudo-geometric factor definitionxocz。
Preferably, in the step S100, S1, S2, S3, S4, S5, S6 are respectively defined as follows:
S1=2·(MLR4a-MLR1a)/(MLR4a+MLR1a) (1)
S2=2·(MLR4a-MLR2a)/(MLR4a+MLR2a) (2)
S3=2·(MLR4a-MLR3a)/(MLR4a+MLR3a) (3)
S4=2·(MLR3a-MLR1a)/(MLR3a+MLR1a) (4)
S5=2·(MLR3a-MLR2a)/(MLR3a+MLR2a) (5)
S6=2·(MLR2a-MLR1a)/(MLR2a+MLR1a) (6)
wherein, MLR1a、MLR2a、MLR3a、MLR4aThe apparent resistivity of the array laterolog provided for four basic well logging curves is represented, and the unit is omega.m; s1, S2, S3, S4, S5 and S6 are dimensionless.
Preferably, in step S200, an initial value r of the radius of the rinsing zone is establishedxoczThe quantitative relation with the resistivity amplitude difference coefficient is as follows:
rxocz=925.14·S1+1004.04·S2-1835.46·S3-1295.69·S4-221.40·S5+70.02·S6+1.39 (7)。
preferably, in step S300, the array laterolog provided for the base log MLR4aApparent resistivity of the earth layer is used as an initial value R of the formation resistivitytczThen defining the obtained impulse based on the pseudo-geometric factorInitial value of resistivity R of washing beltxocz:
Wherein, K1Is MLR1aA corresponding instrument constant; r is0Is the borehole radius; h is1Is MLR1aThe corresponding main current thickness; k1、r0、h1Are all known values.
Preferably, also comprises
Step S400: and comparing the initial value Rxocz of the resistivity of the flushing zone with the actual formation parameters.
The implementation of the invention has the following beneficial effects: according to the resistivity inversion initial value selection method based on the array lateral logging information, the deviation between the initial value and the true value can be effectively reduced without supplementing other logging information, and the resistivity logging inversion accuracy and speed can be improved.
Drawings
The invention will be further described with reference to the accompanying drawings and examples, in which:
FIG. 1 is a schematic diagram of an initial resistivity inversion selection method based on array laterolog data according to the present invention;
FIG. 2 is a schematic representation of a radial step formation model of the present invention;
FIG. 3 is a graphical representation of a plate of the resistivity magnitude difference coefficient of the present invention as a function of the radius of the washband.
Detailed Description
As shown in fig. 1, the resistivity inversion initial value selection method based on the array laterolog data of the present invention solves the problem that the selection of the array laterolog inversion initial value by using the marquit algorithm is not appropriate, and is beneficial to improving the inversion accuracy and the algorithm convergence speed.
The method comprises the following steps:
step S100: forward calculation is carried out based on the radial step formation model, resistivity amplitude difference coefficients S1, S2, S3, S4, S5 and S6 are defined, and a radial step formation model needle is obtainedPaired array laterolog MLR1a、MLR2a、MLR3a、MLR4aAnd (4) correcting the apparent resistivities of the four basic logging curves, well bores and surrounding rocks, and establishing a chart of the variation of the resistivity amplitude difference coefficient along with the radius of a flushing zone. As shown in fig. 2, the model is a radial step stratigraphic model. FIG. 3 is a chart of resistivity magnitude difference coefficient as a function of wash zone radius.
In this step S100, S1, S2, S3, S4, S5, S6 are defined as follows, respectively:
S1=2·(MLR4a-MLR1a)/(MLR4a+MLR1a) (1)
S2=2·(MLR4a-MLR2a)/(MLR4a+MLR2a) (2)
S3=2·(MLR4a-MLR3a)/(MLR4a+MLR3a) (3)
S4=2·(MLR3a-MLR1a)/(MLR3a+MLR1a) (4)
S5=2·(MLR3a-MLR2a)/(MLR3a+MLR2a) (5)
S6=2·(MLR2a-MLR1a)/(MLR2a+MLR1a) (6)
wherein, MLR1a、MLR2a、MLR3a、MLR4aThe apparent resistivity of the array laterolog provided for four basic well logging curves is represented, and the unit is omega.m; s1, S2, S3, S4, S5 and S6 are dimensionless.
Further, MLR1 for radial step formation model (theoretical model) is obtaineda、MLR2a、MLR3a、MLR4aThe apparent resistivity of the base log was measured and borehole, wall rock corrections were made as in table 1.
TABLE 1 theoretical model array laterolog response
Response MLR1 for array laterologa、MLR2a、MLR3a、MLR4aAnd preferably selecting corresponding borehole and surrounding rock correction charts, and performing automatic correction of a numerical simulation computer on the basis of fitting the chart readings into a correction formula to eliminate the influence of the borehole and the surrounding rock on the response of the array lateral logging so as to avoid the influence of other factors on the accuracy degree of resistivity inversion initial value selection.
Step S200: establishing an initial value r of the radius of the flushing zone through linear regression based on a chart of the variation of the resistivity amplitude difference coefficient along with the radius of the flushing zonexoczAnd the quantitative relation with the resistivity amplitude difference coefficient.
In this step S200, an initial value r of the radius of the rinsing zone is establishedxoczThe quantitative relation with the resistivity amplitude difference coefficient is as follows:
rxocz=925.14·S1+1004.04·S2-1835.46·S3-1295.69·S4-221.40·S5+70.02·S6+1.39 (7)。
wherein r isxoczThe initial value of the radius of the rinsing tape, in.
Step S300: MLR4 for basic well logging curve provided by array lateral well loggingaApparent resistivity of the earth layer is used as an initial value R of the formation resistivitytczThen, the initial value R of the resistivity of the flushing zone is obtained based on the pseudo-geometric factor definitionxocz。
In this step S300, the array laterolog provided for the base log MLR4aApparent resistivity of the earth layer is used as an initial value R of the formation resistivitytczThen, the initial value R of the resistivity of the flushing zone is obtained based on the pseudo-geometric factor definitionxocz:
Wherein, K1Is MLR1aA corresponding instrument constant; r is0Is the borehole radius; h is1Is MLR1aThe corresponding main current thickness; k1、r0、h1Are all known values.
In this embodiment, it further comprises
Step S400: initial resistivity value R of the flushing zonexoczCompared to actual formation parameters. It can be understood that the inversion initial value is compared with the actual formation parameter, the effectiveness and the superiority of the resistivity inversion initial value selection method based on the array lateral logging data are verified, and the comparison result is shown in table 2.
TABLE 2 inversion initial value vs. actual formation parameters
By comparing the two inversion initial value selection methods, the resistivity logging initial value selection method based on the array lateral logging data has better effect and can be popularized and applied to the field.
According to the resistivity inversion initial value selection method based on the array lateral logging information, the deviation between the initial value and the true value can be effectively reduced without supplementing other logging information, and the resistivity logging inversion accuracy and speed can be improved.
It is to be understood that the foregoing examples, while indicating the preferred embodiments of the invention, are given by way of illustration and description, and are not to be construed as limiting the scope of the invention; it should be noted that, for those skilled in the art, the above technical features can be freely combined, and several changes and modifications can be made without departing from the concept of the present invention, which all belong to the protection scope of the present invention; therefore, all equivalent changes and modifications made within the scope of the claims of the present invention should be covered by the claims of the present invention.
Claims (5)
1. A resistivity inversion initial value selection method based on array laterolog data is characterized by comprising the following steps:
step S100: forward calculation is carried out based on a radial step stratum model, resistivity amplitude difference coefficients S1, S2, S3, S4, S5 and S6 are defined, and the radial step stratum model is obtained forIn array laterolog MLR1a、MLR2a、MLR3a、MLR4aCorrecting the apparent resistivities of the four basic logging curves, performing borehole and surrounding rock correction, and establishing a chart of variation of resistivity amplitude difference coefficients along with the radius of a flushing zone;
step S200: establishing an initial value r of the radius of the flushing zone through linear regression based on a chart of the variation of the resistivity amplitude difference coefficient along with the radius of the flushing zonexoczA quantitative relationship with a resistivity magnitude difference coefficient;
step S300: MLR4 for basic well logging curve provided by array lateral well loggingaApparent resistivity of the earth layer is used as an initial value R of the formation resistivitytczThen, the initial value R of the resistivity of the flushing zone is obtained based on the pseudo-geometric factor definitionxocz。
2. The method for selecting the initial value of resistivity inversion based on the array laterolog data as claimed in claim 1, wherein in the step S100, S1, S2, S3, S4, S5 and S6 are respectively defined as follows:
S1=2·(MLR4a-MLR1a)/(MLR4a+MLR1a) (1)
S2=2·(MLR4a-MLR2a)/(MLR4a+MLR2a) (2)
S3=2·(MLR4a-MLR3a)/(MLR4a+MLR3a) (3)
S4=2·(MLR3a-MLR1a)/(MLR3a+MLR1a) (4)
S5=2·(MLR3a-MLR2a)/(MLR3a+MLR2a) (5)
S6=2·(MLR2a-MLR1a)/(MLR2a+MLR1a) (6)
wherein, MLR1a、MLR2a、MLR3a、MLR4aThe apparent resistivity of the array laterolog provided for four basic well logging curves is represented, and the unit is omega.m; s1, S2, S3, S4, S5 and S6 are dimensionless.
3. The method of claim 2, wherein the initial value r of the radius of the washzone is established in step S200xoczThe quantitative relation with the resistivity amplitude difference coefficient is as follows:
rxocz=925.14·S1+1004.04·S2-1835.46·S3-1295.69·S4-221.40·S5+70.02·S6+1.39 (7)。
4. the method for selecting initial resistivity inversion values based on array laterolog data as claimed in claim 3, wherein in step S300, the MLR4 for the base well log provided by the array laterolog is appliedaApparent resistivity of the earth layer is used as an initial value R of the formation resistivitytczThen, the initial value R of the resistivity of the flushing zone is obtained based on the pseudo-geometric factor definitionxocz:
Wherein, K1Is MLR1aA corresponding instrument constant; r is0Is the borehole radius; h is1Is MLR1aThe corresponding main current thickness; k1、r0、h1Are all known values.
5. The method of claim 1, further comprising selecting an initial value for resistivity inversion based on the array laterolog data
Step S400: setting the initial resistivity value R of the flushing beltxoczCompared to actual formation parameters.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910786024.7A CN110488367B (en) | 2019-08-23 | 2019-08-23 | Resistivity inversion initial value selection method based on array laterolog data |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910786024.7A CN110488367B (en) | 2019-08-23 | 2019-08-23 | Resistivity inversion initial value selection method based on array laterolog data |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110488367A CN110488367A (en) | 2019-11-22 |
CN110488367B true CN110488367B (en) | 2021-04-13 |
Family
ID=68553495
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910786024.7A Active CN110488367B (en) | 2019-08-23 | 2019-08-23 | Resistivity inversion initial value selection method based on array laterolog data |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110488367B (en) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111305834B (en) * | 2020-02-25 | 2022-09-23 | 北京工业大学 | Three-dimensional inversion initial model construction method based on multi-detection-mode resistivity logging |
CN112253090B (en) * | 2020-10-14 | 2022-10-18 | 中海油田服务股份有限公司 | Data parameter inversion method and device for multi-frequency electrical imaging |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
SU1188690A1 (en) * | 1984-05-14 | 1985-10-30 | Всесоюзный научно-исследовательский институт нефтепромысловой геофизики | Arrangement for well lateral logging |
RU2402047C1 (en) * | 2009-09-08 | 2010-10-20 | Владимир Алексеевич Королев | Measurement procedure at lateral logging of wells |
RU2478223C1 (en) * | 2011-09-28 | 2013-03-27 | Закрытое акционерное общество Научно-производственная фирма "ГИТАС" (ЗАО НПФ "ГИТАС") | Evaluation method of formation resistivity at investigations of wells cased with metal string |
CN102678102B (en) * | 2012-03-31 | 2015-02-25 | 中国石油大学(华东) | Array electric imaging logging based reservoir oil-water identification method and system |
CN103821506B (en) * | 2014-03-04 | 2016-04-20 | 中国石油大学(华东) | Well all resistivity of media cubical arraies imaging measurement method |
CN108729898B (en) * | 2017-04-17 | 2023-11-28 | 中国石油集团长城钻探工程有限公司 | Lateral logging instrument electrode system and lateral logging instrument |
CN109577963A (en) * | 2018-10-19 | 2019-04-05 | 中国石油天然气股份有限公司 | Device for simulating anisotropic stratum array lateral logging response and construction method thereof |
CN109885927B (en) * | 2019-02-18 | 2022-08-02 | 中国石油天然气集团有限公司 | Continuous inversion method for stratum radial resistivity |
-
2019
- 2019-08-23 CN CN201910786024.7A patent/CN110488367B/en active Active
Also Published As
Publication number | Publication date |
---|---|
CN110488367A (en) | 2019-11-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11118441B2 (en) | Optimized geosteering using real-time geological models | |
CN107783187B (en) | Method for establishing three-dimensional velocity field by combining logging velocity and seismic velocity | |
EP2810101B1 (en) | Improving efficiency of pixel-based inversion algorithms | |
CN106154322B (en) | Log curve correction method and apparatus | |
CN108345719B (en) | Post-drilling inversion method for electromagnetic wave logging data while drilling of highly deviated well/horizontal well | |
CN108005646B (en) | Stratum anisotropic resistivity extraction method based on electromagnetic wave logging while drilling data | |
CN110488367B (en) | Resistivity inversion initial value selection method based on array laterolog data | |
CN107218033B (en) | Method for identifying stratum attitude and inversion method of stratum parameters | |
US11841476B2 (en) | Methods and systems of determining parameters characterizing porous media from data gathered by a plurality of different tools | |
Zhou et al. | Interpretation of azimuthal propagation resistivity measurements: Modeling, inversion, application and discussion | |
CN117950065B (en) | Layer thickness correction method for resistivity logging data of horizontal well array | |
CN109885927B (en) | Continuous inversion method for stratum radial resistivity | |
CN108665545B (en) | Logging parameter three-dimensional geological model establishing method | |
CN111965712B (en) | Method for correcting static effect of controllable source audio magnetotelluric method | |
CN112114374B (en) | Multi-density interface inversion method of complex geologic body | |
CN112379462A (en) | Electromagnetic seismic data joint processing method and device | |
Wu et al. | Enhanced reservoir geosteering and geomapping from refined models of ultra-deep LWD resistivity inversions using machine-learning algorithms | |
CN115263287A (en) | High-precision comprehensive exploration method for geothermal resources | |
CN118011509B (en) | Formation fracture effectiveness evaluation method and system | |
CN110513105A (en) | A kind of cubical array induction logging instrument data processing system and method | |
CN115267927B (en) | Multi-boundary curtain type geosteering method based on ant colony-gradient series algorithm | |
CN113468727B (en) | Layer density modeling method based on prior structure and known point double constraint | |
CN114109349B (en) | Method for determining porosity index/saturation index of tight sandstone reservoir | |
Boesing et al. | High-Angle Formation Response Modeling for More Accurate Petrophysics–Challenges and Experiences from Multiple Cases | |
CN118011509A (en) | Formation fracture effectiveness evaluation method and system |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |