CN106094035B - A kind of residual static correction computational methods quickly handled based on unit - Google Patents
A kind of residual static correction computational methods quickly handled based on unit Download PDFInfo
- Publication number
- CN106094035B CN106094035B CN201610682656.5A CN201610682656A CN106094035B CN 106094035 B CN106094035 B CN 106094035B CN 201610682656 A CN201610682656 A CN 201610682656A CN 106094035 B CN106094035 B CN 106094035B
- Authority
- CN
- China
- Prior art keywords
- group
- static correction
- residual static
- subdatas
- value
- 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/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides a kind of residual static correction computational methods quickly handled based on unit.Methods described includes step：Adaptive grouping is carried out to 3D seismic data, obtains n group subdatas；Complete the calculating to the non-linear Value of residual static correction of every group of subdata in n group subdatas；The Value of residual static correction that every group of subdata is calculated in splicing n group subdatas.The present invention fast can calculate non-linear residual statics using unit, can abundant appliance computer memory size, solve the problems, such as that unit processing big data efficiency is low, and splice the lap for the static correction value that every group of data calculate using segmentation interception, meet the white noise characteristic of Value of residual static correction random distribution, the negative effect of splicing tape can be reduced as far as possible, it is ensured that stacking image effect.
Description
Technical field
The present invention relates to seismic prospecting data processing technology field, it relates in particular to a kind of unit quickly handle it is surplus
Remaining static correction computational methods, can adaptively splice Value of residual static correction, applied to oil seismic exploration, can especially be used as field
The effective ways of non-linear Value of residual static correction are quickly asked in in-situ processing based on unit.
Background technology
In recent years, petroleum exploration domain gradually extends to complex area, and these Earthquakes surface conditions are complicated, earthquake money
Expect that signal to noise ratio is low, residual static correction problem is serious, only relies on conventional linear reflective ripple residual static correction and can not obtain and makes us full
The imaging effect of meaning.Non-linear residual static correction mainly has simulated annealing, genetic algorithm etc., possesses more quiet than the linear residue of routine
Correct the stronger ability for solving the problems, such as low SNR data residual static correction.
Computational efficiency is a bottleneck of non-linear residual static correction technology, and it is related to the anti-of prestack common midpoint gather
Re-reading to write, cross-correlation function is repeatedly asked for.Field condition processing uses unit operation, and calculator memory is usually no more than
100G, in face of hundreds of G, or even upper T three-dimensionals trace gather data, if dynamic correction data is placed on into hard disk, reading will expend greatly every time
Measure the time；If dynamic correction data is placed on into internal memory can save read access time, but at present the internal memory of computer can not store it is big
Geological data.Conventional way is calculated using geological data packet, then splicing Value of residual static correction manually, it has not been convenient to Ye Bu sections
Learn.
The content of the invention
It is an object of the invention to solve at least one in deficiencies of the prior art.
For example, an object of the present invention is that providing a kind of more reasonable and science unit quickly calculates non-linear remain
Yu Jing schools amount method, method of the invention be particularly suitable for field condition processing when quickly asked for using unit operation it is non-linear
Residual statics, so as to fully apply single calculator memory size, solve three-dimensional big data Value of residual static correction and quickly ask
Take problem.
The invention provides a kind of residual static correction computational methods quickly handled based on unit.Methods described includes step
Suddenly：Adaptive grouping is carried out to 3D seismic data, obtains n group subdatas；Complete to every group of subdata in n group subdatas
The calculating of non-linear Value of residual static correction；The Value of residual static correction that every group of subdata is calculated in splicing n group subdatas, its
In, include the step of the adaptive grouping：3D seismic data is read, obtains size, the first master of the 3D seismic data
Survey line, main profile total number, statistic computer physical memory size；By main profile bar number since the first main profile it is right in order
The 3D seismic data is grouped, and the main profile bar number of preceding n-1 distribution is k, the main profile bar number of last time distribution
For m, wherein, k=inlineCnt/n+size × β and rounded using a method is entered, m=inlineCnt-(n-1) × k, n=(1+
β) × size/ (phym × α) and use are entered a method and rounded, wherein, β is the coincidence factor of subdata group adjacent to each other, and is met
size×β<Phym × (1- α), α are calculator memory occupancy, and inlineCnt is main profile total number, and size is described three
The size of geological data is tieed up, phym is computer physical memory size；The step of splicing, is carried out in a manner of being segmented interception.
In one exemplary embodiment of the present invention, the 3D seismic data is CRP gather after dynamic correction.
In one exemplary embodiment of the present invention, the α can be selected in the range of 70%~90%.
In one exemplary embodiment of the present invention, the object the step of splicing is phase each other in n group subdatas
The shared Value of residual static correction of the intersection of adjacent subdata group, it may include the step of the splicing：Inspection to the object
The geodetic coordinates of wave point or shot point carries out spatial domain spread, finds out the i-th subdata group adjacent to each other and i+1 subdata group
The min coordinates minX and maximum coordinates maxX and min coordinates minY and maximum on y-coordinate axle on x coordinate axle
Coordinate maxY, demarcated with being used as through (minX, 0.5 × (maxY+minY)) and (maxX, 0.5 × (maxY+minY)) straight line
Line, line of demarcation toward the smaller direction of main profile, take the i-th subdata group calculate obtained by Value of residual static correction, line of demarcation is toward main survey
Line takes i+1 subdata group to calculate resulting Value of residual static correction, so repeatedly, completed in n group subdatas compared with general orientation
The splicing of all adjacent subdata groups, wherein, 1 ＜ i ＜ n.
Compared with prior art, beneficial effects of the present invention include：It is quiet fast non-linear residue can be calculated using unit
School is measured, can abundant appliance computer memory size, solves the problems, such as that unit processing big data efficiency is low, and spelling is intercepted using segmentation
The lap for the static correction value that every group of data calculate is connect, meets the white noise characteristic of Value of residual static correction random distribution, can try one's best
Reduce the negative effect of splicing tape, it is ensured that stacking image works well.
Embodiment
Hereinafter, it is remaining quiet to describe quickly being handled based on unit for the present invention in detail that exemplary embodiment will be combined
Correction calculation method.
In one exemplary embodiment of the present invention, the residual static correction computational methods that are quickly handled based on unit can
In the case of unit, the fast automatic Value of residual static correction for asking for three-dimensional big data.In the present example embodiment, it is based on
The residual static correction computational methods that unit is quickly handled may include following steps：
(1) geological data adaptive grouping
Geological data adaptive grouping step obtains n group subdatas by carrying out adaptive grouping to 3D seismic data.
Specifically, geological data adaptive grouping can be carried out in the following manner.
3D seismic data is read first, and the data demand is CRP gather after correction, obtains size of data
Size, first master line FirstInline, main profile total number inlinCnt, statistics calculate meter physical memory size phym.
When data adaptive is grouped, 2 factors must be considered, first, the occupancy of internal memory, second, between data group and group coincidence factor (for
The splicing of Value of residual static correction lays the foundation below).
N=(1+ β) × size/ (phym × α) (1)
Wherein β is the coincidence factor of adjacent data group, must meet a constraints, i.e. size × β<Phym × (1- α), α
For memory usage.For example, in the present exemplary embodiment, α can be 80%, i.e., calculate the size of data for being stored in internal memory every time
For the 80% of whole internal memory；Certainly, α can also select to determine in the range of 70%~90%, but the invention is not restricted to this.n
Rounded using a method is entered, i.e., whole work area data are divided into n groups by program acquiescence.In actual calculating, data press the bar number of main profile
Packet.Before (n-1) secondary distribution main profile bar number it is identical, be k.
K=inlineCnt/n+size × β, rounded using a method is entered.
The main profile bar number of the data of last time distribution is m.
M=inlineCnt-(n-1) × k (2)
(2) the non-linear residual static correction of grouped data
The non-linear residual static correction step of grouped data is to the non-linear surplus of every group of subdata in n group subdatas
Remaining static correction value is calculated respectively.For example, non-linear residual static correction side can be applied in order to the n groups subdata after packet
Method (for example, simulated annealing, genetic algorithm etc.) calculates Value of residual static correction.
(3) splicing after non-linear residual static correction
Splicing step after non-linear residual static correction obtains the surplus of whole work area by splicing based on segmentation interception technology
Remaining static correction value.Why n group Value of residual static correction is spliced, be derived from static correction value in field of seismic exploration all
It is to be based on earth's surface-consistent principle, i.e., same physical points, which have, only only has a correcting value.
The object of splicing is the coincidence of subdata group adjacent to each other in the n sub- data groups that adaptive grouping step obtains
The shared Value of residual static correction in part.By taking geophone station Value of residual static correction as an example, (it is small that n1, n2 are all higher than 1 to two adjacent groups n1, n2
In n) there is identical geophone station, but static correction value is variant, and using based on segmentation interception technology splicing, i.e., program is according to geophone station
Geodetic coordinates, area of space spread is carried out to it, program finds out data n1 and n2 minimum, maximum x coordinate minX, maxX,
Minimum y-coordinate minY, maxY, calculate line of demarcation, i.e., through (minX, 0.5 × (maxY+minY)) and (maxX, 0.5 ×
(maxY+minY) straight line).Line of demarcation is toward the smaller direction of main profile, the Value of residual static correction for taking n1 to calculate, otherwise takes n2
Static correction value.The geophone station Value of residual static correction that other packets calculate, complete to splice by same procedure.
The connecting method of the static correction value of shot point is identical with above-mentioned geophone station joining method.
Generally speaking, the residual static correction computational methods of the invention quickly handled based on unit can fast apply unit
Calculate non-linear residual statics, can abundant appliance computer memory size, solve the problems, such as that unit processing big data efficiency is low,
And splice the lap for the static correction value that every group of data calculate using segmentation interception, meet Value of residual static correction random distribution
White noise characteristic, the negative effect of splicing tape can be reduced as far as possible, it is ensured that stacking image effect.
Although combined exemplary embodiment describes the present invention above, those of ordinary skill in the art should be clear
Chu, in the case where not departing from spirit and scope by the claims, various modifications can be carried out to above-described embodiment.
Claims (4)
1. a kind of residual static correction computational methods quickly handled based on unit, it is characterised in that methods described includes step：
Adaptive grouping is carried out to 3D seismic data, obtains n group subdatas；
Complete the calculating to the non-linear Value of residual static correction of every group of subdata in n group subdatas；
The Value of residual static correction that every group of subdata is calculated in splicing n group subdatas, wherein,
The step of adaptive grouping, includes：3D seismic data is read, obtains size, the first master of the 3D seismic data
Survey line, main profile total number, statistic computer physical memory size；By main profile bar number since the first main profile it is right in order
The 3D seismic data is grouped, and the main profile bar number of preceding n-1 distribution is k, the main profile bar number of last time distribution
For m, wherein,
K=inlineCnt/n+size × β and use are entered a method and rounded,
M=inlineCnt-(n-1) × k,
N=(1+ β) × size/ (phym × α) and use are entered a method and rounded,
Wherein, β is the coincidence factor of two groups of subdatas adjacent to each other, and meets size × β<Phym × (1- α), α is computer
Memory usage, inlineCnt are main profile total number, and size is the size of the 3D seismic data, and phym is computer
Physical memory size；
The step of splicing, is carried out in a manner of being segmented interception.
2. the residual static correction computational methods according to claim 1 quickly handled based on unit, it is characterised in that described
3D seismic data is CRP gather after dynamic correction.
3. the residual static correction computational methods according to claim 1 quickly handled based on unit, it is characterised in that described
α is selected in the range of 70%~90%.
4. the residual static correction computational methods according to claim 1 quickly handled based on unit, it is characterised in that described
The object of the step of splicing is the shared residual static correction of the intersection of two groups of subdatas adjacent to each other in n group subdatas
Amount, include the step of the splicing：The geodetic coordinates of geophone station or shot point to the object carries out spatial domain spread, finds out
The min coordinates minX and maximum coordinates on x coordinate axle of i-th group of subdata adjacent to each other and i+1 group subdata
MaxX and min coordinates minY and maximum coordinates maxY on y-coordinate axle, with through (minX, 0.5 × (maxY+
MinY)) and the straight line of (maxX, 0.5 × (maxY+minY)) takes i-th as line of demarcation, line of demarcation toward the smaller direction of main profile
Group subdata calculates resulting Value of residual static correction, and line of demarcation, compared with general orientation, takes i+1 group subdata to calculate institute toward main profile
Obtained Value of residual static correction, so repeatedly, the splicing to all two groups of adjacent subdatas in n group subdatas is completed, wherein,
1 ＜ i ＜ n.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610682656.5A CN106094035B (en) | 2016-08-18 | 2016-08-18 | A kind of residual static correction computational methods quickly handled based on unit |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201610682656.5A CN106094035B (en) | 2016-08-18 | 2016-08-18 | A kind of residual static correction computational methods quickly handled based on unit |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106094035A CN106094035A (en) | 2016-11-09 |
| CN106094035B true CN106094035B (en) | 2017-11-14 |
Family
ID=58070437
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201610682656.5A Active CN106094035B (en) | 2016-08-18 | 2016-08-18 | A kind of residual static correction computational methods quickly handled based on unit |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106094035B (en) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108535779B (en) * | 2018-03-30 | 2019-07-23 | 中国石油化工股份有限公司 | A kind of optimization method for seismic exploration data |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102231155A (en) * | 2011-06-03 | 2011-11-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for managing and organizing three-dimensional seismic data |
| CN103675908A (en) * | 2012-09-21 | 2014-03-26 | 中国石油化工股份有限公司 | Wave-equation reverse-time migration imaging method for mass-data graphic processing unit |
| CN105403913A (en) * | 2015-10-30 | 2016-03-16 | 中国石油天然气集团公司 | Pre-stack depth migration method and device |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9182913B2 (en) * | 2011-10-18 | 2015-11-10 | Ubiterra Corporation | Apparatus, system and method for the efficient storage and retrieval of 3-dimensionally organized data in cloud-based computing architectures |
-
2016
- 2016-08-18 CN CN201610682656.5A patent/CN106094035B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102231155A (en) * | 2011-06-03 | 2011-11-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Method for managing and organizing three-dimensional seismic data |
| CN103675908A (en) * | 2012-09-21 | 2014-03-26 | 中国石油化工股份有限公司 | Wave-equation reverse-time migration imaging method for mass-data graphic processing unit |
| CN105403913A (en) * | 2015-10-30 | 2016-03-16 | 中国石油天然气集团公司 | Pre-stack depth migration method and device |
Non-Patent Citations (1)
| Title |
|---|
| Fast residual static correction method using first arrivals in the 3D shot–receiver united domain;Meng Zou et al.;《Exploration Geophysics》;20160420;第A-I页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106094035A (en) | 2016-11-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US10288752B2 (en) | Automatic focus identification method and system for Karst cave reservoir | |
| CN104483705B (en) | A kind of three-dimensional residual static corrections | |
| CN104122587B (en) | A kind of abnormal first arrival recognition methods based on migration in offset domain and system | |
| CN107918156B (en) | Detect the method and device of subsea node acquisition earthquake data polarity | |
| CN106680869B (en) | The detection of micro-seismic event and localization method and device | |
| CN105954797B (en) | The fault recognition method and apparatus of seismic data | |
| CN104181599A (en) | Refracted wave static correction handling method and system based on near-surface layer | |
| CN104656139A (en) | Static correction method and equipment for converted wave | |
| CN103713318A (en) | Three-dimensional stacked seismic data splicing method | |
| CN108845358B (en) | Tomography and the recognition methods of structural anomaly body and device | |
| CN106556877B (en) | A kind of earth magnetism Tonghua method and device | |
| CN105319576A (en) | Concave polygon region observation system automatic full coverage quick layout method | |
| CN105467460A (en) | Method and device for electromagnetic prospecting | |
| CN104502967B (en) | The method and device of quick obtaining seismic prospecting observation system bin information | |
| CN109696705A (en) | Deep-sea OBS uplink and downlink Seismic imaging data correcting method and device | |
| CN106094035B (en) | A kind of residual static correction computational methods quickly handled based on unit | |
| Jiang et al. | 3D seismic geometry quality control and corrections by applying machine learning | |
| CN107329174A (en) | The method and device of well logging information is shown on seismic profile | |
| CN106054252B (en) | A kind of method and device of pre-stack time migration | |
| CN105866839B (en) | Static correction method and static correction device on basis of common-attitude trace gathers | |
| CN102353993B (en) | Method for checking closing error of velocity field during two-dimensional work area speed analysis in real time | |
| CN109782355A (en) | The detection method and device of OBS detection point drift | |
| CN103592627A (en) | Method and device for determining combined center of controlled seismic sources | |
| CN108364349A (en) | Metope scanning modeling method, apparatus, terminal and computer readable storage medium | |
| CN104216016A (en) | Converted-wave dynamic correcting method and system for longitudinal-wave constraint scanning |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20180208 Address after: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee after: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp. Address before: Shuangliu County Huayang Huayang Road in Chengdu city of Sichuan Province in 610213 section of No. 216, Igawa geophysical exploration company of the Ministry of science and technology Patentee before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd. |
|
| TR01 | Transfer of patent right | ||
| TR01 | Transfer of patent right |
Effective date of registration: 20200914 Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No. Co-patentee after: BGP Inc., China National Petroleum Corp. Patentee after: CHINA NATIONAL PETROLEUM Corp. Address before: 072751 Zhuozhou, Baoding, Fan Yang Road West, No. 189 Patentee before: BGP Inc., China National Petroleum Corp. |