TARGET-ORIENTED SEISMIC ACQUISITION METHOD AND APPARATUS, MEDIUM AND DEVICE

Abstract

The present invention relates to a target-oriented seismic acquisition method and apparatus, a medium and a device. The target-oriented seismic acquisition method comprises the steps of: giving parameters of an initial velocity model and a three-dimensional seismic layout aiming to an underground target position; conducting wave field continuation and focusing analysis on the three-dimensional seismic layout, and calculating distribution of seismic energy on the ground in an underground target region; conducting normalization processing on distribution of the seismic energy on the ground, and then conducting level partitioning to obtain a primary energy region and a secondary energy region; adding the number of shot points in the primary energy region to achieve target-oriented acquisition, and obtaining a target-oriented inhomogeneous laying acquired data imaging result. By using the method of the present invention, automatic feedback adjustment on excitation and receiving sites and parameters thereof is achieved.

Claims

1. A target-oriented seismic acquisition method, comprising the steps of: giving parameters of an initial velocity model and a three-dimensional seismic layout aiming to an underground target position; conducting wave field continuation and focusing analysis on the three-dimensional seismic layout, and calculating distribution of seismic energy on the ground in an underground target region; conducting normalization processing on distribution of the seismic energy on the ground, and then conducting level partitioning to obtain a primary energy region and a secondary energy region; adding the number of shot points in the primary energy region to achieve target-oriented acquisition, and obtaining a target-oriented inhomogeneous laying acquired data imaging result; and comparing a conventional acquired data imaging result with the target-oriented inhomogeneous laying acquired data imaging result.

2. The target-oriented seismic acquisition method according to claim 1, wherein the conducting wave field continuation and focusing analysis on the three-dimensional seismic layout, and calculating distribution of the seismic energy on the ground in the underground target region specifically comprises the steps of: conducting large-step continuation by using FFD, and obtaining wave fields in various small layers between large steps by combining with Born-Kirchhoff interpolation at the same time; after completing wave field continuation, calculating the focusing attribute of receivers and the focusing attribute of shot points of the three-dimensional seismic layout, and then quantitatively analyzing an expected resolution and an AVP attribute of the whole three-dimensional seismic layout; and calculating distribution of the seismic energy on the ground in the underground target region.

3. The target-oriented seismic acquisition method according to claim 1, wherein the conducting normalization processing on distribution of the seismic energy on the ground, and then conducting level partitioning specifically comprises the steps of: with an energy drop 1/e as a limit, partitioning the seismic energy into the primary energy region if the energy drop is above the limit and the seismic energy into the secondary energy region if the energy drops to the energy drop ½ below the limit.

4. A target-oriented seismic acquisition apparatus, comprising: a first processing unit for giving the parameters of the initial velocity model and the three-dimensional seismic layout aiming to the underground target position; a second processing unit for conducting wave field continuation and focusing analysis on the three-dimensional seismic layout and calculating distribution of the seismic energy on the ground in the underground target region; a third processing unit for conducting normalization processing on distribution of the seismic energy on the ground and then conducting level partitioning to obtain the primary energy region and the secondary energy region; a fourth processing unit for adding the number of shot points in the primary energy region to achieve target-oriented acquisition and obtaining the target-oriented inhomogeneous laying acquired data imaging result; and a fifth processing unit for comparing the conventional acquired data imaging result with the target-oriented inhomogeneous laying acquired data imaging result.

5. A computer readable storage medium, wherein a computer program is stored in the computer readable storage medium; and when the computer program is executed by a processor, the steps of the target-oriented seismic acquisition method according to claim lis implemented.

Description

BRIEF DESCRIPTION OF DRAWINGS

[0029] FIG. 1 is a flowchart of a target-oriented seismic acquisition method according to an embodiment of the present invention.

[0030] FIG. 2 is a schematic diagram of focal beam analysis in complex media according to the embodiment of the present invention.

[0031] FIG. 3 is a schematic diagram of algorithm combining Fourier finite difference (FFD) large-step with Born-Kirchhoff small-step continuation according to the embodiment of the present invention.

[0032] FIG. 4 shows distribution of seismic energy on the ground in an underground target region according to the embodiment of the present invention.

[0033] FIG. 5 shows schematic diagrams of energy distribution according to the embodiment of the present invention, wherein a is a schematic diagram of energy distribution corresponding to conventional shot points, and b is a schematic diagram of energy distribution corresponding to shot points after being added/adjusted in position.

[0034] FIG. 6 shows imaging section comparison on data obtained after target-oriented acquisition (b) according to the embodiment of the present invention and traditional acquisition (a)-taking a certain target region in Zhongyuan oil field as an example.

[0035] FIG. 7 shows imaging section comparison on data obtained after target-oriented acquisition (b) according to the embodiment of the present invention and traditional acquisition (a)-taking a certain target region in Dagang oil field as an example.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0036] As shown in FIG. 1, the target-oriented seismic acquisition method provided by the present invention comprises the steps of:

[0037] 1) Giving parameters of an initial velocity model and a three-dimensional seismic layout aiming to an underground target position.

[0038] 2) Conducting large-step continuation by using FFD, and obtaining wave fields in various small layers between large steps by combining with Born-Kirchhoff interpolation at the same time (as shown in FIG. 2, it is made every 5 layers). The FFD large-step continuation algorithm belongs to the prior art (referring to the Patent Application No. 201010557332.1, entitled with “Migration Noise Acquisition Method And Apparatus For Three-Dimensional Seismic Layout” in detail), and thus FFD large-step continuation will not be described herein. As the calculated amount is huge by purely using the FFD wave field continuation technology, the present invention proposes the technology combining FFD with Born-Kirchhoff interpolation, and a wave field u(r) is obtained by Born-Kirchhoff interpolation deduced by a Lippmann-Schwinger integral formula:

[00001]$u\left(r\right)={\int }_z^{z+\Delta z}\frac{i{k}_0}{2\pi r}\exp \left(i{k}_0\Delta z\left(n-1\right)\right)\exp \left(i{k}_0\Delta r\right)u_s\left(r\right)dr$

[0039] wherein

r=√{square root over ((x.sub.1−x.sub.2).sup.2+(y.sub.1−y.sub.2).sup.2+Δz.sup.2)}

[0040] In the formula, n is a reflection coefficient; k.sub.0 is a background wave number; i is a pure imaginary number; the position of a target point before continuation is (x.sub.1,y.sub.2); the position of the target point after continuation is (x.sub.2,y.sub.2); an ordinate of the target point before continuation is z; a difference between ordinates of the target point before and after continuation is Δz; and u.sub.s(r) is a wave field before continuation.

[0041] As shown in FIG. 3, the wave field propagates to the receiver at the position zo on the ground from a seismic source grid point at the target point r.sub.t through an inhomogeneous medium v(x,y,z). The focusing property of the observation system at the target point is estimated by inverse continuation of an observed wave field. Z.sub.0 is an ordinate of the ground, and z.sub.m is an ordinate of the target point.

[0042] After wave field continuation is completed, the focusing property of the receivers and the focusing property of the shot points of the three-dimensional seismic layout are calculated (the detailed calculation or analysis process of the focusing property of the receivers and the focusing property of the shot points refers to the Patent Application No. 201710081114.7, entitled with “Quantitative Analysis Evaluation Method For Vertical Resolution Of Three-Dimensional Seismic Layout”). Therefore, an expected resolution and an AVP attribute of the whole observation system in a spatial frequency domain and a Radon domain are obtained respectively, and distribution of ground energy in the underground target region is given (as shown in FIG. 4).

[0043] In this step, the whole continuation based on the three-dimensional seismic layout is completed through Fourier finite difference (FFD) large-step continuation and Born-Kirchhoff operator small-step wave field interpolation; and FFD large-step continuation belongs to the prior art and will not be described herein. Compared with a Kirchhoff integral formula under the homogeneous condition, for a Born-Kirchhoff integral formula under incident waves at a small angle, a slowness disturbance correction term is added; and therefore, an operator can be adapted to a lateral variation of a velocity of a medium. After wave field continuation, focal beam analysis is conducted through the receivers, the shot points and the underground target body to solve a resolution function and an amplitude versus ray-parameter (AVP) function. Focal beam analysis belongs to the prior art (referring to the Patent Application No. 201710081114.7, entitled with “Quantitative Analysis Evaluation Method For Vertical Resolution Of Three-Dimensional Seismic Layout” in detail), which will not be described herein.

