CN105445790A - Amplitude regularization method of blind source seismic wave field - Google Patents
Amplitude regularization method of blind source seismic wave field Download PDFInfo
- Publication number
- CN105445790A CN105445790A CN201410477842.6A CN201410477842A CN105445790A CN 105445790 A CN105445790 A CN 105445790A CN 201410477842 A CN201410477842 A CN 201410477842A CN 105445790 A CN105445790 A CN 105445790A
- Authority
- CN
- China
- Prior art keywords
- amplitude
- seismic
- regularization
- blind
- survey line
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention provides an amplitude regularization method of a blind source seismic wave field, and the amplitude regularization method belongs to the field of seismic exploration and development. The amplitude regularization method comprises the steps of: step 1, reading survey line practical seismic trace data; step 2, preprocessing the survey line practical seismic trace data, and eliminating low-frequency and high-frequency noise; step 3, carrying out amplitude regularization processing to obtain blind source survey line data after amplitude regularization; step 4, performing seismic response recovery and trace set construction, regarding the blind source survey line data after amplitude regularization as input of the seismic response recovery and trace set construction, and further recovering seismic response queried by a receiving point to obtain a recovery result; and step 5, judging whether the recovery result is satisfying, if not, returning to the step 3, and if so, then terminating programs.
Description
Technical field
The invention belongs to seismic prospecting and development field, be specifically related to a kind of amplitude rule method of blind focus earthquake wave field.
Background technology
Blind focus earthquake technology is a kind of brand-new seismic method that development in recent years is got up, excite from neighbourhood noise or subsurface geology activity the seismic event propagated through earth interior at ground table record, based on seismic interference method, therefrom recover the seismic response (Green function) between acceptance point, and then as active source seismic data, carry out follow-up various application.
Blind focus earthquake non-passive, but utilize mankind's activity near the various noise of nature or survey line as focus, therefore, the source type in blind source is various, source strength difference is large.The wave detector dynamic range that it adopts is large, and frequency low (< 1Hz).Compared with active source earthquake (as surface seismic), the feature of blind focus earthquake record is: (1) noise energy is strong, and useful signal is submerged among various noise; (2) capacity volume variance between each seismic trace is large especially; (3) seismic event type is various, comprises ground roll, direct wave, transmitted wave and multiple reflection thereof, various noise etc.
Blind seedbed seismic data processing, as shown in Figure 1, generally includes: (1) pre-processing of the information, as denoising, signal enhancing etc.But due to blind focus earthquake record, and the requirement of blind focus earthquake self-technique, make the conventional preprocess method that directly can not adopt active source seismic data, this is the difficult point in the seismic data processing of blind seedbed.(2) seismic response recovers, and based on seismic interference method, mainly contains cross-correlation, deconvolution at present and does three kinds of recovery algorithms mutually.(3) subsequent applications, as the seismic response between the acceptance point of recovery can be utilized to carry out tomography, thus sets up the rate pattern of underground medium, or carries out migration imaging process etc.
Summary of the invention
The object of the invention is to solve the difficult problem existed in above-mentioned prior art, a kind of amplitude rule method of blind focus earthquake wave field is provided, for the feature of blind focus earthquake record, namely the capacity volume variance between each seismic trace is large, eliminate the amplitude difference between each seismic trace in blind focus earthquake record, give prominence to and enhance effective wave field (as the terrible ripple) information recovered for seismic response, thus significantly improve precision and the resolution of seismic response recovery, for follow-up successful reflection wave imaging lays the foundation.
The present invention is achieved by the following technical solutions:
An amplitude rule method for blind focus earthquake wave field, comprising:
Step one reads survey line actual seismic track data;
Step 2 carries out pre-service to described survey line actual seismic track data, eliminates low frequency and high frequency noise;
Step 3 carries out amplitude regularization process, obtains the blind source survey line data after amplitude regularization;
Step 4 seismic response recovers to build with road collection: using the blind source survey line data after described amplitude regularization as the input that seismic response recovers and road collection builds, and then recovers the seismic response between acceptance point, and be restored result;
Step 5, judges whether described restoration result is satisfied with (as whether having the reflection wave of recovery to occur.), if dissatisfied, then return step 3, if satisfied, then terminator.
Described step 3 is achieved in that
(1) amplitude regularization factor p is selected;
(2) formula (1) is utilized to ask for the amplitude gain of a certain seismic trace, the amplitude after then utilizing formula (2) to obtain this seismic trace regularization:
x′
i＝g(t)x
i(2)
In formula, g (t) represents the amplitude gain of a certain seismic trace of blind focus earthquake record.X
irepresent the amplitude of some sampled points; X '
irepresent the amplitude after amplitude regularization.N represents total number of samples of a certain seismic trace in blind source record; P represents the amplitude regularization factor.
The span of described amplitude regularization factor p is p=1,1.5,2,2.5,3 ...According to amplitude regularization experience, when specifically implementing, usually get p=2.5.
Compared with prior art, the invention has the beneficial effects as follows: real data test in blind source shows: the amplitude regularization techniques of this blind source wave field, compared with the energy normalized technology of routine, effectively can eliminate the amplitude difference in blind source record between each seismic trace, give prominence to and enhance effective wave field (as the terrible ripple) information recovered for seismic response, between follow-up successfully recovery acceptance point, low seismic response provides safeguard.
Accompanying drawing explanation
Fig. 1: the step block diagram of blind seedbed seismic data processing
Fig. 2: blind focus earthquake wave field amplitude regularization treatment scheme
Fig. 3 a: the waterfrac treatment ground monitoring work area in the waterfrac treatment monitoring of somewhere, Sichuan and ground survey line record thereof
Fig. 3 b: the original wave field of the wz29-3 survey line in the waterfrac treatment monitoring of somewhere, Sichuan and ground survey line record thereof
Fig. 4 a: the wave field after bandpass filtering
Fig. 4 b: the wave field after energy normalized (i.e. p=2)
Fig. 4 c: wave field after amplitude regularization p=2.5
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail:
Blind focus earthquake wave field amplitude rule method is the chief component of blind source wave field catoptric imaging technology (as shown in Figure 1), and after blind focus earthquake wave field bandpass filtering, between the wavefield reception point of blind source, seismic response carries out before recovering.
The present invention has expanded the formula of energy normalized conventional in surface-seismic data process, that is:
x′
i＝g(t)x
i(2)
In formula, g (t) represents the amplitude gain of a certain seismic trace of blind focus earthquake record.X
irepresent the amplitude of some sampled points; X '
irepresent the amplitude after amplitude regularization.N represents total number of samples of a certain seismic trace in blind source record.P represents the amplitude regularization factor, and its span is generally p=1,1.5,2,2.5,3 ... Deng.When p equals 2, be energy normalized.Although adopt conventional energy normalized technology, also amplitude difference between each seismic trace can be eliminated to a certain extent, for the seismic response recovered between acceptance point, its effect is poor.
It should be noted that: the amplitude difference in blind source record between each seismic trace is eliminated in the effect of amplitude regularization factor p, its value principle is depending on amplitude difference between each seismic trace, and its value is larger, and effect is more obvious, usually p=2.5 is got, on rare occasion desirable p=3 or larger.On the other hand, whether the value of amplitude regularization factor p is suitable, can judge from the quality of seismic response restoration result.
As shown in Figure 2, the specific implementation step of blind focus earthquake wave field amplitude regularization is as follows:
Step one reads survey line actual seismic track data: the physical record reading in survey line WZ29-3 as shown in Figure 3 b.
Step 2 real data pre-service: as bandpass filtering, fundamental purpose eliminates low frequency and high frequency noise (Fig. 4 a).
Step 3 amplitude regularization process
(1) select amplitude regularization factor p, according to actual treatment experience, usually get p=2.5, on rare occasion desirable p=3 or larger;
(2) formula (1) is utilized to ask for the amplitude gain of a certain seismic trace, the amplitude after then utilizing formula (2) to obtain this seismic trace regularization, as illustrated in fig. 4 c.
Step 4 seismic response recovers to build with road collection: the blind source survey line data after amplitude regularization can be used as seismic response recover and road collection build (can with reference to Wapenaar, K. (2004), RetrievingtheelastodynamicGreen ' sfunctionofanarbitraryinhomogeneousmediumbycrosscorrelat ion, Phys.Rev.Lett., 93 (254301), input 1-4), and then have the seismic response between efficient recovery acceptance point.If be unsatisfied with restoration result, then got back to step 3, again choose the amplitude regularization factor (as desirable p=3).If satisfied to restoration result, then terminator.
Shown in Fig. 3 a is earth's surface, well waterfrac treatment work area, somewhere, Sichuan 29 monitoring survey line schematic diagram, in figure, has 10 surveys line, be respectively WZ29-1, WZ29-2, WZ29-3 ..., WZ29-9, WZ29-10.Every bar survey line 125 road, track pitch 25m, smallest offset distance 300m, maximum offset 3400m, Continuous Observation 48 hours.The raw monitored data of survey line 3 as shown in Figure 3 b, the result after bandpass filtering treatment, as shown in fig. 4 a.On this basis, carried out the process of conventional energy normalized (p=2) and above-mentioned amplitude regularization techniques (p=2.5) respectively, its result as shown in figures 4 b and 4 c.Can see: although conventional energy normalized result can reduce to a certain extent between amplitude difference (as shown in Fig. 4 b arrow), can not give prominence to and strengthen the information of terrible ripple after direct wave; And amplitude regularization (p=2.5) technical finesse result, not only eliminate the amplitude difference between each seismic trace, and give prominence to and enhance the information of terrible ripple after direct wave, as shown in arrow in Fig. 4 c, therefore, to contribute between acceptance point seismic response recover and road collection structure.
Technique scheme is one embodiment of the present invention, for those skilled in the art, on the basis that the invention discloses application process and principle, be easy to make various types of improvement or distortion, and the method be not limited only to described by the above-mentioned embodiment of the present invention, therefore previously described mode is just preferred, and does not have restrictive meaning.
Claims (4)
1. an amplitude rule method for blind focus earthquake wave field, is characterized in that: described method comprises:
Step one reads survey line actual seismic track data;
Step 2 carries out pre-service to described survey line actual seismic track data, eliminates low frequency and high frequency noise;
Step 3 carries out amplitude regularization process, obtains the blind source survey line data after amplitude regularization;
Step 4 seismic response recovers to build with road collection: using the blind source survey line data after described amplitude regularization as the input that seismic response recovers and road collection builds, and then the seismic response that recovery acceptance point is asked, be restored result;
Step 5, judges whether described restoration result is satisfied with, if dissatisfied, then returns step 3, if satisfied, then and terminator.
2. the amplitude rule method of blind focus earthquake wave field according to claim 1, is characterized in that: described step 3 is achieved in that
(1) amplitude regularization factor p is selected;
(2) formula (1) is utilized to ask for the amplitude gain of a certain seismic trace, the amplitude after then utilizing formula (2) to obtain this seismic trace regularization:
x′
i＝g(t)x
i(2)
In formula, g (t) represents the amplitude gain of a certain seismic trace of blind focus earthquake record, x
irepresent the amplitude of some sampled points; X '
irepresent the amplitude after amplitude regularization, N represents total number of samples of a certain seismic trace in blind source record; P represents the amplitude regularization factor.
3. the amplitude rule method of blind focus earthquake wave field according to claim 2, is characterized in that: the span of described amplitude regularization factor p is p=1,1.5,2,2.5,3 ...
4. the amplitude rule method of blind focus earthquake wave field according to claim 2, is characterized in that: described amplitude regularization factor p gets p=2.5.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410477842.6A CN105445790A (en) | 2014-09-18 | 2014-09-18 | Amplitude regularization method of blind source seismic wave field |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410477842.6A CN105445790A (en) | 2014-09-18 | 2014-09-18 | Amplitude regularization method of blind source seismic wave field |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN105445790A true CN105445790A (en) | 2016-03-30 |
Family
ID=55556202
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410477842.6A Pending CN105445790A (en) | 2014-09-18 | 2014-09-18 | Amplitude regularization method of blind source seismic wave field |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN105445790A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107884828A (en) * | 2016-09-30 | 2018-04-06 | 中国石油化工股份有限公司 | It is a kind of in spatial frequency domain based on the theoretical terrible ripple drawing method of Green |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100274493A1 (en) * | 2009-04-10 | 2010-10-28 | Harold Szu | Nowcast of Natural Resource and Calamity Warning |
| CN102866425A (en) * | 2012-09-17 | 2013-01-09 | 中国石油大学(华东) | Blind source seismic signal stable-superposition model based blind separation method |
| CN103364826A (en) * | 2012-03-30 | 2013-10-23 | 中国石油化工股份有限公司 | An earthquake blind source deconvolution method based on independent component analysis |
-
2014
- 2014-09-18 CN CN201410477842.6A patent/CN105445790A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20100274493A1 (en) * | 2009-04-10 | 2010-10-28 | Harold Szu | Nowcast of Natural Resource and Calamity Warning |
| CN103364826A (en) * | 2012-03-30 | 2013-10-23 | 中国石油化工股份有限公司 | An earthquake blind source deconvolution method based on independent component analysis |
| CN102866425A (en) * | 2012-09-17 | 2013-01-09 | 中国石油大学(华东) | Blind source seismic signal stable-superposition model based blind separation method |
Non-Patent Citations (4)
| Title |
|---|
| 张林颖 等: ""盲源分离算法在混合震动信号分离中的应用"", 《计算机技术与应用》 * |
| 张涛 等: ""地震数据的空间采样与能量联合规则化方法及应用"", 《新疆石油地质》 * |
| 张贤达: "《信号分析与处理》", 31 December 2011, 清华大学出版社 * |
| 马建仓 等: "《盲信号处理》", 30 June 2006, 国防工业出版社-北京 * |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN107884828A (en) * | 2016-09-30 | 2018-04-06 | 中国石油化工股份有限公司 | It is a kind of in spatial frequency domain based on the theoretical terrible ripple drawing method of Green |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103091714B (en) | A kind of self-adaptation surface wave attenuation method | |
| CN104007469A (en) | Weak seismic signal reconstruction method based on curvelet transform | |
| CN103399348B (en) | Based on the Denoising of Seismic Data method of Shearlet conversion | |
| CN102879825B (en) | High power pulse noise detection and suppression method for vibroseis seismic data | |
| CN106291684B (en) | A kind of seismic response of blind focus earthquake wave field restores and virtual source trace gather construction method | |
| CN103645497B (en) | Emergence angle based multi-component wave field separation method | |
| CN103869364A (en) | Multiple wave suppression method based on dual parabolic Radon transformation | |
| CN106646637A (en) | Method for removing peak noise in nuclear magnetism signal | |
| CN104614769A (en) | Beam-forming filtering method for suppressing seismic surface waves | |
| CN107179550B (en) | A kind of seismic signal zero phase deconvolution method of data-driven | |
| Wang et al. | Seismic data denoising for complex structure using BM3D and local similarity | |
| CN111458749B (en) | Surface wave and body wave separation method and system applied to passive source seismic exploration | |
| CN107402405B (en) | Quiet phase virtual source trace gather construction method | |
| CN105259579A (en) | A high-amplitude shielding layer rejecting method based on seismic data instantaneous attributes | |
| CN104330826A (en) | A method for removing various noises under the condition of complex surface | |
| CN102736109A (en) | Method for de-noising, correcting and superposing CRP (Common Reflection Point) gather | |
| CN106526678A (en) | Reflection acoustic logging wave field separation method and device | |
| CN104597502A (en) | Novel petroleum seismic exploration data noise reduction method | |
| CN104216010A (en) | Method for increasing quality of seismic data by using harmonic waves of controllable seismic focus | |
| CN105319593A (en) | Combined denoising method based on curvelet transform and singular value decomposition | |
| CN104635264B (en) | The processing method of earthquake data before superposition and equipment | |
| CN105301648A (en) | Method of acquiring common reflection surface stacking parameters | |
| CN105510975A (en) | Method and device for improving signal-to-noise ratio of seismic data | |
| CN103412325A (en) | Method for separating and denoising wave fields of p-waves and s-waves of earthquake | |
| CN102338886A (en) | Polarization filtering method for effectively attenuating surface waves in three-component seismic record |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20160330 |
|
| WD01 | Invention patent application deemed withdrawn after publication |