CN104730581A - Method for locating microseism event point - Google Patents

Abstract

The invention provides a method for locating a microseism event point. The method comprises the steps that a, based on acoustic logging data of a monitoring well and the longitudinal wave first arrival time, monitored by a ground detector, of a preset emitting hole, a longitudinal wave speed model about the included angle between the connecting line of the microseism event point and the position of the ground detector and the vertical direction and the projection coordinates of the microseism event point on the horizontal plane is built; b, the position of the microseism event point to be located is obtained according to the built longitudinal wave speed model and the longitudinal wave first arrival time, monitored by the ground detector, of the microseism event point to be located in an iteration mode. According to the method, the microseism event point can be precisely located in a fracture well area with a complex subsurface structure or severe transverse speed changes.

Description

The method of location micro-seismic event point

Technical field

The present invention relates to microearthquake ground monitoring technical field, more particularly, relate to a kind of method of locating micro-seismic event point.

Background technology

At present, the micro-seismic event point location mode usually adopted in waterfrac treatment microearthquake ground monitoring technical field comprises: the mode of focus spacescan stack power playback, use the modes such as even speed model or horizontal layer rate pattern position.But, the mode of existing location micro-seismic event point can meet the needs of production in comparatively simple, that lateral variation in velocity the is less pressure break wellblock of underground structure, then there is comparatively significantly error in comparatively complicated or that lateral variation in velocity is the comparatively violent pressure break wellblock of underground structure.

Summary of the invention

Exemplary embodiment of the present invention is to provide a kind of method of locating micro-seismic event point, and it can solve in location, the pressure break wellblock micro-seismic event point problem not accurately that underground structure is comparatively complicated or lateral variation in velocity is comparatively violent.

The invention provides a kind of method of locating micro-seismic event point, comprise: the compressional wave first break time of the predetermined perforation that (a) monitors based on Sonic Logging Data and the surface geophone of monitor well, set up the velocity of longitudinal wave model about angle formed by the line of the position of micro-seismic event point and surface geophone and vertical direction and micro-seismic event point projection coordinate in the horizontal plane; The compressional wave first break time iteration of b micro-seismic event point to be positioned that () monitors according to the velocity of longitudinal wave model set up and surface geophone asks for the position of micro-seismic event point to be positioned.

Alternatively, step (a) comprising: (a1) obtains initial compressional wave recurrence speed according to the Sonic Logging Data of monitor well; (a2) the compressional wave first break time of the described predetermined perforation monitored according to initial compressional wave recurrence speed, the position of described predetermined perforation, surface geophone determines the velocity of longitudinal wave model about angle formed by the line of the position of micro-seismic event point and surface geophone and vertical direction and micro-seismic event point projection coordinate in the horizontal plane.

Alternatively, the expression formula of described velocity of longitudinal wave model is:

v＝v ₀(1+k ₁x)(1+k ₂y)(1+k ₃θ+k ₄θ ²)，

Wherein, the velocity of longitudinal wave of v instruction from micro-seismic event point to surface geophone, v ₀for initial compressional wave recurrence speed, (x, y) indicates micro-seismic event point projection coordinate in the horizontal plane, and θ indicates angle formed by the line of the position of micro-seismic event point and described surface geophone and vertical direction, k ₁, k ₂, k ₃, k ₄for coefficient.

Alternatively, k ₁, k ₂, k ₃, k ₄by using the position of described predetermined perforation as (x, y) expression formula that, described predetermined perforation substitutes into described velocity of longitudinal wave model as angle formed by the line of the position of v, described predetermined perforation and described surface geophone and vertical direction as θ to the velocity of longitudinal wave of surface geophone respectively calculates, wherein, described predetermined perforation is to the velocity of longitudinal wave of surface geophone, and the compressional wave first break time of the described predetermined perforation monitored to the distance of described predetermined perforation divided by described surface geophone by described surface geophone calculates.

Alternatively, step (b) comprising: the compressional wave first break time iteration of micro-seismic event point to be positioned monitored according to the velocity of longitudinal wave model set up and surface geophone ask for target function value minimum time iteration stepping Δ, when the iteration stepping Δ asked for is less than predetermined threshold, stop iteration using the coordinate of the coordinate of the current micro-seismic event point determined as micro-seismic event to be positioned point, wherein, objective function is:

$f\left(\mathrm{\Δ}\right)=\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2,$

Wherein, Δ=[Δ x, Δ y, Δ z] ^t, $A_{m,n}=[\frac{x-x_m}{d_m{v}_m}-\frac{x-x_n}{d_n{v}_n},\frac{y-y_m}{d_m{v}_m}-\frac{y-y_n}{d_n{v}_n},\frac{z-z_m}{d_m{v}_m}-\frac{z-z_n}{d_n{v}_n}],$ B _m,n=(T _m-T _n)-(t _m-t _n), (x, y, z) indicates the coordinate of the current micro-seismic event point determined, (x _m, y _m, z _m) indicate the coordinate of m surface geophone, (x _n, y _n, z _n) indicate the coordinate of the n-th surface geophone, d _m, d _nindicate current location to the distance of m, the n-th surface geophone respectively, v _m, v _nindicate respectively by described velocity of longitudinal wave model calculate from the current micro-seismic event determined o'clock to m, the velocity of longitudinal wave of the n-th surface geophone, R is the quantity of surface geophone, T _m, T _nthe compressional wave first break time of the micro-seismic event point to be positioned indicating m surface geophone and the n-th surface geophone to monitor respectively, (x+ Δ x, y+ Δ y, z+ Δ z) coordinate of micro-seismic event point determined by current iteration of instruction, wherein, (x, y during primary iteration, z) be the coordinate of perforation, (x, y, z) during follow-up each iteration is (the x+ Δ x determined by last iteration, y+ Δ y, z+ Δ z).

Alternatively, when the quantity of perforation is more than or equal to 2, objective function is:

$f\left(\mathrm{\Δ}\right)=\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2+\mathrm{\Ω}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=1,j\≠i}{\overset{N}{\mathrm{\Σ}}}{\left|\right|\mathrm{\Δz}-z_i+(z_j-z_i)\·r-z\left|\right|}_2^2,$

Wherein, N is the quantity of perforation, and Ω is weight coefficient, z _i, z _jindicate i-th, the degree of depth of a jth perforation respectively, r be i-th, the slope of the line of the position of a jth perforation.

Alternatively, during primary iteration, the value of Ω is 0, and the value of Ω when determining next iteration at the end of each iteration, during next iteration, the value of Ω meets following formula:

$\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2/\mathrm{\Ω}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=1,j\≠i}{\overset{N}{\mathrm{\Σ}}}{\left|\right|\mathrm{\Δz}-z_i+(z_j-z_i)\·r-z\left|\right|}_2^2\≈1.$

Alternatively, k is calculated ₁, k ₂, k ₃, k ₄time the described predetermined perforation that uses and primary iteration time the perforation that uses be same perforation.

Alternatively, predetermined threshold is 5 meters-10 meters.

Locate the method for micro-seismic event point according to an exemplary embodiment of the present invention, micro-seismic event point can be located accurately in comparatively complicated or that lateral variation in velocity is the comparatively violent pressure break wellblock of underground structure, for microearthquake ground monitoring technology provides strong support in the application in complex structure area, expand the range of application of microearthquake ground monitoring technology.

Part in ensuing description is set forth general plotting of the present invention other in and/or advantage, some will be clearly by describing, or can learn through the enforcement of general plotting of the present invention.

Accompanying drawing explanation

By below in conjunction with exemplarily illustrating the description that the accompanying drawing of embodiment carries out, the above and other object of exemplary embodiment of the present and feature will become apparent, wherein:

Fig. 1 illustrates the process flow diagram of the method for locating micro-seismic event point according to an exemplary embodiment of the present invention;

Fig. 2 illustrates that iteration asks for the schematic diagram of the position of micro-seismic event point to be positioned according to an exemplary embodiment of the present invention.

Embodiment

Now will 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 by referring to accompanying drawing, described embodiment will be described, to explain the present invention.

Fig. 1 illustrates the process flow diagram of the method for locating micro-seismic event point according to an exemplary embodiment of the present invention.

As shown in Figure 1, in step S10, based on the compressional wave first break time of the predetermined perforation that Sonic Logging Data and the surface geophone of monitor well monitor, set up the velocity of longitudinal wave model about angle formed by the line of the position of micro-seismic event point and surface geophone and vertical direction and micro-seismic event point projection coordinate in the horizontal plane.That is, set up utilize angle formed by the line of the position of micro-seismic event point and surface geophone and vertical direction, micro-seismic event point projection coordinate in the horizontal plane describes anisotropic velocity model from micro-seismic event point to the velocity of longitudinal wave of surface geophone jointly.

Particularly, first initial compressional wave recurrence speed can be obtained according to the Sonic Logging Data of monitor well.Should be appreciated that, existing various mode can be used to obtain initial compressional wave recurrence speed according to the Sonic Logging Data of monitor well, do not repeat them here.

Then, the compressional wave first break time of the described predetermined perforation that can monitor according to initial compressional wave recurrence speed, the position of described predetermined perforation, surface geophone determines the velocity of longitudinal wave model about angle formed by the line of the position of micro-seismic event point and surface geophone and vertical direction and micro-seismic event point projection coordinate in the horizontal plane.

Due to longitudinal wave propagation speed, and insensitive for oil gas etc., so vertical wave propagation is insensitive to changes in coordinates, comparatively responsive to angle change.Therefore, preferably, the expression formula of velocity of longitudinal wave model can be:

v＝v ₀(1+k ₁x)(1+k ₂y)(1+k ₃θ+k ₄θ ²) (1)，

Wherein, the velocity of longitudinal wave of v instruction from micro-seismic event point to surface geophone, v ₀for initial compressional wave recurrence speed, (x, y) indicates micro-seismic event point projection coordinate in the horizontal plane, and θ indicates angle formed by the line of the position of micro-seismic event point and this surface geophone and vertical direction, k ₁, k ₂, k ₃, k ₄for coefficient.

The expression formula of above-mentioned velocity of longitudinal wave model embodies the characteristic of longitudinal wave propagation well, thus can calculate the velocity of longitudinal wave of micro-seismic event point to surface geophone exactly by above-mentioned velocity of longitudinal wave model.

Exemplarily, k ₁, k ₂, k ₃, k ₄by using the position of described predetermined perforation as (x, y), described predetermined perforation to the velocity of longitudinal wave of surface geophone as v, the expression formula that substitutes into described velocity of longitudinal wave model respectively as θ of angle formed by the line of the position of described predetermined perforation and described surface geophone and vertical direction calculates.Wherein, described predetermined perforation is to the velocity of longitudinal wave of surface geophone, and the compressional wave first break time of the described predetermined perforation that can be monitored divided by described surface geophone to the distance of described predetermined perforation by described surface geophone calculates.

Particularly, the compressional wave first break time of the described predetermined perforation that the position of known predetermined perforation, surface geophone monitor, the position of surface geophone, formed by the described predetermined perforation position of described predetermined perforation, the compressional wave first break time of described predetermined perforation that monitors according to surface geophone calculated to the line of the position of the velocity of longitudinal wave of surface geophone, described predetermined perforation and described surface geophone and vertical direction, angle substitutes into x, y, v, θ respectively successively to form multiple equation, and solving equations can obtain k ₁, k ₂, k ₃, k ₄value.Such as, the expression formula that angle formed by the horizontal level of described predetermined perforation, described predetermined perforation to the line of the position of the velocity of longitudinal wave of a surface geophone, described predetermined perforation and this surface geophone and vertical direction can be substituted into velocity of longitudinal wave model obtains one about k ₁, k ₂, k ₃, k ₄equation, correspondingly, the expression formula that angle formed by the horizontal level of described predetermined perforation, described predetermined perforation to the line of the position of the velocity of longitudinal wave of another surface geophone, described predetermined perforation and another surface geophone and vertical direction can be substituted into velocity of longitudinal wave model obtains another about k ₁, k ₂, k ₃, k ₄equation, correspondingly, obtain in the manner described above some about k ₁, k ₂, k ₃, k ₄equation.

In addition, should be appreciated that, also by other applicable mode determination coefficient k ₁, k ₂, k ₃, k ₄, this is not restricted.

Those skilled in the art should understand that, formula (1) is only as preferred example, the compressional wave first break time of the perforation that other the Sonic Logging Data based on monitor well and surface geophone monitor determine and utilize angle formed by the line of the position of micro-seismic event point and surface geophone and vertical direction, expression formula that micro-seismic event point projection coordinate in the horizontal plane describes jointly from micro-seismic event point to the anisotropic velocity model of the velocity of longitudinal wave of surface geophone all can be applicable to the present invention, this is not restricted.

In step S20, the compressional wave first break time iteration of micro-seismic event point to be positioned monitored according to the velocity of longitudinal wave model set up and surface geophone asks for the position of micro-seismic event point to be positioned.

Fig. 2 illustrates that iteration asks for the schematic diagram of the position of micro-seismic event point to be positioned according to an exemplary embodiment of the present invention.

As shown in Figure 2, (x, y, z) is the coordinate of the current micro-seismic event point determined, utilizes the current micro-seismic event point determined to the compressional wave first break time t of Different Ground wave detector _m, t _n, the current micro-seismic event point determined is to the distance d of Different Ground wave detector _m, d _n, and the compressional wave first break time T of micro-seismic event point to be positioned that Different Ground wave detector monitors _m, T _n(not shown in Fig. 2) asks for iteration stepping Δ=[Δ x, Δ y, Δ z] ^tby (the x+ Δ x that current iteration calculates, y+ Δ y, z+ Δ z), be (x, y during next iteration, z), until when iteration stepping Δ is less than predetermined threshold, stop iteration, and using the coordinate of the coordinate of the current micro-seismic event point determined as micro-seismic event to be positioned point.Wherein, the current micro-seismic event point determined is to the compressional wave first break time t of Different Ground wave detector _mand t _n, by the velocity of longitudinal wave v from the current micro-seismic event point determined to different surface geophones calculated by velocity of longitudinal wave model _m, v _nand d _m, d _ncalculate.In addition, should be appreciated that, Fig. 2 only illustrates m surface geophone and the n-th surface geophone example as Different Ground wave detector, and Different Ground wave detector is not limited thereto.

Exemplarily, the compressional wave first break time iteration of micro-seismic event point to be positioned monitored according to the velocity of longitudinal wave model set up and surface geophone ask for target function value minimum time iteration stepping Δ, when the iteration stepping Δ asked for is less than predetermined threshold, stop iteration using the coordinate of the coordinate of the current micro-seismic event point determined as micro-seismic event to be positioned point.

Exemplarily, objective function can be:

$f\left(\mathrm{\Δ}\right)=\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2---\left(2\right),$

Wherein, Δ=[Δ x, Δ y, Δ z] ^t, $A_{m,n}=[\frac{x-x_m}{d_m{v}_m}-\frac{x-x_n}{d_n{v}_n},\frac{y-y_m}{d_m{v}_m}-\frac{y-y_n}{d_n{v}_n},\frac{z-z_m}{d_m{v}_m}-\frac{z-z_n}{d_n{v}_n}],$ B _m,n=(T _m-T _n)-(t _m-t _n), (x, y, z) indicates the coordinate of the current micro-seismic event point determined, (x _m, y _m, z _m) indicate the coordinate of m surface geophone, (x _n, y _n, z _n) indicate the coordinate of the n-th surface geophone, d _m, d _nindicate current location to the distance of m, the n-th surface geophone respectively, v _m, v _nindicate respectively by described velocity of longitudinal wave model calculate from the current micro-seismic event determined o'clock to m, the velocity of longitudinal wave of the n-th surface geophone, R is the quantity of surface geophone, T _m, T _nthe compressional wave first break time of the micro-seismic event point to be positioned indicating m surface geophone and the n-th surface geophone to monitor respectively, the coordinate of the micro-seismic event point that (x+ Δ x, y+ Δ y, z+ Δ z) instruction is determined by current iteration.

Here, the coordinate that (x, y, z) during primary iteration is perforation, (x, y, z) during follow-up each iteration is (x+ Δ x, y+ Δ y, the z+ Δ z) that determined by last iteration.

Exemplarily, k is calculated ₁, k ₂, k ₃, k ₄time the predetermined perforation that uses can be same perforation with the perforation used during primary iteration, calculate micro-seismic event point to the applicability of the velocity of longitudinal wave of surface geophone and accuracy to improve by velocity of longitudinal wave model.

Exemplarily, when the quantity of perforation is more than or equal to 2, add macroscopical bound term in objective function, that is, the objective function with the constraint of degree of depth macroscopic view is:

$f\left(\mathrm{\Δ}\right)=\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2+\mathrm{\Ω}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=1,j\≠i}{\overset{N}{\mathrm{\Σ}}}{\left|\right|\mathrm{\Δz}-z_i+(z_j-z_i)\·r-z\left|\right|}_2^2---\left(3\right),$

Wherein, N is the quantity of perforation, and Ω is weight coefficient, z _i, z _jindicate i-th, the degree of depth of a jth perforation respectively, r be i-th, the slope of the line of the position of a jth perforation.

Here, exemplarily, during primary iteration, the value of Ω is 0, and the value of Ω when determining next iteration at the end of each iteration, during next iteration, the value of Ω meets following formula:

$\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2/\mathrm{\Ω}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=1,j\≠i}{\overset{N}{\mathrm{\Σ}}}{\left|\right|\mathrm{\Δz}-z_i+(z_j-z_i)\·r-z\left|\right|}_2^2\≈1---\left(4\right).$

Exemplarily, predetermined threshold can be 5-10, and unit can be rice.Should be understood that and can suitably adjust predetermined threshold and corresponding unit according to actual conditions.

In addition, those skilled in the art should understand that, formula (2) and formula (3) are only as preferred example, the compressional wave first break time of other the micro-seismic event point to be positioned utilizing the current micro-seismic event point determined to monitor to the velocity of longitudinal wave of Different Ground wave detector and distance, Different Ground wave detector is all applicable to the present invention to the objective function asking for iteration stepping, further, the mode be also applicable to by other sets the value of Ω.

In addition, the said method according to exemplary embodiment of the present invention may be implemented as computer program, thus when running this program, realizes said method.In addition, said method also performs by the device for locating micro-seismic event point.

Locate the method for micro-seismic event point according to an exemplary embodiment of the present invention, micro-seismic event point can be located accurately in comparatively complicated or that lateral variation in velocity is the comparatively violent pressure break wellblock of underground structure, for microearthquake ground monitoring technology provides strong support in the application in complex structure area, expand the range of application of microearthquake ground monitoring technology.

Although show and described exemplary embodiments more of the present invention, but those skilled in the art should understand that, when not departing from by the principle of the present invention of claim and its scope of equivalents thereof and spirit, can modify to these embodiments.

Claims (9)

1. locate a method for micro-seismic event point, comprising:

The compressional wave first break time of a predetermined perforation that () monitors based on Sonic Logging Data and the surface geophone of monitor well, sets up the velocity of longitudinal wave model about angle formed by the line of the position of micro-seismic event point and surface geophone and vertical direction and micro-seismic event point projection coordinate in the horizontal plane;

The compressional wave first break time iteration of b micro-seismic event point to be positioned that () monitors according to the velocity of longitudinal wave model set up and surface geophone asks for the position of micro-seismic event point to be positioned.

2. method according to claim 1, wherein, step (a) comprising:

(a1) initial compressional wave recurrence speed is obtained according to the Sonic Logging Data of monitor well;

(a2) the compressional wave first break time of the described predetermined perforation monitored according to initial compressional wave recurrence speed, the position of described predetermined perforation, surface geophone determines the velocity of longitudinal wave model about angle formed by the line of the position of micro-seismic event point and surface geophone and vertical direction and micro-seismic event point projection coordinate in the horizontal plane.

3. method according to claim 2, wherein, the expression formula of described velocity of longitudinal wave model is:

v＝v ₀(1+k ₁x)(1+k ₂y)(1+k ₃θ+k ₄θ ²)，

Wherein, the velocity of longitudinal wave of v instruction from micro-seismic event point to surface geophone, v ₀for initial compressional wave recurrence speed, (x, y) indicates micro-seismic event point projection coordinate in the horizontal plane, and θ indicates angle formed by the line of the position of micro-seismic event point and described surface geophone and vertical direction, k ₁, k ₂, k ₃, k ₄for coefficient.

4. method according to claim 3, wherein, k ₁, k ₂, k ₃, k ₄by using the position of described predetermined perforation as (x, y) expression formula that, described predetermined perforation substitutes into described velocity of longitudinal wave model as angle formed by the line of the position of v, described predetermined perforation and described surface geophone and vertical direction as θ to the velocity of longitudinal wave of surface geophone respectively calculates

Wherein, described predetermined perforation is to the velocity of longitudinal wave of surface geophone, and the compressional wave first break time of the described predetermined perforation monitored to the distance of described predetermined perforation divided by described surface geophone by described surface geophone calculates.

5. method according to claim 4, wherein, step (b) comprising:

The compressional wave first break time iteration of micro-seismic event point to be positioned monitored according to the velocity of longitudinal wave model set up and surface geophone ask for target function value minimum time iteration stepping Δ, when the iteration stepping Δ asked for is less than predetermined threshold, stop iteration using the coordinate of the coordinate of the current micro-seismic event point determined as micro-seismic event to be positioned point

Wherein, objective function is:

$f\left(\mathrm{\Δ}\right)=\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2,$

Wherein, Δ=[Δ x, Δ y, Δ z] ^t, $A_{m,n}=[\frac{x-x_m}{d_m{v}_m}-\frac{x-x_n}{d_n{v}_n},\frac{y-y_m}{d_m{v}_m}-\frac{y-y_n}{d_n{v}_n},\frac{z-z_m}{d_m{v}_m}-\frac{z-z_n}{d_n{v}_n}],$ B _m,n=(T _m-T _n)-(t _m-t _n), (x, y, z) indicates the coordinate of the current micro-seismic event point determined, (x _m, y _m, z _m) indicate the coordinate of m surface geophone, (x _n, y _n, z _n) indicate the coordinate of the n-th surface geophone, d _m, d _nindicate current location to the distance of m, the n-th surface geophone respectively, v _m, v _nindicate respectively by described velocity of longitudinal wave model calculate from the current micro-seismic event determined o'clock to m, the velocity of longitudinal wave of the n-th surface geophone, R is the quantity of surface geophone, T _m, T _nthe compressional wave first break time of the micro-seismic event point to be positioned indicating m surface geophone and the n-th surface geophone to monitor respectively, the coordinate of the micro-seismic event point that (x+ Δ x, y+ Δ y, z+ Δ z) instruction is determined by current iteration,

Wherein, the coordinate that (x, y, z) during primary iteration is perforation, (x, y, z) during follow-up each iteration is (x+ Δ x, y+ Δ y, the z+ Δ z) that determined by last iteration.

6. method according to claim 5, wherein, when the quantity of perforation is more than or equal to 2, objective function is:

$f\left(\mathrm{\Δ}\right)=\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2+\mathrm{\Ω}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=1,j\≠i}{\overset{N}{\mathrm{\Σ}}}{\left|\right|\mathrm{\Δz}-z_i+(z_j-z_i)\·r-z\left|\right|}_2^2,$

Wherein, N is the quantity of perforation, and Ω is weight coefficient, z _i, z _jindicate i-th, the degree of depth of a jth perforation respectively, r be i-th, the slope of the line of the position of a jth perforation.

7. method according to claim 6, wherein, during primary iteration, the value of Ω is 0, and the value of Ω when determining next iteration at the end of each iteration, during next iteration, the value of Ω meets following formula:

$\underset{n=1}{\overset{R}{\mathrm{\Σ}}}\underset{m=1,m\≠n}{\overset{R}{\mathrm{\Σ}}}{\left|\right|A_{m,n}\mathrm{\Δ}-b_{m,n}\left|\right|}_2^2/\mathrm{\Ω}\underset{i=1}{\overset{N}{\mathrm{\Σ}}}\underset{j=1,j\≠i}{\overset{N}{\mathrm{\Σ}}}{\left|\right|\mathrm{\Δz}-z_i+(z_j-z_i)\·r-z\left|\right|}_2^2\≈1.$

8. method according to claim 5, wherein, calculates k ₁, k ₂, k ₃, k ₄time the described predetermined perforation that uses and primary iteration time the perforation that uses be same perforation.

9. method according to claim 5, wherein, predetermined threshold is 5 meters-10 meters.