CN107450103A - A kind of terrible ripple drawing method based on boundary integral inverse operator - Google Patents
A kind of terrible ripple drawing method based on boundary integral inverse operator Download PDFInfo
- Publication number
- CN107450103A CN107450103A CN201710799722.1A CN201710799722A CN107450103A CN 107450103 A CN107450103 A CN 107450103A CN 201710799722 A CN201710799722 A CN 201710799722A CN 107450103 A CN107450103 A CN 107450103A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- terrible
- ripple
- terrible ripple
- 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
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/38—Seismology; Seismic or acoustic prospecting or detecting specially adapted for water-covered areas
Abstract
The present invention relates to a kind of terrible ripple drawing method based on boundary integral inverse operator, it is characterised in that comprises the following steps：1) according to the Boundary integral formulation of seimic wave propagation and earthquake depth data, the boundary integral operator formula of terrible ripple is calculated；2) the terrible ripple boundary integral operator formula obtained according to step 1), using terrible wave operator and the coupled relation of terrible ripple inverse operator is suppressed, obtains suppressing the calculation formula of terrible ripple inverse operator；3) geological data according to the calculation formula of the terrible ripple inverse operator of compacting obtained in step 2) and containing there's something fishy ripple, is calculated the geological data after suppressing terrible ripple.The seismic data that the present invention can obtain according to different acquisition mode, obtain suppressing the geological data after terrible ripple, and computational methods efficiency high, thus the present invention can be widely applied to during marine seismic prospectiong.
Description
Technical field
The present invention relates to the Processing Seismic Data in marine seismic prospectiong field, more particularly to one kind to be based on boundary integral
The terrible ripple drawing method of inverse operator.
Background technology
At sea during seismic prospecting, due to the influence of terrible ripple, marine streamer seismic data collection wants to obtain relatively low
The seismic data of frequency, it is general to use compared with bathyseism source and larger towing cable sinking depth, in order to obtain upper frequency seismic data one
As use shallower focus and smaller towing cable sinking depth.In order to study the geological condition of mid-deep reservoir, the target of seismic prospecting
The target imaging of mid-deep strata should be obtained, while requires high-resolution seismic data.This just needs the influence for eliminating terrible ripple, with up to
To extension seismic band so as to obtain preferably in deep formation at target locations geological information purpose.
With the improvement of earthquake-capturing equipment, seismic processes in recent years, many new earthquakes are have developed both at home and abroad
Acquisition method, try to eliminate the influence of terrible ripple by earthquake-capturing and earthquake information processing.Relative to the place of common seismic data
Flow is managed, one of subject matter existing for broad-band teleseismic information Processing Technology is how to suppress or weaken terrible ripple, so as to obtain more
Wide band geological data.The influence of terrible ripple is put forward to be mitigated or eliminated for different acquisition methods to obtain wideband
Seismic data processing technology, it can obtain to finer seismic data and high-resolution Seismic section, but by
In the change of acquisition mode, acquisition cost is materially increased, while for the terrible wave pressure side of new acquisition mode
Method complex disposal process, difficulty are larger, it is impossible to the seismic data applied to the collection of former wide variety of conventional.
The content of the invention
In view of the above-mentioned problems, it is an object of the invention to provide a kind of offshore earthquake that can be suppressed different acquisition mode and obtain
The terrible ripple drawing method based on boundary integral inverse operator of terrible ripple in data.
To achieve the above object, the present invention takes following technical scheme：A kind of terrible wave pressure based on boundary integral inverse operator
Method processed, it is characterised in that comprise the following steps：1) according to the Boundary integral formulation of seimic wave propagation and earthquake depth data,
The boundary integral operator of terrible ripple is calculated；2) the boundary integral operator of the terrible ripple obtained according to step 1) and terrible wave operator and
The coupled relation of the inverse operator of terrible ripple is suppressed, obtains suppressing the inverse operator of terrible ripple；3) according to the inverse operator of the obtained terrible ripple of compacting
And the geological data containing there's something fishy ripple, the geological data after suppressing terrible ripple is calculated.
In the step 1), the boundary integral operator computational methods of terrible ripple comprise the following steps：1.1) according to the ground of collection
The depth of focus and cable depth information in data are shaken, obtains the geometrical relationship figure of earthquake-capturing, and then obtain being connect on cable
The relation of seismic signal of the seismic signal received with reaching cable wave detector respectively from underground and from sea surface reflection；1.2)
According to the geometrical relationship figure of earthquake-capturing in the border condition of continuity and step 1.1) of free interface, obtain receiving on seawater face
The relation of seismic signal of the seismic signal arrived with reaching seawater face respectively from underground and from sea surface reflection；1.3) according to ripple
The boundary element integral formula of dynamic equation, obtain seismic signal received on cable and to reach seawater face from underground and come
The relation of the Boundary integral formulation of the seismic signal in seawater face is reached from sea surface reflection；1.4) according to step 1.2) and step 1.3)
Obtained relation, seismic signal received on cable is expressed as reaching to the border of the seismic signal in seawater face from underground
Integrated form, i.e., the expression formula of seismic signal and the boundary integral operator of terrible ripple without terrible ripple.
In the step 1.1), received seismic signal u (r) is to reach cable wave detector from underground on cable
Signal ut(r) the signal u of cable wave detector and from sea surface reflection is reachedg(r) sum, i.e.,：
U (r)=ut(r)+ug(r),
R is the position vector where cable in formula.
In the step 1.2), received seismic signal u (r on seawater face0) it is to reach seawater face from underground
Signal ut(r0) with reaching the signal u in seawater face from sea surface reflectiong(r0) sum and equal to 0, i.e.,：
u(r0)=ut(r0)+ug(r0)=0；
In formula, r0For the position vector where seawater face.
In the step 1.3), on cable received seismic signal with to reach seawater face and from sea from underground
The relation that face reflexes to the Boundary integral formulation of the seismic signal up to seawater face is：
In the step 1.3), received seismic signal is expressed as reaching the earthquake in seawater face from underground on cable
The boundary integral form of signal is：
Wherein, wherein A (r, r0) be terrible ripple boundary integral operator formula, i.e.,：
In the step 2), the calculation formula of the inverse operator of the obtained terrible ripple of compacting is：
B(r,r0)=A-(r,r0)；
In formula, A-(r,r0) for terrible ripple boundary integral operator discretization after inverse of a matrix matrix.
In the step 3), according to the calculation formula for the inverse operator for suppressing terrible ripple and the earthquake record containing there's something fishy ripple, meter
Calculate the seismic signal u after suppressing terrible ripplet(r0) formula be：
ut(r0)=B (r, r0)u(r)。
For the present invention due to taking above technical scheme, it has advantages below：1st, the present invention required for seismic data only
It is related to sampling depth data, reduces the requirement to seismic data, effectively reduce acquisition cost.2nd, the present invention provides according to earthquake
Material obtains the boundary integral inverse operator for suppressing terrible ripple, and boundary integral inverse operator then is applied into the geological data containing terrible ripple, phase
Than the method in Wave field analysis in the prior art etc. using optimization algorithm, fast more than 80 times of calculating speed of the present invention, effectively carry
High computational efficiency.3rd, the present invention is according to wave field principle of stacking, it is contemplated that the fluctuation characteristic and amplitude of Seismic Wave Propagation Equations become
Change feature, avoid the ray path based on seimic wave propagation and carry out shadow of the terrible wave pressure processing method to seismic data amplitude
Ring, the accuracy for follow-up reservoir prediction and hydrocarbon indication provides more preferable guarantee.The present invention can be widely applied to sea
During foreign seismic prospecting.
Brief description of the drawings
Fig. 1 is the terrible wave pressure method flow diagram of the invention based on boundary integral inverse operator；
Fig. 2 is the seismic signal schematic diagram of cable reception；
Fig. 3 is the effect contrast figure before and after the geological data application present invention containing there's something fishy ripple.
Embodiment
The present invention is described in detail with reference to the accompanying drawings and examples.
As shown in figure 1, a kind of terrible ripple drawing method based on boundary integral inverse operator provided by the invention, including following step
Suddenly：
1) according to the Boundary integral formulation of seimic wave propagation and the earthquake depth data of collection, the side of terrible ripple is calculated
Boundary's integral operator formula；Wherein, earthquake depth data include the depth of focus and cable depth information；
2) the boundary integral operator formula of the terrible ripple obtained according to step 1), using terrible wave operator and terrible ripple inverse operator is suppressed
Coupled relation, obtain suppressing the calculation formula of the inverse operator of terrible ripple；
3) earthquake record according to the calculation formula of the terrible ripple inverse operator of compacting obtained in step 2) and containing there's something fishy ripple, i.e.,
The geological data after suppressing terrible ripple can be calculated.
Above-mentioned steps 1) in, according to the Boundary integral formulation of seimic wave propagation and the seismic basis data of collection, calculate
The method of the boundary integral operator of terrible ripple, comprises the following steps：
1.1) depth of focus in the geological data of collection and cable depth information, the geometry for obtaining earthquake-capturing close
System's figure, and then seismic signal received on cable is obtained with reaching cable detection respectively from underground and from sea surface reflection
The relation of the seismic signal of device.
In offshore towing rope seismic exploration, in seismic data acquisition process, earthquake ship can record earthquake class and notify breath, earthquake
The depth of focus and cable depth that class notifies when collection can be included in breath, using the depth of focus and cable depth information, can be obtained
To the geometrical relationship figure of earthquake-capturing.
As shown in Fig. 2 seismic signal u (r) received on cable is the signal that cable wave detector is reached from underground
ut(r) the signal u of cable wave detector and from sea surface reflection is reachedg(r) sum, i.e.,：
U (r)=ut(r)+ug(r) (1)
R is the position vector where cable in formula.
1.2) according to the geometrical relationship figure of earthquake-capturing in the border condition of continuity and step 1.1) of free interface, obtain
The pass of seismic signal of the seismic signal received on seawater face with reaching seawater face respectively from underground and from sea surface reflection
System.
It can be seen from the geometrical relationship figure of earthquake-capturing, received seismic signal u (r on seawater face0) it is from ground
The lower signal u for reaching seawater facet(r0) with reaching the signal u in seawater face from sea surface reflectiong(r0) sum, i.e.,：
u(r0)=ut(r0)+ug(r0) (2)
In formula, r0For the position vector where seawater face.
Assuming that seawater face is free interface, according to the border condition of continuity of free interface, we can draw：
u(r0)=ut(r0)+ug(r0)=0 (3)
1.3) according to the boundary element integral formula of wave equation, obtain seismic signal received on cable with from
Underground reaches seawater face and the relation of the Boundary integral formulation of the seismic signal in seawater face is reached from sea surface reflection.
Wherein, the boundary element integral formula of wave equation is：
For two-dimentional Seawater, For border normal derivative.WhereinFor Hunk
That function, k=ω/v are background wave number, and ω is angular frequency, and v is seawater speed.
Formula (1), (3) joint can be obtained：
1.4) relation obtained according to step 1.2) and step 1.3), seismic signal received on cable is expressed as
The boundary integral form of the seismic signal in seawater face is reached from underground, i.e., the border of the seismic signal without terrible ripple and terrible ripple is accumulated
Divide the expression formula of operator.
Formula (3) is substituted into formula (5) to obtain：
Wherein A (r, r0) be terrible ripple boundary integral operator formula, i.e.,：
Above-mentioned steps 2) in, it is right first when the inverse operator for suppressing terrible ripple is calculated according to the boundary integral operator of terrible ripple
Formula (7) carries out discrete integration calculating, by boundary integral operator A (r, the r of terrible ripple0) be converted to matrix form A (r, r0), afterwards
Ask for the matrix A (r, r0) inverse matrix B (r, r0), i.e.,
B(r,r0)=A-(r,r0) (8)
B(r,r0) it is the calculation formula for suppressing terrible ripple inverse operator.
Above-mentioned steps 3) in, according to the calculation formula for the inverse operator for suppressing terrible ripple and the geological data containing there's something fishy ripple, meter
Calculate the seismic signal u after suppressing terrible ripplet(r0) when, it is only necessary to terrible ripple inverse operator and the u of the geological data containing there's something fishy ripple will be suppressed
(r) product is done, i.e.,：
ut(r0)=B (r, r0)u(r) (9)
The effect of terrible ripple drawing method of the present invention based on boundary integral inverse operator is carried out below by specific embodiment
Describe in detail：
As shown in figure 3, one two-dimension earthquake survey line of selection, the two-dimension earthquake survey line seawater mean depth more than 1000m,
During earthquake-capturing, the depth of focus used is 8m, cable depth 13m, and the seismic data has relatively strong terrible ripple, due to ghost
The influence of ripple, the wave group feature of seismic data are covered by terrible ripple, and after carrying out terrible wave pressure system using the present invention, seismic wave groups feature becomes
Clearly, Geologic Structure Feature is more easy to identify.
The various embodiments described above are merely to illustrate the present invention, wherein the structure of each part, connected mode and manufacture craft etc. are all
It can be varied from, every equivalents carried out on the basis of technical solution of the present invention and improvement, should not exclude
Outside protection scope of the present invention.
Claims (8)
1. a kind of terrible ripple drawing method based on boundary integral inverse operator, it is characterised in that comprise the following steps：
1) according to the Boundary integral formulation of seimic wave propagation and earthquake depth data, the boundary integral that terrible ripple is calculated is calculated
Son；
2) coupling of the boundary integral operator of the terrible ripple obtained according to step 1) and terrible wave operator and the inverse operator of the terrible ripple of compacting
Relation, obtain suppressing the inverse operator of terrible ripple；
3) geological data according to the inverse operator of the obtained terrible ripple of compacting and containing there's something fishy ripple, is calculated the ground after suppressing terrible ripple
Shake data.
A kind of 2. terrible ripple drawing method based on boundary integral inverse operator as claimed in claim 1, it is characterised in that：The step
It is rapid 1) in, the boundary integral operator computational methods of terrible ripple comprise the following steps：
1.1) depth of focus in the geological data of collection and cable depth information, obtain the geometrical relationship of earthquake-capturing
Figure, and then seismic signal received on cable is obtained with reaching cable wave detector respectively from underground and from sea surface reflection
Seismic signal relation；
1.2) according to the geometrical relationship figure of earthquake-capturing in the border condition of continuity and step 1.1) of free interface, seawater is obtained
The relation of seismic signal of the seismic signal received on face with reaching seawater face respectively from underground and from sea surface reflection；
1.3) according to the boundary element integral formula of wave equation, seismic signal received on cable is obtained and with from underground
Reach seawater face and the relation of the Boundary integral formulation of the seismic signal in seawater face is reached from sea surface reflection；
1.4) relation obtained according to step 1.2) and step 1.3), seismic signal received on cable is expressed as coming from
Underground reaches the boundary integral form of the seismic signal in seawater face, i.e., the boundary integral of the seismic signal without terrible ripple and terrible ripple is calculated
The expression formula of son.
A kind of 3. terrible ripple drawing method based on boundary integral inverse operator as claimed in claim 2, it is characterised in that：The step
It is rapid 1.1) in, received seismic signal u (r) is the signal u that cable wave detector is reached from underground on cablet(r) with coming
The signal u of cable wave detector is reached from sea surface reflectiong(r) sum, i.e.,：
U (r)=ut(r)+ug(r),
R is the position vector where cable in formula.
A kind of 4. terrible ripple drawing method based on boundary integral inverse operator as claimed in claim 2, it is characterised in that：The step
It is rapid 1.2) in, received seismic signal u (r on seawater face0) it is the signal u that seawater face is reached from undergroundt(r0) with coming from
Sea surface reflection reaches the signal u in seawater faceg(r0) sum and equal to 0, i.e.,：
u(r0)=ut(r0)+ug(r0)=0；
In formula, r0For the position vector where seawater face.
A kind of 5. terrible ripple drawing method based on boundary integral inverse operator as claimed in claim 2, it is characterised in that：The step
It is rapid 1.3) in, on cable received seismic signal with to reach seawater face from underground and reach seawater from sea surface reflection
The relation of the Boundary integral formulation of the seismic signal in face is：
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>+</mo>
<msub>
<mi>u</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mo>&Integral;</mo>
<mi>&Gamma;</mi>
</msub>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mo>+</mo>
<msub>
<mi>u</mi>
<mi>g</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mo>&Integral;</mo>
<mi>&Gamma;</mi>
</msub>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mo>.</mo>
</mrow>
</mrow>
A kind of 6. terrible ripple drawing method based on boundary integral inverse operator as claimed in claim 2, it is characterised in that：The step
It is rapid 1.3) in, received seismic signal is expressed as reaching the boundary integral of the seismic signal in seawater face from underground on cable
Form is：
<mrow>
<mi>u</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<msub>
<mi>u</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<msub>
<mi>u</mi>
<mi>t</mi>
</msub>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>{</mo>
<msub>
<mo>&Integral;</mo>
<mi>&Gamma;</mi>
</msub>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<mi>&Gamma;</mi>
</msub>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mo>}</mo>
<mo>,</mo>
</mrow>
Wherein, wherein A (r, r0) be terrible ripple boundary integral operator formula, i.e.,：
<mrow>
<mi>A</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mo>{</mo>
<msub>
<mo>&Integral;</mo>
<mi>&Gamma;</mi>
</msub>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<mi>r</mi>
<mo>,</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
<mi>d</mi>
<mi>r</mi>
</mrow>
<mo>-</mo>
<msub>
<mo>&Integral;</mo>
<mi>&Gamma;</mi>
</msub>
<mrow>
<mo>&lsqb;</mo>
<mfrac>
<mrow>
<mo>&part;</mo>
<mi>G</mi>
<mrow>
<mo>(</mo>
<msub>
<mi>r</mi>
<mn>0</mn>
</msub>
<mo>,</mo>
<mi>r</mi>
<mo>)</mo>
</mrow>
</mrow>
<mrow>
<mo>&part;</mo>
<mi>n</mi>
</mrow>
</mfrac>
<mo>&rsqb;</mo>
</mrow>
<mo>}</mo>
<mo>.</mo>
</mrow>
A kind of 7. terrible ripple drawing method based on boundary integral inverse operator as claimed in claim 1, it is characterised in that：The step
It is rapid 2) in, the calculation formula of the inverse operator of the obtained terrible ripple of compacting is：
B(r,r0)=A-(r,r0)；
In formula, A-(r,r0) for terrible ripple boundary integral operator discretization after inverse of a matrix matrix.
A kind of 8. terrible ripple drawing method based on boundary integral inverse operator as claimed in claim 1, it is characterised in that：The step
It is rapid 3) in, according to the calculation formula for the inverse operator for suppressing terrible ripple and the earthquake record containing there's something fishy ripple, calculate after suppressing terrible ripple
Seismic signal ut(r0) formula be：
ut(r0)=B (r, r0)u(r)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710799722.1A CN107450103B (en) | 2017-09-07 | 2017-09-07 | Ghost wave compression method based on boundary integral inverse operator |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710799722.1A CN107450103B (en) | 2017-09-07 | 2017-09-07 | Ghost wave compression method based on boundary integral inverse operator |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107450103A true CN107450103A (en) | 2017-12-08 |
| CN107450103B CN107450103B (en) | 2020-04-28 |
Family
ID=60495075
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710799722.1A Active CN107450103B (en) | 2017-09-07 | 2017-09-07 | Ghost wave compression method based on boundary integral inverse operator |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107450103B (en) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114167501A (en) * | 2021-11-04 | 2022-03-11 | 中国海洋大学 | FTX processing method, storage medium and processing terminal for suppressing strong-energy sea ghost waves |
| CN114355450A (en) * | 2022-03-21 | 2022-04-15 | 中国石油大学(华东) | Method, system and equipment for suppressing ghost waves of marine plough type cable full waveform inversion |
| CN114442173A (en) * | 2020-11-02 | 2022-05-06 | 中国石油化工股份有限公司 | Computer program product and method for predicting and eliminating multiples in the beam domain |
| CN116774280A (en) * | 2023-06-25 | 2023-09-19 | 中海石油(中国)有限公司深圳分公司 | Ghost wave pressure formulation quality control method and device, electronic equipment and storage medium |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104536045A (en) * | 2015-01-16 | 2015-04-22 | 中国海洋石油总公司 | Ghost wave compression method based on wavelet processing |
| CN106199699A (en) * | 2016-06-30 | 2016-12-07 | 恒泰艾普（北京）能源科技研究院有限公司 | Utilize the method that transfer matrix method removes ghost ripple |
| CN106526677A (en) * | 2016-10-26 | 2017-03-22 | 中海石油(中国)有限公司 | Marine self-adaptive ghost reflection-suppressing broadband reverse time migration imaging method |
-
2017
- 2017-09-07 CN CN201710799722.1A patent/CN107450103B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104536045A (en) * | 2015-01-16 | 2015-04-22 | 中国海洋石油总公司 | Ghost wave compression method based on wavelet processing |
| CN106199699A (en) * | 2016-06-30 | 2016-12-07 | 恒泰艾普（北京）能源科技研究院有限公司 | Utilize the method that transfer matrix method removes ghost ripple |
| CN106526677A (en) * | 2016-10-26 | 2017-03-22 | 中海石油(中国)有限公司 | Marine self-adaptive ghost reflection-suppressing broadband reverse time migration imaging method |
Non-Patent Citations (9)
| Title |
|---|
| A.B. WEGLEIN ET AL.: "New approaches to deghosting towed-streamer and ocean-bottom pressure measurements", 《SEG INT"L EXPOSITION AND 72ND ANNUAL MEETING》 * |
| JAMES D. MAYHAN ET AL.: "First application of Green’s theorem derived source and receiver deghosting on deep water Gulf of Mexico synthetic (SEAM) and field data", 《SEG LAS VEGAS 2012 ANNUAL MEETING》 * |
| LIN TANG ET AL.: "Using Green’s theorem to satisfy data requirements of multiple removal methods: The impact of acquisition design", 《SEG HOUSTON 2013 ANNUAL MEETING》 * |
| TAO JIE ET AL.: "A deghosting method based on source wavelet processing", 《SEG NEW ORLEANS ANNUAL MEETING》 * |
| 李洪建等: "基于格林函数理论的波场预测和鬼波压制方法", 《地球物理学报》 * |
| 杨金龙等: "基于格林理论的鬼波压制方法及其应用", 《石油物探》 * |
| 王芳芳等: "基于逆散射级数法的鬼波压制方法", 《地球物理学报》 * |
| 管西竹等: "基于波动方程的上下缆地震数据鬼波压制方法研究", 《地球物理学报》 * |
| 陈宝书等: "基于确定性子波处理的鬼波压制方法", 《中国海上油气》 * |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114442173A (en) * | 2020-11-02 | 2022-05-06 | 中国石油化工股份有限公司 | Computer program product and method for predicting and eliminating multiples in the beam domain |
| CN114442173B (en) * | 2020-11-02 | 2023-05-02 | 中国石油化工股份有限公司 | Computer program product and method for predicting and eliminating multiple in beam domain |
| CN114167501A (en) * | 2021-11-04 | 2022-03-11 | 中国海洋大学 | FTX processing method, storage medium and processing terminal for suppressing strong-energy sea ghost waves |
| CN114167501B (en) * | 2021-11-04 | 2023-06-30 | 中国海洋大学 | FTX processing method for suppressing strong energy sea surface ghost wave, storage medium and processing terminal |
| CN114355450A (en) * | 2022-03-21 | 2022-04-15 | 中国石油大学(华东) | Method, system and equipment for suppressing ghost waves of marine plough type cable full waveform inversion |
| US11709286B1 (en) | 2022-03-21 | 2023-07-25 | China University Of Petroleum (East China) | Method, system, and device for full-waveform inversion deghosting of marine variable depth streamer data acquisition |
| CN116774280A (en) * | 2023-06-25 | 2023-09-19 | 中海石油(中国)有限公司深圳分公司 | Ghost wave pressure formulation quality control method and device, electronic equipment and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107450103B (en) | 2020-04-28 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| AU2012202582B2 (en) | Methods and apparatus for seismic exploration using pressure changes caused by sea-surface variations | |
| CN107450103A (en) | A kind of terrible ripple drawing method based on boundary integral inverse operator | |
| US9829592B2 (en) | Seismic imaging with visco-acoustic reverse-time migration using pseudo-analytical method | |
| US9651695B2 (en) | Construction and application of angle gathers from three-dimensional imaging of multiples wavefields | |
| US10935680B2 (en) | Generating geophysical images using directional oriented wavefield imaging | |
| US20130077832A1 (en) | Method for analyzing seismic data | |
| CN1529823A (en) | Method and apparatus for determining orientation of seismic receiver during seismic survey | |
| EP3059615A1 (en) | Amplitude-versus-angle analysis for quantative interpretation | |
| US10459100B2 (en) | Survey techniques using streamers at different depths | |
| KR20180067650A (en) | FWI model domain angular stacks with amplitude preservation | |
| NO20141519A1 (en) | Methods and systems for optimizing the generation of seismic images | |
| CN103645503A (en) | Three-dimensional time domain illumination analysis and amplitude compensation method | |
| CN106896409A (en) | A kind of varying depth cable ghost ripple drawing method based on wave equation boundary values inverting | |
| EP3136130B1 (en) | Wavefield interpolation and regularization in imaging of multiple reflection energy | |
| CN104570116A (en) | Geological marker bed-based time difference analyzing and correcting method | |
| CN110031900A (en) | Diffraction wave imaging method, device, electronic equipment and storage medium | |
| AU2011380936B2 (en) | Seismic imaging systems and methods employing correlation-based stacking | |
| Wu et al. | Side-scan sonar and sub-bottom profiler surveying | |
| CN109975873B (en) | Method and system for removing low-frequency noise by reverse time migration imaging | |
| Xiong et al. | An effective method for submarine pipeline inspection using three-dimensional (3D) models constructed from multisensor data fusion | |
| CN113777653B (en) | Synthetic aperture type shallow seismic section exploration method and system | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
| US11092710B2 (en) | Inversion techniques using streamers at different depths | |
| CN115184987B (en) | Stratum inclination angle information extraction method and device, server and storage medium | |
| CN104375176B (en) | Method for determining in-phase superposition highest cutoff frequency through first arrival time error |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |