CN105487106A - Method for supplementing shotpoint based on energy illumination on target stratum by Gaussian beams - Google Patents

Abstract

The invention discloses a method for supplementing shotpoints based on energy illumination on a target stratum by Gaussian beams. The method comprises the following steps: according to existing geological data, establishing a three-dimensional geologic model of an exploration region, based on Gaussian beams of shotpoints, performing illumination analysis on the geologic model to determine element displacement energy of a target stratum; counting distribution condition of the element displacement energy of the target stratum, to determine an energy weak region in the target stratum; according to an illumination reference value of each shotpoint on the energy weak region, determining a surface shotpoint supplementing region corresponding to the energy weak region in the target stratum; and arranging supplemented shotpoints in original shotpoint grids of the surface shotpoint supplementing region at equal interval. The method can accurately determine the energy weak region in the target stratum, and determine the number of supplemented shotpoints quantitatively, so that the position in the target stratum where illumination energy is weak is compensated, thereby improving signal to noise ratio after seismic data is processed, and improving seismic acquisition data quality.

Description

A kind of benefit big gun method based on the illumination of Gaussian ray bundle zone of interest energy

Technical field

The present invention relates to technical field of geophysical exploration, particularly relate to a kind of benefit big gun method based on the illumination of Gaussian ray bundle zone of interest energy.

Background technology

Earthquake-capturing is the first link of seismic prospecting, is most critical one step determining whole seismic data quality quality, and recording geometry design is the emphasis in earthquake-capturing link.Recording geometry design mainly rationally arranges the distributing position of shot point and geophone station, improves efficiency and the order of accuarcy of field data collection work.

Conventional vision systems design is mainly theoretical based on HORIZONTAL LAYERED MEDIUM WITH HIGH ACCURACY, and recording geometry property parameters all obtains based on HORIZONTAL LAYERED MEDIUM WITH HIGH ACCURACY, does not consider actual architectonic impact.Because the zone of interest in actual tectonic structure is not in the distribution of strict horizontal layer, the seismic data that existing recording geometry obtains there will be the more weak region of signal energy.

In the region that this part energy is more weak, signal to noise ratio (S/N ratio) is lower, and can not do any process to the weak district of energy in conventional vision systems, the seismic data therefore observed is difficult to the requirement reaching seismic processing and explanation.Causing can not for the subsurface geological structure of energy weak district rediscover in follow-up seismic processing process.

Therefore, a kind of can carrying out mending the method for big gun in the region more weak to zone of interest energy is needed badly.

Summary of the invention

The present invention is directed to the deficiencies in the prior art, propose a kind of benefit big gun method based on the illumination of Gaussian ray bundle zone of interest energy, comprise the following steps:

Set up the three-dimensional geological model of survey area according to existing geologic information, based on shot point Gaussian ray bundle, illumination analysis is done to geologic model and determine zone of interest bin deflection energy;

Statistics zone of interest bin deflection energy distribution situation is to draw a circle to approve the weak district of energy in zone of interest;

The earth's surface Bu Bao region corresponding to the weak district of energy in the illumination reference value determination zone of interest in the weak district of energy according to each shot point;

Arrange at original shot point grid equal intervals in Bu Bao region, earth's surface and mend big gun shot point.

According to one embodiment of present invention, describedly based on shot point Gaussian ray bundle, illumination analysis is done to geologic model and determines, in the step of zone of interest bin deflection energy, to comprise further:

Based on shot point Gaussian ray bundle, illumination analysis is done to geologic model, determine the projectile energy of the whole Gaussian ray bundles inciding each zone of interest bin;

Calculate via after zone of interest reflection, the Gauss Ray beam energy of the covering wave detector place grid scope received by each earth's surface wave detector;

According to incide a zone of interest bin and by whole earth's surfaces wave detector receive gaussian line beam energy determination zone of interest bin deflection energy.

According to one embodiment of present invention, describedly based on shot point Gaussian ray bundle, illumination analysis is done to geologic model, determines that the projectile energy of the whole Gaussian ray bundles inciding each zone of interest bin comprises:

Statistics incides whole Gaussian ray bundles of zone of interest;

Statistics covers the Gaussian ray bundle of each zone of interest bin respectively;

Superposition calculation incides the projectile energy of the Gaussian ray bundle of each zone of interest bin, obtains each zone of interest bin projectile energy.

According to one embodiment of present invention, described calculating is via after zone of interest reflection, and the Gauss Ray beam energy of the covering wave detector place grid scope received by each earth's surface wave detector comprises:

Statistics reflexes to all Gaussian ray bundles on earth's surface by zone of interest;

Statistics covers the Gaussian ray bundle of each wave detector place, earth's surface grid scope respectively;

The projectile energy of whole Gaussian ray bundles that superposition calculation is received by each earth's surface wave detector, obtains the received energy of each earth's surface wave detector.

According to one embodiment of present invention, described basis incide a zone of interest bin and by whole earth's surfaces wave detector receive gaussian line beam energy determination zone of interest bin deflection energy comprise:

Calculate the ratio that each zone of interest bin projectile energy accounts for the total projectile energy of shot point Gaussian ray bundle;

Add up the received energy of each earth's surface wave detector, calculate the Gauss Ray beam energy from each shot point reflected by a zone of interest bin that all earth's surfaces wave detector receives;

The Gauss Ray beam energy from all shot points that superposition calculation is reflected by zone of interest bin, obtains the usable reflection energy of zone of interest bin;

The product accounting for the ratio of the total projectile energy of shot point Gaussian ray bundle by the usable reflection energy of zone of interest bin and the projectile energy of zone of interest bin obtains the bin deflection energy of zone of interest bin.

According to one embodiment of present invention, the described benefit big gun region corresponding to the weak district of energy in the illumination reference value determination zone of interest in the weak district of energy according to each shot point comprises:

Statistics incides shot point corresponding to the Gaussian ray bundle in the weak district of zone of interest energy;

Calculate respectively and incide the Gauss Ray beam energy from each shot point in the weak district of zone of interest energy and the ratio of the total projectile energy of the weak Qu of energy, obtain the illumination reference value of each shot point to the weak district of energy;

Select to be more than or equal to the shot point composition Bu Bao region, earth's surface of establishing threshold value in advance to the illumination reference value in the weak district of energy.

According to one embodiment of present invention, describedly arrange and mend big gun shot point and be mending original shot point grid equal intervals in big gun region:

A benefit big gun shot point is set in the centre position of original two adjacent shot points in Bu Bao region, earth's surface.

According to one embodiment of present invention, the described three-dimensional geological model setting up survey area according to existing geologic information comprises:

Existing geologic information is utilized to build mid-deep strata tectonic model and terrain model respectively;

Set the velocity of longitudinal wave on each stratum, shear wave velocity and density, obtain three-dimensional block geologic model.

According to one embodiment of present invention, the bin deflection energy of described zone of interest bin is expressed as:

$\mathrm{ECellRecv}(i,j)=\underset{p=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{ER}\left(p\right)\·\mathrm{PCell}(i,j)$

Wherein, (i, j) horizontal coordinate in zone of interest is represented, for the purpose of p, aspect unit reflexes to the beam numbering on earth's surface, for the purpose of ER (p), in aspect unit, p article of beam reflection arrives ground by all wave detector received energy superposition summations, for the purpose of N, aspect unit reflexes to the sum of the beam on ground, and PCell (i, j) represents that zone of interest bin projectile energy accounts for the ratio of the total projectile energy of shot point beam.

According to one embodiment of present invention, described preset energy threshold value is 60%.

The concept that the present invention proposes zone of interest bin deflection energy really can be used for the effective received energy of underground bin of ground surface imaging for characterizing, thus can the quality of quantitative evaluation recording geometry design, more reasonably design view examining system parameter.It is more integrated that the proposition of bin deflection energy simultaneously makes earthquake-capturing and data be processed into picture, for providing foundation with the Seismic acquisition designs being imaged as object.

The present invention accurately can also draw a circle to approve the weak district of energy of zone of interest, quantitatively determines to mend big gun number, thus the weak position of illumination energy in zone of interest is compensated, improve the signal to noise ratio (S/N ratio) after seism processing, improve earthquake-capturing data quality.

Other features and advantages of the present invention will be set forth in the following description, and partly become apparent from instructions, or understand by implementing the present invention.Object of the present invention and other advantages realize by structure specifically noted in instructions, claims and accompanying drawing and obtain.

Accompanying drawing explanation

Fig. 1 a is the Gaussian ray bundle distribution of amplitudes schematic diagram according to the embodiment of the present invention;

Fig. 1 b is the principle schematic based on Gaussian ray bundle energy supposition according to the embodiment of the present invention;

Fig. 2 is the ray coordinator schematic diagram according to the embodiment of the present invention;

Fig. 3 is the flow chart of steps of the benefit big gun method based on the illumination of Gaussian ray bundle zone of interest energy according to the embodiment of the present invention;

Fig. 4 is the flow chart of steps of the determination zone of interest bin deflection energy according to the embodiment of the present invention;

Fig. 5 a is the Gaussian ray bundle distribution schematic diagram inciding zone of interest bin;

Fig. 5 b is the distribution schematic diagram of the zone of interest reflection wire harness that earth's surface wave detector receives;

Fig. 5 c is zone of interest bin deflection energy beam distribution schematic diagram;

Fig. 6 a is the three-dimensional geological structural representation of multi-thread bundle incident illumination;

Fig. 6 b is the three-dimensional geological structural representation that multi-thread bundle receives illumination;

Fig. 6 c is the three-dimensional geological structural representation of multi-thread bundle skew illumination;

Fig. 7 is the three-dimensional geologic structure model schematic of certain engineering example;

Fig. 8 is the distribution plan that before mending big gun, recording geometry is thrown light on to zone of interest;

Fig. 9 is the distribution plan that after mending big gun, recording geometry is thrown light on to zone of interest;

Figure 10 a is the distribution plan in Bu Bao region, earth's surface;

Figure 10 b is the amplification display result in Bu Bao region, earth's surface.

Embodiment

Build the three-dimensional geological model based on survey area in embodiments of the invention, utilize the lighting engineering based on Gaussian beam beam to obtain the deflection energy distribution of underground zone of interest.The location-appropriate finding zone of interest energy weak by distribution results carries out benefit big gun in relevant position, earth's surface, quantitatively determines to mend big gun number, and the weak position of illumination energy in zone of interest is compensated.After making seism processing, zone of interest energy and signal to noise ratio (S/N ratio) are improved, and finally reach the object improving earthquake-capturing data quality.

Below the lighting engineering based on Gaussian ray bundle is described in detail.

As shown in Figure 1a, Gaussian ray bundle can be regarded as one from the energy pipe of focus centered by ray, and the distribution of amplitudes of beam exponentially decays with the square distance departing from central ray.As shown in Figure 1 b, the wave field at acceptance point R or subsurface imaging point place can be regarded as by many from focal point s0, the superposition of the Gauss Ray beam energy in R point certain limit.

In fig. 2, P be space a bit, its vertical projection point on central ray S is P ', s is the arc length of P ' point to focal point s0, and n represents the distance of P ' point to P point, and namely P point is to the distance of central ray S.In Fig. 2, ray coordinator is defined by vector n (law vector along ray) and t (tangent vector along ray).

Gaussian ray bundle formula is the high-frequency asymptote solution that the wave equation be based upon under ray center coordinate system concentrates near ray.Near central ray, the distribution of high-frequency energy represents with following formula:

$u(s,n,\mathrm{\ω})=\frac{\mathrm{\ψ}}{\sqrt{\mathrm{\ρ}\left(s\right)v\left(s\right)q\left(s\right)}}\exp \±\{\mathrm{i\ω\τ}\left(s\right)+\frac{\mathrm{i\ω}}{2}\frac{p\left(s\right)}{q\left(s\right)}{n}^2\}---\left(1\right)$

In above formula, u represents P-wave displacement, exponential part-number expression forward continuation (just drilling) ,+number expression backward extension.ψ is a certain setting constant value, and ω represents Gaussian ray bundle wave field frequency.(s, n) is the ray coordinates of calculation level under central ray coordinate system.V represents the speed of central ray, ray whilst on tour centered by τ.

(1) in formula, p (s) and q (s) is the complex value kinetic parameter of centrally ray change, meets following ordinary differential system:

$\begin{array}{cc}\frac{\mathrm{dq}}{\mathrm{ds}}=\mathrm{vp}& \frac{\mathrm{dp}}{\mathrm{ds}}=-v^{-2}\frac{{\∂}^{2}v}{\∂n^2}q\end{array}---\left(2\right)$

Formula is above converted, the Gaussian ray bundle expression-form of physical significance can be had more as follows:

$u(s,n,\mathrm{\ω})=A\left(s\right)\exp \{\mathrm{i\ω\τ}\left(s\right)+\frac{\mathrm{i\ω}}{2v\left(s\right)}K\left(s\right)n^2-\frac{n^2}{L^2\left(s\right)}\}---\left(3\right)$

Wherein:

$A\left(s\right)=\frac{\mathrm{\ψ}}{\sqrt{\mathrm{\ρ}\left(s\right)v\left(s\right)q\left(s\right)}}---\left(4\right)$

K(s)＝v(s)Re[p(s)/q(s)](5)

$L\left(s\right)={\left\{\frac{\mathrm{\ω}}{2}\mathrm{Im}\right[p\left(s\right)/q\left(s\right)\left]\right\}}^{-(1/2)}---\left(6\right)$

In formula, A is Gaussian ray bundle amplitude, and K represents the wavefront curvature of beam, and L represents effective half width of ray, and L determines the distribution of Gaussian ray bundle amplitude near central ray, is similar to Gaussian distribution so must be called Gaussian ray bundle due to distribution of amplitudes.Formula above can be utilized like this to carry out just drilling of Gaussian ray bundle.

Embodiments of the invention provide a kind of benefit big gun method based on the illumination of Gaussian ray bundle zone of interest energy.Figure 3 shows that the flow chart of steps of the benefit big gun method of the present embodiment.

First, in step S101, set up the three-dimensional geological model of survey area according to existing geologic information, based on shot point Gaussian ray bundle, illumination analysis is done to geologic model, determine zone of interest bin deflection energy.

Wherein, three-dimensional geological modeling is the basis of illumination and skew in the subsequent step of the present embodiment.The present embodiment carries out modeling mainly for beam, therefore needs to determine three-dimensional block geologic model.

The basic skills of three-dimensional modeling utilizes interpolation tectonic level, tomography curved surface, by cutting, sews up the consistent surface model of building topology, utilize automatic three-dimensional block method for tracing to set up three-dimensional bits body Model.Specifically, section or Depth Domain layer bit data/time domain layer bit data obtain mid-deep strata tectonic model to utilize the Depth Domain in work area to explain, wherein, when time domain data need, dark chi carries out time and depth transfer; Then, GIS elevation or SPS geophone station Coordinate generation real surface model is utilized.And then set the velocity of longitudinal wave on each stratum, shear wave velocity and density, obtain three-dimensional block geologic model.

Utilize Gaussian ray bundle to carry out orthodromic illumination, virtual wave detector can be placed on reflecting interface or earth's surface to calculate the energy received by virtual wave detector, determine zone of interest bin deflection energy.

The detailed process determining zone of interest bin deflection energy is described in step S101 in detail below with reference to Fig. 4.

First, in sub-step S1011, based on shot point Gaussian ray bundle, illumination analysis is done to geologic model, determine the projectile energy of the whole Gaussian ray bundles inciding each zone of interest bin.

Specifically, as shown in Figure 5 a, add up all Gaussian ray bundles inciding zone of interest, thus statistics covers whole Gaussian ray bundles of each zone of interest bin respectively.Superposition calculation incides the projectile energy of the whole Gaussian ray bundles on each zone of interest bin, obtains each zone of interest bin projectile energy ECell (x, y), forms the incident illumination result of the multi beam ray shown in Fig. 6 a.

Next, perform sub-step S1012, calculate via after zone of interest reflection, the Gauss Ray beam energy of the covering wave detector place grid scope received by each earth's surface wave detector.Namely calculate the received energy of wave detector, namely wave detector receives the energy from shot point.

As shown in Figure 5 b, add up all Gaussian ray bundles being reflexed to earth's surface by zone of interest, statistics covers the Gaussian ray bundle of each wave detector place, earth's surface grid scope respectively.The projectile energy of whole Gaussian ray bundles that superposition calculation is received by each earth's surface wave detector, obtains the received energy ERecv (x, y) of each earth's surface wave detector, forms the reception lighting result of the multi beam ray shown in Fig. 6 b.

Subsequently, in step S1013 according to incide a zone of interest bin and by whole earth's surfaces wave detector receive gaussian line beam energy determination zone of interest bin deflection energy, namely reflected by bin and the effective energy that can be used for realizing migration imaging received by earth's surface wave detector.

As shown in Figure 5 c, bin deflection energy refers to that the seismic event of epicenter excitation is by after the reflection of zone of interest bin, is observed the energy summation that multiple wave detector in system receives, represents the effective received energy that really may be used for imaging underground bin.

Specifically, the ratio PCell (i, j) that each zone of interest bin projectile energy accounts for the total projectile energy of shot point beam is first calculated.

Add up the received energy of each earth's surface wave detector, calculate the Gauss Ray beam energy from each shot point reflected by a zone of interest bin that all earth's surfaces wave detector receives, then the Gauss Ray beam energy from all shot points that reflected by a zone of interest bin of superposition calculation, obtains the usable reflection energy of a zone of interest bin.

Finally, the projectile energy of the usable reflection energy and this zone of interest bin that calculate this zone of interest bin accounts for the product of the ratio of the total projectile energy of shot point Gaussian ray bundle, obtains the bin deflection energy of this zone of interest bin:

$\mathrm{ECellRecv}(i,j)=\underset{p=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{ER}\left(p\right)\·\mathrm{PCell}(i,j)---\left(7\right)$

Wherein, (i, j) horizontal coordinate in zone of interest is represented, for the purpose of p, aspect unit reflexes to the beam numbering on ground, for the purpose of ER (p), in aspect unit, p article of beam reflection arrives ground by all wave detector received energy superposition summations, for the purpose of N, aspect unit reflexes to the sum of the beam on ground represent the usable reflection energy of zone of interest bin.

The numerical value of ER (p) can calculate to (6) according to expression formula (1), then can obtain the skew lighting result of the multi beam ray shown in Fig. 6 c according to expression formula (7).Like this, can analyze three-dimensional geological model based on the lighting engineering of Gaussian beam beam, obtain the deflection energy distribution of zone of interest.

Compared with prior art, the concept of zone of interest bin deflection energy is proposed for characterizing the effective received energy of underground bin that really can be used for ground surface imaging in the present embodiment, thus can quantitative evaluation recording geometry design quality, more reasonably design view examining system parameter.It is more integrated that the proposition of bin deflection energy simultaneously makes earthquake-capturing and data be processed into picture, for providing foundation with the Seismic acquisition designs being imaged as object.

Subsequently, perform step S102, statistics zone of interest bin deflection energy distribution situation is to draw a circle to approve the weak district of energy in zone of interest.

Concrete, the zone of interest of first selected seismic prospecting, namely the geological problem solving which zone of interest of underground is mainly wanted in this seismic prospecting, then bin deflection energy statistics is carried out for selected underground zone of interest, according to the actual configuration of geologic province or the threshold value in the exploration object setting weak district of energy, the weak district of energy in delineation zone of interest.

Next, corresponding to the weak district of energy in the illumination reference value determination zone of interest in the weak district of energy according to each shot point in step s 103 benefit big gun region.

First the energy threshold that shot point is irradiated to the weak district of energy is set.Such as can according to 60% of territory, shot point irradiation energy occupied area gross energy as energy threshold.

Then the shot point that the Gaussian ray bundle in the statistics arrival weak district of zone of interest energy is corresponding, calculate respectively and arrive the Gauss Ray beam energy from each shot point in the weak district of zone of interest energy and the ratio of the total projectile energy of the weak Qu of energy, obtain the illumination reference value of each shot point to the weak district of energy.Select the shot point composition Bu Bao region, the earth's surface illumination reference value in the weak district of energy being more than or equal to predetermined threshold value.

Finally, perform step S104, arrange at the original shot point grid equal intervals mending big gun region and mend big gun shot point.

Preferably, a benefit big gun shot point is set in the centre position of original two adjacent shot points in Bu Bao region, earth's surface, encryption shot point is set like this and can makes original shot point mesh refinement one times.Be not limited thereto, also two benefit big gun shot points can be set in the centre position of original two adjacent shot points, or more, specifically can set according to implementation cost.So, by benefit big gun region encryption shot point, the weak position of illumination energy in zone of interest can be compensated, improve the signal to noise ratio (S/N ratio) of seismic data.

case history

Below for DEH area, describe whole benefit big gun process in detail.Recording geometry parameter is as shown in table 1 below.

Table 1

A three-dimensional structure model (as shown in Figure 7) is set up by three depth sections crossing work area.Fig. 8 is the distribution that before mending big gun, recording geometry is thrown light on to zone of interest, and as can be seen from Figure 8, bin deflection energy difference in distribution is very large, therefore chooses position shown in black circle in figure to carry out the demonstration of benefit big gun.Fig. 9 is the distribution of the zone of interest energy after mending big gun, and known black circled positions energy is compensated (compensating big gun number is 100 big guns, with zone of interest ceiling capacity for standard).Figure 10 a and Figure 10 b is the amplification display that big gun number was put and mended in benefit emplacement, ground, and as we can see from the figure, the position of benefit big gun is all the position after original shot point mesh refinement, and original shot point grid is encrypted at double.Comparison diagram 8 and Fig. 9 can find out, the present invention accurately can draw a circle to approve the weak district of energy of zone of interest, quantitatively determines to mend big gun number, thus the weak position of illumination energy in zone of interest is compensated, improve the signal to noise ratio (S/N ratio) after seism processing, improve earthquake-capturing data quality.

Although embodiment disclosed in this invention is as above, the embodiment that described content just adopts for the ease of understanding the present invention, and be not used to limit the present invention.Technician in any the technical field of the invention; under the prerequisite not departing from spirit and scope disclosed in this invention; any amendment and change can be done what implement in form and in details; but scope of patent protection of the present invention, the scope that still must define with appending claims is as the criterion.

Claims (10)

1., based on a benefit big gun method for Gaussian ray bundle zone of interest energy illumination, it is characterized in that, comprise the following steps:

Set up the three-dimensional geological model of survey area according to existing geologic information, based on shot point Gaussian ray bundle, illumination analysis is done to geologic model and determine zone of interest bin deflection energy;

Statistics zone of interest bin deflection energy distribution situation is to draw a circle to approve the weak district of energy in zone of interest;

The earth's surface Bu Bao region corresponding to the weak district of energy in the illumination reference value determination zone of interest in the weak district of energy according to each shot point;

Arrange at original shot point grid equal intervals in Bu Bao region, earth's surface and mend big gun shot point.

2. method according to claim 1, is characterized in that, describedly does illumination analysis based on shot point Gaussian ray bundle to geologic model and determines that the step of zone of interest bin deflection energy comprises:

Based on shot point Gaussian ray bundle, illumination analysis is done to geologic model, determine the projectile energy of the whole Gaussian ray bundles inciding each zone of interest bin;

Calculate via after zone of interest reflection, the Gauss Ray beam energy of the covering wave detector place grid scope received by each earth's surface wave detector;

According to incide a zone of interest bin and by whole earth's surfaces wave detector receive gaussian line beam energy determination zone of interest bin deflection energy.

3. method according to claim 2, is characterized in that, describedly does illumination analysis based on shot point Gaussian ray bundle to geologic model, determines that the projectile energy of the whole Gaussian ray bundles inciding each zone of interest bin comprises:

Statistics incides whole Gaussian ray bundles of zone of interest;

Statistics covers the Gaussian ray bundle of each zone of interest bin respectively;

Superposition calculation incides the projectile energy of the Gaussian ray bundle of each zone of interest bin, obtains each zone of interest bin projectile energy.

4. method according to claim 3, is characterized in that, described calculating is via after zone of interest reflection, and the Gauss Ray beam energy of the covering wave detector place grid scope received by each earth's surface wave detector comprises:

Statistics reflexes to all Gaussian ray bundles on earth's surface by zone of interest;

Statistics covers the Gaussian ray bundle of each wave detector place, earth's surface grid scope respectively;

The projectile energy of whole Gaussian ray bundles that superposition calculation is received by each earth's surface wave detector, obtains the received energy of each earth's surface wave detector.

5. method according to claim 4, is characterized in that, described basis incide a zone of interest bin and by whole earth's surfaces wave detector receive gaussian line beam energy determination zone of interest bin deflection energy comprise:

Calculate the ratio that each zone of interest bin projectile energy accounts for the total projectile energy of shot point Gaussian ray bundle;

Add up the received energy of each earth's surface wave detector, calculate the Gauss Ray beam energy from each shot point reflected by a zone of interest bin that all earth's surfaces wave detector receives;

The Gauss Ray beam energy from all shot points that superposition calculation is reflected by zone of interest bin, obtains the usable reflection energy of zone of interest bin;

The product accounting for the ratio of the total projectile energy of shot point Gaussian ray bundle by the usable reflection energy of zone of interest bin and the projectile energy of zone of interest bin obtains the bin deflection energy of zone of interest bin.

6. the method according to any one of claim 1-5, is characterized in that, the described benefit big gun region corresponding to the weak district of energy in the illumination reference value determination zone of interest in the weak district of energy according to each shot point comprises:

Statistics incides shot point corresponding to the Gaussian ray bundle in the weak district of zone of interest energy;

Calculate respectively and incide the Gauss Ray beam energy from each shot point in the weak district of zone of interest energy and the ratio of the total projectile energy of the weak Qu of energy, obtain the illumination reference value of each shot point to the weak district of energy;

Select the shot point composition Bu Bao region, the earth's surface illumination reference value in the weak district of energy being more than or equal to preset energy threshold value.

7. the method according to any one of claim 1-5, is characterized in that, described original shot point grid equal intervals mending big gun region arranges benefit big gun shot point and is:

A benefit big gun shot point is set in the centre position of original two adjacent shot points in Bu Bao region, earth's surface.

8. the method according to any one of claim 1-5, is characterized in that, the described three-dimensional geological model setting up survey area according to existing geologic information comprises:

Existing geologic information is utilized to build mid-deep strata tectonic model and terrain model respectively;

Set the velocity of longitudinal wave on each stratum, shear wave velocity and density, obtain three-dimensional block geologic model.

9. method according to claim 5, is characterized in that, the bin deflection energy of described zone of interest bin is expressed as:

$\mathrm{ECellRecv}(i,j)=\underset{p=1}{\overset{N}{\mathrm{\Σ}}}\mathrm{ER}\left(p\right)\·\mathrm{PCell}(i,j)$

Wherein, (i, j) horizontal coordinate in zone of interest is represented, for the purpose of p, aspect unit reflexes to the beam numbering on earth's surface, for the purpose of ER (p), in aspect unit, p article of beam reflection arrives ground by all wave detector received energy superposition summations, for the purpose of N, aspect unit reflexes to the sum of the beam on ground, and PCell (i, j) represents that zone of interest bin projectile energy accounts for the ratio of the total projectile energy of shot point beam.

10. method according to claim 6, is characterized in that, described preset energy threshold value is 60%.