CN109164487A - A kind of method based on model foundation average velocity field and finely at the method for figure - Google Patents
A kind of method based on model foundation average velocity field and finely at the method for figure Download PDFInfo
- Publication number
- CN109164487A CN109164487A CN201811045739.9A CN201811045739A CN109164487A CN 109164487 A CN109164487 A CN 109164487A CN 201811045739 A CN201811045739 A CN 201811045739A CN 109164487 A CN109164487 A CN 109164487A
- Authority
- CN
- China
- Prior art keywords
- speed
- velocity
- velocity field
- volcanic
- well point
- 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 claims abstract description 43
- 239000010410 layer Substances 0.000 claims abstract description 67
- 239000011435 rock Substances 0.000 claims abstract description 63
- 239000011229 interlayer Substances 0.000 claims abstract description 47
- 230000000977 initiatory effect Effects 0.000 claims abstract description 27
- 238000012937 correction Methods 0.000 claims abstract description 12
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 230000015572 biosynthetic process Effects 0.000 claims description 16
- 238000012546 transfer Methods 0.000 claims description 7
- 239000000284 extract Substances 0.000 claims description 6
- 238000003786 synthesis reaction Methods 0.000 claims description 4
- 238000013507 mapping Methods 0.000 abstract description 6
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 238000004088 simulation Methods 0.000 description 9
- 238000010276 construction Methods 0.000 description 8
- 238000001228 spectrum Methods 0.000 description 8
- 238000011160 research Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 4
- 238000005553 drilling Methods 0.000 description 3
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 3
- 208000035126 Facies Diseases 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 241001163600 Bathylaco nigricans Species 0.000 description 1
- 241001074085 Scophthalmus aquosus Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000011230 binding agent Substances 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 239000004575 stone Substances 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
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/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/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
- G01V1/50—Analysing data
-
- 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
- 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 present invention relates to geological exploration fields, and in particular to a kind of method based on model foundation average velocity field and finely at the method for figure.The present invention obtains initiation layer velocity field in conjunction with horizon calibration by establishing Stratum of Volcanic Rocks rate pattern and carrying out forward modeling；Well point interlayer speed is obtained according to pit shaft analysis, well point velocity information and the forward modeling, is average velocity field using well point interlayer velocity correction initiation layer velocity field, and by the interlayer rate conversion after correction.The present invention average velocity field that interval velocity model is calculated in combination is more accurate, and the fine figure in obtained Search for Ore Deposits in Volcanic Areas, which is able to reflect down, covers the true structural configuration in stratum, improves the precision of variable speed mapping, eliminates and constructs illusion caused by high speed volcanic rock body.
Description
Technical field
The present invention relates to geological exploration fields, and in particular to a kind of method based on model foundation average velocity field and fine
At the method for figure.
Background technique
With going deep into for oil-gas exploration work, subsurface investigation object becomes increasingly complex.Especially in Search for Ore Deposits in Volcanic Areas, by
In being influenced to cause formation velocity vertically and horizontally to change greatly by volcanic rock, after volcanic rock is developed in superstratum, by high speed volcanic rock
Formation influence is caused to cover seismic reflection lineups under volcanic rock and be arched upward upwards, and the structural configuration shown on seismic profile is not
Really.After the seismic data interpretation to time-domain, only in the case where establishing accurate velocity field, time-depth conversion using variant velocity ability is carried out
Obtain the structural map of reflection subsurface structure form.During time and depth transfer, rationally accurately velocity field is to guarantee structural map essence
True important foundation.Conventional speeds field modeling method mainly has individual well constant speed method, more well interpolation methods, normal-moveout spectrum correction method etc..This
The method for establishing velocity field a bit not can be carried out accurate time and depth transfer, it is difficult to eliminate earthquake and cut open to the Volcanic area being unevenly distributed
The illusion in face.Currently, how to eliminate the high speed of volcanic rock in tectonic ore-forming influences to turn deeply when carries out accurate to interpretation results
It changes, is an extremely important technical problem, there are no highly developed methods so far.
Summary of the invention
The object of the present invention is to provide a kind of method based on model foundation average velocity field and finely at the method for figure, use
Inaccurately lead to structural map and physical presence error with the velocity field for solving the problems, such as that the prior art is established.
To achieve the above object, the present invention provides a kind of sides based on model foundation Search for Ore Deposits in Volcanic Areas average velocity field
Method, comprising the following steps:
Building Stratum of Volcanic Rocks rate pattern simultaneously carries out forward modeling；
Seismic data layer position is carried out according to horizon calibration and the forward modeling to explain, obtains initiation layer velocity field；
Well point interlayer speed is obtained according to pit shaft analysis, well point velocity information and the forward modeling；
Interlayer speed after being corrected is corrected to the initiation layer velocity field using the well point interlayer speed；
Average speed is calculated according to interlayer speed after the correction, average speed is established according to the average speed
?.
Further, Stratum of Volcanic Rocks rate pattern is established according to the layering of well, lithology, speed data in wellblock.
Further, the process for obtaining initiation layer velocity field further include:
Under horizon calibration and forward modeling guidance, the high speed reflection circle of volcanic rock is found out on each borehole logging tool curve
Face, Seismology and Geology mark timing determination volcanic rock reflecting interface；The bottom interface pair of top bottom surface and other stratum to Stratum of Volcanic Rocks
The reflecting interface answered explains, and constructs the top surface structural configuration of target zone；
The root mean sequare velocity data that seismic data of well is handled is arranged, exceptional value is removed, in conjunction with earthquake solution
The control of layer position is released, the root mean sequare velocity between each interval is extracted from mean-square-root velocity field, the root mean sequare velocity of each layer is converted to
Interval velocity obtains the initiation layer velocity field.
Further, the process for obtaining well point interlayer speed includes:
Make synthetic seismogram using the sound wave and density curve of well, find out the speed of well point, according to VSP at well point,
The when depth relationship of synthetic seismogram seeks the well point interlayer speed of corresponding each earthquake interval at well point.
The present invention also provides a kind of fine drawing methods in the Search for Ore Deposits in Volcanic Areas based on model, comprising the following steps:
Building Stratum of Volcanic Rocks rate pattern simultaneously carries out forward modeling；
Seismic data layer position is carried out according to horizon calibration and the forward modeling to explain, obtains initiation layer velocity field；
Well point interlayer speed is obtained according to pit shaft analysis, well point velocity information and the forward modeling；
Interlayer speed after being corrected is corrected to the initiation layer velocity field using the well point interlayer speed；
Average speed is calculated according to interlayer speed after the correction, average velocity field is established according to the average speed
And time and depth transfer is carried out, then it is patterned.
Further, Stratum of Volcanic Rocks rate pattern is established according to the layering of well, lithology, speed data in wellblock.
Further, the process for obtaining initiation layer velocity field further include:
Under horizon calibration and forward modeling guidance, the high speed reflection circle of volcanic rock is found out on each borehole logging tool curve
Face, Seismology and Geology mark timing determination volcanic rock reflecting interface；The bottom interface pair of top bottom surface and other stratum to Stratum of Volcanic Rocks
The reflecting interface answered explains, and constructs the top surface structural configuration of target zone；
The root mean sequare velocity data that the seismic data of well is handled is arranged, removes exceptional value, in combination
Interpretation horizon control is shaken, extracts the root mean sequare velocity between each interval from mean-square-root velocity field, the root mean sequare velocity of each layer is turned
It is changed to interval velocity, obtains the initiation layer velocity field.
Further, the process for obtaining well point interlayer speed includes:
Make synthetic seismogram using the sound wave and density curve of well, find out the speed of well point, according to VSP at well point,
The when depth relationship of synthetic seismogram seeks the well point interlayer speed of corresponding each earthquake interval at well point.
The beneficial effects of the present invention are: by establishing Stratum of Volcanic Rocks rate pattern and carrying out forward modeling, binder course position
Calibration obtains initiation layer velocity field；Well point interlayer is obtained according to pit shaft analysis, well point velocity information and the forward modeling
Speed is average velocity field using well point interlayer velocity correction initiation layer velocity field, and by the interlayer rate conversion after correction.This
The invention average velocity field that interval velocity model is calculated in combination is more accurate, and the fine figure in obtained Search for Ore Deposits in Volcanic Areas can
The true structural configuration in stratum is covered under reflection, improves the precision of variable speed mapping, eliminates construction caused by high speed volcanic rock body
Illusion.
Detailed description of the invention
Fig. 1 is the fine drawing methods flow chart in Search for Ore Deposits in Volcanic Areas based on model；
Fig. 2 is Depth Domain geological model；
Fig. 3 is time-domain forward simulation section；
Fig. 4 was volcano rock pillar mean velocity profile；
Fig. 5 was volcano rock pillar seismic profile；
Fig. 6 is seismic data interpretation section；
Fig. 7 is speed and depth curve at well point；
Fig. 8 is the field visualized figure of comprehensive speed；
Fig. 9 is the T0 such as the reflecting layer Tx figure；
Figure 10 is the reflecting layer Tx Depth Domain structural map.
Specific embodiment
The present invention will be further described in detail with reference to the accompanying drawing.
Stratum real construction form is covered in order to get clear down, using forward modeling technique verifying high speed volcanic rock to interval of interest
It is influenced caused by reflection line-ups；The normal-moveout spectrum velocity information obtained using pit shaft velocity information, seism processing, extract layer
Speed.Interval velocity is composed using well point interval velocity correction rate, and the interval velocity after correction is converted into average velocity field.Actually answer
With the precision shown the method increase variable speed mapping, eliminates and construct illusion caused by high speed volcanic rock body.
One, basic principle
1. rate pattern forward modeling
In Search for Ore Deposits in Volcanic Areas, influence caused by covering stratum reflecting interface lineups to reflect high speed stratum under is used
Forward technology, to study the layering of area's well, lithology, speed data establish Search for Ore Deposits in Volcanic Areas geological model.Foundation well brill chance (or
Person's periphery is appeared) thickness determine the thickness and other formation thicknesses of volcanic rock, filling speed is by well point in Stratigraphic framework
Sound wave curve calculating acquires.By the time section after forward modelling, analyze high speed volcanic rock following stratum can be generated it is bright
Aobvious lineups pull up phenomenon.
2. volcano rock stratum position control methods establish rate pattern
It explains 2.1 seismic data layer position
Under horizon calibration and forward modeling guidance, the high speed reflecting interface of volcanic rock is found out on each borehole logging tool curve,
The quasi- volcanic rock reflecting interface of Seismology and Geology mark timing card.Seismic data is explained, the top surface structural configuration of target zone is constructed.
2.2 use the layer position of seismic data interpretation as the Stratigraphic framework of model, acquire stratum with sound wave curve in research area
Filling speed in screen work obtains time section by forward simulation, analyzes the phase of forward simulation section and actual seismic section
Like property, if similitude difference illustrates that the frame model on stratum is not fine, that is to say, that the layer position of seismic data interpretation it is not fine (or
It is unreasonable), need return to reinterpret, until rationally.
The generation of 2.3 interval velocity fields
2.3.1 establishing initiation layer velocity field
The root mean sequare velocity data on stratum is arranged first, goes exceptional value, with the control of earthquake interpretation horizon from just
Root velocity field extracts the root mean sequare velocity between each interval, and each layer root mean sequare velocity is converted to interval velocity by formula (1), is obtained
To initiation layer velocity field.
In formula, vr,nFor the root mean sequare velocity of n-th layer；vr,n-1For (n-1)th layer of root mean sequare velocity；VnFor the layer speed of n-th layer
Degree；tnFor the reflection interval of n-th layer.
2.3.2 the acquisition of well point interlayer speed
Fine synthesis earthquake record is made using the sound wave and density curve of well, according to VSP, synthetic seismogram at well point
Depth relationship seeks the formation velocity of corresponding each earthquake interval at well point whens equal.
2.3.3 initial interval velocity is corrected using well point interlayer speed, obtains accurate interlayer speed.
Proportionality coefficient is obtained divided by the initial interval velocity of point with well point interval velocity, then two-dimensional mesh is carried out to the proportionality coefficient
It formats, forms proportionality coefficient trend surface, then coefficient trend surface is multiplied with initial interval velocity, the interlayer speed after being corrected
Degree.
3. seeking average velocity field
After obtaining accurate interlayer speed, every aspect can be sought with the average speed of overlying strata using formula (2).
In formula, VAv, nFor earth's surface to the average speed between n-th layer, viFor i-th layer of interval velocity, tiFor i-th layer of ground
Shake reflection interval.
Average velocity field is established, velocity field combination logging speed information and speed spectrum information make full use of well logging vertical
To velocity accuracy height and the accurate advantage of normal-moveout spectrum lateral velocity trend.
4. time and depth transfer
It is multiplied using the speed of velocity field with the reflection interval of target zone, obtains construction depth value.Utilize construction depth value
Make structural map.
Two, implementation steps
Drilling well layering is compared, horizon calibration and on the basis of the Fine structural interpretation of layer position, be guidance with forward simulation,
Under target zone control, initiation layer velocity field is calculated by the normal-moveout spectrum that seism processing obtains, utilizes the well logging of pit shaft
Speed calculates interval velocity at well point, corrects initial interval velocity using well point interval velocity, and the interval velocity after correction is converted to flat
Equal velocity field.Implementation step is shown in Fig. 1.
(1) research area has well Strata Comparison, determines the formation thickness of Stratum of Volcanic Rocks thickness and other intervals.
(2) speed on Stratum of Volcanic Rocks speed and other stratum is calculated by the sound wave curve of well point.
(3) to study the layering of area's well, lithology, speed data establish Search for Ore Deposits in Volcanic Areas geological model.It bores and meets according to well
Thickness determine the thickness and other formation thicknesses of volcanic rock, filling speed is calculated by the sound wave curve of well point in Stratigraphic framework
It acquires.By the time section after forward modelling, the volcanic rock for analyzing high speed can generate on apparent lineups following stratum
Draw phenomenon.
(4) seismic data layer position is explained
Under horizon calibration and forward modeling guidance, the high speed reflecting interface of volcanic rock is found out on each borehole logging tool curve,
The quasi- volcanic rock reflecting interface of Seismology and Geology mark timing card.The bottom interface of top bottom surface and other stratum to Stratum of Volcanic Rocks is corresponding
Reflecting interface explains, and constructs the top surface structural configuration of target zone.
(5) initiation layer velocity field is established
The root mean sequare velocity data obtained to research area's seism processing arranges, and goes exceptional value, uses seismic interpretation
The control of layer position, extracts the root mean sequare velocity between each interval from mean-square-root velocity field, by formula (1) each layer root mean sequare velocity
Interval velocity is converted to, initiation layer velocity field (normal-moveout spectrum interval velocity) is obtained.
(6) acquisition of well point interlayer speed
Fine synthesis earthquake record is made using the sound wave and density curve of well, according to VSP, synthetic seismogram at well point
Depth relationship seeks the formation velocity (well point interlayer speed) of corresponding each earthquake interval at well point whens equal.
(7) initial interval velocity is corrected using well point interlayer speed, obtains accurate interlayer speed.
Proportionality coefficient is obtained divided by the initial interval velocity of point with well point interval velocity, then two-dimensional mesh is carried out to the proportionality coefficient
It formats, forms proportionality coefficient trend surface, then coefficient trend surface is multiplied with initial interval velocity, the interlayer speed after being corrected
Degree.
(8) average velocity field is sought
After obtaining accurate interlayer speed, every aspect can be sought with the average speed of overlying strata using formula (2).
(9) time and depth transfer
It is multiplied using the speed of velocity field with the reflection interval of target zone, obtains construction depth value.Utilize construction depth value
Make structural map.
Three, application example
Research work area area has more than 500 sq-kms to develop more volcanic rock, and the overflow facies volcanic rock of the same period time is not wide
It is general to be distributed within the scope of work area, the different volcanic rock high speed stratum of thickness is formed, there are more than ten mouthfuls of drilling well in area, part well, which bores, meets
Volcanic rock, part well is far from volcanic crater.Drilling well discloses the feature that stratum longitudinal velocity differs greatly, and volcanic vent and thickness are not
Equal volcanic rock cap covers stratigraphic structure under and produces very big influence.Volcanic rock cap and each phase subintrusion phase, the overflow facies Black Warrior
Rock, tufa stone formation velocity are apparently higher than country rock formation velocity up and down, and in same stratum, due to apart from volcanic crater distance
Horizontally there is also certain variations for difference, formation thickness and formation interval velocity.It is anti-that stratum is covered in order to eliminate volcanic rock under
Influence caused by structure is penetrated, to get down the structural configuration for covering stratum clear, provides for well site deployment according to accurate structural map, utilizes
Purpose of horizon control methods establish high-precision average velocity field, carry out variable speed mapping, improve mapping precision.
1. the foundation of geological model
In order to confirm the presence of this pull-up effect, Stratum of Volcanic Rocks rate pattern is made based on well, such as Fig. 2 institute
Show.W1 and W2 well, which bore, in the model meets the different volcanic rock of thickness, and Stratum of Volcanic Rocks thickness is layered according to actual well drilled at well point
, upper and lower stratum is assumed to horizontal layered earth, and the sound wave curve at each stratum filling speed well point is calculated, for letter
Change model, the interlayer speed on each stratum takes the average formation velocity of two mouthfuls of wells.
2. forward modeling
As shown in figure 3, found out by the time section after forward modelling, in the case where not considering prime stratum inclination angle,
The volcanic rock of high speed can generate apparent lineups to following stratum and pull up phenomenon, T3, T4, T5 Bottom surfaces of strata sheet at W1 well
Being horizontal interface, and on the time section of forward modeling, the reflecting layer T3, T4, T5 has apparent pull-up, pulls up 10-20ms, far
At volcanic crater, as Stratum of Volcanic Rocks thickness gradually becomes balanced, obvious office is not caused to underlying strata reflection line-ups
Portion pulls up phenomenon, and forward modeling result and time earthquake profile morphology are more consistent, as shown in Figure 4 and Figure 5.
3. seismic data layer position is explained
Under horizon calibration and forward modeling guidance, the high speed reflecting interface of volcanic rock is found out on each borehole logging tool curve,
The quasi- volcanic rock reflecting interface of Seismology and Geology mark timing card.The bottom interface of top bottom surface and other stratum to Stratum of Volcanic Rocks is corresponding
Reflecting interface explains, and constructs the top surface structural configuration of target zone, as shown in Figure 6.
4. forward simulation corrects interpretation model
Using primary explanation layer position T1, T2, T3, T4 as model framework, the speed data of research on utilization area well is as mould
The interval velocity of type carries out forward simulation, the similitude of forward simulation section and earthquake real data is compared, if forward simulation cuts open
Face and actual seismic section similitude are poor, show that the layer position of seismic data interpretation is problematic, need again to explain layer position, adjust
Whole interpretation scheme, until forward simulation section is similar to actual seismic section.
5. establishing initiation layer velocity field
The root mean sequare velocity data obtained to research area's seism processing arranges, and goes exceptional value, uses seismic interpretation
The control of layer position, extracts the root mean sequare velocity between each interval from mean-square-root velocity field, by formula (1) each layer root mean sequare velocity
Interval velocity is converted to, initiation layer velocity field is obtained.
6. the acquisition of well point interlayer speed
The when depth relation curve at the well point is found out according to VSP, synthetic seismogram at well point, as shown in Figure 7.
According to the individual-layer data of the interval of interest of each well in the area, the interval velocity of well point is sought.Seek the layer speed of well point
The method of degree is as follows:
VnFor the interval velocity of target zone at well point；Hn、Hn+1The depth at the top bottom of target zone respectively at well point；an、an+1Point
It Wei not depth corresponding average speed in target zone top bottom at well point.
7. correcting initial interval velocity using well point interlayer speed, accurate interlayer speed is obtained.
Due to stratigraphic dip and resolution ratio, normal-moveout spectrum interval velocity usually has large error with practical interval velocity, benefit
The normal-moveout spectrum low disadvantage of precision in the longitudinal direction can be made up by being corrected with the interval velocity at well point to this.It is removed with well point interval velocity
With the initial interval velocity of point, proportionality coefficient is obtained, then two-dimensional mesh is carried out to the proportionality coefficient and is formatted, forms proportionality coefficient trend
Then coefficient trend surface is multiplied by face with initial interval velocity, the interlayer speed after being corrected.
8. seeking average velocity field
After obtaining accurate interlayer speed, every aspect can be sought with the average speed of overlying strata, shape using formula (2)
At average velocity field, as shown in Figure 8.
9. time and depth transfer
It is multiplied using the speed of velocity field with the reflection interval of target zone, obtains construction depth value.Utilize construction depth value
Make structural map.
Fig. 9 and Figure 10 are T0 figure and the variable speed structure figures such as the top volcanic rock cap underlying strata T3 interface respectively, wherein A well spacing from
Farther out, B well spacing is closer from volcanic crater at volcanic crater, waits on T0 figure, and B well ratio A well constructs high 20ms, and (the figure after variable speed mapping
10), B well ratio A well constructs low 20m, structural map and equal T0 figure contradiction, but consistent with geology comparison conclusion.Cause this contradiction
The reason of phenomenon be likely to be influenced apart from the closer B well in volcanic crater by volcanic rock cap high speed stratum it is bigger, in time section
On cause the pull-up of lineups, so as to form the high illusion of local structure on time section.The fire based on model is proved simultaneously
The fine drawing methods in the development area Shan Yan have preferable effect to the tectonic ore-forming for solving velocity anomaly geologic body.
Specific embodiment of the present invention is presented above, but the present invention is not limited to described embodiment,
Such as other establish the method for model or the simple transformation of speed calculation formula, the technical solution formed in this way is to above-mentioned
Embodiment, which is finely adjusted, to be formed, and this technical solution still falls in protection scope of the present invention.
Claims (8)
1. a kind of method based on model foundation Search for Ore Deposits in Volcanic Areas average velocity field, which comprises the following steps:
Building Stratum of Volcanic Rocks rate pattern simultaneously carries out forward modeling；
Seismic data layer position is carried out according to horizon calibration and the forward modeling to explain, obtains initiation layer velocity field；
Well point interlayer speed is obtained according to pit shaft analysis, well point velocity information and the forward modeling；
Interlayer speed after being corrected is corrected to the initiation layer velocity field using the well point interlayer speed；
Average speed is calculated according to interlayer speed after the correction, average velocity field is established according to the average speed.
2. a kind of method based on model foundation Search for Ore Deposits in Volcanic Areas average velocity field according to claim 1, feature
It is: Stratum of Volcanic Rocks rate pattern is established according to the layering of well, lithology, speed data in wellblock.
3. a kind of method based on model foundation Search for Ore Deposits in Volcanic Areas average velocity field according to claim 1 or 2, special
Sign is, the process for obtaining initiation layer velocity field further include:
Under horizon calibration and forward modeling guidance, the high speed reflecting interface of volcanic rock is found out on each borehole logging tool curve,
Seismology and Geology mark timing determination volcanic rock reflecting interface；The bottom interface of top bottom surface and other stratum to Stratum of Volcanic Rocks is corresponding
Reflecting interface explains, and constructs the top surface structural configuration of target zone；
The root mean sequare velocity data that seismic data of well is handled is arranged, exceptional value is removed, in conjunction with earthquake interpretation layer
Position control, extracts the root mean sequare velocity between each interval from mean-square-root velocity field, and the root mean sequare velocity of each layer is converted to layer speed
Degree, obtains the initiation layer velocity field.
4. a kind of method based on model foundation Search for Ore Deposits in Volcanic Areas average velocity field according to claim 3, feature
It is, the process for obtaining well point interlayer speed includes:
Synthetic seismogram is made using the sound wave and density curve of well, finds out the speed of well point, according to VSP, synthesis at well point
The when depth relationship of earthquake record seeks the well point interlayer speed of corresponding each earthquake interval at well point.
5. a kind of fine drawing methods in Search for Ore Deposits in Volcanic Areas based on model, which comprises the following steps:
Building Stratum of Volcanic Rocks rate pattern simultaneously carries out forward modeling；
Seismic data layer position is carried out according to horizon calibration and the forward modeling to explain, obtains initiation layer velocity field；
Well point interlayer speed is obtained according to pit shaft analysis, well point velocity information and the forward modeling；
Interlayer speed after being corrected is corrected to the initiation layer velocity field using the well point interlayer speed；
Average speed is calculated according to interlayer speed after the correction, average velocity field is established according to the average speed and is gone forward side by side
Row time and depth transfer, is then patterned.
6. the fine drawing methods in a kind of Search for Ore Deposits in Volcanic Areas based on model according to claim 5, it is characterised in that: root
Stratum of Volcanic Rocks rate pattern is established according to the layering of well, lithology, speed data in wellblock.
7. the fine drawing methods in a kind of Search for Ore Deposits in Volcanic Areas based on model according to claim 5 or 6, feature exist
In the process for obtaining initiation layer velocity field further include:
Under horizon calibration and forward modeling guidance, the high speed reflecting interface of volcanic rock is found out on each borehole logging tool curve,
Seismology and Geology mark timing determination volcanic rock reflecting interface；The bottom interface of top bottom surface and other stratum to Stratum of Volcanic Rocks is corresponding
Reflecting interface explains, and constructs the top surface structural configuration of target zone；
The root mean sequare velocity data that the seismic data of well is handled is arranged, exceptional value is removed, in conjunction with earthquake solution
The control of layer position is released, the root mean sequare velocity between each interval is extracted from mean-square-root velocity field, the root mean sequare velocity of each layer is converted to
Interval velocity obtains the initiation layer velocity field.
8. the fine drawing methods in a kind of Search for Ore Deposits in Volcanic Areas based on model according to claim 7, which is characterized in that institute
The process for stating to obtain well point interlayer speed includes:
Synthetic seismogram is made using the sound wave and density curve of well, finds out the speed of well point, according to VSP, synthesis at well point
The when depth relationship of earthquake record seeks the well point interlayer speed of corresponding each earthquake interval at well point.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811045739.9A CN109164487A (en) | 2018-09-07 | 2018-09-07 | A kind of method based on model foundation average velocity field and finely at the method for figure |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811045739.9A CN109164487A (en) | 2018-09-07 | 2018-09-07 | A kind of method based on model foundation average velocity field and finely at the method for figure |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN109164487A true CN109164487A (en) | 2019-01-08 |
Family
ID=64894403
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811045739.9A Pending CN109164487A (en) | 2018-09-07 | 2018-09-07 | A kind of method based on model foundation average velocity field and finely at the method for figure |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109164487A (en) |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109884700A (en) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | Multi-information fusion seismic velocity modeling method |
| CN111060969A (en) * | 2019-12-25 | 2020-04-24 | 恒泰艾普（北京）能源科技研究院有限公司 | Well control Q compensation method |
| CN111308549A (en) * | 2020-03-09 | 2020-06-19 | 中国石油化工股份有限公司 | Variable-speed mapping method based on model inversion |
| CN111815769A (en) * | 2019-04-12 | 2020-10-23 | 中国石油天然气股份有限公司 | Modeling method, computing device and storage medium for thrust-driven tectonic belt structure |
| CN112649893A (en) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | Thin reservoir oriented multi-data multi-parameter fusion modeling method and system |
| CN113658016A (en) * | 2021-08-17 | 2021-11-16 | 淮北矿业（集团）有限责任公司 | Fine time-depth conversion method and device applied to few-control-point horizon of coal mine |
| CN115453618A (en) * | 2022-09-05 | 2022-12-09 | 中海石油国际能源服务(北京)有限公司 | Method for establishing forward model and structural depth model of overburden stratum of paste-salt-containing layer |
| US11561313B2 (en) * | 2019-07-03 | 2023-01-24 | Saudi Arabian Oil Company | Creating seismic depth grids using horizontal wells |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106125133A (en) * | 2016-07-04 | 2016-11-16 | 中国海洋石油总公司 | A kind of based on the fine velocity modeling method under the constraint of gas cloud district |
| CN107728206A (en) * | 2017-09-14 | 2018-02-23 | 中国石油大学(华东) | A kind of velocity field modeling method |
-
2018
- 2018-09-07 CN CN201811045739.9A patent/CN109164487A/en active Pending
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106125133A (en) * | 2016-07-04 | 2016-11-16 | 中国海洋石油总公司 | A kind of based on the fine velocity modeling method under the constraint of gas cloud district |
| CN107728206A (en) * | 2017-09-14 | 2018-02-23 | 中国石油大学(华东) | A kind of velocity field modeling method |
Non-Patent Citations (2)
| Title |
|---|
| 朱广生 等编: "《勘探地震学教程》", 31 December 2005 * |
| 杜伟 等: "变速成图技术在K1地区的应用", 《新疆石油天然气》 * |
Cited By (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109884700A (en) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | Multi-information fusion seismic velocity modeling method |
| CN109884700B (en) * | 2019-03-20 | 2021-02-26 | 中国石油化工股份有限公司 | Multi-information fusion seismic velocity modeling method |
| CN111815769A (en) * | 2019-04-12 | 2020-10-23 | 中国石油天然气股份有限公司 | Modeling method, computing device and storage medium for thrust-driven tectonic belt structure |
| CN111815769B (en) * | 2019-04-12 | 2024-04-30 | 中国石油天然气股份有限公司 | Modeling method, computing device and storage medium for thrust covered zone construction |
| US11561313B2 (en) * | 2019-07-03 | 2023-01-24 | Saudi Arabian Oil Company | Creating seismic depth grids using horizontal wells |
| CN112649893A (en) * | 2019-10-10 | 2021-04-13 | 中国石油化工股份有限公司 | Thin reservoir oriented multi-data multi-parameter fusion modeling method and system |
| CN112649893B (en) * | 2019-10-10 | 2024-04-09 | 中国石油化工股份有限公司 | Multi-data and multi-parameter fusion modeling method and system for thin reservoir |
| CN111060969A (en) * | 2019-12-25 | 2020-04-24 | 恒泰艾普（北京）能源科技研究院有限公司 | Well control Q compensation method |
| CN111308549A (en) * | 2020-03-09 | 2020-06-19 | 中国石油化工股份有限公司 | Variable-speed mapping method based on model inversion |
| CN111308549B (en) * | 2020-03-09 | 2022-01-11 | 中国石油化工股份有限公司 | Variable-speed mapping method based on model inversion |
| CN113658016A (en) * | 2021-08-17 | 2021-11-16 | 淮北矿业（集团）有限责任公司 | Fine time-depth conversion method and device applied to few-control-point horizon of coal mine |
| CN115453618A (en) * | 2022-09-05 | 2022-12-09 | 中海石油国际能源服务(北京)有限公司 | Method for establishing forward model and structural depth model of overburden stratum of paste-salt-containing layer |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109164487A (en) | A kind of method based on model foundation average velocity field and finely at the method for figure | |
| US11313994B2 (en) | Geophysical deep learning | |
| CN106597533B (en) | A kind of Depth Domain velocity modeling method for mountain front seism processing | |
| CN105334535B (en) | A kind of method for identifying the thin hidden lithologic deposit of reservoir | |
| CN103527184B (en) | A kind of Forecasting Methodology of dolostone reservoirs and system | |
| US20170051598A1 (en) | System For Hydraulic Fracturing Design And Optimization In Naturally Fractured Reservoirs | |
| CN107121699A (en) | A kind of sedimentary facies identification method under earthquake phase control | |
| CN105510993A (en) | Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method | |
| CN109738945A (en) | A method of structural map is directly generated using pre-stack depth migration achievement | |
| CN102466815A (en) | Triassic clastic rock petroleum reservoir identification method | |
| CN104570067A (en) | Phase-controlled earthquake inversion method in geophysical exploration | |
| CN109839660A (en) | A method of velocity depth model is established using prestack trace gather data | |
| CN109541685A (en) | A kind of river channel sand recognition methods | |
| CN111983683B (en) | Prediction method and system for lake-facies limestone reservoir under low-well condition | |
| CN105842733A (en) | Shale reservoir earthquake identification method | |
| Khatiwada et al. | A window into the Proterozoic: Integrating 3D seismic, gravity, and magnetic data to image subbasement structures in the southeast Fort Worth basin | |
| CN106033127A (en) | Ground stress orientation seismic prediction method based on shear wave speed variation rate | |
| Nwaezeapu et al. | Sequence stratigraphic approach to hydrocarbon exploration: A case study of Chiadu field at eastern onshore Niger Delta basin, Nigeria | |
| Singleton | Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas | |
| CN109581499A (en) | A method of structural map is generated using anisotropy pre-stack depth migration | |
| CN109839659A (en) | The method that a kind of pair of pre-stack depth migration section is iterated optimization | |
| CN110941029B (en) | Speed modeling method related to geological capping | |
| Curkan | Reservoir characterization of channel-belt strata, McMurray Formation, northeastern Alberta | |
| Hillman | Integrating Full-Bore Formation Micro-Imager (FMI) Data for Niobrara Reservoir Characterization, Postle Area, Wattenberg Field, Colorado, USA | |
| Ahmed et al. | Seismic facies modelling of potwar basin using seismic and well log data |
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: 20190108 |
|
| RJ01 | Rejection of invention patent application after publication |