CN103698807B - Scalariform two-dimensional wide-band observation system design method - Google Patents

Scalariform two-dimensional wide-band observation system design method Download PDF

Info

Publication number
CN103698807B
CN103698807B CN201210366027.3A CN201210366027A CN103698807B CN 103698807 B CN103698807 B CN 103698807B CN 201210366027 A CN201210366027 A CN 201210366027A CN 103698807 B CN103698807 B CN 103698807B
Authority
CN
China
Prior art keywords
spacing
line
wave
width
velocity
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
Application number
CN201210366027.3A
Other languages
Chinese (zh)
Other versions
CN103698807A (en
Inventor
宁宏晓
章多荣
李国顺
尹吴海
吴永国
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201210366027.3A priority Critical patent/CN103698807B/en
Publication of CN103698807A publication Critical patent/CN103698807A/en
Application granted granted Critical
Publication of CN103698807B publication Critical patent/CN103698807B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Abstract

The present invention is the Scalariform two-dimensional wide-band observation system design method improving seismic prospecting imaging section signal to noise ratio. Apply known earthquake and well data, and investigate apparent velocity, the wavelength of exploratory area scattered noise, determine the minimum of recording geometry and the element interval of maximum width, maximum reception line-spacing, reception line number and receiver pattern, centered by acceptance point, wave-detector is transversely become elementary cell with longitudinal Symmetric Composite respectively, transverse combination base receives line-spacing apart from equaling, and array length equals distance in length and breadth, by elementary cell according to the coupling combination of wave-detector vertical and horizontal. Present invention achieves Spatial continual pressure to make an uproar, thus really realize the horizontal pressure of associating between recording geometry reception line and make an uproar, also making up spatial sampling deficiency avoids alias noise simultaneously, is conducive to raising starting material quality and imaging section signal to noise ratio.

Description

Scalariform two-dimensional wide-band observation system design method
Technical field
The present invention relates to seismic data acquisition parameter designing, mainly a kind of Scalariform two-dimensional wide-band observation system design method.
Background technology
Complicated mountain region is all generally that tectonic movement strong deformation region, earth's surface and underground structure are all very complicated, grows polytype strong energy noise, and starting material signal to noise ratio is extremely low. Seismic wave field is extremely complicated, and data imaging is very difficult. From at the end of last century, Mountainous Seismic Exploration Technology has had and has developed on a large scale very much, but in some extremely complicated areas, seismic data quality still achieves no breakthrough, and seriously constrains the oil-gas exploration and development progress in relevant area.
Existing two-dimensional observational system designs, proving bin size according to physical model and testing data and guarantee that the requirement of geologic objective lateral resolution, enough spreadss guarantee that the sampling to steep dip geologic body, suitable covering number of times guarantee that stack velocity and residual static correction are asked for, such recording geometry is adapted to the region that starting material has certain signal to noise ratio. In low signal-to-noise ratio district, conventional two-dimensional observational system does not meet the requirement of data imaging, so having developed the laterally big combination technique based on conventional two-dimensional observational system, object is that compacting side noise improves starting material quality, such observation mode is adapted to side noise direction and extremely by force and combines the time difference relatively zonule, and the region violent at topographic relief is limited. So for the low signal-to-noise ratio district that landform acutely rises and falls, develop wide line observation technology, but conventional wide line collection method, general only consideration increases the horizontal line number that receives and improves covering number of times, and consider deficiency to pressing continuously between the lines to make an uproar, cause wide line laterally to press ability of making an uproar greatly to reduce.
Summary of the invention
The object of the invention is to provide a kind of Scalariform two-dimensional wide-band observation system design method adapting to the violent region of topographic relief.
The present invention is realized by following step:
1) earthquake and well data is gathered; or apply known seismic section, velocity spectrum and well data; extracting zone of interest double-pass reflection time, the stratigraphic dip time difference and stack velocity or buried depth of strata information, the frequency protected in conjunction with geological tasks determines bin size, spreads.
2) applying the known seismic interpretation section in exploratory area and structural map, velocity spectrum or well data, extract stratigraphic dip and interval velocity information, the frequency protected in conjunction with geological tasks determines maximum width.
Step 2) described in maximum width be protect the significant wave of the most short-sighted wavelength can in-phase stacking and the horizontal width that limits.
3) apparent wavelength of single big gun material computation strong refraction interference wave that the apparent velocity of investigation exploratory area scattered noise, apparent cycle, group number and intensity or application have obtained, it is determined that minimum width.
Minimum width described in step 3) is the horizontal width suppressed the strongest scattered noise and limit.
4) by the greatest combined base distance of top layer model velocity, the restriction of the thickness data calculation combination time difference, it is determined that maximum reception line-spacing.
Maximum reception line-spacing described in step 4) is the greatest combined base distance that protection receiver pattern does not suppress significant wave and limit. Not distorted according to significant wave, the calculating significant wave described in step 4) should not distort.
5) expectation designed according to multiple stacking Statistical Effect covers number of times and determines to receive line number;
Determination described in step 5) receives line number: reception line number is multiplied by reception line-spacing after subtracting 1 should equal recording geometry width, receives line-spacing and is not more than maximum reception line-spacing, and recording geometry width is not less than minimum width and is not more than maximum width;
6) according to the Coherence radii of the random noise in work area, it is determined that the element interval of receiver pattern;
7) centered by acceptance point, transversely with longitudinal Symmetric Composite respectively, wave-detector being become elementary cell, transverse combination base receives line-spacing apart from equaling, and array length equals distance in length and breadth;
8) by elementary cell according to the coupling combination of wave-detector vertical and horizontal.
Coupling combination described in step 8) receives line to be connected in the horizontal with the elementary cell received between line, and the elementary cell of same reception line is connected in the vertical, forms scalariform two-dimensional observational system.
Instant invention overcomes conventional two-dimensional observational system pressure make an uproar characteristic discontinuity, receive and combine pressure between line and make an uproar combination time difference that alias noise that ability, spatial sampling deficiency cause and laterally combining greatly causes to the deficiency of the aspects such as the disadvantageous effect of significant wave and static correction and intensity of illumination are inadequate, achieve Spatial continual pressure to make an uproar, thus really realize receiving the horizontal pressure of associating between line and make an uproar, the present invention simultaneously also makes up spatial sampling deficiency and avoids alias noise, is conducive to improving starting material quality and imaging section signal to noise ratio. The present invention is applicable to that scattering noise is grown, original signal to noise ratio is low, it is necessary to strengthens laterally pressure and makes an uproar ability and the high area covering number of times.
Accompanying drawing explanation
Fig. 1 schematic diagram of the present invention.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in detail.
The concrete implementation step of the present invention is as follows:
1) earthquake and well data is gathered; or apply known seismic section, velocity spectrum and well data; extracting zone of interest double-pass reflection time, the stratigraphic dip time difference and stack velocity or buried depth of strata information, the frequency protected in conjunction with geological tasks determines bin size, spreads.
Bin size��stratigraphic dip time difference, (unit: rice/millisecond) was divided by the protection frequencies of 2 times (unit: hertz)
Spreads �� equals buried depth of strata
2) applying the known seismic interpretation section in exploratory area and structural map, velocity spectrum or well data, extract stratigraphic dip and interval velocity information, the frequency protected in conjunction with geological tasks determines maximum width.
Step 2) described in maximum width be protect the significant wave of the most short-sighted wavelength can in-phase stacking and the horizontal width that limits.
Maximum width Wy is the twice of the horizontal width of CDP bin.
If the horizontal width of CDP bin is Ly, destination layer apparent dip in the horizontal is ��, and the reflection time arrived at first in same CDP bin is t1, and the reflection time finally arrived is t2, then have:
T2-t1=2 �� Ly �� tg ��/Vi (Vi is interval velocity)
Reflection wave in CDP bin width can in-phase stacking, then have:
2Ly �� tg ��/Vi��T/4=1/4Fp (T is for protecting frequency cycle, Fp for protection frequency)
That is, maximum width Wy=2 �� Ly��Vi/4Fptg ��
3) apparent wavelength of single big gun material computation strong refraction interference wave that the apparent velocity of investigation exploratory area scattered noise, apparent cycle, group number and intensity or application have obtained, it is determined that minimum width.
Minimum width described in step 3) is the horizontal width suppressed the strongest scattered noise and limit.
Minimum width is more than or equal to the maximum apparent wavelength of strong scattered noise.
4) by the greatest combined base distance of top layer model velocity, the restriction of the thickness data calculation combination time difference, it is determined that maximum reception line-spacing.
Maximum reception line-spacing described in step 4) is the greatest combined base distance that protection receiver pattern does not suppress significant wave and limit.
T / 2 ≥ ( x + L / 2 ) 2 + 4 h 0 2 - ( x - L / 2 ) 2 + 4 h 0 2 v + ΔH v H
In upper formula: T is significant wave period, L is array length, h0For significant wave reflects depth of stratum, V is formation velocity, and �� H is the combination discrepancy in elevation, VHFor the formation velocity in the combination discrepancy in elevation, x is offset distance.
5) expectation designed according to multiple stacking Statistical Effect covers number of times and determines to receive line number;
RLN=FoldexX��2SI
In upper formula: FoldexFor expecting to cover number of times, RLN is for receiving line number, and X is spreads, and SI is shotpoint spacing
Described determination receives line number: reception line number is multiplied by reception line-spacing after subtracting 1 should equal recording geometry width, receives line-spacing and is not more than maximum reception line-spacing, and recording geometry width is not less than minimum width and is not more than maximum width;
6) according to the Coherence radii of the random noise in work area, it is determined that the element interval of receiver pattern;
7) centered by acceptance point, transversely with longitudinal Symmetric Composite respectively, wave-detector being become elementary cell, transverse combination base receives line-spacing apart from equaling, and array length equals distance in length and breadth;
8) by elementary cell according to the coupling combination of wave-detector vertical and horizontal.
Described coupling combination receives line to be connected in the horizontal with the elementary cell received between line, and the elementary cell of same reception line is connected in the vertical, forms scalariform two-dimensional observational system.
The present invention implements in the following ways further in concrete embodiment:
Zone of fracture, inverse nappe, the low wave resistance geologic body etc. covered are explored, and will increase the intensity of illumination to geologic body, receiving the horizontal big gun line number of increase between line, should be formed " Scalariform two-dimensional wide-band observation system " that many big guns are multi-thread.
Earth's surface is smooth, the exploration of buried hill, underground or weak reflection geologic body, suitable " Scalariform two-dimensional wide-band observation system " adopting many big guns two-wire.
High-dip structure district, when acceptance point spacing does not meet high-dip structure imaging, it is possible to docking take-up carries out part road apart from the changing of the relative positions, forms " Scalariform two-dimensional wide-band observation system " of distortion.

Claims (3)

1. a Scalariform two-dimensional wide-band observation system design method, feature is realized by following step:
1) earthquake and well data is gathered, or apply known seismic section, velocity spectrum and well data, extracting zone of interest double-pass reflection time, the stratigraphic dip time difference and stack velocity and buried depth of strata information, the frequency protected in conjunction with geological tasks determines bin size, spreads;
Bin size��stratigraphic dip time difference divided by the protection frequency of 2 times, spreads �� buried depth of strata;
2) applying the known seismic section in exploratory area, structural map, velocity spectrum and well data, extract stratigraphic dip and interval velocity information, the frequency protected in conjunction with geological tasks determines maximum width;
3) apparent wavelength of the apparent velocity of investigation exploratory area scattered noise, apparent cycle, group number and Strength co-mputation strong refraction interference wave, it is determined that minimum width, or apply the apparent wavelength of the single big gun material computation strong refraction interference wave obtained, it is determined that minimum width;
Described minimum width is the horizontal width suppressed the strongest scattered noise and limit;
4) by the greatest combined base distance of top layer model velocity, the restriction of the thickness data calculation combination time difference, it is determined that maximum reception line-spacing;
Described maximum reception line-spacing is the greatest combined base distance that protection receiver pattern does not suppress significant wave and limit, and determines that significant wave does not distort by following formula:
T / 2 ≥ ( x + L / 2 ) 2 + 4 h 0 2 - ( x - L / 2 ) + 4 h 0 2 v + Δ H v H
In upper formula: T is significant wave period, L is array length, h0For significant wave reflects depth of stratum, v is formation velocity, and �� H is the combination discrepancy in elevation, vHFor the formation velocity in the combination discrepancy in elevation, x is offset distance;
5) expectation designed according to multiple stacking Statistical Effect covers number of times and determines to receive line number, and determines to receive line number by following formula:
RLN=Foldex/X��2SI
In upper formula: FoldexFor expecting to cover number of times, RLN is for receiving line number, and X is spreads, and SI is shotpoint spacing;
6) according to the Coherence radii of the random noise in work area, it is determined that the element interval of receiver pattern;
7) centered by acceptance point, transversely with longitudinal Symmetric Composite respectively, wave-detector being become elementary cell, transverse combination base receives line-spacing apart from equaling, and longitudinal combination base is apart from equaling distance;
8) by elementary cell according to the coupling combination of wave-detector vertical and horizontal;
Described coupling combination receives line to be connected in the horizontal with the elementary cell received between line, and the elementary cell of same reception line is connected in the vertical, forms scalariform two-dimensional observational system.
2. method according to claim 1, feature is step 2) described in maximum width be protect the significant wave of the most short-sighted wavelength can in-phase stacking and the horizontal width that limits.
3. method according to claim 1, feature is step 5) described in determination receive line number be: receive line number subtract 1 after be multiplied by receive line-spacing should equal recording geometry width, receiving line-spacing and be not more than maximum reception line-spacing, recording geometry width is not less than minimum width and is not more than maximum width.
CN201210366027.3A 2012-09-28 2012-09-28 Scalariform two-dimensional wide-band observation system design method Active CN103698807B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201210366027.3A CN103698807B (en) 2012-09-28 2012-09-28 Scalariform two-dimensional wide-band observation system design method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201210366027.3A CN103698807B (en) 2012-09-28 2012-09-28 Scalariform two-dimensional wide-band observation system design method

