CN1130571C - Three-D earthquake prediction and observation system restoring technology - Google Patents
Three-D earthquake prediction and observation system restoring technology Download PDFInfo
- Publication number
- CN1130571C CN1130571C CN 01119969 CN01119969A CN1130571C CN 1130571 C CN1130571 C CN 1130571C CN 01119969 CN01119969 CN 01119969 CN 01119969 A CN01119969 A CN 01119969A CN 1130571 C CN1130571 C CN 1130571C
- Authority
- CN
- China
- Prior art keywords
- coordinate
- shot
- big gun
- arrival
- shot 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.)
- Expired - Lifetime
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The present invention relates to a seismic data processing method which is the restoring technology of a three-dimensional seismic prediction and observation system. The present invention determines X coordinates of a shot point, Y coordinates of the shot point and X coordinates of first channels of rows and columns receiving each shot for conventional three-dimensional seismic field collecting parameters, and other auxiliary parameters are directly determined by single shot recordation; once three key parameters of the X coordinates and the Y coordinates of the shot point and the X coordinate of the first channels of the rows and the columns are determined, the shot point, the arrangement position relation and the position of the shot point in the field can be determined, the three-dimensional original prediction can be restored and can be combined with the single shot recordation to obtain a corresponding observation system, and a position distribution picture of common reflection points can be further obtained; a position geographic coordinate of the common reflection points can be obtained by comparing a position distribution picture of whole zone common reflection points treated in the past, the restoration of the prediction and observation system can be realized, and a restoration results is completely correct after practically treated and validated.
Description
The invention belongs to the petroleum and natural gas field of seismic exploration, relate to the seismic data disposal route, is a kind of three-D earthquake prediction and observation system restoring technology.
The resulting source book of 3-D seismics field acquisition mainly comprises open-air grandfather tape and class's newspaper and the recording geometry (geographic coordinate of all common reflection point positions, exploratory area, the geographic coordinate of shot point and detection), these source book have been arranged, just can carry out seismic data and handle, explain, do structural map, the make location, drilling well and oil recovery, therefore these source book are the most basic raw data of whole seismic survey work, along with the fast development of seismic exploration technique, the 3-D seismics source book will be handled once again every 3-5; Because class's newspaper and recording geometry data volume are big, use frequently, thereby cause it to lose or lose phenomenon happening occasionally, and be inevitable.In the past,, to carry out the processing again of seismic data again, and just must carry out field design and construction again, and not only will spend great amount of manpower and material resources and fund, and will influence the process of seismic prospecting in case generation three-D earthquake prediction and recording geometry are lost.
The objective of the invention is to utilize the single big gun information that is write down on the three-dimensional open-air grandfather tape of only depositing, in conjunction with Seismic Data Processing Technique and handled achievement in the past, three-D earthquake prediction and the observation system restoring technology realized by Computer Processing.
The concrete technical scheme of the present invention is achieved in that
Adopt original field acquisition recording tape and combination processing achievement in the past, adopt following steps:
1. according to single shot record first arrival, arrangement number and reception channel number, determine the stereo observing system type that to recover;
2. on computer screen, show single shot record, according to the Y coordinate grouping of first arrival with every bunch shot point by big gun;
3. arrange at the adjacent nearly shot points of two of every big gun, read the first break time of the identical seismic trace of 5 groups of above X coordinates, utilize direct wave whilst on tour relational expression:
Wherein: h is a well depth, and v is an apparent velocity, and the two is constant, t
1, t
2Be respectively the first arrival take-off time of recording channel, x be along X to geophone offset, y be along Y to geophone offset;
Calculate the value of every group of x, y respectively, judge the actual minimum geophone offset of every big gun, and the Y coordinate of all big guns of same bunch is divided into groups again by reality with this;
4. according to shot point Y coordinate and X to minimum geophone offset, the X coordinate of all shot points is fixed on the same coordinate, field layout is carried out the definition first time, and from the first supreme first arrival apparent velocity of measuring this district of arbitrary normal single shot record, and with this speed the result after defining is carried out the linear moving school of first arrival and even up, obtain the linear moving school of first arrival section, be used for proofreading and correct the Y coordinate of shot point and X to minimum geophone offset;
5. calculate the amplitude root mean square sum curve of each seismic trace by big gun, contrast the amplitude root mean square sum curve of the same arrangement reception channel of adjacent big gun, obtain the relative position coordinates between each big gun arrangement;
6. contrast a seismic trace amplitude root mean square sum curve that receives arranging line that repeats between the adjacent wire harness, determine the relative position coordinates between each wire harness, finish the unified geometric definition of the whole district;
7. utilize new common reflection point location map to compare, obtain geographic coordinate and the shot point of the whole district and the geographic coordinate of geophone station of common reflection point position with handle the common reflection point location map that achievement obtains in the past.
Step shows single shot record by big gun in 2. on computer screen, be that per 100 big guns are selected 1 big gun and shown.
The present invention can directly recover original class newspaper and the recording geometry that institute loses or loses, and reverting to power is 100%, can avoid carrying out once more open-air earthquake construction and gather, and reduces man power and material's waste, guarantees carrying out smoothly of seismic data work of treatment.
Description of drawings of the present invention is as follows:
Fig. 1 is a process flow diagram of the present invention;
Fig. 2 is the three-dimensional original single big gun figure of the embodiment of the invention;
Fig. 3 is that conventional four lines six of the embodiment of the invention are concocted the recording geometry synoptic diagram;
Fig. 4 is that the open-air shot point of the embodiment of the invention excites synoptic diagram;
Fig. 5 is the single big gun figure of embodiment of the invention Y1 or Y2;
Fig. 6 is the single big gun figure of embodiment of the invention Y3;
Fig. 7 is the single big gun figure of embodiment of the invention Y4;
Fig. 8 is the single big gun figure of embodiment of the invention Y5 or Y6;
Fig. 9 is the moving school of a first arrival sectional view before embodiment of the invention Y coordinate is adjusted;
Figure 10 is the moving school of a first arrival sectional view after the adjustment of embodiment of the invention Y coordinate;
Figure 11 is the moving school of a first arrival sectional view before embodiment of the invention X coordinate is adjusted;
Figure 12 is the moving school of a first arrival sectional view after the adjustment of embodiment of the invention X coordinate;
Figure 13 is a seismic trace amplitude root mean square curve comparison diagram in the embodiment of the invention wire harness;
Figure 14 is a seismic trace amplitude root mean square curve comparison diagram between embodiment of the invention wire harness;
Figure 15 is that embodiment of the invention shot point, geophone station recover result map;
Figure 16 embodiment of the invention common reflection point recovers result map;
The three-dimensional class of Figure 17 embodiment of the invention newspaper recovers result map;
Figure 18 embodiment of the invention stereo observing system recovers result map.
The embodiment of the invention is as follows:
In 3-D seismics data processing procedure again, for various reasons, 3-D seismics original class newspaper and recording geometry are all lost, and have only original field to be with a cover, and the achievement data of in the past handling (achievement band and common reflection point location map), can not carry out thereby cause dealing with the work again.
Technology contents of the present invention and concrete implementation step are: for the parameter of conventional 3-D seismics field acquisition, determine that the X coordinate of shot point, Y coordinate and every big gun of shot point receive the first road X coordinate of arranging, these three parameters are the key parameters that constitute three-dimensional class newspaper and recording geometry, other auxiliary parameter is (as arranging number, the road number of each arrangement, the window's position etc. during split shooting) all can directly determine by single shot record.And shot point X, Y coordinate and arrange three key parameters of X coordinate in first road in case determine, the position relation and the residing in the open air position of shot point of shot point and arrangement can be determined, three-dimensional original class newspaper is restored, class's count off certificate has been arranged, obtain corresponding recording geometry in conjunction with single shot record again, and then obtain the common reflection point location map, by comparing with whole district's common reflection point location map of handling in the past, just can obtain the geographic coordinate of common reflection point position, realize recovery class's newspaper and recording geometry.
1. per 100 big guns of every bunch are selected a big gun list big gun and are shown on computer screen, first arrival, arrangement number and reception channel number according to single big gun, determine the stereo observing system type in this district, judge the stereo observing system (seeing accompanying drawing 3-accompanying drawing 4) that this district receives for conventional four lines, six processes of preparing Chinese medicine 240 roads by accompanying drawing 2.
2. show all single big guns of every bunch by big gun on computer screen, respectively arrange the difference of first break time according to single big gun, the Y coordinate of every bunch shot point is divided into four groups, i.e. the big gun of first order first arrival minimum is Y
1Or Y
2One group, the big gun of second order first arrival minimum is Y
3One group, the 3rd big gun of arranging the first arrival minimum is Y
4One group, the 4th big gun of arranging the first arrival minimum is Y
5Or Y
6One group (seeing accompanying drawing 5, accompanying drawing 6, accompanying drawing 7, accompanying drawing 8)
3. concoct the direct wave whilst on tour relational expression of recording geometry according to three-dimensional four lines six of routine:
Wherein h is a well depth, and v is an apparent velocity, and the two is constant, t
1, t
2Be respectively the first arrival take-off time of recording channel, x be along X to geophone offset, y be along Y to geophone offset.
Arrange at the adjacent nearly shot points of two of every big gun, read the first arrival take-off time (t of 10 groups of X coordinates identical (as the 1st road and the 61st road is a group, and the 5th road and the 65th road are a group etc.) seismic trace
1, t
2), separate top system of equations, obtain 10 groups of x, y, because the variation of direct wave speed, blow out the deviation of position and the error of first arrival take-off time, cause the on duty and different of these 10 groups of x, y, and according to the routine four lines six big gun pencil recording geometrys of field acquisition, its shot point Y coordinate and X distribute to offset distance all its fixing rule, i.e. Y
1, Y
2, Y
3, Y
4, Y
5, Y
6Value have only two kinds may, promptly can only be 50m or 150m, and X also can only be the integral multiple relation of 50m to offset distance, therefore after calculating 10 groups of x, y, it is approximate objectively to the 50m integral multiple that x, y are averaged respectively again, the X that so just can obtain every big gun is to the grouping Y of the minimum geophone offset of reality and every big gun coordinate, thereby the shot point Y coordinate of every bunch is sub-divided into six groups from four groups.
4. obtained every big gun shot point grouping Y coordinate and X to minimum geophone offset after, the X coordinate of all shot points is fixed on the same coordinate figure 100, three-dimensional field layout to all big guns carries out the definition first time, and from the first supreme direct wave apparent velocity (V=2000m/s) of measuring this district of arbitrary normal big gun record, result after using this speed to definition carries out the linear moving school of single big gun first arrival and evens up, obtain the linear moving school of first arrival section, be used for proofreading and correct the Y coordinate of shot point and X to minimum geophone offset, if the value of Y and X is accurate, the first arrival of six arrangements of every big gun should be evened up, and on the same standard time, otherwise one of following four kinds of situations will appear: 1) respectively arrange the first arrival section time and all reduce, and it is big that nearly arrangement reduces, and far away the arrangement lacking of minimizing illustrates that shot point Y coordinate should change Y2 into or change Y5 into from Y6 from Y1; 2) respectively arrange the first arrival section time and all increase, and closely arrange increase many, the far away arrangement lacking of increasing illustrates that shot point Y coordinate should change Y1 into or change Y6 into from Y5 from Y2; 3) respectively arrange the first arrival section time and all reduce, and the amount that reduces is the same, illustrates that X should reduce to geophone offset, and first break time whenever reduces 25ms, X should reduce 50m to geophone offset; 4) respectively arranges the first arrival section time and all increase, and the amount that increases is the same, illustrate X to geophone offset should increase, and the every increase of first break time 25ms, X should increase 50m to geophone offset.(seeing accompanying drawing 9, accompanying drawing 10, accompanying drawing 11, accompanying drawing 12)
Through after this step, just can determine exactly the Y coordinate of each big gun and X thereof to minimum geophone offset.
5. go out the amplitude root mean square sum curve of each record seismic trace on computers by the big gun statistical computation, contrast the same arrangement amplitude root mean square sum curve of the adjacent big gun of same bunch, obtain the relative position relation between each big gun arrangement, on X coordinate stuck-at-00 position with the first bombard of every bunch, just can determine the X coordinate of every big gun shot point successively and receive the Shou Dao position of arranging, as accompanying drawing 13, contrast the amplitude root mean square sum curve of the 10th big gun and the 11st big gun, it is 150m that the sp location that just can obtain the sp location of the 10 big gun and the 11st big gun differs 3 roads, and the rest may be inferred for each big gun contrast of other of same bunch.
6. obtained the X of single bunch, after the first road coordinate that Y coordinate and reception are arranged, utilize an arrangement that repeats between the adjacent wire harness to receive line again, statistical computation goes out the amplitude root mean square sum curve of the different big guns of same arrangement, contrast the amplitude root mean square sum curve of the identical big gun of adjacent wire harness number same arrangement, just can obtain the relative position relation between adjacent two bunch, as accompanying drawing 14, the amplitude root mean square sum curve that contrasts the 7th bunch and the 8th bunch first bombard as can be known, it is 150m that two bunch positions differ 3 roads, and the like, can obtain the relative position relation between all wire harness of the whole district, X with the whole district's the 1st bunch the 1st big gun, the Y coordinate is fixed on (X=50, Y=50) on, just can obtain all shot points of the whole district and receive the X that arranges, the Y coordinate, thus the unified field layout definition of the whole district finished.
7. by whole district field layout definition result after reunification, obtain shot point and the geophone station distribution plan and the common reflection point location map of the whole district on computers, with new shot point and geophone station distribution plan and common reflection point location map with handled resulting common reflection point location map contrast in the past, just can obtain the geographic coordinate of common reflection point position, define the geographic coordinate that just can release all shot points of the whole district and geophone station according to new field layout, realize recovery fully thus three-dimensional original class newspaper and recording geometry.As accompanying drawing 15 is exactly the old and new's shot point and the contrast of geophone station distribution plan, and accompanying drawing 16 is exactly the contrast of the old and new's common reflection point location map.
Use the present invention, successfully recovered above-mentioned area three-dimensional 14 bunch 8821 big guns, area is 100.26km
2Original class newspaper of 3-D seismics and recording geometry, through the actual treatment checking, restoration result is entirely true.
Claims (1)
1, a kind of three-D earthquake prediction and observation system restoring technology adopt original field acquisition recording tape and combination processing achievement in the past, it is characterized in that adopting following steps:
1. according to single shot record first arrival, arrangement number and reception channel number, determine the stereo observing system type that to recover;
2. pursue big gun and show single shot record on computer screen, according to the Y coordinate grouping of first arrival with every bunch shot point, the Y coordinate branch of every bunch shot point can be more than four groups or four groups, shows single shot record by big gun on computer screen, is that per 100 big guns are selected the demonstration of 1 big gun;
3. arrange at the adjacent nearly shot points of two of every big gun, read the first break time of the identical seismic trace of 5 groups of above X coordinates, utilize direct wave whilst on tour relational expression:
Wherein: h is a well depth, and v is an apparent velocity, and the two is constant, t
1, t
2Be respectively the first arrival take-off time of recording channel, x be along X to geophone offset, y be along Y to geophone offset;
Calculate the value of every group of x, y respectively, judge the actual minimum geophone offset of every big gun, and the Y coordinate of all big guns of same bunch is divided into groups again by reality with this;
4. according to shot point Y coordinate and X to minimum geophone offset, the X coordinate of all shot points is fixed on the same coordinate, field layout is carried out the definition first time, and from the first supreme first arrival apparent velocity of measuring this district of arbitrary normal single shot record, and with this speed the result after defining is carried out the linear moving school of first arrival and even up, obtain the linear moving school of first arrival section, be used for proofreading and correct the Y coordinate of shot point and X to minimum geophone offset;
5. calculate the amplitude root mean square sum curve of each seismic trace by big gun, contrast the amplitude root mean square sum curve of the same arrangement reception channel of adjacent big gun, obtain the relative position coordinates between each big gun arrangement;
6. contrast a seismic trace amplitude root mean square sum curve that receives arranging line that repeats between the adjacent wire harness, determine the relative position coordinates between each wire harness, finish the unified geometric definition of the whole district;
7. utilize new common reflection point location map to compare, obtain geographic coordinate and the shot point of the whole district and the geographic coordinate of geophone station of common reflection point position with handle the common reflection point location map that achievement obtains in the past.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 01119969 CN1130571C (en) | 2001-07-05 | 2001-07-05 | Three-D earthquake prediction and observation system restoring technology |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN 01119969 CN1130571C (en) | 2001-07-05 | 2001-07-05 | Three-D earthquake prediction and observation system restoring technology |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN1395114A CN1395114A (en) | 2003-02-05 |
| CN1130571C true CN1130571C (en) | 2003-12-10 |
Family
ID=4663833
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN 01119969 Expired - Lifetime CN1130571C (en) | 2001-07-05 | 2001-07-05 | Three-D earthquake prediction and observation system restoring technology |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN1130571C (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102062869B (en) * | 2009-11-17 | 2012-11-14 | 中国石油化工股份有限公司 | Method for designing seismic exploration and observation system |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US7561493B2 (en) * | 2003-05-30 | 2009-07-14 | Fairfield Industries, Inc. | Method and apparatus for land based seismic data acquisition |
| CN1310040C (en) * | 2004-06-25 | 2007-04-11 | 中国石油化工股份有限公司 | Correcting method and system for gun point or demodulator probe in earthquake prospection |
| CN102254321B (en) * | 2011-05-19 | 2014-03-26 | 中国石油集团川庆钻探工程有限公司 | Method for automatically identifying reversed polarity channel based on first-motion wave |
| CN111352152B (en) * | 2018-12-21 | 2022-11-04 | 中国石油天然气集团有限公司 | Rapid quality control method and device for seismic data observation system |
| CN111965729B (en) * | 2019-05-20 | 2023-04-07 | 中国石油天然气集团有限公司 | Real-time monitoring method, system and device for vibroseis combination center |
| CN112180428B (en) * | 2019-07-03 | 2023-12-26 | 中国石油天然气集团有限公司 | Push-pull observation system receiving relation generation method and device |
-
2001
- 2001-07-05 CN CN 01119969 patent/CN1130571C/en not_active Expired - Lifetime
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102062869B (en) * | 2009-11-17 | 2012-11-14 | 中国石油化工股份有限公司 | Method for designing seismic exploration and observation system |
Also Published As
| Publication number | Publication date |
|---|---|
| CN1395114A (en) | 2003-02-05 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU612552B2 (en) | Model-based depth processing of seismic data | |
| CN1170170C (en) | Spectral decomposition for seismic interpretation | |
| CN1313841C (en) | Layer stripping converted reflected waveforms for dipping fractures | |
| CN101031821A (en) | Offeset-azimuth binning for migration and velocity analysis | |
| CN102128627B (en) | Difference analysis method and device for electronic map | |
| CN110687602A (en) | Shallow seismic multi-wave combined exploration method | |
| CN1130571C (en) | Three-D earthquake prediction and observation system restoring technology | |
| CN110824563A (en) | Reservoir lithology prediction method based on Xgboost algorithm | |
| CN108828668A (en) | A kind of pre-stack time migration data processing method and device | |
| CN1130572C (en) | Method for sorting seismic data | |
| CN105652341B (en) | Rapid quantitative ancient landform restoration method considering ancient structure influence | |
| US20220221604A1 (en) | Method of analysing seismic data | |
| Burns | Analysis of geological events | |
| CN104267432B (en) | Converted wave common conversion point gather high-precision extracting method based on regularization | |
| McEvilly | The earthquake sequence of November 1964 near Corralitos, California | |
| US6249746B1 (en) | Automated seismic isochron analysis | |
| CN113534259B (en) | Real-time prestack time migration imaging method for efficient collection of controllable seismic source | |
| CA1240026A (en) | Method for interpretation of seismic records to yield indications of gaseous hydrocarbons | |
| McKoy et al. | Tight gas reservoir simulation: Modeling discrete irregular strata-bound fracture networks and network flow, including dynamic recharge from the matrix | |
| CN1481511A (en) | System and mtehod for analyzing and imaging enhanced three-D volume data set using one or more attributes | |
| US4063216A (en) | Displaying seismic sections in isometric view | |
| CN1015730B (en) | The method in computing formation dip and orientation | |
| CN111767662B (en) | Method and device for correcting azimuth information of single sand body of oil reservoir | |
| CN103278852A (en) | Method for forecasting oil gas by seismic data volume waveform structure characteristic model | |
| CN1310040C (en) | Correcting method and system for gun point or demodulator probe in earthquake prospection |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| C14 | Grant of patent or utility model | ||
| GR01 | Patent grant | ||
| CX01 | Expiry of patent term |
Granted publication date: 20031210 |
|
| CX01 | Expiry of patent term |