CN110579806B - Rapid well seismic calibration method under well pattern tight condition - Google Patents
Rapid well seismic calibration method under well pattern tight condition Download PDFInfo
- Publication number
- CN110579806B CN110579806B CN201910856307.4A CN201910856307A CN110579806B CN 110579806 B CN110579806 B CN 110579806B CN 201910856307 A CN201910856307 A CN 201910856307A CN 110579806 B CN110579806 B CN 110579806B
- Authority
- CN
- China
- Prior art keywords
- well
- seismic
- horizon
- time domain
- depth
- 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 27
- 238000005553 drilling Methods 0.000 claims abstract description 10
- 239000003550 marker Substances 0.000 claims description 10
- 238000005259 measurement Methods 0.000 claims description 7
- 230000000694 effects Effects 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000002349 favourable effect Effects 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/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)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (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 rapid well-seismic calibration method under a dense well pattern condition. The method mainly solves the problems that the prior well-to-seismic calibration method is low in speed and can influence the subsequent development and application efficiency under the condition of a dense well pattern. The method is characterized in that: 1) evenly selecting 1/10-1/4 wells of the total number of wells in the research area, and making single-well depth domain synthetic seismic records; 2) adjusting the single well depth domain synthetic seismic record, and establishing a corresponding relation between time and depth of the single well; 3) converting the geological horizon of the single well from a depth domain to a time domain, and establishing a time domain geological horizon model; 4) and automatically matching the residual drilling horizon of the research area to the time domain geological horizon model, thereby realizing the rapid well-to-seismic calibration work of all drilling wells. The well seismic calibration method can improve the well seismic calibration speed on the premise of ensuring the well seismic calibration precision.
Description
Technical Field
The invention relates to the technical field of oil field development, in particular to a rapid well seismic calibration method under a dense well pattern condition.
Background
The dense well pattern referred in the patent of the invention refers to that the number of drilled wells is more than 50 holes/square kilometer, and at present, many old oil fields in China, including Daqing oil field, Liaohe oil field, Jilin oil field, etc., all meet the conditions of the dense well pattern. At present, in the field of oil development, the combination of well drilling data and seismic data for oil field development and application becomes a common technology, but because the well drilling data belongs to a depth domain and the seismic data belongs to a time domain, well seismic calibration is required to be carried out in order to classify the well drilling data and the seismic data into the same domain.
The method is characterized in that a single well synthetic seismic record is made, and the single well synthetic seismic record is matched with an actually measured seismic channel beside a well in a coaxial mode by adopting an artificial method. And secondly, the actual measurement seismic marker axis is automatically matched with the geological horizon, the method is low in precision, the correlation between the synthetic record of the well and the seismic waveform is low under the condition of a dense well pattern, and the subsequent development and application effects are influenced.
Disclosure of Invention
The invention provides a rapid well-seismic calibration method under the condition of a dense well pattern, aiming at overcoming the problems that the prior well-seismic calibration method in the background art is low in speed and the subsequent development and application effects are influenced due to large workload under the condition of a dense well pattern. The method for rapidly calibrating the well earthquake under the condition of the tight well pattern can improve the calibration speed of the well earthquake on the premise of ensuring the calibration precision of the well earthquake.
The invention can solve the problems by the following technical scheme: a fast well-seismic calibration method under the condition of a close well pattern comprises the following steps:
1) evenly selecting 1/10-1/4 well drills of the total number of wells in the research area, carrying out actual measurement on the selected well drills to obtain the product of a sound wave velocity well logging curve and a density well logging curve of a depth domain, and making a single well depth domain synthetic seismic record;
2) adjusting the single-well depth domain synthetic seismic record manufactured in the step (1), performing seismic measurement on a research area to obtain three-dimensional seismic data of a time domain, and moving the single-well depth domain synthetic seismic record up and down to enable the waveform of the single-well synthetic record of the depth domain to be matched with the waveform of the time domain seismic data near a well, so that the corresponding relation between the time and the depth of the single well is established;
3) converting the geological horizon of the single well from the depth domain to the time domain by utilizing the corresponding relation between the time and the depth of the single well established in the step (2), and establishing a time domain geological horizon model by adopting a seismic marker axis (the strongest waveform of seismic data and better continuity) constraint interpolation algorithm;
4) and (4) automatically matching (moving) the residual drilled well positions 9/10-3/4 in the research area to the time domain geological position model obtained in the step (3), so that the rapid well-seismic calibration work of all drilled wells is realized.
In the step (1), the single-well depth domain synthetic seismic record is obtained by multiplying the acoustic velocity logging curve and the density logging curve to obtain a longitudinal wave impedance curve of a depth domain, and then performing convolution operation on the longitudinal wave impedance curve and the wavelets;
in the step (3), the time domain geological horizon model is established in a concrete implementation mode by using the time domain horizon of a single well at a well point, and trend constraint is carried out among wells by using a seismic marker axis, so that a space horizon model is established; the specific method comprises the following steps:
known seismic marker axes and geologic horizons of different wells (A, B, C, D);
calculating the distance di (i =1,2,3,4) between the geological horizon and the seismic surface;
interpolating di to obtain dx at any position in the known surface range;
and iv, adding the distance dx to the axial surface of the known seismic mark to obtain the time domain geological horizon model.
The step (4) is specifically realized by the following steps:
(1) generating synthetic seismic records of 9/10-3/4 remaining boreholes in the study area;
(2) recording the time value of the corresponding horizon of each well, and calculating to obtain the time difference value delta t between the well horizon and the horizon model;
(3) and automatically moving the synthetic seismic record by a computer according to the size of delta t to correct the well point horizon to the corresponding horizon model, thereby quickly realizing the well seismic calibration work.
Compared with the background technology, the invention has the following beneficial effects: the invention provides a rapid well-seismic calibration method under the condition of a dense well pattern, which can realize rapid matching of well synthetic records and seismic waveforms under the condition of the dense well pattern, improves the speed by 4-10 times, ensures the well-seismic combined development and application effect, and has important significance for ensuring the stable yield of an oil field.
Description of the drawings:
FIG. 1 is a flow chart of the method for fast well seismic calibration of a tight well pattern according to the invention;
FIG. 2 is a diagram of the single well synthetic seismic record making and calibration with the well-side seismic traces in the example;
FIG. 3 is a schematic diagram of a time domain geological horizon model established by a seismic marker axis constrained interpolation algorithm in the embodiment;
FIG. 4 is an example time domain geological horizon model;
FIG. 5 is a matching graph of well horizons and horizon models for the 9/10 scale well numbers in an example.
The specific implementation mode is as follows:
the invention will be further described with reference to the following drawings and specific embodiments:
example 1
Taking the second north-west block of the Daqing oil field as an example, the second north-west block has 2000 wells in total, as shown in fig. 1, the research method for fast well-seismic calibration of the dense well pattern is divided into four steps:
uniformly selecting 1/10 well-numbered drilling wells (2000 wells are selected for the two northwest blocks, and 200 wells are selected), and actually measuring 200 wells to obtain a sound wave curve and a density curve of a depth domain; multiplying the acoustic wave curve obtained in the first step by the density curve to obtain a longitudinal wave impedance curve of a depth domain, and performing convolution operation on the longitudinal wave impedance curve and the wavelets to obtain a single-well depth domain synthetic seismic record of 200 wells, as shown in fig. 2;
secondly, adjusting the single-well depth domain synthetic seismic record manufactured in the step (1), performing seismic measurement on a research area to obtain three-dimensional seismic data of a time domain, moving the single-well depth domain synthetic seismic record up and down to enable the waveform of the single-well synthetic record of the depth domain to be matched with the waveform of the time domain seismic data near a well, and accordingly establishing the corresponding relation between the time and the depth of the single well;
thirdly, converting the geological horizon of the single well from a depth domain to a time domain by utilizing the corresponding relation between time and depth of the single well established in the second step, establishing a time domain geological horizon model by adopting a seismic marking axis (the strongest waveform of seismic data and better continuity) constraint interpolation algorithm, wherein the time domain geological horizon model is shown as a figure 4, the time domain geological horizon model is specifically realized by utilizing the time domain horizon of the single well at a well point and utilizing the seismic marking axis to carry out trend constraint among wells, so that a space horizon model is established; as shown in fig. 3, the specific method includes:
known seismic marker axes and geologic horizons of different wells (A, B, C, D);
calculating the distance di (i =1,2,3,4) between the geological horizon and the seismic surface;
interpolating di to obtain dx at any position in the known surface range;
and iv, adding the distance dx to the axial surface of the known seismic mark to obtain the time domain geological horizon model.
And fourthly, automatically matching (moving) the residual 9/10 well drilling horizons of the research area to the time domain geological horizon model obtained in the step (3) through a computer (as shown in figure 5), so that the rapid well seismic calibration work of all well drilling is realized. And automatically correcting the well point horizon to the corresponding horizon model according to the longitudinal relation between the geological horizon of the well and the horizon model.
The specific implementation process is as follows: first, synthetic seismic records of the remaining 9/10 drilled wells of the study area are generated; secondly, recording the time value of the corresponding horizon of each well, and calculating to obtain the time difference value delta t between the well horizon and the horizon model; and finally, automatically moving the synthetic seismic record by a computer according to the size of delta t to correct the well point horizon to the corresponding horizon model, thereby quickly realizing the well seismic calibration work.
The well-seismic calibration method can realize rapid well-seismic calibration under the condition of a dense well pattern, improves the application efficiency of well-seismic combination development, and is favorable for ensuring the stable yield of an oil field.
Compared with the prior art that the actual measurement seismic marker axis and the geological horizon are used for automatic matching (or in a manual mode), the well seismic calibration method has the advantages that the precision is guaranteed, and meanwhile the speed is improved by 4-10 times.
Claims (2)
1. A fast well-seismic calibration method under the condition of a close well pattern comprises the following steps:
(1) evenly selecting 1/10-1/4 well drills of the total number of wells in the research area, carrying out actual measurement on the selected well drills to obtain the product of a sound wave velocity well logging curve and a density well logging curve of a depth domain, and making a single well depth domain synthetic seismic record;
(2) adjusting the single-well depth domain synthetic seismic record manufactured in the step (1), performing seismic measurement on a research area to obtain three-dimensional seismic data of a time domain, and moving the single-well depth domain synthetic seismic record up and down to enable the waveform of the single-well depth domain synthetic seismic record to be matched with the waveform of the time domain seismic data near a well, so that the corresponding relation between the time and the depth of the single well is established;
(3) converting the geological horizon of the single well from the depth domain to the time domain by utilizing the corresponding relation between the time and the depth of the single well established in the step (2), and establishing a time domain geological horizon model by adopting a seismic marker axis constraint interpolation algorithm;
(4) automatically matching the residual well drilling horizons 9/10-3/4 in the research area to the time domain geological horizon model obtained in the step (3), thereby realizing the rapid well seismic calibration work of all well drilling;
the time domain geological horizon model is established in the step (3) in a mode that: the method is characterized in that time domain horizons of single wells are utilized at well points, trend constraint is carried out among wells by utilizing seismic marker axes, and the method specifically comprises the following steps:
knowing the seismic marker axis and geological horizons of different wells;
ii, calculating the distance di between the geological horizon and the seismic surface;
interpolating di to obtain dx at any position in the known surface range;
iv, adding the distance dx to the axial surface of the known seismic mark to obtain a time domain geological horizon model;
the step (4) is specifically realized by the following steps:
(1) generating synthetic seismic records of 9/10-3/4 remaining boreholes in the study area;
(2) recording the time value of the corresponding horizon of each well, and calculating to obtain the time difference value delta t between the well horizon and the horizon model;
(3) and automatically moving the synthetic seismic record by a computer according to the size of delta t to correct the well point horizon to the corresponding horizon model, thereby quickly realizing the well seismic calibration work.
2. The method for rapid well seismic calibration under the condition of a tight well pattern according to claim 1, characterized in that: and (2) in the step (1), the single-well depth domain synthetic seismic record is obtained by multiplying the acoustic velocity logging curve and the density logging curve to obtain a longitudinal wave impedance curve of the depth domain, and performing convolution operation on the longitudinal wave impedance curve and the wavelets.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910856307.4A CN110579806B (en) | 2019-09-11 | 2019-09-11 | Rapid well seismic calibration method under well pattern tight condition |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910856307.4A CN110579806B (en) | 2019-09-11 | 2019-09-11 | Rapid well seismic calibration method under well pattern tight condition |
Publications (2)
Publication Number | Publication Date |
---|---|
CN110579806A CN110579806A (en) | 2019-12-17 |
CN110579806B true CN110579806B (en) | 2022-03-22 |
Family
ID=68812844
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910856307.4A Active CN110579806B (en) | 2019-09-11 | 2019-09-11 | Rapid well seismic calibration method under well pattern tight condition |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN110579806B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111413737A (en) * | 2020-05-20 | 2020-07-14 | 大庆油田有限责任公司 | Method for explaining seismic horizon to deposition unit in close well pattern area |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2846175A2 (en) * | 2013-09-06 | 2015-03-11 | Services Petroliers Schlumberger | Seismic survey analysis |
CN105759312A (en) * | 2016-03-01 | 2016-07-13 | 中国海洋石油总公司 | Reservoir characteristic method well-seismic calibration method |
CN106324681A (en) * | 2016-08-24 | 2017-01-11 | 中国石油化工股份有限公司 | Time-frequency continuous wavelet transform well-to-seismic calibration method |
CN107783187A (en) * | 2016-08-24 | 2018-03-09 | 中国石油化工股份有限公司 | It is a kind of that logging speed and seismic velocity are combined to the method for establishing 3D velocity field |
CN107831542A (en) * | 2017-07-24 | 2018-03-23 | 中国石油化工股份有限公司 | DDW high accuracy depths domain well shake matching process |
CN108254785A (en) * | 2017-12-28 | 2018-07-06 | 中国石油天然气股份有限公司 | Reservoir Determination Method and Device |
KR101893800B1 (en) * | 2017-08-09 | 2018-09-04 | 제주대학교 산학협력단 | Method of sedimentary environment interpretation through electrofacies construction |
CN108802812A (en) * | 2017-04-28 | 2018-11-13 | 中国石油天然气股份有限公司 | Well-seismic fusion stratum lithology inversion method |
CN109541685A (en) * | 2017-09-22 | 2019-03-29 | 中国石油天然气股份有限公司 | River channel sand body identification method |
CN109932749A (en) * | 2019-03-05 | 2019-06-25 | 中国石油天然气股份有限公司 | Well-seismic calibration method and device |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7826973B2 (en) * | 2007-06-15 | 2010-11-02 | Chevron U.S.A. Inc. | Optimizing seismic processing and amplitude inversion utilizing statistical comparisons of seismic to well control data |
-
2019
- 2019-09-11 CN CN201910856307.4A patent/CN110579806B/en active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP2846175A2 (en) * | 2013-09-06 | 2015-03-11 | Services Petroliers Schlumberger | Seismic survey analysis |
CN105759312A (en) * | 2016-03-01 | 2016-07-13 | 中国海洋石油总公司 | Reservoir characteristic method well-seismic calibration method |
CN106324681A (en) * | 2016-08-24 | 2017-01-11 | 中国石油化工股份有限公司 | Time-frequency continuous wavelet transform well-to-seismic calibration method |
CN107783187A (en) * | 2016-08-24 | 2018-03-09 | 中国石油化工股份有限公司 | It is a kind of that logging speed and seismic velocity are combined to the method for establishing 3D velocity field |
CN108802812A (en) * | 2017-04-28 | 2018-11-13 | 中国石油天然气股份有限公司 | Well-seismic fusion stratum lithology inversion method |
CN107831542A (en) * | 2017-07-24 | 2018-03-23 | 中国石油化工股份有限公司 | DDW high accuracy depths domain well shake matching process |
KR101893800B1 (en) * | 2017-08-09 | 2018-09-04 | 제주대학교 산학협력단 | Method of sedimentary environment interpretation through electrofacies construction |
CN109541685A (en) * | 2017-09-22 | 2019-03-29 | 中国石油天然气股份有限公司 | River channel sand body identification method |
CN108254785A (en) * | 2017-12-28 | 2018-07-06 | 中国石油天然气股份有限公司 | Reservoir Determination Method and Device |
CN109932749A (en) * | 2019-03-05 | 2019-06-25 | 中国石油天然气股份有限公司 | Well-seismic calibration method and device |
Non-Patent Citations (2)
Title |
---|
密井网条件下井震联合低级序断层识别方法;张昕 等;《石油地球物理勘探》;20120630;第47卷(第3期);第462-468页 * |
油田开发后期井间砂体识别;徐立恒 等;《石油地球物理勘探》;20190430;第54卷(第2期);第390-397页 * |
Also Published As
Publication number | Publication date |
---|---|
CN110579806A (en) | 2019-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106353792B (en) | Method suitable for positioning micro-seismic source of hydraulic fracturing | |
CN106405651B (en) | Full waveform inversion initial velocity model construction method based on logging matching | |
CN107783187B (en) | Method for establishing three-dimensional velocity field by combining logging velocity and seismic velocity | |
CN109339778B (en) | Acoustic logging method for quantitatively evaluating perforation penetration depth | |
CN102943669B (en) | High-accuracy depth domain stochastic simulation inversion method suitable for designing track of horizontal well | |
CN107132578A (en) | A kind of microseism ground monitoring velocity model corrections algorithm | |
CN111257946B (en) | Geophysical drilling guiding method and method for updating stratum seismic velocity | |
CN104122581A (en) | Poststack acoustic wave impedance inversion method | |
CN110531410A (en) | A kind of least square reverse-time migration gradient Preconditioning method based on through wave field | |
CN110579806B (en) | Rapid well seismic calibration method under well pattern tight condition | |
CN104991272A (en) | Earthquake speed disturbance modeling method for well-free earthquake reversion | |
CN107506556B (en) | Method for determining sound wave longitudinal wave velocity value of fresh complete rock mass | |
CN111488328B (en) | Automatic drawing method for uranium ore geological drilling core geophysical prospecting and recording curve | |
CN101451432A (en) | High precision digital sonic wave variable density scale method | |
CN112558180A (en) | Method for rapidly detecting seismic horizon calibration accuracy by using horizontal isochronous surface | |
CN116203633A (en) | Tunnel surrounding rock integrity evaluation method based on natural source vibration imaging | |
CN103091728B (en) | The secondary sound wave analogy method of waterfrac treatment and induced fracture monitoring technique evaluating method | |
CN113658335B (en) | Method for creating complex geological model of karst development area | |
CN109343122A (en) | A kind of reservoir prediction sample extending method towards deep learning | |
CN114510808B (en) | Fine characterization method for numerical simulation of seam hole attribute parameters by using acid pressure of seam finding hole | |
CN109061737B (en) | Reservoir prediction method and device based on synthetic seismic record | |
CN110764136B (en) | Combined positioning method for time-lapse linear combination and nonlinear combination of anisotropic longitudinal and transverse waves | |
CN112882097A (en) | Calibration method for highly deviated well and horizontal well | |
CN2851573Y (en) | Impulse neutron logging tool | |
CN110609328A (en) | Oil saturation prediction method based on well-seismic combination |
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 | ||
TA01 | Transfer of patent application right |
Effective date of registration: 20211231 Address after: 163453 Heilongjiang Province, Daqing City Ranghulu District No. 233 South Central Avenue Applicant after: Daqing Oilfield Co.,Ltd. Applicant after: PetroChina Company Limited Address before: 163453 Heilongjiang Province, Daqing City Ranghulu District No. 233 South Central Avenue Applicant before: Daqing Oilfield Co.,Ltd. |
|
TA01 | Transfer of patent application right | ||
GR01 | Patent grant | ||
GR01 | Patent grant |