CN1517720A - Space correction method for 2D seisimc overlapping shift profile drawing - Google Patents

Abstract

A space correction method for 2D seismic overlapping shift profile drawing includes such steps as choosing a reflection interface of underground object layer, taking a CDP point on the interface, drawing a normal plane perpendicular to the interface by a measuring line and said CDP point, drawing another normal plane by another measuring line and said CDP point, intersecting between said two normal planes, and deriving and calculating the vertical distance from the CDP point to ground surface and the distance from the point of the CDP point reflected from two normal plant to ground surface and the real position of the CDP point on ground surface. Its advantages are high efficiency and high precision.

Description

2D seisimc overlapping shift profile becomes the empty calibration method of figure

Technical field:

The present invention relates to seismic data interpretation method in the oil seismic exploration, belong to 2D seisimc overlapping shift profile one-tenth and scheme empty calibration method.

Background technology:

2D seisimc overlapping shift profile is employed aborning widely in recent years, and the fault interpretation precision is guaranteed, and the ratio of precision of structural map had had more significantly raising in the past like this.But 2D seisimc overlapping shift profile just playbacks earthquake reflected wave along the empty school of two-dimentional line direction, only is along the empty school playback of two-dimentional line direction, and does not do the playback of empty school at vertical line direction, and two-dimentional migrated stack section is still plane normal, the T that two survey line intersection point places read ₀Time is unequal, can not be closed, and bring difficulty for the making structural map.

At present, among the interpre(ta)tive system Zycor of introduction, Landmark or the Geoquest, become between empty school of figure and the interactive interpretation layer bit data transmission speed slow, can not widespread use; And add up again by the strata division playback, operation is loaded down with trivial details, and is inapplicable; Present 2D seisimc overlapping shift profile one-tenth is schemed empty calibration method mainly based on manual empty school, and manual empty calibration method is only to T such as 2D seisimc overlapping shift profiles ₀Figure is along the empty school of vertical main profile direction, along the vertical not empty school of cross-track direction.This method gained structural map mapping precision is low, labour intensity big, work efficiency is also low.

Summary of the invention:

Scheme the problem that empty calibration method mapping precision is low, labour intensity big, work efficiency is also low in order to overcome existing two-dimensional earthquake migrated stack section one-tenth, the invention provides a kind of 2D seisimc overlapping shift profile one-tenth and scheme empty calibration method, this method is to finish on seismic data is finished dealing with the basis that is offset, can improve degree of depth structural map mapping precision, realize the automatic empty school of computing machine, reduce the personnel's of explaining labour intensity, improved work efficiency.

The technical solution adopted for the present invention to solve the technical problems is: utilize computing machine to finish automatically by following step:

A, S are ground level, and P is underground zone of interest reflecting interface, and the reflecting interface inclination angle is ;

The I plane was the plane normal of survey line 1 perpendicular to P, handed over plane P in 1 ' line, and 1 ' line apparent dip is  ₁

The II plane was the plane normal of survey line 2 perpendicular to P, handed over plane P in 2 ' line, and apparent dip is  ₂, 1 ' line and 2 ' line meet at the D point, and this point is exactly a CDP point;

1 line and 2 lines are two seismic line that intersect arbitrarily on the ground, meet at the O point, and angle is ω, ω ₁And ω ₂Be respectively the angle of 1,2 lines and stratum updip direction;

The C point is respectively the projection to 2 lines and 1 line in D o'clock for the subpoint of D point to ground, A, B point, and their distance is respectively h ₂, h ₁, the distance that the D point is ordered to O is h ₀, L ₁, L ₂Be respectively the distance that the C point is ordered to B, A, H is the distance that the D point is ordered to C;

Can derive according to earthquake geometry:

h ₁＝h ₀cos ₁ 1-3-1

h ₂＝h ₀cos ₂ 1-3-2

l ₁＝h ₀sin ₁ 1-3-3

l ₂＝h ₀sin ₂ 1-3-4

L ₁＝h ₁sin′ ₂ 1-3-5

L ₂＝h ₂sin′ ₁ 1-3-6

H＝h ₁cos′ ₂ 1-3-7

H＝h ₂cos′ ₁ 1-3-8

In formula 1-3-1,1-3-2,1-3-3,1-3-4,1-3-5,1-3-6,1-3-7,1-3-8, h ₁, h ₂,  ₁,  ₂All can calculate  ' at migrated section ₁Be ∠ CDA, the  ' among the right-angle triangle CDA ₂Be the ∠ CDB among the right-angle triangle CDB.

Can derive in conjunction with Fig. 1 and top formula:

$\sin {{\mathrm{\φ}}^{\′}}_2=\frac{{\sin \mathrm{\φ}}_2}{{\cos \mathrm{\φ}}_1}\frac{{\sin \mathrm{\ω}}_1}{\cos (\mathrm{\ω}-{\mathrm{\ω}}_1)}---1-3-9$

$\sin {{\mathrm{\φ}}^{\′}}_1=\frac{{\sin \mathrm{\φ}}_1}{{\cos \mathrm{\φ}}_2}\frac{{\sin \mathrm{\ω}}_2}{\cos (\mathrm{\ω}-{\mathrm{\ω}}_2)}---1-3-10$

When the vertical survey line 2 of survey line 1, i.e. ω=90 °:

$\sin {{\mathrm{\φ}}^{\′}}_2=\frac{{\sin \mathrm{\φ}}_2}{{\cos \mathrm{\φ}}_1}---1-3-11$

$\sin {{\mathrm{\φ}}^{\′}}_1=\frac{{\sin \mathrm{\φ}}_1}{{\cos \mathrm{\φ}}_2}---1-3-12$

When the main profile interconnection is any angle, can push away ω ₁, ω ₂Formula:

${\mathrm{tg\ω}}_1=\frac{{\sin \mathrm{\φ}}_2-{\sin \mathrm{\φ}}_1\cos \mathrm{\ω}}{{\sin \mathrm{\φ}}_1\sin \mathrm{\ω}}---1-3-13$

${\mathrm{tg\ω}}_2=\frac{{\sin \mathrm{\φ}}_1-{\sin \mathrm{\φ}}_2\cos \mathrm{\ω}}{\sin {\mathrm{\φ}}_2\sin \mathrm{\ω}}---1-3-14$

ω is the angle of two surveys line in the formula, can directly calculate.

Utilize formula 1-3-9,1-3-10, and the result of calculation of formula 1-3-13,1-3-14 be updated to formula 1-3-5,1-3-6,1-3-7,1-3-8, in can calculate the side-play amount L in sky school ₁, L ₂And the vertical depth H value behind the empty school, promptly finish in empty school; Can draw structural map with the H value.

The invention has the beneficial effects as follows: because this 2D seisimc overlapping shift profile one-tenth schemes the data creating structural map that empty calibration method has utilized main profile and cross-track simultaneously, lost data point not; What draw is the true depth structural map, and structural feature is true, and high some position is accurate; The automatic empty school that uses a computer explains that personnel labor intensity is little, efficient is high, mapping precision is high.

Description of drawings:

The invention will be further described below in conjunction with accompanying drawing.

Fig. 1 is self excitation and self receiving raypath spatial relationship figure of the present invention;

Fig. 2 is fragrant 148 wellblock T ₂Layer is schemed depth value behind the empty calibration method sky school with this 2D seisimc overlapping shift profile one-tenth and is projected to scatter diagram on the survey grid;

Fig. 3 is the degree of depth structural map that generates automatically after Fig. 2 is playbacked to main profile and contact side line on computers;

Fig. 4 is the tame defensive wall T of Xu of pine the Liao Dynasty ₄Layer t ₀Scatter diagram;

Fig. 5 is the degree of depth scatter diagram that Fig. 4 is obtained after with the empty school of this method.

Embodiment:

Fig. 1 is self excitation and self receiving raypath spatial relationship figure, in accompanying drawing 1,

A, to establish D be ground level, and P is underground zone of interest reflecting interface, and the reflecting interface inclination angle is ;

The I plane was the plane normal of survey line 1 perpendicular to P, handed over plane P in 1 ' line, and 1 ' line apparent dip is  ₁

The II plane was the plane normal of survey line 2 perpendicular to P, handed over plane P in 2 ' line, and apparent dip is  ₂, 1 ' line and 2 ' line meet at the D point, and this point is exactly a CDP point;

1 line and 2 lines are two seismic line that intersect arbitrarily on the ground, meet at the O point, and angle is ω, ω ₁And ω ₂Be respectively the angle of 1,2 lines and stratum updip direction;

The C point is respectively the projection to 2 lines and 1 line in D o'clock for the subpoint of D point to ground level, A, B point, and their distance is respectively h ₂, h ₁, the distance that the D point is ordered to O is h ₀, L ₁, L ₂Be respectively the distance that the C point is ordered to B, A, H is the distance that the D point is ordered to C;

B, can derive according to earthquake geometry:

In right-angle triangle OBD, can obtain (angle OBD is the right angle):

h ₁＝h ₀cos ₁ 1-3-1

l ₁＝h ₀sin ₁ 1-3-3

In right-angle triangle OAD, can obtain (angle OAD is the right angle):

h ₂＝h ₀cos ₂ 1-3-2

l ₂＝h ₀sin ₂ 1-3-4

In right-angle triangle BCD, can obtain (angle BCD is the right angle):

L ₁＝h ₁sin′ ₂ 1-3-5

H＝h ₁cos′ ₂ 1-3-7

In right-angle triangle ACD, can obtain (angle ACD is the right angle):

L ₂＝h ₂sin′ ₁ 1-3-6

H＝h ₂cos′ ₁ 1-3-8

Above-mentioned various in, h ₁, h ₂,  ₁,  ₂All can calculate  ' at migrated section ₁Be the ∠ CDA among the right-angle triangle CDA,  ' ₂Be the ∠ CDB among the right-angle triangle CDB.

Can derive in conjunction with Fig. 1 and top formula:

$\sin {{\mathrm{\φ}}^{\′}}_2=\frac{{\sin \mathrm{\φ}}_2}{{\cos \mathrm{\φ}}_1}\frac{{\sin \mathrm{\ω}}_1}{\cos (\mathrm{\ω}-{\mathrm{\ω}}_1)}---1-3-9$

$\sin {{\mathrm{\φ}}^{\′}}_1=\frac{{\sin \mathrm{\φ}}_1}{{\cos \mathrm{\φ}}_2}\frac{{\sin \mathrm{\ω}}_2}{\cos (\mathrm{\ω}-{\mathrm{\ω}}_2)}---1-3-10$

When the vertical survey line 2 of survey line 1, promptly during ω=90 °:

$\sin {{\mathrm{\φ}}^{\′}}_2=\frac{{\sin \mathrm{\φ}}_2}{{\cos \mathrm{\φ}}_1}---1-3-11$

$\sin {{\mathrm{\φ}}^{\′}}_1=\frac{{\sin \mathrm{\φ}}_1}{{\cos \mathrm{\φ}}_2}---1-3-12$

When two main profiles 1 and 2 are any angle, can push away ω ₁, ω ₂Formula:

${\mathrm{tg\ω}}_1=\frac{{\sin \mathrm{\φ}}_2-{\sin \mathrm{\φ}}_1\cos \mathrm{\ω}}{{\sin \mathrm{\φ}}_2\sin \mathrm{\ω}}---1-3-13$

${\mathrm{tg\ω}}_2=\frac{{\sin \mathrm{\φ}}_1-{\sin \mathrm{\φ}}_2\cos \mathrm{\ω}}{\sin {\mathrm{\φ}}_1\sin \mathrm{\ω}}---1-3-14$

ω is the angle of two surveys line in the formula, can directly calculate.

Utilize formula 1-3-9,1-3-10, and the result of calculation of formula 1-3-13,1-3-14 be updated to formula 1-3-5,1-3-6,1-3-7,1-3-8, in can calculate the side-play amount L in sky school ₁, L ₂And the vertical depth H value behind the empty school, promptly finish in empty school; In formula 1-3-7 and 1-3-8 two formulas, the H value of calculating gained is very approaching, can draw structural map with the H value.

This method is finished automatically by computing machine, reduces the personnel's of explaining labour intensity, has improved work efficiency, the mapping precision height.

Embodiment 1:

Fig. 2 is fragrant 148 wellblock T ₂Layer is schemed depth value behind the empty school of empty calibration method with this 2D seisimc overlapping shift profile one-tenth and is projected to scatter diagram on the survey grid, "+" expression CDP point among the figure, " ° " represent breakpoint; As can be seen from the figure, each CDP point has all reached underground real reflection spot position behind the empty school; Fig. 3 is to be according to the automatic degree of depth structural map that generates after on computers main profile and contact side line all being playbacked with this method with Fig. 2.Just finished with 5 minutes time of less than on computers, and method by hand needs several days time just can finish, thereby improved work efficiency greatly, the structural map precision also is improved simultaneously.

Embodiment 2:

Fig. 4 is the tame defensive wall T of Xu of pine the Liao Dynasty ₄Layer t ₀Scatter diagram; Fig. 5 is the tame defensive wall T of Xu ₄The degree of depth scatter diagram that layer obtains behind the empty school of this method, "+" expression CDP point among the figure, " ° " the expression breakpoint.Can find that from Fig. 4 and Fig. 5 each the empty school side of CDP point among Fig. 5 is to all being along the playback of stratum updip direction, by the skew principle, each CDP point of playback back has all reached underground real reflection spot position.

This shows that this 2D seisimc overlapping shift profile becomes the empty calibration method of figure, utilized the data creating structural map of main profile and cross-track simultaneously, not the lost data point; What draw is the true depth structural map, and structural feature is true, and high some position is accurate; The automatic empty school that uses a computer explains that personnel labor intensity is little, efficient is high, mapping precision is high.

Claims (1)

1, a kind of 2D seisimc overlapping shift profile becomes the empty calibration method of figure, it is characterized in that: utilize computing machine to finish automatically by following step:

A, to establish S be ground level, and P is underground zone of interest reflecting interface, and the reflecting interface inclination angle is ;

The I plane was the plane normal of survey line 1 perpendicular to P, handed over plane P in 1 ' line, and 1 ' line apparent dip is  ₁

The II plane was the plane normal of survey line 2 perpendicular to P, handed over plane P in 2 ' line, and apparent dip is  ₂, 1 ' line and 2 ' line meet at the D point, and this point is exactly a CDP point;

1 line and 2 lines are two seismic line that intersect arbitrarily on the ground, meet at the O point, and angle is ω, ω ₁And ω ₂Be respectively the angle of 1,2 lines and stratum updip direction;

The C point is respectively the projection to 2 lines and 1 line in D o'clock for the subpoint of D point to ground level, A, B point, and their distance is respectively h ₂, h ₁, the distance that the D point is ordered to O is h ₀, L ₁, L ₂Be respectively the distance that the C point is ordered to B, A, H is the distance that the D point is ordered to C.

B, can derive according to earthquake geometry:

h ₁＝h ₀cos ₁ 1-3-1

h ₂＝h ₀cos ₂ 1-3-2

l ₁＝h ₀sin ₁ 1-3-3

l ₂＝h ₀sin ₂ 1-3-4

L ₁＝h ₁sin′ ₂ 1-3-5

L ₂＝h ₂sin′ ₁ 1-3-6

H＝h ₁cos′ ₂ 1-3-7

H＝h ₂cos′ ₁ 1-3-8

In formula 1-3-1,1-3-2,1-3-3,1-3-4,1-3-5,1-3-6,1-3-7,1-3-8, h ₁, h ₂,  ₁,  ₂All can calculate  ' at migrated section ₁Be the ∠ CDA among the right-angle triangle CDA,  ' ₂Be the ∠ CDB among the right-angle triangle CDB.

Can derive in conjunction with Fig. 1 and top formula:

$\sin {{\mathrm{\φ}}^{\′}}_2=\frac{{\sin \mathrm{\φ}}_2}{{\cos \mathrm{\φ}}_1}\frac{{\sin \mathrm{\ω}}_1}{\cos (\mathrm{\ω}-{\mathrm{\ω}}_1)}---1-3-9$

$\sin {{\mathrm{\φ}}^{\′}}_1=\frac{{\sin \mathrm{\φ}}_1}{{\cos \mathrm{\φ}}_2}\frac{{\sin \mathrm{\ω}}_2}{\cos (\mathrm{\ω}-{\mathrm{\ω}}_2)}---1-3-10$

When the vertical survey line 2 of survey line 1, i.e. ω=90 °:

$\sin {{\mathrm{\φ}}^{\′}}_2=\frac{{\sin \mathrm{\φ}}_2}{{\cos \mathrm{\φ}}_1}---1-3-11$

$\sin {{\mathrm{\φ}}^{\′}}_1=\frac{{\sin \mathrm{\φ}}_1}{{\cos \mathrm{\φ}}_2}---1-3-12$

When the main profile interconnection is any angle, can push away ω ₁, ω ₂Formula:

${\mathrm{tg\ω}}_1=\frac{{\sin \mathrm{\φ}}_2-{\sin \mathrm{\φ}}_1\cos \mathrm{\ω}}{{\sin \mathrm{\φ}}_1\sin \mathrm{\ω}}---1-3-13$

${\mathrm{tg\ω}}_2=\frac{{\sin \mathrm{\φ}}_1-{\sin \mathrm{\φ}}_2\cos \mathrm{\ω}}{\sin {\mathrm{\φ}}_2\sin \mathrm{\ω}}---1-3-14$

ω is the angle of two surveys line in the formula, can directly calculate.

Utilize formula 1-3-9,1-3-10, and the result of calculation of formula 1-3-13,1-3-14 be updated to formula 1-3-5,1-3-6,1-3-7,1-3-8, in can calculate the side-play amount L in sky school ₁, L ₂And the vertical depth H value behind the empty school, promptly finish in empty school; Can draw structural map with the H value.

Images