CN115898368A - Horizontal well shale oil reservoir anisotropic resistivity extraction method - Google Patents

Horizontal well shale oil reservoir anisotropic resistivity extraction method Download PDF

Info

Publication number
CN115898368A
CN115898368A CN202211407726.8A CN202211407726A CN115898368A CN 115898368 A CN115898368 A CN 115898368A CN 202211407726 A CN202211407726 A CN 202211407726A CN 115898368 A CN115898368 A CN 115898368A
Authority
CN
China
Prior art keywords
resistivity
stratum
model
array
anisotropic
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.)
Pending
Application number
CN202211407726.8A
Other languages
Chinese (zh)
Inventor
赵伟娜
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Qingdao University of Technology
Original Assignee
Qingdao University of Technology
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Qingdao University of Technology filed Critical Qingdao University of Technology
Priority to CN202211407726.8A priority Critical patent/CN115898368A/en
Publication of CN115898368A publication Critical patent/CN115898368A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention discloses a horizontal well shale oil reservoir anisotropic resistivity extraction method, which comprises the following steps: acquiring adjacent well information, N sub-arrays of horizontal well array induction logging and a synthesized apparent resistivity curve; dividing horizontal well stratum interfaces based on the short source distance subarray curve and the actually measured GR curve; establishing an explanation model and constructing a layered isotropic processing model M 1 And a homogeneous anisotropy processing model M 2 (ii) a Inverting the ith sub-array curve to obtain a formation resistivity profile, and taking the profile as an initial value of a next long-source-distance sub-array; sequentially inverting the apparent resistivity curves of the sub arrays to obtain corresponding formation resistivity profiles; for the jth stratum, model M is adopted 2 Inverting the anisotropic resistivity R of the current layer hj And R vj (ii) a And circulating to obtain the anisotropic resistivity of each stratum. The method can accurately extract the anisotropic resistivity information of the shale oil reservoir, effectively eliminate the influence of well deviation, surrounding rocks and anisotropy on logging response, and ensure the calculation precision of the oil saturation of the reservoir.

Description

Horizontal well shale oil reservoir anisotropic resistivity extraction method
Technical Field
The invention relates to the technical field of petroleum exploration and development, belongs to the field of electrical logging methods, and particularly relates to a method for extracting anisotropic resistivity of a horizontal well shale oil reservoir.
Background
The array induction logging is widely applied to oil and gas resource exploration, has remarkable characteristics for shale oil reservoirs, has small invasion abnormality and is widely applied to oil-based mud. And the separation of curves in the inclined well and the horizontal well causes the inaccuracy of saturation parameters, and the synthetic curve can not reflect the real information of the stratum under the influence of the layer thickness. However, the simulation of the array induction logging response in the horizontal well usually adopts a three-dimensional processing method, the calculation precision is low, the speed is slow, and meanwhile, the simulation in the horizontal well is seriously influenced by various factors such as instrument inclination, mud filtrate invasion and the like, and the real stratum condition cannot be accurately reflected. At present, the accurate extraction method for the resistivity of the horizontal well is less, and a new extraction method for the anisotropic resistivity of the horizontal well needs to be researched and developed urgently for effectively measuring the anisotropic resistivity information of the shale oil reservoir in the horizontal well.
Disclosure of Invention
In order to solve the technical problem, the invention discloses a horizontal well shale oil reservoir anisotropic resistivity extraction method which is based on an array induction logging original signal and aims to provide accurate stratum information for shale oil reservoir evaluation.
In order to realize the purpose, the invention adopts the following technical scheme:
a horizontal well shale oil reservoir anisotropic resistivity extraction method comprises the following steps:
s1. obtaining adjacent well stratum layer thickness, stratum horizontal resistivity and natural GR curve information, and actually measuring N sub-arrays of horizontal well array induction logging and a synthesized apparent resistivity curve, wherein the horizontal well array induction logging sub-array comprises a short source distance sub-array and a long source distance sub-array;
s2, dividing horizontal well stratum interfaces based on the short source distance subarray curve and an actually measured horizontal well GR curve;
s3. establishes an interpretation model and a layered isotropic processing model M 1 And a homogeneous anisotropy processing model M 2
s4. inverses the ith sub-array curve to obtain the formation resistivity profile Rts i And taking the result as the initial value of the next long-source-distance sub array;
s5. sequentially inverts apparent resistivity curves of the sub-arrays to respectively obtain corresponding formation resistivity profiles Rts i ,i=1,…,N;
s6. model M for layered isotropic processing 1 And (5) combining the inversion result of the step s5 with the jth stratum, and adopting a homogeneous anisotropy processing model M 2 Inverting the anisotropic resistivity R of the current layer hj And R vj
s7. sequentially loop through step s6 to process the layered isotropic processing model M 1 And obtaining the anisotropic resistivity of each stratum.
Optionally, step s3 comprises the steps of:
s3.1, aiming at the electrical characteristics of a shale oil reservoir and the drilling environment, considering the influences of a well hole, stratum folds, layer thickness and well deviation, and establishing an actual stratum model;
s3.2, neglecting the influence of stratum folds, and after borehole correction is carried out on the N sub-array signals in the step s1, establishing a corresponding interpretation model;
and S3.3, respectively considering the layer thickness and the anisotropic influence, and equating the stratum to a stratum with the layer thickness influence and without the anisotropic influence and a stratum with the anisotropic influence and without the layer thickness influence, and further sequentially constructing a layered isotropic processing model M 1 And a homogeneous anisotropy processing model M 2
S3.4 layered Isotropic treatment model M 1 Including the parameter Rts i Theta and H j Wherein Rts is i For equivalent isotropic resistivity, H j The thickness of the stratum is shown, theta is well deviation, i represents the ith sub-array, j represents the jth stratum, and the layered isotropic processing model M1 is used for eliminating the influence of the layer thickness;
s3.5. homogeneous Anisotropic Process model M 2 Comprising a parameter R hj Theta and R vj Wherein R is hj For horizontal resistivity, R vj Homogeneous anisotropy processing model M for vertical resistivity, theta for well deviation, j for jth formation 2 For obtaining anisotropy information.
Optionally, step s4 comprises the steps of:
s4.1, providing an inversion initial value by using the reference resistivity of the adjacent well, the actual measurement subarray curve and the synthetic curve of the horizontal well;
s4.2. For the first subarray curve, combining the inversion initial value provided in step s4.1 and the layered isotropic processing model M constructed in step s3 1 And obtaining formation resistivity profile Rts by inversion i And taking the result as the initial value of the next long-source-distance sub-array;
s4.3, for the ith sub-array curve, carrying out inversion on the previous sub-array to obtain a formation resistivity profile Rts i-1 As an initial value of inversion, obtaining formation resistivity profile Rts by inversion i And the result is used as the initial value of the next long source distance sub-array.
Optionally, step s6 comprises the steps of:
s6.1. homogeneous anisotropy processing model M established in step s3 2 Including a parameter R h And R v Wherein R is h For horizontal resistivity, R v For vertical resistivity, the model ignores layer thickness effects;
s6.2. model M for lamellar Isotropic processing 1 Taking the N resistivity profiles of the current layer obtained in the step s5 as inversion input information of the jth stratum, and combining the homogeneous anisotropic processing model M constructed in the step s3 2 And (5) inverting the interpretation result of the step s5 to obtain the anisotropic resistivity R of the current layer hj And R vj
S6.3. Model M for lamellar Isotropic treatment 1 And (4) performing layer-by-layer inversion on all the layers, and repeating the process to obtain the anisotropic resistivity of each stratum.
The method has the beneficial effects that the method for extracting the anisotropic resistivity of the shale oil reservoir of the horizontal well based on the array induction logging original signal is provided, the information of adjacent wells is fully utilized, the position of a layer interface is determined, and an initial value is provided for subsequent processing; two equivalent processing strategies are adopted for a complex three-dimensional stratum, and a plurality of groups of initial values are provided for each equivalent model, so that the speed and the precision of the inversion process are greatly improved; in the later explanation evaluation, the oil saturation, the mobile oil saturation and the like of the reservoir can be accurately calculated through the anisotropic inversion result of the resistivity of each layer of stratum, and reliable parameters are provided for the reservoir evaluation.
The method provided by the invention can effectively eliminate the influence of well deviation, surrounding rocks and anisotropy on logging response, and effectively extract the anisotropic resistivity information of the horizontal well shale oil reservoir.
Drawings
FIG. 1 is a flow chart of a horizontal well shale oil reservoir anisotropic resistivity extraction method based on an array induction logging original signal in the invention;
FIG. 2 is a schematic diagram of an equivalent interpretation model of the present invention;
FIG. 3 is a schematic view of a layered isotropic processing model according to the present invention;
FIG. 4 is a schematic view of a homogeneous anisotropic process model according to the present invention;
FIG. 5 is a 65 and 80 triple layer earth model of the present invention wherein (a) is a 65 triple layer earth model and (b) is an 80 triple layer earth model;
FIG. 6 is the raw response of each subarray of the induction log of the model array of FIG. 5, wherein (a) is the raw response of each subarray of the induction log of the 65-degree formation model array, and (b) is the raw response of each subarray of the induction log of the 80-degree formation model array;
fig. 7 is an inversion result of the lamellar isotropic processing model corresponding to each sub-array of the present invention, in which (a) is an inversion result of the 65 ° lamellar isotropic processing model, and (b) is an inversion result of the 80 ° lamellar isotropic processing model;
FIG. 8 is the inversion result of the homogeneous anisotropy processing model of the invention, wherein (a) is the inversion result of the 65 ° homogeneous anisotropy processing model, and (b) is the inversion result of the 80 ° homogeneous anisotropy processing model.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
A horizontal well shale oil reservoir anisotropic resistivity extraction method is shown in figure 1 and comprises the following steps:
s1. obtains the stratum layer thickness, the stratum horizontal resistivity and the natural GR curve information of an adjacent well, and carries out induction logging on N sub-arrays (SGM) of an actually measured horizontal well array 0 、SGM 1 ……SGM N-1 ) And a synthesized apparent resistivity curve, wherein the horizontal well array induction logging subarray comprises a short source distance subarray and a long source distance subarray;
s2, dividing horizontal well stratum interfaces based on the short source distance subarray curve and an actually measured horizontal well GR curve;
s3. establishes an explanation model and a layered isotropic processing model M 1 And a homogeneous anisotropy processing model M 2
Step s3 specifically includes the following steps:
s3.1, aiming at the electrical characteristics of the shale oil reservoir and the drilling environment, considering the influences of a well hole, stratum folds, layer thickness and well deviation, and establishing an actual stratum model;
s3.2. As shown in FIG. 2, because the detection range is limited, neglecting the influence of stratum folds, after borehole correction is carried out on the N sub-array signals in the step s1, establishing a corresponding interpretation model;
and S3.3, respectively considering the layer thickness and the anisotropic influence, and equating the stratum to a stratum with the layer thickness influence and without the anisotropic influence and a stratum with the anisotropic influence and without the layer thickness influence, and further sequentially constructing a layered isotropic processing model M 1 And a homogeneous anisotropy processing model M 2
S3.4. As shown in FIG. 3, layered Isotropic Process model M 1 Including the parameter Rts i Theta and H j Wherein Rts is i For equivalent isotropic resistivity, H j The stratum thickness is shown, theta is well deviation, i represents the ith sub-array, j represents the jth stratum, and the layered isotropic processing model M1 is used for eliminating the layer thickness influence;
s3.5. homogeneous Anisotropic Process model M, as shown in FIG. 4 2 Comprising a parameter R hj Theta and R vj Wherein R is hj Is horizontal resistivity, R vj Homogeneous anisotropy processing model M for vertical resistivity, theta for well deviation, j for jth formation 2 For obtaining anisotropy information.
s4. inverses the ith sub-array curve to obtain the formation resistivity profile Rts i And using the result as the next long source spacerAn initial value of the array;
step s4 specifically includes the following steps:
s4.1, taking the stratum model shown in fig. 5 as an example, providing an inversion initial value for data processing by using the reference resistivity of the adjacent well, the actual measurement subarray curve and the synthetic curve of the horizontal well;
s4.2. For the first subarray curve, as shown in fig. 6, combining the initial inversion value provided in step s4.1 with the layered isotropic processing model M constructed in step s3 1 And obtaining formation resistivity profile Rts by inversion i And taking the result as the initial value of the next long-source-distance sub-array;
s4.3, for the ith sub-array curve, performing inversion on the previous sub-array to obtain a formation resistivity profile Rts i-1 As an initial value of inversion, obtaining formation resistivity profile Rts by inversion i And the result is used as the initial value of the next long source range sub-array.
s5. as shown in FIG. 7, sequentially inverting apparent resistivity curves of each sub-array to obtain corresponding formation resistivity profile Rts i ,i=1,…,N;
s6. model M for layered isotropic processing 1 And (5) processing the model M by adopting homogeneous anisotropy in combination with the inversion result of the step s5 for the jth stratum 2 Inverting the anisotropic resistivity R of the current layer hj And R vj
Step s6 specifically includes the following steps:
s6.1. homogeneous anisotropy processing model M established in step s3 2 Including the parameter R h And R v Wherein R is h For horizontal resistivity, R v For vertical resistivity, the model ignores layer thickness effects;
s6.2. As shown in FIG. 8, model M is processed isotropically to a layer 1 In the j-th stratum, taking the N resistivity profiles of the current stratum obtained in the step s5 as inversion input information, and combining the homogeneous anisotropy processing model M constructed in the step s3 2 And (5) inverting the interpretation result of the step s5 to obtain the anisotropic resistivity R of the current layer hj And R vj
S6.3. LaminatingIsotropic processing model M 1 And (4) performing layer-by-layer inversion on all the layers, and repeating the process to obtain the anisotropic resistivity of each stratum.
s7. sequentially loop through step s6 to process the layered isotropic processing model M 1 And obtaining the anisotropic resistivity of each stratum.
It is to be understood that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make modifications, alterations, additions or substitutions within the spirit and scope of the present invention.

Claims (4)

1. A horizontal well shale oil reservoir anisotropic resistivity extraction method is characterized by comprising the following steps:
s1. obtaining the stratum layer thickness, the stratum horizontal resistivity and the natural GR curve information of an adjacent well, and actually measuring N sub-arrays of the horizontal well array induction logging and a synthesized apparent resistivity curve, wherein the horizontal well array induction logging sub-array comprises a short source distance sub-array and a long source distance sub-array;
s2, dividing horizontal well stratum interfaces based on the short source distance subarray curve and an actually measured horizontal well GR curve;
s3. establishes an explanation model and a layered isotropic processing model M 1 And a homogeneous anisotropy processing model M 2
s4. inverses the ith sub-array curve to obtain the formation resistivity profile Rts i And taking the result as an initial value of a next long source distance sub array;
s5. sequentially inverts apparent resistivity curves of the sub-arrays to respectively obtain corresponding formation resistivity profiles Rts i ,i=1,…,N;
s6. model M for layered isotropic processing 1 And (5) combining the inversion result of the step s5 with the jth stratum, and adopting a homogeneous anisotropy processing model M 2 Inverting the anisotropic resistivity R of the current layer hj And R vj
s7. sequentially loops step s6 to process the layered isotropic spotPhysical model M 1 And obtaining the anisotropic resistivity of each stratum.
2. The method for extracting the anisotropic resistivity of the horizontal well shale oil reservoir according to claim 1, wherein the step s3 comprises the following steps:
s3.1, aiming at the electrical characteristics of a shale oil reservoir and the drilling environment, considering the influences of a well hole, stratum folds, layer thickness and well deviation, and establishing an actual stratum model;
s3.2, neglecting stratum fold influence, and establishing a corresponding interpretation model after borehole correction is carried out on the N sub-array signals in the step s 1;
and S3.3, respectively considering the layer thickness and the anisotropic influence, and equating the stratum to a stratum with the layer thickness influence and without the anisotropic influence and a stratum with the anisotropic influence and without the layer thickness influence, and further sequentially constructing a layered isotropic processing model M 1 And a homogeneous anisotropy processing model M 2
S3.4 layered Isotropic treatment model M 1 Including the parameter Rts i Theta and H j Wherein Rts is i For equivalent isotropic resistivity, H j The stratum thickness is shown, theta is well deviation, i represents the ith sub-array, j represents the jth stratum, and the layered isotropic processing model M1 is used for eliminating the layer thickness influence;
s3.5. homogeneous Anisotropic Process model M 2 Comprising a parameter R hj Theta and R vj Wherein R is hj For horizontal resistivity, R vj Homogeneous anisotropy processing model M for vertical resistivity, theta for well deviation, j for jth formation 2 For obtaining anisotropy information.
3. The method for extracting the anisotropic resistivity of the horizontal well shale oil reservoir according to claim 1, wherein the step s4 comprises the following steps:
s4.1, providing an inversion initial value by utilizing the reference resistivity of the adjacent well, the actually measured subarray curve of the horizontal well and the synthetic curve;
s4.2. For the first subarray curve, combine the stepss4.1 and the layered isotropic processing model M constructed in step s3 1 And obtaining formation resistivity profile Rts by inversion i And taking the result as the initial value of the next long-source-distance sub-array;
s4.3, for the ith sub-array curve, performing inversion on the previous sub-array to obtain a formation resistivity profile Rts i-1 As an initial value of inversion, obtaining formation resistivity profile Rts by inversion i And the result is used as the initial value of the next long source range sub-array.
4. The method for extracting the anisotropic resistivity of the horizontal well shale oil reservoir according to claim 1, wherein the step s6 comprises the following steps:
s6.1. homogeneous anisotropy processing model M established in step s3 2 Including a parameter R h And R v Wherein R is h For horizontal resistivity, R v For vertical resistivity, the model ignores layer thickness effects;
s6.2. model M for layer isotropy 1 Taking the N resistivity profiles of the current layer obtained in the step s5 as inversion input information of the jth stratum, and combining the homogeneous anisotropic processing model M constructed in the step s3 2 And (5) inverting the interpretation result of the step s5 to obtain the anisotropic resistivity R of the current layer hj And R vj
S6.3. Model M is treated isotropically in layers 1 And (4) performing layer-by-layer inversion on all the layers, and repeating the process to obtain the anisotropic resistivity of each stratum.
CN202211407726.8A 2022-11-10 2022-11-10 Horizontal well shale oil reservoir anisotropic resistivity extraction method Pending CN115898368A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211407726.8A CN115898368A (en) 2022-11-10 2022-11-10 Horizontal well shale oil reservoir anisotropic resistivity extraction method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211407726.8A CN115898368A (en) 2022-11-10 2022-11-10 Horizontal well shale oil reservoir anisotropic resistivity extraction method

Publications (1)

Publication Number Publication Date
CN115898368A true CN115898368A (en) 2023-04-04

Family

ID=86483488

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211407726.8A Pending CN115898368A (en) 2022-11-10 2022-11-10 Horizontal well shale oil reservoir anisotropic resistivity extraction method

Country Status (1)

Country Link
CN (1) CN115898368A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118504294A (en) * 2024-07-18 2024-08-16 中国石油大学(华东) Method for rapidly calculating response of all-anisotropic medium of coaxial coil instrument

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118504294A (en) * 2024-07-18 2024-08-16 中国石油大学(华东) Method for rapidly calculating response of all-anisotropic medium of coaxial coil instrument
CN118504294B (en) * 2024-07-18 2024-10-18 中国石油大学(华东) Method for rapidly calculating response of all-anisotropic medium of coaxial coil instrument

Similar Documents

Publication Publication Date Title
CN110795878B (en) Tunnel water inflow prediction method
CN107783187B (en) Method for establishing three-dimensional velocity field by combining logging velocity and seismic velocity
CN102373923B (en) Reservoir stratum identification method
CN103233727B (en) Inversion method of stratum shear wave velocity radial sections
CN108005646B (en) Stratum anisotropic resistivity extraction method based on electromagnetic wave logging while drilling data
CN107329171A (en) Depth domain reservoir stratum seismic inversion method and device
CN104749617A (en) Multi-scale fractured reservoir forward model establishing method
CN102606138B (en) Method for correcting dielectric constant of electromagnetic wave resistivity logger while drilling by means of phase-amplitude method
CN107784159B (en) Method for determining anisotropic coefficient of reservoir resistivity
CN110596757B (en) Method for correcting longitudinal wave and transverse wave velocities of shale formation
CN105259581B (en) Seismic data time-depth conversion method
CN110501757B (en) Point-line-surface three-dimensional point-by-point well depth design method
CN105114063A (en) Inclined-well dual laterolog curve correction method
CN115898368A (en) Horizontal well shale oil reservoir anisotropic resistivity extraction method
CN102877828A (en) CT (Computed Tomography) imaging method of three-dimensional multi-well combined well land
CN110532507B (en) Method for improving drilling rate of compact oil reservoir of horizontal well
CN103790579B (en) Method and device for determining distance between drill bit and stratum interface in while-drilling geological steering
CN110488367B (en) Resistivity inversion initial value selection method based on array laterolog data
CN110500089B (en) Shale gas horizontal well bedding crack evaluation method based on array induction logging data
CN106908838A (en) The method for building target area stratum inelastic attenuation quality factor three-dimensional model
CN109991663B (en) Work area seismic velocity body correction method and device
CN115434696A (en) Lithology identification database construction method and lithology identification method
CN106646668B (en) A kind of method for building up of radar logging standard well model
CN110687591B (en) Method for determining physical property parameters of coal bed and surrounding rock based on waveform matching of prior data
CN109061737B (en) Reservoir prediction method and device based on synthetic seismic record

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