Publications (2)

Publication Number Publication Date
CN103698807A CN103698807A (en) 2014-04-02
CN103698807B true CN103698807B (en) 2016-06-08

Family

ID=50360398

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201210366027.3A Active CN103698807B (en) 2012-09-28 2012-09-28 Scalariform two-dimensional wide-band observation system design method

Country Status (1)

Country Link
CN (1) CN103698807B (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104199083A (en) * 2014-08-29 2014-12-10 中国石油天然气集团公司 Method and device for determining wide-line seismological observation system
CN105093268B (en) * 2015-06-30 2018-04-03 中国石油天然气集团公司 A kind of method and device based on wave field Lighting Design observation system parameter
CN106501841B (en) * 2016-09-19 2018-08-14 中国石油天然气集团公司 A kind of optimum design method and device of two dimension wode line seismic observation system
CN108845349B (en) * 2018-06-14 2020-10-16 中国海洋大学 Energy-based arrangement width design method
CN111983673B (en) * 2019-05-21 2023-08-22 中国石油天然气集团有限公司 Method and device for determining receiving line distance of three-dimensional seismic observation system
CN110174697A (en) * 2019-05-31 2019-08-27 中国煤炭地质总局物测队 A kind of the change sight simulation system and its analogy method of 3-d seismic exploration
CN112698394B (en) * 2019-10-22 2022-11-04 中国石油天然气集团有限公司 Method and device for determining size of design surface element of two-dimensional observation system
CN111025390B (en) * 2019-12-13 2022-07-08 中国海洋石油集团有限公司 Seismic data acquisition method and device
CN111538078A (en) * 2020-05-08 2020-08-14 中国石油天然气集团有限公司 Observation mode determination method and device of two-dimensional wide-line seismic observation system

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3885225A (en) * 1972-07-21 1975-05-20 Seiscom Delta Inc Broad line seismic profiling
CN101105536A (en) * 2006-07-12 2008-01-16 中国石油集团东方地球物理勘探有限责任公司 Complicated earth surface banding earthquake collection method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US3885225A (en) * 1972-07-21 1975-05-20 Seiscom Delta Inc Broad line seismic profiling
CN101105536A (en) * 2006-07-12 2008-01-16 中国石油集团东方地球物理勘探有限责任公司 Complicated earth surface banding earthquake collection method

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
《准噶尔盆地东部阜康断裂带宽线采集技术及其应用效果》;顾雯等;《天然气勘探与开发》;20111031;第34卷(第4期);第40—43页 *
《姜山洼陷宽线地震采集技术应用与效果》;张怀榜;《内蒙古石油化工》;20091231(第12期);第104—105页 *
《宽线观测大组合接收技术在阜康断裂带的应用》;朱鹏宇等;《勘探地球物理进展》;20101031;第33卷(第5期);第359-362页 *
《库车坳陷复杂山地地震采集技术》;彭才等;《西部探矿工程》;20091231(第11期);第136-138页 *
《柴达木盆地复杂山地采集技术及效果》;胡杰等;《中国石油勘探》;20061231(第6期);第52-57页 *

Also Published As

Publication number Publication date
CN103698807A (en) 2014-04-02

Similar Documents

Publication Publication Date Title
CN103698807B (en) Scalariform two-dimensional wide-band observation system design method
CN107526101B (en) A kind of acquisition and processing method obtaining earthquake reflected wave
CN108037526B (en) Reverse-time migration method based on all-wave wave field VSP/RVSP seismic data
CN101915939B (en) Surface wave pressing method
CN101581790B (en) Seismic sensor array device and data collecting method thereof
CN102073064B (en) Method for improving velocity spectrum resolution by using phase information
CN101819279B (en) High-precision three-dimensional static correction based on variogram fitting reconfiguration
CN102590862A (en) Prestack time migration method for compensating absorptive attenuation
CN102928878A (en) Amplitude balance quantitative evaluation method for three-dimensional earthquake observing system
CN104570116A (en) Geological marker bed-based time difference analyzing and correcting method
CN109669223A (en) The method for obtaining shallow earth's surface acoustic logging and synthetic seismogram using micro logging
CN102103215B (en) Method for suppressing surface waves of three-dimensional high-density seismic prospecting records before stack
CN102565852B (en) Angle domain pre-stack offset data processing method aiming to detect oil-gas-bearing property of reservoir
CN104007466B (en) The reservoir that a kind of no restriction from borehole data prestack inversion based on P-wave amplitude realizes and fluid prediction method
Wang et al. High-resolution crustal velocity imaging using ambient noise recordings from a high-density seismic array: An example from the Shangrao section of the Xinjiang basin, China
CN109856680A (en) A kind of Coastal beach area pull-type seismic reflection survey method
CN104459787A (en) Speed analysis method through seismic record of vertical receiving array
CN102830424B (en) A kind of receiver pattern calculation method of parameters
CN101900830B (en) Seismic scattering P-P wave imaging method
CN201489120U (en) Array device of earthquake sensors
CN104076395B (en) Based on minute surface Energy extraction and the formation method of filtering combination
Jiang et al. Distributed acoustic sensing for shallow structure imaging using mechanical noise: A case study in Guangzhou, China
CN104991271A (en) Prestack depth migration method based on VTI medium bidirectional wave field superposition
CN104155690A (en) Ellipsoid expansion-based three-dimensional seismic data stacking velocity calculation method
CN109212606A (en) Method for marine seismic data earth surface consistency vibration amplitude compensation processing 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
C14 Grant of patent or utility model
GR01 Patent grant