CN107422384B - A method of vertical seismic profile data is imitated using well-log information - Google Patents
A method of vertical seismic profile data is imitated using well-log information Download PDFInfo
- Publication number
- CN107422384B CN107422384B CN201710166760.3A CN201710166760A CN107422384B CN 107422384 B CN107422384 B CN 107422384B CN 201710166760 A CN201710166760 A CN 201710166760A CN 107422384 B CN107422384 B CN 107422384B
- Authority
- CN
- China
- Prior art keywords
- squared
- wave
- data
- well
- curve
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
- G01V1/44—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
- G01V1/48—Processing data
Abstract
The present invention relates to a kind of methods for imitating vertical seismic profile data using well-log information.This method comprises: 1) in Conventional Logs gamma ray curve and bulk density curve carry out environmental correction；2) lithology layering is carried out using the gamma ray curve after environmental correction；3) unified with nature gamma is layered result, electrode couple sound wave measuring well curve carries out squared processing, and utilize squared compressional wave time difference, squared shear wave slowness, cumulative time, interval velocity is calculated in bulk density after combining environmental correction, the layer imitative vertical seismic profile data such as wave impedance and reflection coefficient, provides more parameters for the processing and explanation of seismic data.This method is on the basis of being layered using natural gamma to stratum, electrode couple Sonic Logging Data is handled, numerous parameters such as cumulative time, interval velocity, layer wave impedance, reflection coefficient under identical lithology units are calculated, the processing and explanation that can be applied to seismic data, improve the precision of explanation results.
Description
Technical field
The invention belongs to petroleum geology exploration fields, and in particular to a kind of to imitate vertical seismic profile data using well-log information
Method.
Background technique
During oil-gas exploration and development, the spatial distribution description of sand body is most important with quantitative assessment.With to oil-containing
The identification of gas reservoir sand body and the raising for describing required precision, obtaining high accuracy three-dimensional seismic data seems most important.One
Since straight, the application of seismic data is based on qualitative description and analysis, and in processing and explanation, the processing of relationship is earthquake when deep
The key of datum processing.Relationship is to realize sound logging and ground when the direct depth that vertical seismic profiling (VSP) (VSP) data provides
The valuable data that face seismic data is correctly demarcated.Since Partial Block does not carry out VSP well logging, seismic data is handled at this time
There are bigger difficulty for explanation.
Conventional Logs are to survey in prospect pit well logging, evaluation borehole logging tool and exploitation borehole logging tool in oil-gas exploration and development at present
The well-log information measured, i.e., so-called " nine " logging program --- natural potential, natural gamma, three lithology curve of hole diameter, shallow,
In, deep three resistivity curves and interval transit time, neutron, density three porosity curve.
In recent years, can tentatively have been met using conventional logging interval travel time data to seismic data carry out processing explanation want
It asks.When conventional logging interval travel time data is applied to surface-seismic data, domestic Shaanxi Province extends Petroleum Group research institute
Guo Chao, burnt top-notch player's routine acoustic velocity logging and VSP logging speed carried out comparative study；The research of Sinopec petroleum exploration and development
The Zhang Weihong etc. of institute's Nanjing petroleum exploration research has also carried out the research of a large amount of theory and actual applications, proposes to utilize conventional logging
Data and seismic data joint application and lithologic interpretation；How the Chen Wenchao of Xi'an Communications University is to by sound logging and surface seismic
Data carries out matching and has carried out the research of practical application, proposes sound logging and surface-seismic data carrying out matched method.
But average speed and the interval velocity etc. that longitudinal wave can only be finally sought using conventional logging interval travel time data, cannot seek shear wave
Average speed, interval velocity and by p-and s-wave velocity data converse Lai the parameters such as speed ratio, Poisson's ratio, Young's modulus, often
Rule well logging sonic differential time data is that the processing of seismic data explains that the parameter provided is less, is had to the raising of earthquake data precision
Limit.
Summary of the invention
The purpose of the present invention is when lacking VSP data, provide a kind of imitative vertical seismic profile data of utilization well-log information
Method, provide more parameters for the processing and explanation of seismic data, and then improve the identification of oil-bearing reservoir sand body and description
Precision.
In order to achieve the goal above, the technical scheme adopted by the invention is that:
A method of vertical seismic profile data is imitated using well-log information, comprising the following steps:
1) in Conventional Logs gamma ray curve and bulk density curve carry out environmental correction；
2) lithology layering is carried out using the gamma ray curve after environmental correction, obtains natural gamma layering result；
Compressional wave time difference, shear wave slowness, P-S wave velocity ratio curve are obtained by dipole acoustic log data；When using longitudinal wave
Poisson's ratio, Young's modulus curve is calculated in bulk density after difference, the correction of shear wave slowness combining environmental；
3) layering of unified with nature gamma is as a result, to compressional wave time difference, shear wave slowness, P-S wave velocity ratio, Poisson's ratio, Young mould
Amount, bulk density curve carry out squared processing, obtain squared compressional wave time difference, squared shear wave slowness, squared
P-S wave velocity ratio, squared Poisson's ratio, squared Young's modulus, squared bulk density；
4) when accumulative under identical lithology units is calculated using squared compressional wave time difference, squared shear wave slowness
Between, interval velocity, by interval velocity combine squared bulk density layer wave impedance and reflection coefficient is calculated；It imitates to obtain the final product vertically
Shake sectional data.
In step 4), VSP data further can be imitated using gained and carry out the determination of seismic reflection layer position, target interval oil-containing
The prediction etc. of gas.For example, gained cumulative time, layer wave impedance and reflection coefficient to be applied to the processing and explanation of seismic data
In.The geologic horizon of seismic reflection layer is determined using interval velocity and reflection coefficient.Utilize squared compressional wave time difference and squared
Shear wave slowness and different lithology interval transit time standard comparing, the formation lithology variation near the integrated interpretation borehole wall.Utilize side
The oil-gas possibility of the when squared Poisson's ratio prediction target interval of the p-and s-wave velocity of wave.
The method provided by the invention for imitating vertical seismic profile data using well-log information is to utilize natural gamma over the ground
On the basis of layer is layered, electrode couple Sonic Logging Data is handled, when accumulative under identical lithology units is calculated
Between, numerous parameters such as interval velocity, layer wave impedance, reflection coefficient, so as to preferably be applied to the processing and solution of seismic data
It releases, improves the precision of explanation results.
The method provided by the invention for imitating vertical seismic profile data using well-log information, in the condition for lacking VSP data
Under, dipole acoustic log data and Conventional Logs are subjected to integrated treatment, provide numerous geologic parameters, it not only can be with
For the processing and explanation of seismic data, it can be further used for determining geology (epoch, lithology) layer position of seismic reflection layer, determine
Property distinguish the borehole wall near formation lithology variation and prediction target interval oil-gas possibility；This method is from the ground of strata division lithology
Matter rule is set out, and gained prediction conclusion is consistent with actual test and assay conclusion, and prospect of the application is good, and suitable popularization is answered
With.
Detailed description of the invention
Fig. 1 is that J1 borehole logging tool curve imitates VSP figure；
Fig. 2 is that Z1 borehole logging tool curve imitates VSP figure.
Specific embodiment
The present invention is further explained in the light of specific embodiments.
Embodiment 1
The method for imitating vertical seismic profile data using well-log information of the present embodiment is carried out for J1 well, including following
Step:
1) in Conventional Logs gamma ray curve and bulk density curve carry out environmental correction；
2) lithology layering is carried out using the gamma ray curve after environmental correction, obtains natural gamma layering result；
Compressional wave time difference, shear wave slowness, P-S wave velocity ratio curve are obtained by dipole acoustic log data；When using longitudinal wave
Poisson's ratio, Young's modulus curve is calculated in bulk density after difference, the correction of shear wave slowness combining environmental；
3) layering of unified with nature gamma is as a result, to compressional wave time difference, shear wave slowness, P-S wave velocity ratio, Poisson's ratio, Young mould
Amount, bulk density curve carry out squared processing, obtain squared compressional wave time difference, squared shear wave slowness, squared
P-S wave velocity ratio, squared Poisson's ratio, squared Young's modulus, squared bulk density；
4) when accumulative under identical lithology units is calculated using squared compressional wave time difference, squared shear wave slowness
Between, interval velocity, by interval velocity combine squared bulk density layer wave impedance and reflection coefficient is calculated；It imitates to obtain the final product vertically
Shake sectional data.
Fig. 1 is that the log of J1 well imitates VSP figure.It can intuitively find out from figure and be sought using natural gamma point lithology
Parameters, these parameters such as cumulative time, interval velocity, layer wave impedance, reflection coefficient can provide ginseng for the processing of seismic data
It examines；The geologic horizon of seismic reflection layer can be determined using interval velocity and reflection coefficient.Utilize squared compressional wave time difference and side
It is comprehensive can to change progress to the formation lithology near the borehole wall for the shear wave slowness of wave and the interval transit time standard comparing of different lithology
It closes and explains.
The oil-gas possibility of J1 well can be predicted using the when squared Poisson's ratio of squared p-and s-wave velocity.Fig. 1 is real
Stratum in color frame has lesser P-S wave velocity ratio and Poisson's ratio, and inferring has preferable gas-bearing property, the test knot of the well
By the deduction for also demonstrating gas-bearing property, the laminated examination daily gas of the well 3061.5-3064m, 3082-3085m two reach 25045 sides/
It.
Embodiment 2
The method for imitating vertical seismic profile data using well-log information of the present embodiment carries out, specific method for Z1 well
It is same as Example 1.
Fig. 2 is that Z1 borehole logging tool curve imitates VSP figure.Also there is relatively small P-S wave velocity ratio from the grey box of Fig. 2
And Poisson's ratio, infer that there is certain gas-bearing property, the live resolution gas of well assay from top to bottom be respectively 0.48 side/
Ton, 0.45 side/ton, 4.75 sides/ton, it was confirmed that above-mentioned deduction.
Claims (5)
1. a kind of method for imitating vertical seismic profile data using well-log information, which comprises the following steps:
1) in Conventional Logs gamma ray curve and bulk density curve carry out environmental correction；
2) lithology layering is carried out using the gamma ray curve after environmental correction, obtains natural gamma layering result；
Compressional wave time difference, shear wave slowness, P-S wave velocity ratio curve are obtained by dipole acoustic log data；Utilize compressional wave time difference, cross
Poisson's ratio, Young's modulus curve is calculated in bulk density after the correction of wave time difference combining environmental；
3) unified with nature gamma layering as a result, to compressional wave time difference, shear wave slowness, P-S wave velocity ratio, Poisson's ratio, Young's modulus,
Bulk density curve carries out squared processing, obtain squared compressional wave time difference, squared shear wave slowness, it is squared in length and breadth
Wave velocity than, squared Poisson's ratio, squared Young's modulus, squared bulk density；
4) be calculated using squared compressional wave time difference, squared shear wave slowness cumulative time under identical lithology units,
Interval velocity combines squared bulk density that layer wave impedance and reflection coefficient is calculated by interval velocity；Up to imitative vertical seismic
Sectional data.
2. the method for imitating vertical seismic profile data using well-log information as described in claim 1, which is characterized in that step 4)
In, gained cumulative time, layer wave impedance and reflection coefficient are applied in the processing and explanation of seismic data.
3. the method for imitating vertical seismic profile data using well-log information as described in claim 1, which is characterized in that step 4)
In, the geologic horizon of seismic reflection layer is determined using interval velocity and reflection coefficient.
4. the method for imitating vertical seismic profile data using well-log information as described in claim 1, which is characterized in that step 4)
In, using the interval transit time standard comparing of squared compressional wave time difference and squared shear wave slowness and different lithology, comprehensive solution
Release the formation lithology variation near the borehole wall.
5. the method for imitating vertical seismic profile data using well-log information as described in claim 1, which is characterized in that step 4)
In, utilize the oil-gas possibility of the when squared Poisson's ratio prediction target interval of squared p-and s-wave velocity.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710166760.3A CN107422384B (en) | 2017-03-20 | 2017-03-20 | A method of vertical seismic profile data is imitated using well-log information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710166760.3A CN107422384B (en) | 2017-03-20 | 2017-03-20 | A method of vertical seismic profile data is imitated using well-log information |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107422384A CN107422384A (en) | 2017-12-01 |
| CN107422384B true CN107422384B (en) | 2019-03-01 |
Family
ID=60423525
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710166760.3A Active CN107422384B (en) | 2017-03-20 | 2017-03-20 | A method of vertical seismic profile data is imitated using well-log information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107422384B (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108104805B (en) * | 2017-12-13 | 2021-02-05 | 西南石油大学 | Method for acquiring transverse wave time difference profile by using well logging and core data |
| CN111077589B (en) * | 2018-10-19 | 2022-06-03 | 中国石油天然气股份有限公司 | Method and device for determining longitudinal-transverse wave velocity ratio of dry rock |
| CN112485844B (en) * | 2019-09-12 | 2023-10-27 | 中国石油天然气股份有限公司 | Volcanic lithology lithofacies distribution prediction method and prediction device |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5654938A (en) * | 1996-05-31 | 1997-08-05 | Western Atlas International, Inc. | Method for identifying alteration of earth formations using dipole acoustic logging |
| CN102650701B (en) * | 2012-04-27 | 2014-02-05 | 中国石油天然气股份有限公司 | Reservoir wave impedance prediction method based on dipole wave |
| CN104656142A (en) * | 2013-11-19 | 2015-05-27 | 中国石油天然气集团公司 | Seismic horizon calibration method utilizing vertical seismic profiling (VSP) and well-logging combination |
| CN105370270A (en) * | 2015-11-06 | 2016-03-02 | 中石化石油工程技术服务有限公司 | Method for determining gas saturation of shale gas reservoir by longitudinal-transverse wave time difference of dipole acoustic waves |
-
2017
- 2017-03-20 CN CN201710166760.3A patent/CN107422384B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US5654938A (en) * | 1996-05-31 | 1997-08-05 | Western Atlas International, Inc. | Method for identifying alteration of earth formations using dipole acoustic logging |
| CN102650701B (en) * | 2012-04-27 | 2014-02-05 | 中国石油天然气股份有限公司 | Reservoir wave impedance prediction method based on dipole wave |
| CN104656142A (en) * | 2013-11-19 | 2015-05-27 | 中国石油天然气集团公司 | Seismic horizon calibration method utilizing vertical seismic profiling (VSP) and well-logging combination |
| CN105370270A (en) * | 2015-11-06 | 2016-03-02 | 中石化石油工程技术服务有限公司 | Method for determining gas saturation of shale gas reservoir by longitudinal-transverse wave time difference of dipole acoustic waves |
Non-Patent Citations (1)
| Title |
|---|
| 垂直地震剖面与地面地震联合成像方法;郭旭光 等;《新疆石油地质》;20070831;第28卷(第4期);第507-509页 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107422384A (en) | 2017-12-01 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN107121699B (en) | A kind of sedimentary facies identification method under earthquake phase control | |
| Willis et al. | Spatial orientation and distribution of reservoir fractures from scattered seismic energy | |
| US6928367B2 (en) | Reservoir fracture characterization | |
| US6941227B2 (en) | Frequency-dependent processing and interpretation (FDPI) of seismic data for identifying, imaging and monitoring fluid-saturated underground reservoirs | |
| US20100232259A1 (en) | Fracture clusters identification | |
| CN104898161B (en) | Effective sandstone predicting method based on logging response simulator | |
| CN108363100A (en) | Based on the coal bed gas dessert area Seismology and Geology recognition methods for sealing unit and rough set up for safekeeping | |
| CN103116184B (en) | A kind of method calculating shear wave statics in near-surface investigation | |
| WO2017035104A1 (en) | Velocity model seismic static correction | |
| Yang et al. | Integrated application of 3D seismic and microseismic data in the development of tight gas reservoirs | |
| CN104516018A (en) | Porosity inversion method under lithological constraint in geophysical exploration | |
| CN105093306A (en) | Method for automatic interpretation and thickness calculation of reservoir in geophysical exploration | |
| CN104005760A (en) | Azimuthal anisotropic elastic impedance based crack detection method | |
| CN107422384B (en) | A method of vertical seismic profile data is imitated using well-log information | |
| CN104975851B (en) | For amplitude with the reservoir model optimization method of geophone offset variation road set analysis | |
| CN107728205B (en) | A kind of Formation pressure prediction method | |
| CN101315428B (en) | High resolution prospecting method of large-ventage clastic rock deposition stratum based on horizontal wave velocity | |
| CN106501872B (en) | A kind of calculation method and device of fracture reservoir ground stress characteristics | |
| Wang et al. | Multi-scale natural fracture prediction in continental shale oil reservoirs: a case study of the Fengcheng Formation in the Mahu Sag, Junggar Basin, China | |
| CN104280773B (en) | Using the time-frequency spectrum changed with geophone offset cross figure predict thickness of thin layer method | |
| CN107132574A (en) | A kind of Forecasting Methodology of marine bed clastic rock lithology combination | |
| Parra et al. | Detection of guided waves between gas wells for reservoir characterization | |
| CN109283577A (en) | A kind of seismic layer labeling method | |
| CN113514884A (en) | Compact sandstone reservoir prediction method | |
| Naji et al. | Prediction of sonic shear wave using artificial neural network |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20220207 Address after: 100028 Chaoyang District, Beijing Hui Xin Street 6, Twelfth level. Patentee after: SINOPEC OILFIELD SERVICE Corp. Patentee after: Sinopec Jingwei Co.,Ltd. Patentee after: SINOPEC NORTH CHINA PETROLEUM ENGINEERING Co.,Ltd. Patentee after: North China Measurement and control company of Sinopec Jingwei Co.,Ltd. Address before: 100728 9 Jishi Kou Road, Chaoyang District, Beijing Patentee before: SINOPEC OILFIELD SERVICE Corp. Patentee before: Sinopec North China Petroleum Engineering Co., Ltd |
|
| TR01 | Transfer of patent right |