Disclosure of Invention
The invention provides a calculation method of an original dip angle of a deposition structure, which is used for solving the problem that the calculation method of the original dip angle of the deposition structure is not available in the prior art; the invention also provides a method for calculating the original tendency of the deposition structure, which is used for solving the problem that the method for calculating the original tendency of the deposition structure does not exist in the prior art; correspondingly, the invention also provides a calculation system of the original dip angle of the deposition structure, which is used for solving the problem that no calculation method of the original dip angle of the deposition structure exists in the prior art; the invention also provides a system for calculating the original tendency of the deposition structure, which is used for solving the problem that no method for calculating the original tendency of the deposition structure exists in the prior art.
In order to achieve the above purpose, the technical scheme provided by the invention is as follows:
a method for calculating an original dip angle of a deposition structure, comprising the steps of:
(1) Acquiring a normal line of a reference surface and a normal line of a deposition structural surface;
(2) And obtaining an intersection line of a plane where the normal of the reference surface and the normal of the deposition structural surface are positioned and a horizontal plane, and calculating an included angle between the intersection line and the normal of the reference surface, wherein the included angle is the original inclination angle of the deposition structural surface.
Further, let the original dip of the deposited structure be x=arccoso (|k) 1 I) and a value range of [0, pi/2 ]]Then:
K 1 =sin(r 1 )sin(r)[1/sin 2 (r 1 )+1/sin 2 (r)-cot 2 (r 1 )-cot 2 (r)+2cos(y 1 -y)cot(r 1 )cot(r)]/2
wherein r is 1 Is the normal L of the reference plane 1 Inclination angle, y 1 The inclination of the normal line of the reference surface is r, the inclination angle of the normal line of the original structural surface is y, and the inclination of the normal line of the original structural surface is y.
A method of calculating an original propensity of a deposit formation, comprising the steps of:
(1) Acquiring an included angle between a normal line of a reference surface and a normal line of a deposition structural surface;
(2) Obtaining the original tendency of the deposition structure according to the included angle between the normal of the reference surface and the normal of the deposition structure surface:
when the included angle between the normal line of the reference surface and the normal line of the deposition structure surface is 0 or pi/2, judging that the original trend of the deposition structure is no trend;
when the included angle between the normal line of the reference surface and the normal line of the deposition structural surface is between 0 and pi/2 and r 1 At r.ltoreq., the original tendency of the deposit structure is judged as y 1 +m;
When the angle between the normal line of the reference surface and the normal line of the deposition structural surface is 0 toBetween pi/2 and r 1 >When the intersection point of the normal line of the deposition structural surface and the horizontal plane is judged to be positioned on the minor arc, the original trend of the deposition structural surface is y 1 -m+π;
When the included angle between the normal line of the reference surface and the normal line of the deposition structural surface is between 0 and pi/2 and r 1 >When the intersection point of the normal line of the deposition structural surface and the horizontal plane is judged to be positioned on the major arc, the original trend of the deposition structural surface is y 1 +m;
When the included angle between the normal line of the reference surface and the normal line of the deposition structure surface is greater than pi/2, the original tendency of the deposition structure is judged to be y 1 +m+π;
m=arcsin(K 2 )
K 2 =sin(y 1 -y)·cos(r)/sin(P)
Set S 1 =tan(r)/tan(r 1 ),Then
When C 1 ≥y 1 -y or 2 pi-C 1 ≤y 1 When y, judging that the intersection point of the normal line of the deposition structural surface and the horizontal plane is positioned on the minor arc;
when C 1 <y 1 -y and 2 pi-C 1 >y 1 -y, determining that the intersection point of the normal line of the deposition structural surface and the horizontal plane is located on the major arc;
wherein r is 1 Is the normal L of the reference plane 1 Inclination angle, y 1 The inclination of the normal line of the reference surface is represented by r, the inclination angle of the normal line of the original construction surface is represented by y, the inclination of the normal line of the original construction surface is represented by P, and the included angle between the normal line of the reference surface and the normal line of the deposition construction surface is represented by P.
A computing system for deposit build primitive tilt comprising a processor and a memory having stored thereon a computer program for execution on the processor; the processor, when executing the computer program, implements the steps of:
(1) Acquiring a normal line of a reference surface and a normal line of a deposition structural surface;
(2) And obtaining an intersection line of a plane where the normal of the reference surface and the normal of the deposition structural surface are positioned and a horizontal plane, and calculating an included angle between the intersection line and the normal of the reference surface, wherein the included angle is the original inclination angle of the deposition structural surface.
Further, let the original dip of the deposited structure be x=arccoso (|k) 1 I) and a value range of [0, pi/2 ]]Then:
K 1 =sin(r 1 )sin(r)[1/sin 2 (r 1 )+1/sin 2 (r)-cot 2 (r 1 )-cot 2 (r)+2cos(y 1 -y)cot(r 1 )cot(r)]/2
wherein r is 1 Is the normal L of the reference plane 1 Inclination angle, y 1 The inclination of the normal line of the reference surface is r, the inclination angle of the normal line of the original structural surface is y, and the inclination of the normal line of the original structural surface is y.
A computing system for deposit build primitive trends comprising a processor and a memory having stored thereon a computer program for execution on the processor; the processor, when executing the computer program, implements the steps of:
(1) Acquiring an included angle between a normal line of a reference surface and a normal line of a deposition structural surface;
(2) Obtaining the original tendency of the deposition structure according to the included angle between the normal of the reference surface and the normal of the deposition structure surface:
when the included angle between the normal line of the reference surface and the normal line of the deposition structure surface is 0 or pi/2, judging that the original trend of the deposition structure is no trend;
when the included angle between the normal line of the reference surface and the normal line of the deposition structural surface is between 0 and pi/2 and r 1 At r.ltoreq., the original tendency of the deposit structure is judged as y 1 +m;
When the included angle between the normal line of the reference surface and the normal line of the deposition structural surface is between 0 and pi/2 and r 1 >When the intersection point of the normal line of the deposition structural surface and the horizontal plane is judged to be positioned on the minor arc, the original trend of the deposition structural surface is y 1 -m+π;
When the included angle between the normal line of the reference surface and the normal line of the deposition structural surface is between 0 and pi/2 and r 1 >When the intersection point of the normal line of the deposition structural surface and the horizontal plane is judged to be positioned on the major arc, the original trend of the deposition structural surface is y 1 +m;
When the reference isWhen the included angle between the surface normal and the normal of the deposition structure surface is larger than pi/2, the original tendency of the deposition structure is judged to be y 1 +m+π;
m=arcsin(K 2 )
K 2 =sin(y 1 -y)·cos(r)/sin(P)
Set S 1 =tan(r)/tan(r 1 ),Then
When C 1 ≥y 1 -y or 2 pi-C 1 ≤y 1 When y, judging that the intersection point of the normal line of the deposition structural surface and the horizontal plane is positioned on the minor arc;
when C 1 <y 1 -y and 2 pi-C 1 >y 1 -y, determining that the intersection point of the normal line of the deposition structural surface and the horizontal plane is located on the major arc;
wherein r is 1 Is the normal L of the reference plane 1 Inclination angle, y 1 The inclination of the normal line of the reference surface is represented by r, the inclination angle of the normal line of the original construction surface is represented by y, the inclination of the normal line of the original construction surface is represented by P, and the included angle between the normal line of the reference surface and the normal line of the deposition construction surface is represented by P.
According to the technical scheme provided by the invention, the original inclination angle and the original tendency of the deposition structure are calculated by utilizing the relation among the reference surface, the horizontal surface and the deposition structure surface, so that a basis is provided for predicting the development trend of the sand body and finding out an advantageous oil and gas reservoir.
Detailed Description
The invention provides a calculation method of an original dip angle of a deposition structure, which is used for solving the problem that the calculation method of the original dip angle of the deposition structure is not available in the prior art; the invention also provides a method for calculating the original tendency of the deposition structure, which is used for solving the problem that the method for calculating the original tendency of the deposition structure does not exist in the prior art; correspondingly, the invention also provides a calculation system of the original dip angle of the deposition structure, which is used for solving the problem that no calculation method of the original dip angle of the deposition structure exists in the prior art; the invention also provides a system for calculating the original tendency of the deposition structure, which is used for solving the problem that no method for calculating the original tendency of the deposition structure exists in the prior art.
In order to achieve the above purpose, the technical scheme provided by the invention is as follows:
a method for calculating an original dip angle of a deposition structure, comprising the steps of:
(1) Acquiring a normal line of a reference surface and a normal line of a deposition structural surface;
(2) And obtaining an intersection line of a plane where the normal of the reference surface and the normal of the deposition structural surface are positioned and a horizontal plane, and calculating an included angle between the intersection line and the normal of the reference surface, wherein the included angle is the original inclination angle of the deposition structural surface.
Embodiments of the present invention will be further described with reference to the accompanying drawings.
Method embodiment:
the embodiment provides a calculation method for recovering the original occurrence of a sedimentary structure, which comprises a calculation method for calculating the original inclination angle of the sedimentary structure and a calculation method for calculating the original trend of the sedimentary structure, and is used for providing a basis for judging the ancient water flow direction of a sedimentary basin and predicting the development trend of a sand body so as to find an advantageous oil and gas reservoir.
In this embodiment, the method for calculating the original inclination angle of the deposition structure is as follows:
acquiring a normal line of a reference surface and a normal line of a deposition structural surface;
and obtaining an intersection line of a plane where the normal of the reference surface and the normal of the deposition structure surface are located and a horizontal plane, and calculating an included angle between the intersection line and the normal of the reference surface, wherein the included angle is the original inclination angle of the deposition structure.
As shown in fig. 1, the horizontal plane is α, and the normal line of the reference plane is L 1 A normal line of a deposition structural surface is L 2 The method comprises the steps of carrying out a first treatment on the surface of the Normal L of reference plane 1 And a deposition structural surface normal L 2 Intersecting at the point A, making a vertical line AD of the horizontal plane alpha from the point A, and intersecting the vertical line AD with the horizontal plane alpha at a point D; normal L of reference plane 1 Intersecting the horizontal plane with point B, depositing the normal L of the structural surface 2 Intersecting with the horizontal plane alpha at a point C, the straight line AB is the normal L of the reference plane 1 And a deposition structural surface normal L 2 The included angle BAC is the original inclination angle of the inclined deposition structure surface of the deposition structure.
The calculation method of the original inclination angle of the deposition structural surface comprises the following steps:
let L be 1 Is at an inclination angle r 1 Tendency is y 1 ;L 2 The inclination angle of (a) is r, the tendency is y, the length of a line segment AD is 1, and the included angle BAC is solved, wherein the steps are as follows:
BD=AD·cot(r 1 )=cot(r 1 )
CD=AD·cot(r)=cot(r)
∠BDC=y 1 -y
BC 2 =BD 2 +DC 2 -2BD·DC·cos∠BDC
=cot 2 (r 1 )+cot 2 (r)-2cot(r 1 )·cot(r)·cos(y 1 -y)
AB=AD/sin(r 1 )=1/sin(r 1 )
AC=AD/sin(r)=1/sin(r)
cos∠BAC=(AB 2 +AC 2 -BC 2 )/(2AB·AC)
=sin(r 1 )sin(r)[1/sin 2 (r 1 )+1/sin 2 (r)-cot 2 (r 1 )-cot 2 (r)+2cos(y 1 -y)cot(r 1 )cot(r)]/2。
let K be 1 =cos +.BAC, then since the threshold of tilt angle is [0, pi/2]The original dip of the deposited structure is x=arccoss (K 1 )。
In calculating the original propensity for deposit formation, two conditions need to be assumed: firstly, the formation is changed only in the vertical direction by the construction movement, and is not rotated in the horizontal direction; secondly, the tilt of the deposition build is reset from the minimum angular orientation to horizontal, i.e. it is assumed that the build tilt angle does not exceed 90 °.
In this embodiment, the method for calculating the original tendency of the deposition structure is as follows:
after the deposition structure is reset, the normal L of the reference surface 1 In a vertical state, the reset of the deposition structure is carried out from the direction of the minimum angle, and the complementary angle of the normal line of the deposition structure and the horizontal plane alpha is the minimum reset angle, so the normal line L of the reference surface 1 Is carried out in a plane defined by the angle ABD. Reference surface normal L 1 Is a directional line of (2)No rotation takes place during the resetting, which is directed to the same before and after the resetting of the deposited structure, so in this embodiment +.>As a reference line.
Thus, the reference surface normal L is to be determined 1 The inclination after reset requires first determining a line of inclination, i.e. a vertical projection line (or a straight line parallel thereto) on the horizontal plane alpha after reset, and then determining the line of inclination and the line of inclinationAn included angle between the two.
As shown in fig. 1, a perpendicular CE to the AB is drawn through the point C, such that the perpendicular CE is perpendicular to the AB;
when the deposition structure is reset, the line AB is perpendicular to the horizontal plane α, and the CE is parallel to the horizontal plane α, i.e. vector, since the CE is perpendicular to the line ABFor depositing the structural surface normal L 2 Is a directional line of (a).
The passing point C is a straight line CF perpendicular to the straight line BD at F, and is connected with the point E and the point F; CF is in horizontal plane α, CF is perpendicular to AD, and CF is perpendicular to BD, so CF is perpendicular to plane ABD; the straight line AB is in the plane ABD, and CF can be obtained to be perpendicular to the straight line AB; the line AB is perpendicular to the line CE, and thus perpendicular to the plane EFC, and perpendicular to the line EF.
When the deposition structure is reset, the normal L of the reference surface 1 Perpendicular to the horizontal plane alpha, the plane ABD intersects with the horizontal plane alpha at BD, so the +.FEC is the pointing line andan included angle between the two.
FC=CD·sin(y 1 -y)=cot(r)·sin(y 1 -y)
CE=AC·sin∠BAC=(1/sin(r))·sin(P)
sin∠FEC=CF/CE
=cot(r)·sin(y 1 -y)·sin(r)/sin(P)
=sin(y 1 -y)·cos(r)/sin(P)
Let K 2 =sin∠FEC
=cot(r)·sin(y 1 -y)·sin(r)/sin(P)
=sin(y 1 -y)·cos(r)/sin(P)
∠FEC=arcsin(K 2 )
Let m= fec=arcsin (K 2 )
Where P is the angle between the normal of the reference surface and the normal of the deposition build surface, and p=arccose (K 1 ) According to the difference of the included angle P between the normal of the reference surface and the normal of the deposition structural surface, the original tendency is classified and calculated under different conditions:
when p=0, the reference surface normal L 1 Normal L to the deposition build surface 2 Coincidence, so that the deposition structure is horizontal after reset, and no tendency exists;
when p=pi/2, referenceSurface normal L 1 Normal L to the deposition build surface 2 Perpendicular to A because of the normal L to the reference surface after reset of the deposited structure 1 Perpendicular to the horizontal plane alpha, thus depositing the normal L of the construction surface after reset 2 Parallel to the horizontal plane alpha, i.e. the deposition structure is vertical to the horizontal plane after reset;
when P takes on a value between 0 and pi/2, if r 1 R.ltoreq.r, as shown in FIG. 2, is provided with a circle O 1 Is the normal L of the reference plane 1 Trajectory of possible intersection with horizontal plane α, circle O 2 Build plane L for deposition 2 A trajectory of possible intersections with the horizontal plane α.
From r 1 Less than or equal to r can obtain cotr 1 Not less than cotr, i.e. circle O 1 Radius of equal to or greater than circle O 2 Radius of (2);
as can be seen from the figure, the perpendicular CF passing through C and BD is always the same as the perpendicular CF regardless of the change of the included angle PIntersecting; and because r 1 Not greater than pi/2, so that the reference surface normal L is made by F 1 Is always intersected at a straight line +>On, the beginning edge of the +.FEC +.>Always pointing to the normal L of the reference plane 1 The direction of opening of the final edge is changed from sin (y 1 -y) when sin (y 1 -y) taking the positive, azimuthal increasing direction, when sin (y) 1 -y) when negative, the azimuth increasing direction is reversed, when sin (y) 1 -y) =0, +fec=0.
Normal L of deposition structural surface after reset 2 The azimuth size is: y is 1 +pi+m, thus the original trend of the deposit structure is y 1 +m。
When r is 1 When r, as shown in FIG. 3, circle O 1 Radius smaller than circle O 2 Radius of (2); at this time, the normal L of the deposition structural surface is firstly judged 2 And the horizontal planeThe intersection point of alpha is in a minor arc or a major arc, and the judging method comprises the following steps:
set S 1 =tan(r)/tan(r 1 ),
When C 1 ≥y 1 -y or 2 pi-C 1 ≤y 1 At y, the normal L of the formation surface is deposited 2 The intersection point with the horizontal plane alpha is in a minor arcAnd (3) upper part.
When C 1 <y 1 -y and 2 pi-C 1 >y 1 At y, the normal L of the formation surface is deposited 2 The intersection point with the horizontal plane alpha is in a major arcAnd (3) upper part.
When depositing the normal L of the structural surface 2 The intersection point with the horizontal plane alpha is in a minor arcAs shown in FIG. 3, the point C is CF and +.>Crossing the inversion line of F, and crossing F as FE is perpendicular to +.>Is at E, ">I.e. normal L of the deposition structural surface 2 Is a directional line of (a). It can be seen that the beginning of the +.FEC at this time +.>Along the normal L of the reference plane 1 Orientation of inclination, so that the normal L of the deposition structural surface after the reset 2 Tendency is y 1 M, i.e. of deposited configurationThe original trend is y 1 -m+π。
When depositing the normal L of the structural surface 2 Intersection with horizontal plane alpha is located in major arc(including G, G) 1 Point) and r as shown in fig. 4 1 The same is true for r or less, the normal L of the deposition structural surface 2 Tendency after reset is y 1 +pi+m, the original trend of the deposit structure is y 1 +m。
When P > pi/2, as shown in FIG. 5, after resetWill intersect the horizontal plane alpha, i.e. the normal L to the deposition build surface 2 Inversion occurs. C is adopted as CF to cross BD at F; f is used as FE to be crossed with BA to E; a connection C, E; because of the normal L of the deposition build surface 2 Reverse rotation occurs, so->The reverse direction of (1) is the normal L of the deposition structural surface after reset 2 Is a directional line of (2); />Reverse->Namely the direction of the starting edge of the angle m; thereby get->And L is equal to 1 Is consistent in the direction of inclination; so that the normal L of the structural surface is deposited 2 Tendency after reset is y 1 +m, the original tendency of the deposit structure is y 1 +m+π。
If the result interval of the original trend Y is not within the [0,2 pi ] interval, performing a cyclic operation of adding or subtracting 2 pi until the result is within the [0,2 pi ].
In summary, the original propensity for deposition build up is:
when P is 0 or pi/2, the original propensity of the deposited structure is no propensity;
when P is between 0 and pi/2, and r 1 At r.ltoreq.the original tendency of the deposit structure is y 1 +m;
When P is between 0 and pi/2, and r 1 >r, and deposit the normal L of the structural surface 2 Intersection point with horizontal plane alpha is located in minor arcOn the upper side, the original tendency of the deposition structural surface is y 1-m+pi;
when P is between 0 and pi/2, and r 1 >r, and deposit the normal L of the structural surface 2 Intersection with horizontal plane alpha is located in major arcOn the upper surface, the original tendency of the deposition structural surface is y 1 +m;
When P is greater than pi/2, the original trend of the deposited structure is y1+m+pi.
System example 1:
the embodiment provides a computing system for a deposition build primitive tilt angle, comprising a processor and a memory, the memory having stored thereon a computer program for execution on the processor; the processor, when executing the computer program, implements the method for calculating the original dip angle of the deposition build in the above method embodiment.
System example 2:
the present embodiment provides a computing system for a raw dip of a deposition build, comprising a processor and a memory having stored thereon a computer program for execution on the processor; the processor, when executing the computer program, implements the method of calculating the original propensity of the deposit formation in the above-described method embodiments.