[0044] 3) Conducting normalization processing on distribution of the seismic energy on the ground, and then conducting level partitioning. With an energy drop 1/e as a limit, the seismic energy is partitioned into the primary energy region if the energy drop is above the limit; and the seismic energy is partitioned into the secondary energy region if the energy drops to the energy drop ½ below the limit.

[0045] 4) In addition to the conventionally designed shot points and receivers, adding the number of the shot points in the primary energy region to achieve target-oriented acquisition (as shown in FIG. 5), wherein the design solution serves as a basis for improving the design on the three-dimensional seismic layout.

[0046] Firstly, the positions of the shot points and the receivers are given to form a design method for underground target-oriented inhomogeneous laying seismic acquisition, i.e. “adding the number of the shot points in the primary energy region” described in the step 4).

[0047] A result obtained herein is data for seismic acquisition. The result is that after the design of the adjusted shot points and receivers, data acquisition can be conducted; and imaging results shown on the right side of FIG. 6 and on the right side of FIG. 7 may be obtained by seismic processing and imaging of the acquired data Imaging diagrams on the left side of FIG. 6 and on the left side of FIG. 7 show that the shot points and the receivers are not adjusted to form a homogeneously laid observation system. Therefore, results on the left side of FIG. 6 and on the left side of FIG. 7 may be obtained by seismic processing and imaging of the acquired data.

[0048] 5) Applying the above analysis method to an actual region, and comparing a conventional acquired data imaging result with the target-oriented inhomogeneous laying acquired data imaging result (as shown in FIG. 6 and FIG. 7). From FIG. 6 and FIG. 7, a target and a fracture cannot be recognized by using the traditional method; whereas by employing the method of the present invention, the target features and the fracture features may be recognized clearly.

[0049] A second aspect of the present invention provides a target-oriented seismic acquisition apparatus, comprising:

[0050] a first processing unit for giving the parameters of the initial velocity model and the three-dimensional seismic layout aiming to the underground target position;

[0051] a second processing unit for conducting wave field continuation and focusing analysis on the three-dimensional seismic layout and calculating distribution of the seismic energy on the ground in the underground target region;

[0052] a third processing unit for conducting normalization processing on distribution of the seismic energy on the ground and then conducting level partitioning to obtain the primary energy region and the secondary energy region;

[0053] a fourth processing unit for adding the number of shot points in the primary energy region to achieve target-oriented acquisition and obtaining the target-oriented inhomogeneous laying acquired data imaging features; and

[0054] a fifth processing unit for comparing the conventional acquired data imaging result with the target-oriented inhomogeneous laying acquired data imaging features.

[0055] A third aspect of the present invention provides a computer readable storage medium. A computer program is stored in the computer readable storage medium; and when the computer program is executed by a processor, the steps of the above target-oriented seismic acquisition method are implemented.

[0056] A fourth aspect of the present invention provides a computer device, comprising a memory, a processor and computer programs which are stored in the memory and may run on the processor. When the processor executes the computer programs, the steps of the above target-oriented seismic acquisition method are implemented.

[0057] The present invention makes descriptions according to the method and the device (system) of the specific implementation and a flowchart and/or a block diagram of a computer program product. It should be understood that each flow and/or block in the flowchart and/or the block diagram as well as combination of flows and/or blocks in the flowchart and/or the block diagram may be implemented by instructions of the computer programs. These instructions of the computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor or other programmable data processing devices to produce a machine; and therefore, an apparatus for achieving specified functions in one or more flows in the flowchart or one or more blocks in the block diagram is produced through the instructions executed by the processor of the general-purpose computer or other programmable data processing devices.

[0058] These instructions of the computer programs may further be stored in a computer readable memory which can guide the computer or other programmable data processing devices to work in a specific mode; and therefore, an article of manufacture comprising an instruction apparatus is produced by the instructions stored in the computer readable memory. The instruction apparatus achieves the specified functions in one or more flows in the flowchart or one or more blocks in the block diagram.

[0059] These instructions of the computer programs may further be loaded into the computer or other programmable data processing devices; and therefore, a series of operation steps are executed on the computer or other programmable devices to generate processing implemented by the computer, and then the instructions executed on the computer or other programmable devices provide the steps for achieving the specified functions in one or more flows in the flowchart or one or more blocks in the block diagram.