CN106094035A - A kind of residual static correction computational methods quickly processed based on unit - Google Patents
A kind of residual static correction computational methods quickly processed based on unit Download PDFInfo
- Publication number
- CN106094035A CN106094035A CN201610682656.5A CN201610682656A CN106094035A CN 106094035 A CN106094035 A CN 106094035A CN 201610682656 A CN201610682656 A CN 201610682656A CN 106094035 A CN106094035 A CN 106094035A
- Authority
- CN
- China
- Prior art keywords
- subdata
- static correction
- group
- residual static
- 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.)
- Granted
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
Abstract
The invention provides a kind of residual static correction computational methods quickly processed based on unit.Described method includes step: 3D seismic data is carried out adaptive grouping, obtains n group subdata；Complete the calculating to the non-linear Value of residual static correction often organizing subdata in n group subdata；Splicing n group subdata often organizes the calculated Value of residual static correction of subdata.The present invention can fast apply unit to calculate non-linear residual statics, can abundant appliance computer memory size, solve unit and process the problem that big data efficiency is low, and use segmentation to intercept the lap of the splicing often static correction value that group data calculate, meet the white noise characteristic of Value of residual static correction random distribution, the negative effect that splicing brings can be reduced, it is ensured that stacking image effect as far as possible.
Description
Technical field
The present invention relates to seismic prospecting data processing technology field, it relates in particular to a kind of unit quickly process surplus
Remaining static correction computational methods, it is possible to self adaptation splicing Value of residual static correction, are applied to oil seismic exploration, especially can be 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, and earthquake provides
Material signal to noise ratio is low, and residual static correction problem is serious, only relies on the linear reflective ripple residual static correction of routine and can not obtain and make us full
The imaging effect of meaning.Non-linear residual static correction mainly has simulated annealing, genetic algorithm etc., possesses more quiet than conventional linear residue
Correct the ability of higher solution low SNR data residual static correction problem.
Computational efficiency is a bottleneck of non-linear residual static correction technology, and it relates to the anti-of prestack common midpoint gather
Re-reading write, repeatedly asking for of cross-correlation function.Field condition processes and uses unit operation, and calculator memory is usually no more than
100G, in the face of hundreds of G, even goes up T three-dimensional road collection data, if dynamic(al) correction data are placed on hard disk, reads big for consuming every time
The amount time；If dynamic(al) correction data being placed on internal memory can save the reading time, but the internal memory of computer can not store big at present
Geological data.Conventional way uses geological data packet to calculate, more manually splices Value of residual static correction, it has not been convenient to Ye Bu section
Learn.
Summary of the invention
It is an object of the invention to solve at least one in deficiencies of the prior art.
Such as, an object of the present invention is to provide a kind of more reasonable and science unit quickly to calculate non-linear surplus
The method of Yu Jing school amount, the method for the present invention is particularly suitable for when field condition processes using unit operation quickly to ask for non-linear
Residual statics, thus can fully apply single calculator memory size, solve three-dimensional big data redundancy static correction value and quickly ask
Take a difficult problem.
The invention provides a kind of residual static correction computational methods quickly processed based on unit.Described method includes step
Rapid: 3D seismic data is carried out adaptive grouping, obtain n group subdata；Complete the often group subdata in n group subdata
The calculating of non-linear Value of residual static correction；Splicing n group subdata often organizes the calculated Value of residual static correction of subdata, its
In, the step of described adaptive grouping includes: read 3D seismic data, obtains the size of this 3D seismic data, the first master
Survey line, main profile total number, statistic computer physical memory size；It is the most right to start from the first main profile by main profile bar number
Described 3D seismic data is grouped, and the main profile bar number of front n-1 distribution is k, the main profile bar number of last distribution
For m, wherein, k=inlineCnt/n+size × β and using rounds into a method, m=inlineCnt (n-1) × k, n=(1+
β) × size/ (phym × α) and using rounds into a method, and wherein, β is the coincidence factor of subdata group adjacent one another are, and meets
Size × β < phym × (1-α), α are calculator memory occupancy, and inlineCnt is main profile total number, and size is described three
The size of dimension geological data, phym is computer physical memory size；The step of described splicing is carried out in the way of segmentation intercepting.
In one exemplary embodiment of the present invention, described 3D seismic data is CRP gather after dynamic(al) correction.
In one exemplary embodiment of the present invention, described α can select in the range of 70%～90%.
In one exemplary embodiment of the present invention, the step of described splicing to as if n group subdata in phase each other
The Value of residual static correction that the intersection of adjacent subdata group is total, the step of described splicing comprises the steps that the inspection to described object
The geodetic coordinates of wave point or shot point carries out spatial domain spread, finds out the i-th subdata group adjacent one another are and i+1 subdata group
Min coordinates minX on x coordinate axle and maximum coordinates maxX and min coordinates minY on y-coordinate axle and maximum
Coordinate maxY, with the straight line through (minX, 0.5 × (maxY+minY)) and (maxX, 0.5 × (maxX+minX)) as boundary
Line, demarcation line forgets the less direction of main profile, takes the Value of residual static correction obtained by the i-th subdata group calculates, and main survey is forgotten in demarcation line
Line relatively general orientation, takes the Value of residual static correction obtained by i+1 subdata group calculates, the most repeatedly, completes in n group subdata
The splicing of all adjacent subdata groups, wherein, 1 ＜ i ＜ n.
Compared with prior art, the beneficial effect comprise that can fast apply unit calculate non-linear residue quiet
School is measured, can abundant appliance computer memory size, solve unit and process the problem that big data efficiency is low, and use segmentation to intercept spelling
Connect the lap of the often static correction value that group data calculate, meet the white noise characteristic of Value of residual static correction random distribution, can try one's best
Reduce the negative effect that splicing brings, it is ensured that stacking image is respond well.
Detailed description of the invention
Hereinafter, quiet to the residue quickly processed based on unit describing the present invention in detail by combining exemplary embodiment
Correction calculation method.
In one exemplary embodiment of the present invention, the residual static correction computational methods quickly processed based on unit, can
In the case of being applied to unit, the fast automatic Value of residual static correction asking for three-dimensional big data.In the present example embodiment, based on
The residual static correction computational methods that unit quickly processes can comprise the following steps that
(1) geological data adaptive grouping
Geological data adaptive grouping step, by 3D seismic data is carried out adaptive grouping, obtains n group subdata.
Specifically, geological data adaptive grouping can be carried out in the following manner.
First reading 3D seismic data, this data demand is CRP gather after dynamic(al) correction, obtains size of data
Size, Article 1 master line FirstInline, main profile total number inlinCnt, statistical computation meter physical memory size phym.
During data adaptive packet, must consider 2 factors, one is the occupancy of internal memory, and two is that the coincidence factor between data set and group (is
The splicing of Value of residual static correction below lays the foundation).
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.Such as, in this exemplary embodiment, α can be 80%, calculates the size of data being stored in internal memory i.e. every time
For whole internal memory 80%；Certainly, α can also select to determine in the range of 70%～90%, but the invention is not restricted to this.n
Use and round into a method, i.e. whole work area data are divided into n group by program acquiescence.In Practical Calculation, data press the bar number of main profile
Packet.Before the main profile bar number of (n-1) secondary distribution identical, for k.
K=inlineCnt/n+size × β, uses and rounds into a method.
The main profile bar number of the data of last 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 and remains the non-linear of subdata of often organizing in n group subdata
Remaining static correction value calculates respectively.Such as, the n group subdata after packet can be applied non-linear residual static correction side in order
Method (such as, 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 is spliced by intercepting technology based on segmentation, obtains the surplus of whole work area
Remaining static correction value.Why to splice n group Value of residual static correction, the static correction value being derived from field of seismic exploration is all
Being based on earth's surface-consistent principle, the most same physical points has only only one of which correcting value.
Splicing to as if the n that obtains of adaptive grouping step sub-data set in the coincidence of subdata group adjacent one another are
The Value of residual static correction that part is total.As a example by geophone station Value of residual static correction, (it is little that n1, n2 are all higher than 1 for two adjacent groups n1, n2
In n) there being identical geophone station, but static correction value is variant, uses and intercepts technology splicing based on segmentation, and i.e. program is according to geophone station
Geodetic coordinates, it is carried out area of space spread, program finds out data n1 and the minimum of n2, maximum x coordinate minX, maxX,
Minimum y-coordinate minY, maxY, calculate demarcation line, i.e. through (minX, 0.5 × (maxY+minY)) and (maxX, 0.5 ×
(maxX+minX) straight line).Demarcation line, toward the less direction of main profile, takes the Value of residual static correction that n1 calculates, otherwise takes n2
Static correction value.The geophone station Value of residual static correction that other packets calculate, completes splicing 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 quickly processed based on unit of the present invention can fast apply unit
Calculate non-linear residual statics, can abundant appliance computer memory size, solve unit and process the problem that big data efficiency is low,
And use segmentation to intercept the lap of the splicing often static correction value that group data calculate, meet Value of residual static correction random distribution
White noise characteristic, can reduce the negative effect that splicing brings, it is ensured that stacking image effect as far as possible.
Although describing the present invention already in connection with exemplary embodiment above, but those of ordinary skill in the art should be clear
Chu, in the case of without departing from spirit and scope by the claims, can carry out various amendment to above-described embodiment.
Claims (4)
1. the residual static correction computational methods quickly processed based on unit, it is characterised in that described method includes step:
3D seismic data is carried out adaptive grouping, obtains n group subdata；
Complete the calculating to the non-linear Value of residual static correction often organizing subdata in n group subdata；
Splice and n group subdata is often organized the calculated Value of residual static correction of subdata, wherein,
The step of described adaptive grouping includes: read 3D seismic data, obtains the size of this 3D seismic data, the first master
Survey line, main profile total number, statistic computer physical memory size；It is the most right to start from the first main profile by main profile bar number
Described 3D seismic data is grouped, and the main profile bar number of front n-1 distribution is k, the main profile bar number of last distribution
For m, wherein,
K=inlineCnt/n+size × β and using rounds into a method,
M=inlineCnt (n-1) × k,
N=(1+ β) × size/ (phym × α) and using rounds into a method,
Wherein, β is the coincidence factor of subdata group adjacent one another are, and meets size × β < phym × (1-α), in α is computer
Depositing occupancy, inlineCnt is main profile total number, and size is the size of described 3D seismic data, and phym is computer thing
Reason memory size；
The step of described splicing is carried out in the way of segmentation intercepting.
The residual static correction computational methods quickly processed based on unit the most according to claim 1, it is characterised in that described
3D seismic data is CRP gather after dynamic(al) correction.
The residual static correction computational methods quickly processed based on unit the most according to claim 1, it is characterised in that described
α selects in the range of 70%～90%.
The residual static correction computational methods quickly processed based on unit the most according to claim 1, it is characterised in that described
The Value of residual static correction total to the intersection liking subdata group adjacent one another are in n group subdata of the step of splicing,
The step of described splicing includes: the geophone station of described object or the geodetic coordinates of shot point are carried out spatial domain spread, finds out that
Min coordinates minX on x coordinate axle of this i-th adjacent subdata group and i+1 subdata group and maximum coordinates maxX,
And min coordinates minY on y-coordinate axle and maximum coordinates maxY, with through (minX, 0.5 × (maxY+minY)) and
The straight line of (maxX, 0.5 × (maxX+minX)) is as demarcation line, and the less direction of main profile is forgotten in demarcation line, takes the i-th subdata group
Value of residual static correction obtained by calculating, demarcation line is forgotten main profile relatively general orientation, is taken remaining obtained by i+1 subdata group calculates
Remaining static correction value, the most repeatedly, completes the splicing to adjacent subdata groups all in n group subdata, 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 true CN106094035A (en) | 2016-11-09 |
| CN106094035B 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) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108535779A (en) * | 2018-03-30 | 2018-09-14 | 中国石油化工股份有限公司 | A kind of optimization method for seismic exploration data |
Citations (4)
| 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 |
| US20130097170A1 (en) * | 2011-10-18 | 2013-04-18 | Ubiterra Corporation | Apparatus, system and method for the efficient storage and retrieval of 3-dimensionally organized data in cloud-based computing architectures |
| 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 |
-
2016
- 2016-08-18 CN CN201610682656.5A patent/CN106094035B/en active Active
Patent Citations (4)
| 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 |
| US20130097170A1 (en) * | 2011-10-18 | 2013-04-18 | Ubiterra Corporation | Apparatus, system and method for the efficient storage and retrieval of 3-dimensionally organized data in cloud-based computing architectures |
| 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 |
|---|
| MENG ZOU ET AL.: "Fast residual static correction method using first arrivals in the 3D shot–receiver united domain", 《EXPLORATION GEOPHYSICS》 * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108535779A (en) * | 2018-03-30 | 2018-09-14 | 中国石油化工股份有限公司 | A kind of optimization method for seismic exploration data |
| CN108535779B (en) * | 2018-03-30 | 2019-07-23 | 中国石油化工股份有限公司 | A kind of optimization method for seismic exploration data |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106094035B (en) | 2017-11-14 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104181599B (en) | A kind of refractive wave static correction processing method based near surface formation and system | |
| CN103713318B (en) | The method that three-dimensional post-stack seismic data is spliced | |
| CN104122587B (en) | A kind of abnormal first arrival recognition methods based on migration in offset domain and system | |
| CN102944896B (en) | The modelling static correcting method of surface survey data | |
| CN104483705B (en) | A kind of three-dimensional residual static corrections | |
| US9329291B2 (en) | Non-hyperbolic correction of seismic data | |
| CN107783186B (en) | A kind of near-surface velocity modeling method and device | |
| CN105954797B (en) | The fault recognition method and apparatus of seismic data | |
| CN101984366A (en) | Fluctuating surface pre-stack time migration method and device based on common aperture surface | |
| CN104656139A (en) | Static correction method and equipment for converted wave | |
| CN106680869A (en) | Microearthquake event detecting and positioning method and apparatus | |
| CN104199103A (en) | Static correction method and static correction device for converted wave | |
| CN105319576A (en) | Concave polygon region observation system automatic full coverage quick layout method | |
| CN107817525A (en) | A kind of near surface constant speed fill method and device | |
| CN109696705A (en) | Deep-sea OBS uplink and downlink Seismic imaging data correcting method and device | |
| CN104502967B (en) | The method and device of quick obtaining seismic prospecting observation system bin information | |
| Jiang et al. | 3D seismic geometry quality control and corrections by applying machine learning | |
| CN106054252B (en) | A kind of method and device of pre-stack time migration | |
| CN108845350A (en) | The method and device of inverting two-dimension speed model | |
| CN105242315A (en) | Geologic structure imagining method and device | |
| CN106094035B (en) | A kind of residual static correction computational methods quickly handled based on unit | |
| EP3612863B1 (en) | Post-stack kirchhoff depth de-migration method for tilted transverse isotropic (tti) and heterogeneous media based on ray tracing on migrated data | |
| CN105717538B (en) | Relief surface earthquake data offset reference plane conversion method and device | |
| CN110376644B (en) | Pre-stack seismic data splicing method and device | |
| CN111337973B (en) | Seismic data reconstruction method and system |
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 |
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. |
|
| TR01 | Transfer of patent right |