CN110673209B - Well-seismic calibration method - Google Patents
Well-seismic calibration method Download PDFInfo
- Publication number
- CN110673209B CN110673209B CN201910970415.4A CN201910970415A CN110673209B CN 110673209 B CN110673209 B CN 110673209B CN 201910970415 A CN201910970415 A CN 201910970415A CN 110673209 B CN110673209 B CN 110673209B
- Authority
- CN
- China
- Prior art keywords
- standard layer
- value
- seismic
- well
- speed
- 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
- 238000000034 method Methods 0.000 title claims abstract description 40
- 239000010410 layer Substances 0.000 claims abstract description 92
- 238000006243 chemical reaction Methods 0.000 claims abstract description 13
- 239000011229 interlayer Substances 0.000 claims abstract description 6
- 238000005553 drilling Methods 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 238000003908 quality control method Methods 0.000 claims description 2
- 238000004062 sedimentation Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013213 extrapolation 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
-
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a well-seismic calibration method, which comprises the steps of firstly, preparing seismic data and logging data; secondly, selecting a top and bottom standard layer based on the seismic interpretation scheme; thirdly, calibrating and controlling the top standard layer; fourthly, calibrating and controlling the quality of the bottom standard layer; fifthly, establishing a speed field based on the calibration result of the top and bottom standard layers; and sixthly, generating a synthetic seismic record and carrying out interlayer calibration of the well seismic. The method takes the top and bottom reference layers as constraint, adopts sound wave velocity to replace sound wave conversion to synthesize the seismic recording velocity between layers, eliminates space errors by means of average velocity fitting, trend surface low frequency effect removal and the like, and obtains top and bottom and standard interlayer precise velocity.
Description
The technical field is as follows:
the invention relates to the technical field of seismic data interpretation in seismic exploration, in particular to a well seismic calibration method.
Background art:
at present, the joint construction and interpretation and reservoir prediction by using seismic and logging data become the indispensable choice for fine reservoir description, but the longitudinal scale of the logging data is depth, the seismic profile is time scale, and the two data description modes are different, so that the two data cannot be directly applied in a combined manner, and therefore, the well seismic calibration between the logging data and the seismic data must be firstly carried out for developing the fine reservoir description.
The existing well-seismic calibration method is calibrated based on the correlation between synthetic seismic records and well-side seismic trace waveforms, and mainly has two problems, namely, the difference between well information and seismic data in time and waveform exists, the well logging data are measured in a depth domain, the velocity is the well-side vertical velocity, the seismic data are measured in a time domain, and the velocity is the velocity of a root-mean-square conversion layer with an incidence angle, so that a closing time error exists between the well information after time-depth conversion and ground seismic data, and the well-seismic calibration precision is not high. Secondly, the stratum lateral heterogeneity is strong, the well calibration information is difficult to extrapolate, the calibration precision of the well position is high, but when the stratum lateral heterogeneity is strong, the calibration precision of the well-free position can not be reliably guaranteed when the well calibration information is extrapolated.
The invention content is as follows:
the invention aims to provide a well-seismic calibration method which is used for solving the problem that the existing well-seismic calibration method is low in precision.
The technical scheme adopted by the invention for solving the technical problems is as follows: the well-seismic calibration method comprises the following steps:
preparing seismic imaging data, logging interpretation data, seismic interpretation standard horizon data and logging interpretation standard horizon data according to the required work area range and the logging number;
selecting a top standard layer and a bottom standard layer based on the seismic interpretation scheme, and selecting stratums with strong continuity of in-phase axes and stable sedimentation characteristics from the interpretation scheme as the top standard layer and the bottom standard layer respectively;
step three, calibrating and controlling the top standard layer:
a, performing linear fitting on all acoustic logging speed curves in a work area, and solving a fitting speed value v of each well at the top standard layer position according to the fitting curves;
b, according to the fitting speed value, a time value t1 of each well at the position of the top standard layer is calculated by using a conventional time-depth conversion method, a difference value t3 exists between the time value and the seismic interpretation time t2 of the top standard layer, namely t2-t1= t3, a trend surface of all the difference values in the work area is calculated by using a least square method, and the numerical value on the trend surface is recorded as t4;
c, correcting t1 in the step b by using the value t4 of the trend surface, wherein the time value of each well at the position of the top standard layer after correction is t5= t1+ t4;
d, knowing that the drilling platform and the seismic data reference surface have static correction time difference t6, performing residual time difference correction on t5 in the step c to obtain a final calibration time value t7= t5+ t6 of the position of the top standard layer of each well;
step four, bottom standard layer calibration and quality control:
(1) leveling the top standard layer, taking the leveled top standard layer as a time starting zero point of speed fitting, performing linear fitting on the acoustic logging speeds between all top and bottom layers in the work area, and calculating a fitting speed value w of each well at the position of the bottom standard layer according to a fitting curve;
(2) according to the fitting speed value, a time value t11 of each well at the position of the bottom standard layer is obtained by using a conventional time-depth conversion method, a difference value t33 exists between the time value and the seismic interpretation time t22 of the bottom standard layer, namely t22-t11= t33, a trend surface of all the difference values in the work area is obtained by using a least square method, and the value on the trend surface is recorded as t 44;
(3) using the value t44 of the trend surface to correct t11 in the step (2), wherein the time value of each well at the bottom standard layer position after correction is t55= t11+ t 44;
establishing a speed field based on the calibration result of the top and bottom standard layers;
and step six, generating a synthetic seismic record based on the velocity field in the step five and a conventional time-depth conversion method, and carrying out interlayer calibration of the well seismic by using the synthetic seismic record.
The fifth step in the scheme is specifically as follows:
a, calculating a difference value between a fitting speed value v corresponding to each well in a top standard layer and a speed of the standard layer, calculating a trend surface of the difference values corresponding to all wells in a work area by using a least square method, and recording a numerical value on the trend surface as vv;
b, correcting the fitting speed v corresponding to the top standard layer by using the numerical value vv on the trend surface to obtain the corrected top standard layer speed v1;
c, calculating the difference value between the fitting speed value w corresponding to each well in the bottom standard layer and the speed of the standard layer, calculating the trend surface of the difference values corresponding to all wells in the work area by using a least square method, and recording the numerical value on the trend surface as ww;
d, correcting the fitting speed w corresponding to the bottom standard layer by using the numerical value ww on the trend surface to obtain a corrected bottom standard layer speed w1;
and e, establishing a velocity field between the top standard layer and the bottom standard layer by using v1 and w1 and adopting a three-dimensional interpolation smoothing method.
The invention has the following beneficial effects:
1. the invention can synthesize the seismic recording speed by using the top and bottom reference layers as constraint and replacing acoustic wave conversion with acoustic wave speed between layers, eliminate space errors by means of mean speed fitting, trend surface low-frequency effect removal and the like, and obtain the well seismic calibration method of accurate speed between the top and bottom and standard layers.
2. The method can solve the problem of mismatching of well information and seismic data in time and waveform, reduces the influence of transverse heterogeneity on well calibration speed extrapolation, and has important application value on oil and gas resource exploration.
Drawings
FIG. 1 is a flow chart of the present invention.
Detailed Description
The invention is further illustrated below:
the well-seismic calibration method comprises the following steps:
firstly, preparing seismic imaging data, well logging interpretation data, seismic interpretation standard horizon data and well logging interpretation standard horizon data according to the required work area range and the well logging quantity.
And secondly, selecting stratums with strong in-phase axis continuity and stable deposition characteristics from the explanation scheme as a top standard layer and a bottom standard layer respectively.
Thirdly, calibrating and controlling the top standard layer:
the method comprises the following specific steps:
a, performing linear fitting on all acoustic logging speed curves in a work area, and solving a fitting speed value v of each well at the top standard layer position according to the fitting curves;
b, according to the fitting speed value, a time value t1 of each well at the position of the top standard layer is calculated by using a conventional time-depth conversion method, a difference value t3 exists between the time value and the seismic interpretation time t2 of the top standard layer, namely t2-t1= t3, a trend surface of all the difference values in the work area is calculated by using a least square method, and the numerical value on the trend surface is recorded as t4;
c, correcting t1 in the step b by using the value t4 of the trend surface, wherein the time value of each well at the position of the top standard layer after correction is t5= t1+ t4;
d, knowing that the drilling platform and the seismic data datum plane have static correction difference t6, performing residual moveout correction on t5 in the step c to obtain a final calibration time value t7= t5+ t6 of the position of the top standard layer of each well.
Fourthly, calibrating and controlling the bottom standard layer:
the method comprises the following specific steps:
a, leveling a top standard layer, taking the leveled top standard layer as a time initial zero point of speed fitting, performing linear fitting on all interlayer acoustic logging speeds in a work area, and calculating a fitting speed value w of each well at the position of the bottom standard layer according to a fitting curve;
b, according to the fitting speed value, a time value t11 of each well at the position of the bottom standard layer is calculated by using a conventional time-depth conversion method, a difference value t33 exists between the time value and the seismic interpretation time t22 of the bottom standard layer, namely t22-t11= t33, a trend surface of all the difference values in the work area is calculated by using a least square method, and the numerical value on the trend surface is recorded as t 44;
c, using the value t44 of the trend surface to correct the t11 in the step b, wherein the time value of each well at the bottom standard layer position after correction is t55= t11+ t 44.
Fifthly, establishing a speed field based on the calibration result of the top and bottom standard layers:
the method comprises the following specific steps:
a, calculating a difference value between a fitting speed value v corresponding to each well in a top standard layer and a speed of the standard layer, calculating a trend surface of the difference values corresponding to all wells in a work area by using a least square method, and recording a numerical value on the trend surface as vv;
c, calculating the difference value between the fitting speed value w of each well corresponding to the bottom standard layer and the speed of the standard layer, and calculating the trend surface of the corresponding difference values of all wells in the work area by using a least square method, wherein the numerical value on the trend surface is marked as ww;
d, correcting the fitting speed w corresponding to the bottom standard layer by using the numerical value ww on the trend surface to obtain a corrected bottom standard layer speed w1;
and e, establishing a velocity field between the top standard layer and the bottom standard layer by using v1 and w1 and adopting a three-dimensional interpolation smoothing method.
And sixthly, generating a synthetic seismic record based on the velocity field in the step five and a conventional time-depth conversion method, and carrying out interlayer calibration of the well seismic by using the synthetic seismic record.
Claims (1)
1. A well-seismic calibration method is characterized by comprising the following steps:
preparing seismic imaging data, logging interpretation data, seismic interpretation standard horizon data and logging interpretation standard horizon data according to the required work area range and the logging number;
selecting a top standard layer and a bottom standard layer based on the seismic interpretation scheme, and selecting stratums with strong continuity of in-phase axes and stable sedimentation characteristics from the interpretation scheme as the top standard layer and the bottom standard layer respectively;
step three, calibrating and controlling the top standard layer:
a, performing linear fitting on all acoustic logging speed curves in a work area, and solving a fitting speed value v of each well at the top standard layer position according to the fitting curves;
b, according to the fitting speed value, a time value t1 of each well at the position of the top standard layer is calculated by using a conventional time-depth conversion method, a difference value t3 exists between the time value and the seismic interpretation time t2 of the top standard layer, namely t2-t1= t3, a trend surface of all the difference values in the work area is calculated by using a least square method, and the numerical value on the trend surface is recorded as t4;
c, correcting t1 in the step b by using the value t4 of the trend surface, wherein the time value of each well at the position of the top standard layer after correction is t5= t1+ t4;
d, knowing that the drilling platform and the seismic data reference surface have static correction time difference t6, performing residual time difference correction on t5 in the step c to obtain a final calibration time value t7= t5+ t6 of the position of the top standard layer of each well;
step four, bottom standard layer calibration and quality control:
(1) leveling the top standard layer, taking the leveled top standard layer as a time starting zero point of speed fitting, performing linear fitting on the acoustic logging speeds between all top and bottom layers in the work area, and calculating a fitting speed value w of each well at the position of the bottom standard layer according to a fitting curve;
(2) according to the fitting speed value, a time value t11 of each well at the position of the bottom standard layer is obtained by using a conventional time-depth conversion method, a difference value t33 exists between the time value and the seismic interpretation time t22 of the bottom standard layer, namely t22-t11= t33, a trend surface of all the difference values in the work area is obtained by using a least square method, and the value on the trend surface is recorded as t 44;
(3) using the value t44 of the trend surface to correct t11 in the step (2), wherein the time value of each well at the bottom standard layer position after correction is t55= t11+ t 44;
step five, establishing a speed field based on the calibration result of the top and bottom standard layers:
a, calculating a difference value between a fitting speed value v corresponding to each well in a top standard layer and a speed of the standard layer, calculating a trend surface of the difference values corresponding to all wells in a work area by using a least square method, and recording a numerical value on the trend surface as vv;
b, correcting the fitting speed v corresponding to the top standard layer by using the numerical value vv on the trend surface to obtain the corrected top standard layer speed v1;
c, calculating the difference value between the fitting speed value w corresponding to each well in the bottom standard layer and the speed of the standard layer, calculating the trend surface of the difference values corresponding to all wells in the work area by using a least square method, and recording the numerical value on the trend surface as ww;
d, correcting the fitting speed w corresponding to the bottom standard layer by using the numerical value ww on the trend surface to obtain a corrected bottom standard layer speed w1;
e, establishing a velocity field between the top standard layer and the bottom standard layer by using v1 and w1 and adopting a three-dimensional interpolation smoothing method;
and step six, generating a synthetic seismic record based on the velocity field in the step five and a conventional time-depth conversion method, and carrying out interlayer calibration of the well seismic by using the synthetic seismic record.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910970415.4A CN110673209B (en) | 2019-10-13 | 2019-10-13 | Well-seismic calibration method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910970415.4A CN110673209B (en) | 2019-10-13 | 2019-10-13 | Well-seismic calibration method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110673209A CN110673209A (en) | 2020-01-10 |
| CN110673209B true CN110673209B (en) | 2021-06-04 |
Family
ID=69082196
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910970415.4A Active CN110673209B (en) | 2019-10-13 | 2019-10-13 | Well-seismic calibration method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110673209B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112147717B (en) * | 2020-09-27 | 2021-11-16 | 北京中恒利华石油技术研究所 | Pre-stack and post-stack well-seismic combination time depth calibration method |
| CN112269212A (en) * | 2020-10-20 | 2021-01-26 | 中国石油天然气集团有限公司 | Method, device, equipment and medium for determining seismic interpretation horizon of small logging layering |
| CN117930359B (en) * | 2024-03-21 | 2024-06-14 | 东北石油大学三亚海洋油气研究院 | Method for correcting fault shadows in post-stack seismic data by establishing high-precision velocity fields |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010005762A2 (en) * | 2008-07-09 | 2010-01-14 | Chevron U.S.A. Inc. | Well log normalization |
| CN104047597A (en) * | 2014-01-26 | 2014-09-17 | 中国石油大学(华东) | Fat gas mud shale stratum well log standardizing method |
| CN104714250A (en) * | 2014-11-07 | 2015-06-17 | 中国石油化工股份有限公司 | Practical internal substratum automatic interpretation method |
| CN105759312A (en) * | 2016-03-01 | 2016-07-13 | 中国海洋石油总公司 | Reservoir characteristic method well-seismic calibration method |
| CN106569280A (en) * | 2016-10-28 | 2017-04-19 | 中国石油化工股份有限公司 | Calibration method for improving seismic interpretation precision of fault terrace structures |
| WO2017142719A1 (en) * | 2016-02-16 | 2017-08-24 | Schlumberger Technology Corporation | Calibrating seismic data using measurements made during drilling operations |
| CN108121009A (en) * | 2016-11-30 | 2018-06-05 | 中国石油化工股份有限公司 | Development late stage complex fault block velocity-variant depth mapping method |
| CN109932749A (en) * | 2019-03-05 | 2019-06-25 | 中国石油天然气股份有限公司 | A kind of well shake scaling method, device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8107317B2 (en) * | 2006-12-28 | 2012-01-31 | Schlumberger Technology Corporation | Technique and system for performing a cross well survey |
-
2019
- 2019-10-13 CN CN201910970415.4A patent/CN110673209B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2010005762A2 (en) * | 2008-07-09 | 2010-01-14 | Chevron U.S.A. Inc. | Well log normalization |
| CN104047597A (en) * | 2014-01-26 | 2014-09-17 | 中国石油大学(华东) | Fat gas mud shale stratum well log standardizing method |
| CN104714250A (en) * | 2014-11-07 | 2015-06-17 | 中国石油化工股份有限公司 | Practical internal substratum automatic interpretation method |
| WO2017142719A1 (en) * | 2016-02-16 | 2017-08-24 | Schlumberger Technology Corporation | Calibrating seismic data using measurements made during drilling operations |
| CN105759312A (en) * | 2016-03-01 | 2016-07-13 | 中国海洋石油总公司 | Reservoir characteristic method well-seismic calibration method |
| CN106569280A (en) * | 2016-10-28 | 2017-04-19 | 中国石油化工股份有限公司 | Calibration method for improving seismic interpretation precision of fault terrace structures |
| CN108121009A (en) * | 2016-11-30 | 2018-06-05 | 中国石油化工股份有限公司 | Development late stage complex fault block velocity-variant depth mapping method |
| CN109932749A (en) * | 2019-03-05 | 2019-06-25 | 中国石油天然气股份有限公司 | A kind of well shake scaling method, device |
Non-Patent Citations (1)
| Title |
|---|
| 井震结合的测井数据标准化处理技术——以克拉玛依油田八区为例;张胜 等;《石油物探》;20120930;第51卷(第5期);第531-538页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110673209A (en) | 2020-01-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110673209B (en) | Well-seismic calibration method | |
| US7840394B2 (en) | Method for generating a 3D earth model | |
| US8670288B2 (en) | Velocity model for well time-depth conversion | |
| CN109738945B (en) | Method for directly generating construction diagram by using prestack depth migration result | |
| CN108663713B (en) | Method for establishing depth domain structure model | |
| CN112505753B (en) | Seismic-driven high-precision sequence grid model construction method | |
| WO2016135505A1 (en) | A method of producing and utilising high resolution impedance logs derived from vsp data for use in assessing an oilfield subterranean formation | |
| CN105223617A (en) | Seismic properties is utilized to control modeling method and the system of impedance initial value | |
| US8547793B2 (en) | Correction of velocity cubes for seismic depth modeling | |
| US8599644B2 (en) | Velocity models for a single well and for a set of wells | |
| CN111257946B (en) | Geophysical drilling guiding method and method for updating stratum seismic velocity | |
| US6073079A (en) | Method of maintaining a borehole within a multidimensional target zone during drilling | |
| CN110927796A (en) | Method for improving time-depth conversion precision of seismic data | |
| CN111257932B (en) | Method for updating stratum seismic velocity | |
| CN111456726B (en) | Acoustic wave time difference correction method and device | |
| CN113945978B (en) | Depth domain seismic data optimization method and device | |
| CN112529981B (en) | Method for mapping contour lines of deep-area coal seam floor | |
| CN116500679A (en) | Well shock calibration method and device and related equipment | |
| CN116840909A (en) | Combined depth domain modeling method for hidden mountain inner curtain | |
| CN116931089B (en) | Time depth calibration method and device based on thickness partition, electronic equipment and medium | |
| CN113589375A (en) | VSP layer velocity inversion method based on inclined layer constraint travel time calculation | |
| CN112305595A (en) | Method for analyzing geologic body structure based on refracted wave and storage medium | |
| CN111722276A (en) | Seismic guiding method and system for rock drilling | |
| CN112394391B (en) | Seismic facies quantitative characterization method of Rake wavelet fractional order derivative | |
| Rowbotham et al. | The implications of anisotropy for seismic impedance 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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |