CN104155703A - Method and device for evaluating three-dimensional observing system - Google Patents

Abstract

The invention provides a method and device for evaluating a three-dimensional observing system. The method includes the steps of (a) establishing a three-dimensional geologic model of a work area; (b) carrying out simulation blasting on the three-dimensional geologic model according to the current to-be-evaluated three-dimensional observing system, and collecting seismic data of CRP surface elements in a full coverage area of the three-dimensional geologic model; (c) determining energy attributes, geophone pair azimuth angle attributes and geophone distance attributes of the current to-be-evaluated three-dimensional observing system according to the seismic data; (d) determining the evaluation coefficient of the current to-be-evaluated three-dimensional observing system according to the energy attributes, the geophone pair azimuth angle attributes and the geophone distance attributes of the current to-be-evaluated three-dimensional observing system. By means of the method and device, evaluation is carried out based on the seismic data of the CRP surface elements, and factors referring to the performance of the three-dimensional observing system are considered more comprehensively, so that the evaluation result is more accurate.

Description

Evaluate the method and apparatus of stereo observing system

Technical field

All things considered of the present invention relates to geophysical exploration technology field, more particularly, relates to a kind of method and apparatus of evaluating stereo observing system.

Background technology

When carrying out field earthquake data acquisition for Given task, the stereo observing system of better performances can meet the requirement to image formation quality, suppresses preferably the various random disturbance in field, and reduces the financial cost of data acquisition.That is to say, the quality of stereo observing system has determined the quality of the geological data that gathers, therefore, how to choose science, effectively stereo observing system is one of the gordian technique that improves the quality of the geological data gathering.

About how evaluating stereo observing system, all have conducted intensive studies both at home and abroad, the factor that current evaluation of programme stresses is different, and the metewand of employing is also different.

For example, the criterion that Berkhout proposed three arguments in 1998 is for evaluating stereo observing system, respectively: (1) signal argument (resolution after amplitude accuracy and data skew); (2) noise argument (noise compacting rate); (3) economic argument (balance between shot point and geophone station stacking fold).

In addition, by underground geologic bodies is thrown light on evaluate stereo observing system aspect, 1997, Berkhout introduced the burnt CFP imaging concept of copolymerization and focused beam acts concept, realizes the migration before stack of seismic wave field with double focusing; Volker has also done similar research at 1998 to 2002 and Veldhuizen in 2003 to 2006, propose the method by image quality analyzing three-dimensional recording geometry; The Li Wanwan of China takes the method for Wave equation forward modeling, the zone of interest of complex area is thrown light on, and then establish acquisition parameter.But said method calculated amount is large, practicality is not strong, and result of calculation relies on the accuracy of model.

Based on common midpoint (CMP) bin attributive analysis, evaluate stereo observing system aspect, the people such as Zhao Diandong, Lv Gonghe uses high-precision three-dimensional seismic data acquisition method on Shengli Oil Field block, from demonstration of acquisition parameters, attributive analysis, to forward modeling, has obtained good result; Yin Wuhai, the people such as Zhen Wensheng, by multiple technologies such as Wave equation forward modeling, illumination simulation and stereo observing system CMP bin attributive analysises, have solved the contrary fracture lower wall zone of interest imaging problem of covering effectively.Because said method is evaluated stereo observing system by hypothesis HORIZONTAL LAYERED MEDIUM WITH HIGH ACCURACY, so work as surface relief, when underground structure is complicated, common midpoint (CMP) bin and common reflection point (CRP) bin do not overlap, and the error of evaluation result is larger.

The scheme of above stereo observing system being evaluated is mainly to consider from factors such as geologic body illumination, amplitude accuracy, noise compacting rate and CMP bin attributes, and accuracy rate is on the low side.

Summary of the invention

Exemplary embodiment of the present invention is to provide a kind of method and apparatus of evaluating stereo observing system, and described method and apparatus has higher accuracy.

According to an aspect of the present invention, provide a kind of method of evaluating stereo observing system, described method comprises: the three-dimensional geological model that (a) creates work area; (b) on described three-dimensional geological model, according to current stereo observing system to be evaluated, simulate and blow out, gather the geological data of each common reflection point (CRP) bin in the full overlay area of described three-dimensional geological model, wherein, the geological data of each CRP bin comprises: the energy intensity receiving at CRP bin and the big gun that occurs to reflect in CRP bin are examined right position angle and geophone offset; (c) according to described geological data, determine that the energy properties of current stereo observing system to be evaluated, big gun examine position angle attribute and geophone offset attribute; (d) according to the energy properties of current stereo observing system to be evaluated, big gun inspection, position angle attribute and geophone offset attribute are determined to the evaluation coefficient of current stereo observing system to be evaluated, wherein, described evaluation coefficient is used to indicate the performance of stereo observing system.

In described method, described energy properties can represent by energy intensity and energy even sex index, and described big gun inspection can represent by position angle index similarity position angle attribute, and described geophone offset attribute can be represented by geophone offset evenness index.

In the step (d) of described method, can determine by following objective function the evaluation coefficient of current stereo observing system to be evaluated:

P＝E/E _max+F _min/F+C _min/C+V _sta/V，

Wherein, the evaluation coefficient of the current stereo observing system to be evaluated of P indication, the energy intensity of the current stereo observing system to be evaluated of E indication, the energy even sex index of the current stereo observing system to be evaluated of F indication, the position angle index similarity of the current stereo observing system to be evaluated of C indication, the geophone offset evenness index of the current stereo observing system to be evaluated of V indication, E _maxfor the maximal value in the energy intensity value of all stereo observing systems to be evaluated, F _minfor the minimum value in the energy even sex index of all stereo observing systems to be evaluated, C _minfor the minimum value in the position angle index similarity of all stereo observing systems to be evaluated, V _staideal value for geophone offset evenness index.

In the step (c) of described method, can determine described energy intensity E by following calculating formula:

$E=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l,$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

In the step (c) of described method, can determine described energy even sex index F by following calculating formula:

$F=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{(Q_l-\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l)}^2,$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

In the step (c) of described method, can determine described position angle index similarity C by following calculating formula:

$C=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\underset{\mathrm{\α}=1}{\overset{360}{\mathrm{\Σ}}}\frac{|H_l\left(\mathrm{\α}\right)-M\left(\mathrm{\α}\right)|}{H_l\left(\mathrm{\α}\right)+M\left(\mathrm{\α}\right)}\right),$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, α indicates angle, H _l(α) big gun that indication occurs to reflect and position angle is α in l CRP bin is examined right number,

In the step (c) of described method, can determine described geophone offset evenness index V by following calculating formula:

$V=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\frac{1}{N_l-1}\underset{i=1}{\overset{N_l-1}{\mathrm{\Σ}}}(1+|1-\frac{S_l(i+1)-S_l\left(i\right)}{S_{l0}}\left|\right)\right),$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, N _lfor there is the big gun of reflection in l CRP bin, examine right number, i indicates all big guns inspections that reflection occurs in same CRP bin sequence number when sorting from small to large by the size of geophone offset, S _l(i) i the big gun that indication occurs to reflect in l CRP bin examined right geophone offset, S _l(i+1) i+1 the big gun that indication occurs to reflect in l CRP bin examined right geophone offset,

According to a further aspect in the invention, provide a kind of equipment of evaluating stereo observing system, described equipment comprises: three-dimensional geological model creating unit, the three-dimensional geological model in establishment work area; Seismic data acquisition unit, on described three-dimensional geological model, according to current stereo observing system to be evaluated, simulate and blow out, gather the geological data of each common reflection point (CRP) bin in the full overlay area of described three-dimensional geological model, wherein, the geological data of each CRP bin comprises: the energy intensity receiving at CRP bin and the big gun that occurs to reflect in CRP bin are examined right position angle and geophone offset; Attribute determining unit, determines that according to described geological data the energy properties of current stereo observing system to be evaluated, big gun examine position angle attribute and geophone offset attribute; Evaluation coefficient determining unit, according to the energy properties of current stereo observing system to be evaluated, big gun inspection, position angle attribute and geophone offset attribute are determined to the evaluation coefficient of current stereo observing system to be evaluated, wherein, described evaluation coefficient is used to indicate the performance of stereo observing system.

In described equipment, described energy properties can represent by energy intensity and energy even sex index, and described big gun inspection can represent by position angle index similarity position angle attribute, and described geophone offset attribute can be represented by geophone offset evenness index.

In described equipment, evaluation coefficient determining unit can be determined by following objective function the evaluation coefficient of current stereo observing system to be evaluated:

P＝E/E _max+F _min/F+C _min/C+V _sta/V，

Wherein, the evaluation coefficient of the current stereo observing system to be evaluated of P indication, the energy intensity of the current stereo observing system to be evaluated of E indication, the energy even sex index of the current stereo observing system to be evaluated of F indication, the position angle index similarity of the current stereo observing system to be evaluated of C indication, the geophone offset evenness index of the current stereo observing system to be evaluated of V indication, E _maxfor the maximal value in the energy intensity value of all stereo observing systems to be evaluated, F _minfor the minimum value in the energy even sex index of all stereo observing systems to be evaluated, C _minfor the minimum value in the position angle index similarity of all stereo observing systems to be evaluated, V _staideal value for geophone offset evenness index.

In described equipment, attribute determining unit can be determined described energy intensity E by following calculating formula:

$E=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l,$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

In described equipment, attribute determining unit can be determined described energy even sex index F by following calculating formula:

$F=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{(Q_l-\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l)}^2,$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

In described equipment, attribute determining unit can be determined described position angle index similarity C by following calculating formula:

$C=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\underset{\mathrm{\α}=1}{\overset{360}{\mathrm{\Σ}}}\frac{|H_l\left(\mathrm{\α}\right)-M\left(\mathrm{\α}\right)|}{H_l\left(\mathrm{\α}\right)+M\left(\mathrm{\α}\right)}\right),$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, α indicates angle, H _l(α) big gun that indication occurs to reflect and position angle is α in l CRP bin is examined right number,

In described equipment, attribute determining unit can be determined described geophone offset evenness index V by following calculating formula:

$V=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\frac{1}{N_l-1}\underset{i=1}{\overset{N_l-1}{\mathrm{\Σ}}}(1+|1-\frac{S_l(i+1)-S_l\left(i\right)}{S_{l0}}\left|\right)\right),$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, N _lfor there is the big gun of reflection in l CRP bin, examine right number, i indicates all big guns inspections that reflection occurs in same CRP bin sequence number when sorting from small to large by the size of geophone offset, S _l(i) i the big gun that indication occurs to reflect in l CRP bin examined right geophone offset, S _l(i+1) i+1 the big gun that indication occurs to reflect in l CRP bin examined right geophone offset,

In evaluating according to an exemplary embodiment of the present invention the method and apparatus of stereo observing system, geological data based on common reflection point (CRP) bin is evaluated, thereby can reflect more exactly the distribution of the effective reflection of buried target layer, be applicable to thus have the work area of relief surface and complex structure, effectively avoid the error of the geological data of CMP bin, and evaluation result is more accurate.In addition, according in the method and apparatus of the evaluation stereo observing system of exemplary embodiment of the present invention, from energy properties, the big gun inspection of stereo observing system, position angle attribute and geophone offset attribute three aspects: are carried out to comprehensive evaluation to stereo observing system, thereby considered more all sidedly to relate to the factor of stereo observing system performance, made evaluation result more accurate.

Accompanying drawing explanation

The description of being undertaken by the accompanying drawing below in conjunction with embodiment is exemplarily shown, above and other object and the feature of exemplary embodiment of the present will become apparent, wherein:

Fig. 1 illustrates the process flow diagram of the method for evaluating according to an exemplary embodiment of the present invention stereo observing system;

Fig. 2 illustrates the example of the three-dimensional geological model in work area according to an exemplary embodiment of the present invention;

Fig. 3 A to Fig. 3 D illustrates the distribution plan of the energy intensity receiving in full overlay area according to an exemplary embodiment of the present invention;

Fig. 4 A to Fig. 4 D illustrate according to an exemplary embodiment of the present invention the big gun that reflection occurs in each CRP bin examine right position angle towards distribution plan;

Fig. 5 A to Fig. 5 D illustrates the distribution plan that the big gun that occurs according to an exemplary embodiment of the present invention to reflect in each CRP bin is examined right geophone offset size;

Fig. 6 illustrates the block diagram of the equipment of evaluating according to an exemplary embodiment of the present invention stereo observing system.

Embodiment

Now will be in detail with reference to embodiments of the invention, the example of described embodiment is shown in the drawings, and wherein, identical label refers to identical parts all the time.Below will be by described embodiment is described with reference to accompanying drawing, to explain the present invention.

Fig. 1 illustrates the process flow diagram of the method for evaluating according to an exemplary embodiment of the present invention stereo observing system.Here, described method can realize by computer program, also can realize by evaluating the equipment of stereo observing system.

With reference to Fig. 1, at step S10, create the three-dimensional geological model in work area.Particularly, can utilize the existing real data in work area to create three-dimensional geological model.Fig. 2 illustrates the example of the three-dimensional geological model in work area according to an exemplary embodiment of the present invention.In the three-dimensional geological model shown in Fig. 2, gray scale has in various degree represented different stratum, as example, as can be seen from Figure 2, and the complicated geological in the work area that described three-dimensional geological model is corresponding, surface relief is uneven, and there is fracture on part stratum.

At step S20, on described three-dimensional geological model, according to current stereo observing system to be evaluated, simulate and blow out, gather the geological data of each common reflection point (CRP) bin in the full overlay area of described three-dimensional geological model.Here, the geological data of each CRP bin of collection can comprise: the energy intensity receiving at CRP bin and the big gun that occurs to reflect in CRP bin are examined right position angle and geophone offset.

As example, the stereo observing system is here candidate's stereo observing system of choosing for exploration task and work area geological condition.By the stereo observing system of choosing is evaluated, can contribute to determine the stereo observing system in applicable work area.Table 1 illustrates the parameter list of four candidate's stereo observing systems according to an exemplary embodiment of the present invention.As example, four stereo observing systems in table 1 are candidate's stereo observing systems of choosing for the work area shown in Fig. 2.

The parameter list of table 1 stereo observing system

Stereo observing system	8L10S460R	16L5S460R	16L10S460R	20L10S600R
					Mode of excitation	Centre excites	Centre excites	Centre excites	Centre excites
Big gun line mode	Perpendicular array	Perpendicular array	Perpendicular array	Perpendicular array
					Receive line number/bar	8	16	16	20

Shot point number/	10	5	10	10
					Every receives Xian road number	460	460	460	600
Total reception channel number	3680	7360	7360	12000
					Track pitch/m	30	30	30	30
Receive line-spacing/m	600	300	300	300
					Shotpoint spacing/m	60	60	60	60
Big gun array pitch/m	300	300	300	300
					Maximum non-advance/m	2370	2370	2550	3150
Longitudinal maximum offset/m	6885	6885	6885	8985
					Maximum offset/m	7281	7281	7342	9521
Longitudinal degree of covering/time	23	23	23	30
					Lateral fold number/time	4	8	8	10
Total degree of covering/time	92	184	184	300
					Bin size/m.m	15×30	15×30	15×30	15×30
Wire harness rolling number/bar	1	1	2	2
					Wire harness rolling wire distance/m	600	300	600	600

Particularly, the big gun that the big gun that can gather energy intensity that each CRP bin receives at step S20, reflection occurs in each CRP bin is examined right position angle and reflection occurred in each CRP bin is examined right geophone offset.

Fig. 3 A to Fig. 3 D illustrates the distribution plan of the energy intensity receiving in full overlay area according to an exemplary embodiment of the present invention.Particularly, Fig. 3 A to Fig. 3 D shows on the three-dimensional geological model shown in Fig. 2 respectively according to stereo observing system 8L10S460R, 16L5S460R, 16L10S460R, 20L10S600R in table 1 and simulates while blowing out, the distribution plan of the energy intensity that each CRP bin in full overlay area collecting receives.With reference to Fig. 3 A to Fig. 3 D, the plane that x axle and y axle form represents the 2 d plane picture of full overlay area, and wherein, gray scale more deeply feels that to be shown in the energy intensity that relevant range receives larger.Full overlay area is comprised of a plurality of CRP bins, because the area of CRP bin is too little with respect to the area of full overlay area, therefore in Fig. 3 A to Fig. 3 D, specifically do not mark each CRP bin, yet, should be understood that the energy intensity in can the full overlay area based on shown in Fig. 3 A to Fig. 3 D distributes to obtain the energy intensity receiving at each CRP bin.

Fig. 4 A to Fig. 4 D illustrate according to an exemplary embodiment of the present invention the big gun that reflection occurs in each CRP bin examine right position angle towards distribution plan.Particularly, Fig. 4 A to Fig. 4 D shows on the three-dimensional geological model shown in Fig. 2 respectively according to stereo observing system 8L10S460R, 16L5S460R, 16L10S460R, 20L10S600R in table 1 and simulates while blowing out, the big gun that reflection occurs in each CRP bin in full overlay area collecting examine right position angle towards distribution plan.With reference to Fig. 4 A to Fig. 4 D, a little rectangle represents a CRP bin, because indication range is limited, only shows part CRP bin, wherein, every line in similar " rice " font pattern indicates each big gun inspection centering to blow out a little and the projection in the plane of CRP bin of the line of geophone station.

Fig. 5 A to Fig. 5 D illustrates the distribution plan that the big gun that occurs according to an exemplary embodiment of the present invention to reflect in each CRP bin is examined right geophone offset size.Particularly, Fig. 5 A to Fig. 5 D shows on the three-dimensional geological model shown in Fig. 2 respectively according to stereo observing system 8L10S460R, 16L5S460R, 16L10S460R, 20L10S600R in table 1 and simulates while blowing out, and the big gun that reflection occurs in each CRP bin in full overlay area collecting is examined the distribution plan of right geophone offset size.With reference to Fig. 5 A to Fig. 5 D, a little rectangle represents a CRP bin, because indication range is limited, only shows part CRP bin, and wherein, the vertical line of arranging from small to large in histogram indicates each big gun to examine right geophone offset.

Referring again to Fig. 1, at step S30, according to the geological data gathering at step S20, determine that the energy properties of current stereo observing system to be evaluated, big gun examine position angle attribute and geophone offset attribute.

Particularly, can current stereo observing system to be evaluated be analyzed according to described geological data, to obtain the parameter that is used for respectively representing above-mentioned three attributes.For example, for energy properties, can analyze from aspects such as energy intensity, energy intensity distributed area, energy intensity distributing homogeneities; For big gun inspection, to position angle attribute, can analyze from aspects such as azimuthal distribution interval, azimuthal distribution homogeneity, azimuthal distribution similarities; For geophone offset attribute, can analyze from aspects such as geophone offset distributed area, geophone offset distributing homogeneities.As example, above-mentioned analysis can be qualitative analysis, quantitative test or both combinations.

Preferably, for energy properties, can analyze from energy intensity, energy intensity distributing homogeneity two aspects, to obtain for representing energy intensity and the energy even sex index of energy properties; For big gun inspection, to position angle attribute, can analyze from azimuthal distribution similarity aspect, to obtain for representing the position angle index similarity of big gun inspection to position angle attribute; For geophone offset attribute, can analyze from geophone offset distributing homogeneity aspect, to obtain for representing the geophone offset evenness index of geophone offset attribute.Below how detailed description is determined to energy intensity, energy even sex index, position angle index similarity and geophone offset evenness index based on geological data.

As example, at step S30, can determine above-mentioned energy intensity by calculating formula (1):

$E=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l---\left(1\right)$

Wherein, the energy intensity of the current stereo observing system to be evaluated of E indication, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

Here, energy intensity is for reflecting the height of the energy receiving in the full overlay area of described three-dimensional geological model, the energy intensity E definite by calculating formula (1) is larger, represent that the energy receiving in the full overlay area of described three-dimensional geological model is higher, in the situation that other attributes are identical, the performance of stereo observing system is better.

In above-mentioned calculating formula (1), the mean value of the energy intensity that each CRP bin is received is as the energy intensity of stereo observing system, but those skilled in the art can understand, also can determine by other computing method the energy intensity of stereo observing system, for example, the energy intensity that all CRP bins in the full overlay area of described three-dimensional geological model can be received is added summation, and the energy intensity using the summing value obtaining as stereo observing system.

As example, at step S30, can determine energy even sex index by calculating formula (2):

$F=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{(Q_l-\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l)}^2---\left(2\right)$

Wherein, the energy even sex index of the current stereo observing system to be evaluated of F indication, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

Here, energy even sex index is for reflecting the situation that is uniformly distributed of the energy intensity that in the full overlay area of described three-dimensional geological model, each bin receives, the energy even sex index F definite by calculating formula (2) is less, represent that the energy intensity that in the full overlay area of described three-dimensional geological model, each bin receives distributes more evenly, in the situation that other attributes are identical, the performance of stereo observing system is better.

In above-mentioned calculating formula (2), the variance of the energy intensity that each CRP bin is received is as energy even sex index, but those skilled in the art can understand, and also can determine energy even sex index by other computing method.

As example, at step S30, can determine position angle index similarity by calculating formula (3):

$C=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\underset{\mathrm{\α}=1}{\overset{360}{\mathrm{\Σ}}}\frac{|H_l\left(\mathrm{\α}\right)-M\left(\mathrm{\α}\right)|}{H_l\left(\mathrm{\α}\right)+M\left(\mathrm{\α}\right)}\right)---\left(3\right)$

Wherein, the position angle index similarity of the current stereo observing system to be evaluated of C indication, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, α indicates angle, H _l(α) big gun that indication occurs to reflect and position angle is α in l CRP bin is examined right number, $M\left(\mathrm{\α}\right)=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{H}_l\left(a\right).$

In above-mentioned calculating formula (3), M (α) examines the mean value of right number for there is big gun that reflection and position angle are α in CRP bin.By the big gun that reflection and position, position angle α occur in each CRP bin being examined to right number and M (α) compares, the big gun that obtains in each CRP bin occurring reflection and position, position angle α is examined the similarity of right number and M (α), calculate again the mean value of described similarity, obtain position angle index similarity.

Here, the big gun that reflection occurs for being reflected in each CRP bin in the full overlay area of three-dimensional geological model position angle index similarity is examined the similarity degree of right azimuthal distribution.The position angle index similarity C definite by calculating formula (3) is lower, be illustrated in reflection occurs each CRP bin big gun examine right azimuthal distribution must be more similar, in the situation that other attributes are identical, the performance of stereo observing system is better, correspondingly, utilize such stereo observing system can obtain the comparatively desirable geological data of signal to noise ratio (S/N ratio) and resolution.

Under actual geologic condition, particularly in western part of China Complex Mountain structure situation, due to big gun, to examine right azimuthal distribution extremely inhomogeneous, the big gun inspection of different stereo observing systems is differing very little to position angle aspect distributing homogeneity, be therefore difficult to azimuthal distributing homogeneity be evaluated to the stereo observing system in the work area of complex geological condition by definite big gun inspection.According to exemplary embodiment of the present invention, by using position angle index similarity as evaluating one of factor of stereo observing system, more can reflect different stereo observing systems in big gun inspection to the difference on the attribute of position angle.

Here, adopted the computing method of above-mentioned calculating formula (3) to carry out computer azimuth angle index similarity, but those skilled in the art can understand, the additive method of similarity that also can be by computational data is determined position angle index similarity.

As example, in step S30, can determine geophone offset evenness index by calculating formula (4):

$V=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\frac{1}{N_l-1}\underset{i=1}{\overset{N_l-1}{\mathrm{\Σ}}}(1+|1-\frac{S_l(i+1)-S_l\left(i\right)}{S_{l0}}\left|\right)\right)---\left(4\right)$

Wherein, the geophone offset evenness index of the current stereo observing system to be evaluated of V indication, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, N _lfor there is the big gun of reflection in l CRP bin, examine right number, i indicates all big guns inspections that reflection occurs in same CRP bin sequence number when sorting from small to large by the size of geophone offset, S _l(i) i the big gun that indication occurs to reflect in l CRP bin examined right geophone offset, S _l(i+1) i+1 the big gun that indication occurs to reflect in l CRP bin examined right geophone offset,

In above-mentioned calculating formula (4), at desirable geophone offset, divide and plant, S _l(i+1)-S _l(i) should level off to S _l0, correspondingly, geophone offset evenness index V should level off to its ideal value, is 1.

Here, the big gun that reflection occurs for being reflected in each CRP bin in the full overlay area of three-dimensional geological model geophone offset evenness index is examined the situation that is uniformly distributed of right geophone offset size.The geophone offset evenness index definite by calculating formula (4) more approaches ideal value, the big gun that reflection occurs each CRP bin in the full overlay area of three-dimensional geological model is examined the size distribution of right geophone offset must be more even, in the situation that other attributes are identical, the performance of stereo observing system is better.

Here, adopted the computing method of above-mentioned calculating formula (4) to calculate geophone offset evenness index, but those skilled in the art can understand, also can determine geophone offset evenness index by calculating the inhomogeneity additive method of difference between adjacent geophone offset.

At step S40, according to the energy properties of current stereo observing system to be evaluated, big gun inspection, position angle attribute and geophone offset attribute are determined to the evaluation coefficient of current stereo observing system to be evaluated, wherein, described evaluation coefficient is used to indicate the performance of stereo observing system.

Those skilled in the art can understand, can consider according to the needs such as the geology characteristic in exploration task and work area because usually designing the funtcional relationship between evaluation coefficient and above-mentioned three attributes, thereby definite evaluation coefficient.

As example, in energy properties, by energy intensity and energy even sex index, represented, big gun inspection is represented by position angle index similarity position angle attribute, in the situation that geophone offset attribute is represented by geophone offset evenness index, can determine by following objective function (5) evaluation coefficient of current stereo observing system to be evaluated:

P＝E/E _max+F _min/F+C _min/C+V _sta/V (5)

Wherein, the evaluation coefficient of the current stereo observing system to be evaluated of P indication, the energy intensity of the current stereo observing system to be evaluated of E indication, the energy even sex index of the current stereo observing system to be evaluated of F indication, the position angle index similarity of the current stereo observing system to be evaluated of C indication, the geophone offset evenness index of the current stereo observing system to be evaluated of V indication, E _maxfor the maximal value in the energy intensity value of all stereo observing systems to be evaluated, F _minfor the minimum value in the energy even sex index of all stereo observing systems to be evaluated, C _minfor the minimum value in the position angle index similarity of all stereo observing systems to be evaluated, V _staideal value for geophone offset evenness index.In geophone offset evenness index be definite by calculating formula (4) in the situation that, the ideal value V of geophone offset evenness index _stabe 1.

When calculating formula (1) to (4) is obtained corresponding energy intensity, energy even sex index, position angle index similarity, geophone offset evenness index as described above, the maximal value by the definite evaluation coefficient P of above-mentioned objective function (5) is 4.In order to make evaluation coefficient can indicate more intuitively the characteristic of stereo observing system, the result that objective function (5) can be calculated is multiplied by 25 as evaluation coefficient, and as following objective function (6), like this, the maximal value of the evaluation coefficient calculating is 100.

P＝(E/E _max+F _min/F+C _min/C+V _sta/V)×25 (6)

Table 2 is listed the evaluating list of four stereo observing systems shown in the table 1 drawing according to above computing method, and evaluating comprises: energy intensity, energy even sex index, position angle index similarity, geophone offset evenness index and evaluation coefficient.Wherein, the energy intensity of stereo observing system 20L10S600R is maximum, the energy even sex index of stereo observing system 8L10S460R is minimum, the index similarity of stereo observing system 8L10S460R is minimum, the geophone offset evenness index of stereo observing system 20L10S600R approaches ideal value most, the evaluation coefficient of stereo observing system 8L10S460R is the highest, and performance is best.

The evaluating list of table 2 stereo observing system

Fig. 6 illustrates the block diagram of the equipment of evaluating according to an exemplary embodiment of the present invention stereo observing system.

As shown in Figure 6, the equipment of evaluating according to an exemplary embodiment of the present invention stereo observing system comprises: three-dimensional geological model creating unit 10, seismic data acquisition unit 20, attribute determining unit 30 and evaluation coefficient determining unit 40.These unit can be realized by common hardware processors such as digital signal processor, field programmable gate arrays, also can realize by dedicated hardware processors such as special chips, also can partly or entirely by computer program, with software mode, realize, for example, be implemented as and be arranged in computing machine for evaluating the modules of the software of stereo observing system.

Particularly, three-dimensional geological model creating unit 10 creates the three-dimensional geological model in work area.Particularly, can utilize the existing real data in work area to create three-dimensional geological model.As example, the three-dimensional geological model in the work area of exemplary embodiment of the present can be as shown in Figure 2.

Seismic data acquisition unit 20 is simulated and is blown out according to current stereo observing system to be evaluated on described three-dimensional geological model, gathers the geological data of each common reflection point (CRP) bin in the full overlay area of described three-dimensional geological model.Here, the geological data of each CRP bin of collection comprises: the energy intensity receiving at CRP bin and the big gun that occurs to reflect in CRP bin are examined right position angle and geophone offset.

As example, the stereo observing system is here candidate's stereo observing system of choosing for exploration task and work area geological condition.By the stereo observing system of choosing is evaluated, can contribute to determine the stereo observing system in applicable work area.As example, four stereo observing systems as shown in table 1 are candidate's stereo observing systems of choosing for the work area shown in Fig. 2.

Particularly, the big gun that seismic data acquisition unit 20 can gather energy intensity that each CRP bin receives, big gun that reflection occurs in each CRP bin is examined right position angle and reflection occurred in each CRP bin is examined right geophone offset.As example, on the three-dimensional geological model shown in Fig. 2, according to stereo observing system 8L10S460R, 16L5S460R, 16L10S460R, 20L10S600R in table 1, simulate while blowing out respectively, the energy intensity receiving at each CRP bin can obtain from the distribution plan of the energy intensity receiving in full overlay area shown in Fig. 3 A to Fig. 3 D, at the big gun that reflection occurs in each CRP bin, examining right azimuthal distribution can be as shown in Fig. 4 A to Fig. 4 D, and the distribution of examining right geophone offset size at the big gun that reflection occurs in each CRP bin can be as shown in Fig. 5 A to Fig. 5 D.

The geological data that attribute determining unit 30 can gather according to data determining unit 20 determines that the energy properties of current stereo observing system to be evaluated, big gun examine position angle attribute and geophone offset attribute.

Particularly, attribute determining unit 30 can be analyzed current stereo observing system to be evaluated according to described geological data, to obtain the parameter that is used for respectively representing above-mentioned three attributes.For example, for energy properties, can analyze from aspects such as energy intensity, energy intensity distributed area, energy intensity distributing homogeneities; For big gun inspection, to position angle attribute, can analyze from aspects such as azimuthal distribution interval, azimuthal distribution homogeneity, azimuthal distribution similarities; For geophone offset attribute, can analyze from aspects such as geophone offset distributed area, geophone offset distributing homogeneities.As example, above-mentioned analysis can be qualitative analysis, quantitative test or both combinations.

Preferably, for energy properties, can analyze from energy intensity, energy intensity distributing homogeneity two aspects, to obtain for representing energy intensity and the energy even sex index of energy properties; For big gun inspection, to position angle attribute, can analyze from azimuthal distribution similarity aspect, to obtain for representing the position angle index similarity of big gun inspection to position angle attribute; For geophone offset attribute, can analyze from geophone offset distributing homogeneity aspect, to obtain for representing the geophone offset evenness index of geophone offset attribute.Below how detailed description is determined to energy intensity, energy even sex index, position angle index similarity and geophone offset evenness index based on geological data.

As example, attribute determining unit 30 can be determined above-mentioned energy intensity by above-mentioned calculating formula (1).Here, energy intensity is for reflecting the height of the energy receiving in the full overlay area of described three-dimensional geological model, the energy intensity E definite by above-mentioned calculating formula (1) is larger, represent that the energy receiving in the full overlay area of described three-dimensional geological model is higher, in the situation that other attributes are identical, the performance of stereo observing system is better.

In above-mentioned calculating formula (1), the mean value of the energy intensity that each CRP bin is received is as the energy intensity of stereo observing system, but those skilled in the art can understand, attribute determining unit 30 also can be determined by other computing method the energy intensity of stereo observing system, for example, the energy intensity that all CRP bins in the full overlay area of described three-dimensional geological model can be received is added summation, and the energy intensity using the summing value obtaining as stereo observing system.

As example, attribute determining unit 30 can be determined energy even sex index by above-mentioned calculating formula (2).Here, energy even sex index is for reflecting the situation that is uniformly distributed of the energy intensity that in the full overlay area of described three-dimensional geological model, each bin receives, the energy even sex index F definite by above-mentioned calculating formula (2) is less, represent that the energy intensity that in the full overlay area of described three-dimensional geological model, each bin receives distributes more evenly, in the situation that other attributes are identical, the performance of stereo observing system is better.

In above-mentioned calculating formula (2), the variance of the energy intensity that each CRP bin is received is as energy even sex index, but those skilled in the art can understand, and attribute determining unit 30 also can be determined energy even sex index by other computing method.

As example, attribute determining unit 30 can be determined position angle index similarity by above-mentioned calculating formula (3).Here, the big gun that reflection occurs for being reflected in each CRP bin in the full overlay area of three-dimensional geological model position angle index similarity is examined the similarity degree of right azimuthal distribution.The position angle index similarity C definite by above-mentioned calculating formula (3) is lower, be illustrated in reflection occurs each CRP bin big gun examine right azimuthal distribution must be more similar, in the situation that other attributes are identical, the performance of stereo observing system is better, correspondingly, utilize such stereo observing system can obtain the comparatively desirable geological data of signal to noise ratio (S/N ratio) and resolution.

Under actual geologic condition, particularly in western part of China Complex Mountain structure situation, due to big gun, to examine right azimuthal distribution extremely inhomogeneous, the big gun inspection of different stereo observing systems is differing very little to position angle aspect distributing homogeneity, be therefore difficult to azimuthal distributing homogeneity be evaluated to the stereo observing system in the work area of complex geological condition by definite big gun inspection.According to exemplary embodiment of the present invention, by using position angle index similarity as evaluating one of factor of stereo observing system, more can reflect different stereo observing systems in big gun inspection to the difference on the attribute of position angle.

Here, attribute determining unit 30 has adopted the computing method of above-mentioned calculating formula (3) to carry out computer azimuth angle index similarity, but those skilled in the art can understand, the additive method of the similarity that attribute determining unit 30 also can be by computational data is determined position angle index similarity.

As example, attribute determining unit 30 can be determined geophone offset evenness index by above-mentioned calculating formula (4).Here, the big gun that reflection occurs for being reflected in each CRP bin in the full overlay area of three-dimensional geological model geophone offset evenness index is examined the situation that is uniformly distributed of right geophone offset size.The geophone offset evenness index definite by above-mentioned calculating formula (4) more approaches ideal value, the big gun that reflection occurs each CRP bin in the full overlay area of three-dimensional geological model is examined the size distribution of right geophone offset must be more even, in the situation that other attributes are identical, the performance of stereo observing system is better.

Here, attribute determining unit 30 has adopted the computing method of above-mentioned calculating formula (4) to calculate geophone offset evenness index, but those skilled in the art can understand, attribute determining unit 30 also can be determined geophone offset evenness index by calculating the inhomogeneity additive method of difference between adjacent geophone offset.

Evaluation coefficient determining unit 40 can be determined the evaluation coefficient of current stereo observing system to be evaluated according to the energy properties of current stereo observing system to be evaluated, big gun inspection to position angle attribute and geophone offset attribute, wherein, described evaluation coefficient is used to indicate the performance of stereo observing system.

Those skilled in the art can understand, can consider according to the needs such as the geology characteristic in exploration task and work area because usually designing the funtcional relationship between evaluation coefficient and above-mentioned three attributes, thereby definite evaluation coefficient.

As example, in energy properties, by energy intensity and energy even sex index, represented, big gun inspection is represented by position angle index similarity position angle attribute, in the situation that geophone offset attribute is represented by geophone offset evenness index, evaluation coefficient determining unit 40 can be determined by above-mentioned objective function (5) or (6) evaluation coefficient of current stereo observing system to be evaluated.

Above-mentioned table 2 is listed the evaluating list of four stereo observing systems shown in the table 1 drawing according to above computing method, and evaluating comprises: energy intensity, energy even sex index, position angle index similarity, geophone offset evenness index and evaluation coefficient.Wherein, the energy intensity of stereo observing system 20L10S600R is maximum, the energy even sex index of stereo observing system 8L10S460R is minimum, the index similarity of stereo observing system 8L10S460R is minimum, the geophone offset evenness index of stereo observing system 20L10S600R approaches ideal value most, the evaluation coefficient of stereo observing system 8L10S460R is the highest, and performance is best.

In evaluating according to an exemplary embodiment of the present invention the method and apparatus of stereo observing system, geological data based on common reflection point (CRP) bin is evaluated, thereby can reflect more exactly the distribution of the effective reflection of buried target layer, be applicable to thus have the work area of relief surface and complex structure, effectively avoid the error of the geological data of CMP bin, and evaluation result is more accurate.

In addition, according in the method and apparatus of the evaluation stereo observing system of exemplary embodiment of the present invention, from energy properties, the big gun inspection of stereo observing system, position angle attribute and geophone offset attribute three aspects: are carried out to comprehensive evaluation to stereo observing system, thereby considered more all sidedly to relate to the factor of stereo observing system performance, made evaluation result more accurate.

It should be noted that above each embodiment of the present invention is only exemplary, the present invention is not limited to this.Those skilled in the art should understand that: without departing from the principles and spirit of the present invention, can change these embodiments, wherein, scope of the present invention limits in claim and equivalent thereof.

Claims (14)

1. a method of evaluating stereo observing system, is characterized in that, comprising:

(a) create the three-dimensional geological model in work area;

(b) on described three-dimensional geological model, according to current stereo observing system to be evaluated, simulate and blow out, gather the geological data of each common reflection point (CRP) bin in the full overlay area of described three-dimensional geological model, wherein, the geological data of each CRP bin comprises: the energy intensity receiving at CRP bin and the big gun that occurs to reflect in CRP bin are examined right position angle and geophone offset;

(c) according to described geological data, determine that the energy properties of current stereo observing system to be evaluated, big gun examine position angle attribute and geophone offset attribute;

(d) according to the energy properties of current stereo observing system to be evaluated, big gun inspection, position angle attribute and geophone offset attribute are determined to the evaluation coefficient of current stereo observing system to be evaluated, wherein, described evaluation coefficient is used to indicate the performance of stereo observing system.

2. method according to claim 1, it is characterized in that, described energy properties represents by energy intensity and energy even sex index, and described big gun inspection represents by position angle index similarity position angle attribute, and described geophone offset attribute is represented by geophone offset evenness index.

3. method according to claim 2, is characterized in that, in step (d), determines the evaluation coefficient of current stereo observing system to be evaluated by following objective function:

P＝E/E _max+F _min/F+C _min/C+V _sta/V，

Wherein, the evaluation coefficient of the current stereo observing system to be evaluated of P indication, the energy intensity of the current stereo observing system to be evaluated of E indication, the energy even sex index of the current stereo observing system to be evaluated of F indication, the position angle index similarity of the current stereo observing system to be evaluated of C indication, the geophone offset evenness index of the current stereo observing system to be evaluated of V indication, E _maxfor the maximal value in the energy intensity value of all stereo observing systems to be evaluated, F _minfor the minimum value in the energy even sex index of all stereo observing systems to be evaluated, C _minfor the minimum value in the position angle index similarity of all stereo observing systems to be evaluated, V _staideal value for geophone offset evenness index.

4. method according to claim 3, is characterized in that, in step (c), by following calculating formula, determines described energy intensity E:

$E=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l,$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

5. method according to claim 3, is characterized in that, in step (c), by following calculating formula, determines described energy even sex index F:

$F=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{(Q_l-\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l)}^2,$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

6. method according to claim 3, is characterized in that, in step (c), by following calculating formula, determines described position angle index similarity C:

$C=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\underset{\mathrm{\α}=1}{\overset{360}{\mathrm{\Σ}}}\frac{|H_l\left(\mathrm{\α}\right)-M\left(\mathrm{\α}\right)|}{H_l\left(\mathrm{\α}\right)+M\left(\mathrm{\α}\right)}\right),$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, α indicates angle, H _l(α) big gun that indication occurs to reflect and position angle is α in l CRP bin is examined right number,

7. method according to claim 3, is characterized in that, in step (c), by following calculating formula, determines described geophone offset evenness index V:

$V=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\frac{1}{N_l-1}\underset{i=1}{\overset{N_l-1}{\mathrm{\Σ}}}(1+|1-\frac{S_l(i+1)-S_l\left(i\right)}{S_{l0}}\left|\right)\right),$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, N _lfor there is the big gun of reflection in l CRP bin, examine right number, i indicates all big guns inspections that reflection occurs in same CRP bin sequence number when sorting from small to large by the size of geophone offset, S _l(i) i the big gun that indication occurs to reflect in l CRP bin examined right geophone offset, S _l(i+1) i+1 the big gun that indication occurs to reflect in l CRP bin examined right geophone offset,

8. an equipment of evaluating stereo observing system, is characterized in that, comprising:

Three-dimensional geological model creating unit, the three-dimensional geological model in establishment work area;

Seismic data acquisition unit, on described three-dimensional geological model, according to current stereo observing system to be evaluated, simulate and blow out, gather the geological data of each common reflection point (CRP) bin in the full overlay area of described three-dimensional geological model, wherein, the geological data of each CRP bin comprises: the energy intensity receiving at CRP bin and the big gun that occurs to reflect in CRP bin are examined right position angle and geophone offset;

Attribute determining unit, determines that according to described geological data the energy properties of current stereo observing system to be evaluated, big gun examine position angle attribute and geophone offset attribute;

Evaluation coefficient determining unit, according to the energy properties of current stereo observing system to be evaluated, big gun inspection, position angle attribute and geophone offset attribute are determined to the evaluation coefficient of current stereo observing system to be evaluated, wherein, described evaluation coefficient is used to indicate the performance of stereo observing system.

9. equipment according to claim 8, it is characterized in that, described energy properties represents by energy intensity and energy even sex index, and described big gun inspection represents by position angle index similarity position angle attribute, and described geophone offset attribute is represented by geophone offset evenness index.

10. equipment according to claim 9, is characterized in that, evaluation coefficient determining unit is determined the evaluation coefficient of current stereo observing system to be evaluated by following objective function:

P＝E/E _max+F _min/F+C _min/C+V _sta/V，

Wherein, the evaluation coefficient of the current stereo observing system to be evaluated of P indication, the energy intensity of the current stereo observing system to be evaluated of E indication, the energy even sex index of the current stereo observing system to be evaluated of F indication, the position angle index similarity of the current stereo observing system to be evaluated of C indication, the geophone offset evenness index of the current stereo observing system to be evaluated of V indication, E _maxfor the maximal value in the energy intensity value of all stereo observing systems to be evaluated, F _minfor the minimum value in the energy even sex index of all stereo observing systems to be evaluated, C _minfor the minimum value in the position angle index similarity of all stereo observing systems to be evaluated, V _staideal value for geophone offset evenness index.

11. equipment according to claim 10, is characterized in that, attribute determining unit is determined described energy intensity E by following calculating formula:

$E=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l,$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

12. equipment according to claim 10, is characterized in that, attribute determining unit is determined described energy even sex index F by following calculating formula:

$F=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{(Q_l-\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}{Q}_l)}^2,$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, Q _lthe energy intensity that indication receives at l CRP bin.

13. equipment according to claim 10, is characterized in that, attribute determining unit is determined described position angle index similarity C by following calculating formula:

$C=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\underset{\mathrm{\α}=1}{\overset{360}{\mathrm{\Σ}}}\frac{|H_l\left(\mathrm{\α}\right)-M\left(\mathrm{\α}\right)|}{H_l\left(\mathrm{\α}\right)+M\left(\mathrm{\α}\right)}\right),$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, α indicates angle, H _l(α) big gun that indication occurs to reflect and position angle is α in l CRP bin is examined right number,

14. equipment according to claim 10, is characterized in that, attribute determining unit is determined described geophone offset evenness index V by following calculating formula:

$V=\frac{1}{K}\underset{l=1}{\overset{K}{\mathrm{\Σ}}}\left(\frac{1}{N_l-1}\underset{i=1}{\overset{N_l-1}{\mathrm{\Σ}}}(1+|1-\frac{S_l(i+1)-S_l\left(i\right)}{S_{l0}}\left|\right)\right),$

Wherein, the quantity of the CRP bin in the full overlay area that K is described three-dimensional geological model, the numbering of l indication CRP bin, N _lfor there is the big gun of reflection in l CRP bin, examine right number, i indicates all big guns inspections that reflection occurs in same CRP bin sequence number when sorting from small to large by the size of geophone offset, S _l(i) i the big gun that indication occurs to reflect in l CRP bin examined right geophone offset, S _l(i+1) i+1 the big gun that indication occurs to reflect in l CRP bin examined right geophone offset,