CN113740912A - Full-stratum quality factor Q body building method - Google Patents
Full-stratum quality factor Q body building method Download PDFInfo
- Publication number
- CN113740912A CN113740912A CN202111127727.2A CN202111127727A CN113740912A CN 113740912 A CN113740912 A CN 113740912A CN 202111127727 A CN202111127727 A CN 202111127727A CN 113740912 A CN113740912 A CN 113740912A
- Authority
- CN
- China
- Prior art keywords
- quality factor
- vsp
- velocity
- layer velocity
- layer
- 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
Links
- 238000000034 method Methods 0.000 title description 10
- 239000010410 layer Substances 0.000 claims abstract description 40
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 239000002344 surface layer Substances 0.000 claims abstract description 9
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 6
- 238000010276 construction Methods 0.000 claims abstract description 5
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000010521 absorption reaction Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000003475 lamination Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/303—Analysis for determining velocity profiles or travel times
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6222—Velocity; travel time
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
The invention provides a full formation quality factor Q construction method, which comprises the following steps: s1, inputting the quality factor Q of the micro loggingwellLayer velocity VwellSuperficial velocity body VsVSP log quality factor QvspLayer velocity VvspMedium-depth layer-by-layer velocity body Vdeep(ii) a S2, counting the micro logging zone velocity V of the step S1wellAnd quality factor QwellUsing the relational expression (A), the surface layer velocity body V of step S1sConversion into surface quality factor Qs(ii) a S3, counting the VSP log layer velocity V of the step S1vspAnd quality factor QvspThe relation of (A) is used to obtain the intermediate layer velocity volume V of step S1deepConversion into surface quality factor Qdeep. The invention obtains the statistical relationship between the quality factor Q and the layer velocity through multi-point micro logging and multi VSP, and converts the velocity body into the Q body of the whole stratum according to the statistical relationship, thereby being fast, simple and effective.
Description
Technical Field
The invention relates to the field of geophysical exploration, in particular to a method for building a full-stratum quality factor Q body, which is a method for building the full-stratum quality factor Q body by utilizing a stratum velocity body of a micro-logging and zero-deviation vertical seismic profile.
Background
The underground medium has the absorption and attenuation effects on seismic waves, and the resolution of seismic data is seriously reduced. The quality factor Q is used as the representation of the stratum absorption attenuation attribute, and is beneficial to improving the resolution ratio of seismic data.
Zhang Jianfeng et al, the patent Q value field modeling method based on ground-received reflection seismic data, which utilizes the ground-received reflection seismic data to realize a quality factor Q volume building method, is applied to seismic exploration to improve the resolution of seismic imaging. The article "research on application technology based on seismic data Q field establishment and compensation processing" in Sharp et al, finds a quality factor Q corresponding to a whole journey through a seismic data time-sharing window, and establishes a space-variant Q body.
The existing methods are all deep Q body building methods. According to the invention, a surface Q body is established by using multi-point micro logging, and a middle-deep Q body is established by using multi-point VSP, so that the whole stratum Q body is built.
Disclosure of Invention
The invention provides a method for converting a surface layer velocity body into a surface layer quality factor Q body by utilizing multipoint micro-logging and counting the corresponding relation between the layer velocity and the quality factor Q; carrying out multi-point VSP logging, counting the corresponding relation between the interval velocity and the quality factor Q, and converting the intermediate-deep layer velocity body into an intermediate-deep layer quality factor Q body; and establishing a whole stratum quality factor body.
The method aims to utilize the corresponding relation between the speed of a multi-point micro logging and multi-point VSP logging statistical analysis layer and a quality factor to be popularized to a three-dimensional space to establish a full-stratum quality factor Q body.
A full formation quality factor Q construction method comprises the following steps:
s1, inputting the micro-logging productQuality factor QwellLayer velocity VwellSuperficial velocity body VsVSP log quality factor QvspLayer velocity VvspMedium-depth layer-by-layer velocity body Vdeep。
S2, counting the micro logging zone velocity V of the step S1wellAnd quality factor QwellUsing the relational expression (A), the surface layer velocity body V of step S1sConversion into surface quality factor Qs。
S21, using the layer velocity V of the micro logging input in the step S1wellQuality factor QwellAnd counting coefficients a and b on the micro logging point by adopting a formula in a power exponent form:
s22, using the coefficients a and b calculated in step S21, the surface layer velocity V input in step S1sConversion into surface quality factor body Qs:
S3, counting the VSP log layer velocity V of the step S1vspAnd quality factor QvspThe relation of (A) is used to obtain the intermediate layer velocity volume V of step S1deepConversion into surface quality factor Qdeep。
S31, layer velocity V of VSP input in step S1vspQuality factor QvspAnd counting coefficients c and d on the VSP logging point by adopting a formula in a power exponent form:
s32, using the coefficients c and d obtained by calculation in step S31, the medium-depth layer velocity body V input in step S1deepConversion into surface quality factor body Qdeep:
The invention obtains the statistical relationship between the quality factor Q and the layer velocity through multi-point micro logging and multi VSP, and converts the velocity body into the Q body of the whole stratum according to the statistical relationship, thereby being fast, simple and effective.
Drawings
FIG. 1 is a figure of merit for an example micro-log; the abscissa is Q (unit: dimensionless); the ordinate is the depth (unit: meter).
FIG. 2 is a interval velocity of an example micro-log; the abscissa is the velocity (unit: m/s); the ordinate is the depth (unit: meter).
FIG. 3 is an example table layer velocity volume; the abscissa is the line number (unit: none); the ordinate is the depth (unit: meter).
FIG. 4 is a figure of merit for an example VSP log; the abscissa is Q (unit: dimensionless); the ordinate is the depth (unit: meter).
FIG. 5 is interval velocity of an example VSP log; the abscissa is the velocity (unit: m/s); the ordinate is the depth (unit: meter).
FIG. 6 is a deep layer velocity volume of an embodiment; the abscissa is the line number (unit: none); the ordinate is the depth (unit: meter).
FIG. 7 is a statistical relationship between quality factor Q and velocity for microlog according to an embodiment; the abscissa is the velocity (unit: km/s); the ordinate is the Q value (unit: dimensionless).
FIG. 8 is an example surface quality factor Qsome; the abscissa line number (unit: none); the ordinate is the depth (unit: meter).
FIG. 9 is a statistical relationship between the deep quality factor Q and the velocity in the embodiment; the abscissa is the velocity (unit: km/s); the ordinate is the Q value (unit: dimensionless).
FIG. 10 is a deep quality factor Qsome of the examples; the abscissa is the line number (unit: none); the ordinate is the depth (unit: meter).
Detailed Description
The specific technical scheme of the invention is described by combining the embodiment.
S1, inputting the quality factor Q of the micro loggingwellLayer velocity VwellSuperficial velocity body VsVSP log quality factor QvspLayer velocity VvspMedium-depth layer-by-layer velocity body Vdeep。
Inputting a micro-logging quality factor, as shown in FIG. 1;
input micro-log formation velocity, fig. 2;
the input superficial layer velocity body, as shown in fig. 3;
the quality factor of the input VSP log, fig. 4;
input VSP log interval velocity, fig. 5;
and inputting a medium-depth layer lamination velocity body, as shown in figure 6.
S2, counting the micro logging zone velocity V of the step S1wellAnd quality factor QwellUsing the relational expression (A), the surface layer velocity body V of step S1sConversion into surface quality factor Qs。
S21, using the layer velocity V of the micro logging input in the step S1wellQuality factor QwellAnd counting coefficients a and b on the micro logging point by adopting a formula in a power exponent form:
s22, using the coefficients a and b calculated in step S21, the surface layer velocity V input in step S1sConversion into surface quality factor body Qs:
The statistical micro-log quality factor Q is related to velocity, as shown in fig. 7.
The surface quality factor Q-body was established as shown in fig. 8.
S3, statistics of VSP logging in step S1Layer velocity VvspAnd quality factor QvspUsing the relation of (1) to obtain the intermediate-depth layer velocity body VdeepConversion into surface quality factor Qdeep。
S31, layer velocity V of VSP input in step S1vspQuality factor QvspAnd counting coefficients c and d on the VSP logging point by adopting a formula in a power exponent form:
s32 using the coefficients c and d calculated in step S31 to obtain the intermediate layer velocity volume V input in step S1deepConversion into surface quality factor body Qdeep:
The statistical relationship between the quality factor Q of the middle and deep layers and the velocity is shown in fig. 9.
The established Q body of the quality factor of the middle-deep layer is shown in figure 10.
Claims (3)
1. A full formation quality factor Q construction method is characterized by comprising the following steps:
s1, inputting the quality factor Q of the micro loggingwellLayer velocity VwellSuperficial velocity body VsVSP log quality factor QvspLayer velocity VvspMedium-depth layer-by-layer velocity body Vdeep;
S2, counting the micro logging zone velocity V of the step S1wellAnd quality factor QwellUsing the relational expression (A), the surface layer velocity body V of step S1sConversion into surface quality factor Qs;
S3, counting the VSP log layer velocity V of the step S1vspAnd quality factor QvspThe relation of (A) is used to obtain the intermediate layer velocity volume V of step S1deepConversion into surface quality factor Qdeep。
2. The full formation quality factor Q construction method according to claim 1, wherein step S2 comprises the following sub-steps:
s21, using the layer velocity V of the micro logging input in the step S1wellQuality factor QwellAnd counting coefficients a and b on the micro logging point by adopting a formula in a power exponent form:
s22, using the coefficients a and b calculated in step S21, the surface layer velocity V input in step S1sConversion into surface quality factor body Qs:
Qs=a×Vs b。
3. The full formation quality factor Q construction method according to claim 1, wherein step S3 comprises the following sub-steps:
s31, layer velocity V of VSP input in step S1vspQuality factor QvspAnd counting coefficients c and d on the VSP logging point by adopting a formula in a power exponent form:
s32, using the coefficients c and d obtained by calculation in step S31, the medium-depth layer velocity body V input in step S1deepConversion into surface quality factor body Qdeep:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111127727.2A CN113740912A (en) | 2021-09-26 | 2021-09-26 | Full-stratum quality factor Q body building method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111127727.2A CN113740912A (en) | 2021-09-26 | 2021-09-26 | Full-stratum quality factor Q body building method |
Publications (1)
Publication Number | Publication Date |
---|---|
CN113740912A true CN113740912A (en) | 2021-12-03 |
Family
ID=78741179
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111127727.2A Pending CN113740912A (en) | 2021-09-26 | 2021-09-26 | Full-stratum quality factor Q body building method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113740912A (en) |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104502977A (en) * | 2014-12-22 | 2015-04-08 | 中国石油天然气集团公司 | Well-control amplitude-preservation high-resolution seismic data processing method |
CN104635268A (en) * | 2015-03-09 | 2015-05-20 | 成都晶石石油科技有限公司 | Method for calculating quality factor under seismic data constraint |
CN105607124A (en) * | 2016-03-09 | 2016-05-25 | 蒋立 | Seismic-wave near-surface stratum quality factor compensation method and device |
CN106908838A (en) * | 2017-03-15 | 2017-06-30 | 徐诗薇 | The method for building target area stratum inelastic attenuation quality factor three-dimensional model |
CN109884707A (en) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | Near surface is layered time-depth curve static correcting method |
-
2021
- 2021-09-26 CN CN202111127727.2A patent/CN113740912A/en active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104502977A (en) * | 2014-12-22 | 2015-04-08 | 中国石油天然气集团公司 | Well-control amplitude-preservation high-resolution seismic data processing method |
CN104635268A (en) * | 2015-03-09 | 2015-05-20 | 成都晶石石油科技有限公司 | Method for calculating quality factor under seismic data constraint |
CN105607124A (en) * | 2016-03-09 | 2016-05-25 | 蒋立 | Seismic-wave near-surface stratum quality factor compensation method and device |
CN106908838A (en) * | 2017-03-15 | 2017-06-30 | 徐诗薇 | The method for building target area stratum inelastic attenuation quality factor three-dimensional model |
CN109884707A (en) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | Near surface is layered time-depth curve static correcting method |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109884707B (en) | Static correction method for near-surface layered time-depth curve | |
CN106869911B (en) | Evaluation method for describing compressibility of shale reservoir | |
CN107783187B (en) | Method for establishing three-dimensional velocity field by combining logging velocity and seismic velocity | |
CN111257926B (en) | Method for predicting ancient valley uranium reservoir by using old seismic data | |
CN106556866B (en) | A method of using rock physics template to reservoir quantitative interpretation | |
WO2017024523A1 (en) | Inversion method for ray elastic parameter | |
CN109884710B (en) | Micro-logging tomography method aiming at excitation well depth design | |
CN103675907A (en) | AVO inversion hydrocarbon detection method based on petrographic constraints | |
CN102288997A (en) | Technology for inverting near-surface quality factor three-dimensional data volume by using microlog | |
CN105259581B (en) | Seismic data time-depth conversion method | |
CN109521474A (en) | It is a kind of three-dimensional dual control under prestack Inversion of geostatistics | |
CN103399346A (en) | Well-to-seismic integration initial wave impedance modeling method | |
CN109444959A (en) | Full range high-precision interval velocity field method for building up | |
CN110824563A (en) | Reservoir lithology prediction method based on Xgboost algorithm | |
CN104570127A (en) | Method for simultaneously inverting porosity and shaliness by using seismic wave velocity | |
KR101415196B1 (en) | Method for selecting fracking intervals of horizontal drilling zone in sweet spot range using resistivity and density logging data in shale play | |
CN101625419A (en) | Method for improving precision of static correction by uphole time | |
CN113740912A (en) | Full-stratum quality factor Q body building method | |
CN107422384B (en) | A method of vertical seismic profile data is imitated using well-log information | |
CN108708715A (en) | The survey logging method of fast prediction shale gas-bearing formation formation fracture pressure gradient | |
CN104484573A (en) | Novel calculating method for rigidity coefficient of stratum | |
CN112255682B (en) | Q body modeling technical method based on VSP data | |
CN102998700A (en) | Surface consistency compensation method utilizing vertical seismic section data | |
CN104459820B (en) | The method for building up of multiple dimensioned pressure field | |
CN104459771A (en) | Reservoir gas-bearing semi-quantitative recognition method based on frequency division AVO inversion |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20211203 |
|
RJ01 | Rejection of invention patent application after publication |