CN105607147B - A kind of method and system of inverting shale gas reservoir resistivity - Google Patents
A kind of method and system of inverting shale gas reservoir resistivity Download PDFInfo
- Publication number
- CN105607147B CN105607147B CN201510967702.1A CN201510967702A CN105607147B CN 105607147 B CN105607147 B CN 105607147B CN 201510967702 A CN201510967702 A CN 201510967702A CN 105607147 B CN105607147 B CN 105607147B
- Authority
- CN
- China
- Prior art keywords
- apparent resistivity
- inverting
- target area
- resistivity
- shale gas
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a kind of method and system of inverting shale gas reservoir resistivity.This method includes：The apparent resistivity of the target area is obtained using wide area electromagnetic method；Acquire the seismic data of target area；Constraint inverting is carried out to the apparent resistivity according to the seismic data, obtains the changes in distribution rule of the shale gas reservoir resistivity of target area.The invention avoids multi-solution caused by inverting is carried out in the prior art with earthquake or electromagnetism single method, the precision and accuracy of inverting shale gas reservoir resistivity are improved.
Description
Technical field
A kind of method in field more particularly to inverting shale gas reservoir resistivity is reconnoitred the present invention relates to geophysics and is
System.
Background technology
Shale gas is preservation in rich organic matter mud shale and its interlayer, with absorption and free state for main existing way
Unconventional gas, ingredient belongs to one kind based on methane, with " coal bed gas ", " tight gas ".The formation and enrichment of shale gas
There is itself uniqueness, is often distributed in that thickness in basin is larger, in widespread shale hydrocarbon source rock formation.
North America region mainly uses earthquake and drilling technology to carry out shale gas exploration, and China carries out marine facies shale gas and surveys at present
The method of spy is divided into seismic prospecting and electromagnetic prospecting substantially, mainly based on earthquake or the single inversion interpretation of electromagnetism, reliability
Poor, being faced with inversion result multi-solution etc., some are challenged, and cause shale gas exploitation inconvenient.
Invention content
The technical problem to be solved by the present invention is to：The problem of existing shale gas reservoirs exploration method poor reliability.
In order to solve the above technical problems, one aspect of the present invention proposes a kind of method of inverting shale gas reservoir resistivity,
This method includes：
The apparent resistivity of the target area is obtained using wide area electromagnetic method；
Acquire the seismic data of target area；
Constraint inverting is carried out to the apparent resistivity according to the seismic data, obtains the shale gas reservoir resistance of target area
The changes in distribution rule of rate.
Optionally, described that constraint inverting is carried out to the apparent resistivity according to the seismic data, including：
According to the buried depth and distribution of shale gas reservoir in the seismic data acquisition target area；
Constraint inverting is carried out to the apparent resistivity according to the buried depth of shale gas reservoir in the target area and distribution.
Optionally, the buried depth and distribution according to shale gas reservoir in the seismic data acquisition target area, packet
It includes：
The pretreatment of denoising is carried out to the seismic data；
Prestack static correction is carried out to pretreated seismic data；
Prestack time and depth migration are carried out to the data after prestack static correction；
The buried depth of shale gas reservoir and distribution from acquisition target area in the seismic data after prestack time and depth migration
Range.
Optionally, it is described according to the buried depth of shale gas reservoir in the target area and distribution to the apparent resistivity into
Row constraint inverting, including：
The polarization mode of the apparent resistivity is identified；
The apparent resistivity is filtered by electromagnetic equation；
The recovery impedance phase of corresponding polarization mode is carried out to the apparent resistivity；
The one-dimensional continuous media inverting of corresponding polarization mode is carried out to the apparent resistivity；
It is based on one-dimensional continuous media inverting as a result, carrying out two-dimentional continuous media inversion imaging to the apparent resistivity；
In two-dimentional continuous media inversion result in conjunction with the target area in shale gas reservoir buried depth and distribution pair
Inverting section is finely adjusted.
Optionally, before the recovery impedance phase for carrying out corresponding polarization mode to the apparent resistivity, further include：
The data interpolating of corresponding polarization mode is carried out to the apparent resistivity.
Optionally, the polarization mode to the apparent resistivity is identified, including：
The polarization mode of the TE and/or TM of the apparent resistivity are identified.
Optionally, the apparent resistivity that the target area is obtained using wide area electromagnetic method, including：
The apparent resistance of the target area is obtained using horizontal current source wide area electromagnetic method or perpendicular magnetic field source wide area electromagnetic method
Rate.
On the other hand, the present invention proposes a kind of system of inverting shale gas reservoir resistivity, which includes：
Apparent resistivity acquiring unit, the apparent resistivity for obtaining the target area using wide area electromagnetic method；
Seismic data acquisition unit, the seismic data for acquiring target area；
Shale gas reservoir resistivity acquiring unit, it is anti-for to the apparent resistivity constrain according to the seismic data
It drills, obtains the changes in distribution rule of the shale gas reservoir resistivity of target area.
Optionally, the shale gas reservoir resistivity acquiring unit is further used for according to the seismic data acquisition mesh
Mark the buried depth and distribution of shale gas reservoir in area；According to the buried depth and distribution pair of shale gas reservoir in the target area
The apparent resistivity carries out constraint inverting.
The method and system of inverting shale gas reservoir resistivity provided by the invention are using multidimensional delineation shale gas reservoir
Method avoids and carries out multi-solution caused by inverting in the prior art with earthquake or electromagnetism single method, improves inverting page
The precision and accuracy of rock gas reservoir resistivity.
Description of the drawings
The features and advantages of the present invention can be more clearly understood by reference to attached drawing, attached drawing is schematically without that should manage
Solution is carries out any restrictions to the present invention, in the accompanying drawings：
Fig. 1 shows the schematic diagram of the method for the inverting shale gas reservoir resistivity of one embodiment of the invention；
Fig. 2 shows the flow diagrams of the method for the inverting shale gas reservoir resistivity of another embodiment of the present invention；
Fig. 3 shows the buried depth of the acquisition shale gas reservoir of one embodiment of the invention and the flow signal of distribution
Figure；
Fig. 4 shows the flow diagram that constraint inverting is carried out to apparent resistivity of one embodiment of the invention；
Fig. 5 shows the structural schematic diagram of the system of the inverting shale gas reservoir resistivity of one embodiment of the invention；
Fig. 6 shows the seismic profile schematic diagram of the target area of one embodiment of the invention；
Fig. 7 shows the diagrammatic cross-section of the inverting shale gas reservoir resistivity of one embodiment of the invention.
Specific implementation mode
Below in conjunction with attached drawing, embodiments of the present invention is described in detail.
Fig. 1 shows the schematic diagram of the method for the inverting shale gas reservoir resistivity of one embodiment of the invention.Such as Fig. 1 institutes
Show, this method includes：
S1：The apparent resistivity of the target area is obtained using wide area electromagnetic method；
S2：Acquire the seismic data of target area；
S3：Constraint inverting is carried out to the apparent resistivity according to the seismic data, obtains the shale gas reservoir of target area
The changes in distribution rule of resistivity.
The method of the inverting shale gas reservoir resistivity of the present embodiment is avoided using the method for multidimensional delineation shale gas reservoir
Multi-solution caused by inverting is carried out with earthquake or electromagnetism single method in the prior art, improves inverting shale gas reservoir electricity
The precision and accuracy of resistance rate.
Fig. 2 shows the flow diagrams of the method for the inverting shale gas reservoir resistivity of another embodiment of the present invention.
As shown in Fig. 2, this method includes：
S21：The apparent resistivity of the target area is obtained using wide area electromagnetic method；
S22：Acquire the seismic data of target area；
S23：According to the buried depth and distribution of shale gas reservoir in the seismic data acquisition target area；
S24：The apparent resistivity is constrained according to the buried depth of shale gas reservoir in the target area and distribution
Inverting.
Fig. 3 shows the buried depth of the acquisition shale gas reservoir of one embodiment of the invention and the flow signal of distribution
Figure.As shown in figure 3, according to the buried depth and distribution of shale gas reservoir in the seismic data acquisition target area, including：
S31：The pretreatment of denoising is carried out to the seismic data；
S32：Prestack static correction is carried out to pretreated seismic data；
S33：Prestack time and depth migration are carried out to the data after prestack static correction；
S34：From in the seismic data after prestack time and depth migration obtain target area in shale gas reservoir buried depth and
Distribution.
Fig. 6 shows that the seismic profile schematic diagram of the target area of one embodiment of the invention, this figure are according to target area
Seismic data obtains the section that target area stratum includes the buried depth of shale gas reservoir and the distribution in target area.
Fig. 4 shows the flow diagram that constraint inverting is carried out to apparent resistivity of one embodiment of the invention.Such as Fig. 4 institutes
Show, it is described that constraint inverting is carried out to the apparent resistivity according to the buried depth of shale gas reservoir in the target area and distribution,
Including：
S41：The polarization mode of the apparent resistivity is identified；
S42：The apparent resistivity is filtered by electromagnetic equation；
S43：The recovery impedance phase of corresponding polarization mode is carried out to the apparent resistivity；
S44：The one-dimensional continuous media inverting of corresponding polarization mode is carried out to the apparent resistivity；
S45：It is based on one-dimensional continuous media inverting as a result, to the apparent resistivity carry out two-dimentional continuous media inverting at
Picture；
S46：The shale gas reservoir interface correlation letter obtained from seismic profile is added in two-dimentional continuous media inversion result
Breath, is finely adjusted inverting section, and error of fitting is made to reach minimum.
In practical applications, it is loaded into target area in integral system GME_3DI is explained in weight magnetoelectricity 3-d inversion imaging
Coordinate and wide area apparent resistivity information.When being filtered by electromagnetic equation to apparent resistivity, the preferable situation of the quality of data
Selection retains more information compared with small parameter；
In a kind of optional embodiment, in the recovery impedance for carrying out corresponding polarization mode to the apparent resistivity
Before phase, further include：
The data interpolating of corresponding polarization mode is carried out to the apparent resistivity.
In practical applications, the data interpolating for carrying out corresponding polarization mode, mainly for a small amount of frequency data of a small amount of measuring point
Second-rate situation.
The polarization mode of the apparent resistivity is identified, including：
The polarization mode of the TE and/or TM of the apparent resistivity are identified.
Further, the apparent resistivity that the target area is obtained using wide area electromagnetic method, including：
The apparent resistance of the target area is obtained using horizontal current source wide area electromagnetic method or perpendicular magnetic field source wide area electromagnetic method
Rate.
Fig. 5 shows the structural schematic diagram of the system of the inverting shale gas reservoir resistivity of one embodiment of the invention.Such as
Shown in Fig. 5, which includes：
Apparent resistivity acquiring unit 51, the apparent resistivity for obtaining the target area using wide area electromagnetic method；
Seismic data acquisition unit 52, the seismic data for acquiring target area；
Shale gas reservoir resistivity acquiring unit 53, for being constrained the apparent resistivity according to the seismic data
Inverting obtains the changes in distribution rule of the shale gas reservoir resistivity of target area.
In a kind of optional embodiment, the shale gas reservoir resistivity acquiring unit is further used for according to institute
State the buried depth and distribution of shale gas reservoir in seismic data acquisition target area；According to shale gas reservoir in the target area
Buried depth and distribution carry out constraint inverting to the apparent resistivity.
The system of inverting shale gas reservoir resistivity described in the present embodiment can be used for executing above method embodiment,
Principle is similar with technique effect, and details are not described herein again.
In a specific application, Fig. 7 shows the inverting shale gas reservoir resistivity of one embodiment of the invention
Diagrammatic cross-section.The target area stratum that seismic profile is obtained includes the buried depth of shale gas reservoir and the distribution model in target area
It encloses as underground given data, constraint inverting is carried out to the wide area electromagnetic method data of target area, includes to obtain target area stratum
The response accuracy of the changes in distribution rule of shale gas reservoir resistivity, main purpose layer depth and thickness is high.
The method and system of inverting shale gas reservoir resistivity provided by the invention are using multidimensional delineation shale gas reservoir
Method avoids and carries out multi-solution caused by inverting in the prior art with earthquake or electromagnetism single method, improves inverting page
The precision and accuracy of rock gas reservoir resistivity.
Although the embodiments of the invention are described in conjunction with the attached drawings, but those skilled in the art can not depart from this hair
Various modifications and variations are made in the case of bright spirit and scope, such modifications and variations are each fallen within by appended claims
Within limited range.
Claims (3)
1. a kind of method of inverting shale gas reservoir resistivity, which is characterized in that including：
The apparent resistivity of target area is obtained using wide area electromagnetic method；
Acquire the seismic data of target area；
Constraint inverting is carried out to the apparent resistivity according to the seismic data, obtains the shale gas reservoir resistivity of target area
Changes in distribution rule；
It is described that constraint inverting is carried out to the apparent resistivity according to the seismic data, including：
According to the buried depth and distribution of shale gas reservoir in the seismic data acquisition target area；
Constraint inverting is carried out to the apparent resistivity according to the buried depth of shale gas reservoir in the target area and distribution；
It is described that constraint inverting is carried out to the apparent resistivity according to the buried depth of shale gas reservoir in the target area and distribution,
Including：
The polarization mode of the apparent resistivity is identified；
Apparent resistivity after identification is filtered by electromagnetic equation；
The recovery impedance phase of corresponding polarization mode is carried out to filtered apparent resistivity；
Apparent resistivity to restoring impedance phase carries out the one-dimensional continuous media inverting of corresponding polarization mode；
It is based on one-dimensional continuous media inverting as a result, carrying out two-dimentional continuous media to the apparent resistivity after one-dimensional continuous media inverting
Inversion imaging；
In two-dimentional continuous media inversion result in conjunction with the target area in shale gas reservoir buried depth and distribution to inverting
Section is finely adjusted；
The polarization mode to the apparent resistivity is identified, including：
The polarization mode of the TE and/or TM of the apparent resistivity are identified；
The apparent resistivity that target area is obtained using wide area electromagnetic method, including：
The apparent resistivity of the target area is obtained using horizontal current source wide area electromagnetic method or perpendicular magnetic field source wide area electromagnetic method.
2. the method for inverting shale gas reservoir resistivity according to claim 1, which is characterized in that regarded to described described
Before resistivity carries out the recovery impedance phase of corresponding polarization mode, further include：
The data interpolating of corresponding polarization mode is carried out to the apparent resistivity.
3. a kind of system of inverting shale gas reservoir resistivity, which is characterized in that including：
Apparent resistivity acquiring unit, the apparent resistivity for obtaining target area using wide area electromagnetic method；
Seismic data acquisition unit, the seismic data for acquiring target area；
Shale gas reservoir resistivity acquiring unit, for carrying out constraint inverting to the apparent resistivity according to the seismic data,
Obtain the changes in distribution rule of the shale gas reservoir resistivity of target area；
The shale gas reservoir resistivity acquiring unit is further used for：
According to the buried depth and distribution of shale gas reservoir in the seismic data acquisition target area；
Constraint inverting is carried out to the apparent resistivity according to the buried depth of shale gas reservoir in the target area and distribution；
It is described that constraint inverting is carried out to the apparent resistivity according to the buried depth of shale gas reservoir in the target area and distribution,
Including：
The polarization mode of the apparent resistivity is identified；
Apparent resistivity after identification is filtered by electromagnetic equation；
The recovery impedance phase of corresponding polarization mode is carried out to filtered apparent resistivity；
Apparent resistivity to restoring impedance phase carries out the one-dimensional continuous media inverting of corresponding polarization mode；
It is based on one-dimensional continuous media inverting as a result, carrying out two-dimentional continuous media to the apparent resistivity after one-dimensional continuous media inverting
Inversion imaging；
In two-dimentional continuous media inversion result in conjunction with the target area in shale gas reservoir buried depth and distribution to inverting
Section is finely adjusted；
The polarization mode to the apparent resistivity is identified, including：
The polarization mode of the TE and/or TM of the apparent resistivity are identified；
The apparent resistivity that target area is obtained using wide area electromagnetic method, including：
The apparent resistivity of the target area is obtained using horizontal current source wide area electromagnetic method or perpendicular magnetic field source wide area electromagnetic method.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510967702.1A CN105607147B (en) | 2015-12-21 | 2015-12-21 | A kind of method and system of inverting shale gas reservoir resistivity |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510967702.1A CN105607147B (en) | 2015-12-21 | 2015-12-21 | A kind of method and system of inverting shale gas reservoir resistivity |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN105607147A CN105607147A (en) | 2016-05-25 |
| CN105607147B true CN105607147B (en) | 2018-07-13 |
Family
ID=55987222
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510967702.1A Active CN105607147B (en) | 2015-12-21 | 2015-12-21 | A kind of method and system of inverting shale gas reservoir resistivity |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105607147B (en) |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106873041B (en) * | 2017-04-19 | 2019-05-31 | 中南大学 | A method of apparent resistivity is obtained by any level electric field component |
| CN108663414A (en) * | 2018-03-22 | 2018-10-16 | 长江大学 | A kind of method and system of determining fracturing fluid saturation degree range |
| CN110824567B (en) * | 2019-09-30 | 2020-12-22 | 中南大学 | Device for quickly positioning leakage point of ditch and control method thereof |
| CN111580181B (en) * | 2020-04-22 | 2021-07-20 | 中国矿业大学(北京) | Water guide collapse column identification method based on multi-field multi-feature information fusion |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103995301A (en) * | 2014-05-07 | 2014-08-20 | 中国石油天然气集团公司 | Method and device for evaluating total organic carbon content in shale gas reservoir |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2007068259A1 (en) * | 2005-12-15 | 2007-06-21 | Zakrytoe Aktsionernoe Obschestvo 'emmet' | Method and device for marine electrical exploration of oil and gas fields |
| US9322910B2 (en) * | 2011-07-15 | 2016-04-26 | Technoimaging, Llc | Method of real time subsurface imaging using electromagnetic data acquired from moving platforms |
| CN105116452A (en) * | 2015-08-24 | 2015-12-02 | 中国石油天然气集团公司 | Method and device of determining resistivity and polarizability of geological abnormal body |
-
2015
- 2015-12-21 CN CN201510967702.1A patent/CN105607147B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103995301A (en) * | 2014-05-07 | 2014-08-20 | 中国石油天然气集团公司 | Method and device for evaluating total organic carbon content in shale gas reservoir |
Non-Patent Citations (1)
| Title |
|---|
| 广域电磁法在湘西北页岩气探测中的应用;何继善等;《石油地球物理勘探》;20141031;第49卷(第5期);第1006、1009-1010页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN105607147A (en) | 2016-05-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Daley et al. | Time-lapse crosswell seismic and VSP monitoring of injected CO 2 in a brine aquifer | |
| Balch et al. | Vertical seismic profiling: technique, applications, and case histories | |
| WO2014143166A1 (en) | Methods of characterizing earth formations using physiochemical model | |
| Li | 4D seismic monitoring of CO 2 flood in a thin fractured carbonate reservoir | |
| Bergmann et al. | Combination of seismic reflection and constrained resistivity inversion with an application to 4D imaging of the CO 2 storage site, Ketzin, Germany | |
| CN105607147B (en) | A kind of method and system of inverting shale gas reservoir resistivity | |
| Jia et al. | Gas hydrate saturation and distribution in the Kumano Forearc Basin of the Nankai Trough | |
| Kim et al. | Lithofacies modeling by multipoint statistics and economic evaluation by NPV volume for the early Cretaceous Wabiskaw Member in Athabasca oilsands area, Canada | |
| Ponce et al. | Gravity and magnetic expression of the San Leandro gabbro with implications for the geometry and evolution of the Hayward fault zone, northern California | |
| Eisner et al. | Effective VTI anisotropy for consistent monitoring of microseismic events | |
| Unsworth et al. | Hydrogeologic assessment of the Amchitka Island nuclear test site (Alaska) with magnetotellurics | |
| Galli et al. | Resistivity modeling of array laterolog tools: An application in an offshore Norway clastic reservoir | |
| Saibi et al. | Audio‐Magnetotelluric Survey for Groundwater Investigation in the Al‐Jaww Plain in Eastern Abu Dhabi, Al‐Ain, United Arab Emirates | |
| Buddo et al. | Joint interpretation of seismic and TEM data from the Kovykta gas-condensate field, East Siberia | |
| Liberty et al. | Regolith and host rock influences on CO2 leakage: Active source seismic profiling across the Little Grand Wash fault, Utah | |
| Tinker et al. | Integrated reservoir characterization of a carbonate ramp reservoir, South Dagger Draw field, New Mexico: Seismic data are only part of the story | |
| Etienam et al. | History matching of Reservoirs by Updating Fault Properties Using 4D Seismic Results and Ensemble Kalman Filter | |
| Inks et al. | Prospecting for gas hydrate accumulations using 2D and 3D seismic data, Milne Point, North Slope, Alaska | |
| Samantray et al. | Interpretation and application of borehole image logs in a new generation of reservoir models for a cluster of fields in southern Oman | |
| Robinson et al. | 4D time-lapse monitoring of Chirag Field | |
| Sultan et al. | Structural Dip Estimation from Crossbedding on Borehole Images | |
| Fang et al. | Applying integrated seismic technology to complex foothill areas of foreland basins in China | |
| Xu et al. | Integration Adds Value to Deepwater Oil Exploration and Development: A Case Study of North Gulf of Mexico From Seismic and Well Log to Simulation | |
| Rijfkogel | Integrated analysis of seismic attributes and well-logs in reservoir characterization: seismic-facies classification and reservoir facies mapping | |
| Sanyal et al. | Integrated Formation Evaluation Using High Resolution Sedimentary and Facies Analysis an Aid for Field Development: A Case Study from North Kuwait |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |