Content of the invention
It is an object of the invention to provide a kind of orientation log Numerical Methods, to solve orientation curve 0 (360)
The problem that degree nearby swings and periodically spins;Present invention also offers a kind of orientation log numeric value display method,
To solve the problems, such as that the orientation curve of current display produces orientation rotation less than side.
The present invention provides a kind of orientation log Numerical Methods for solving above-mentioned technical problem, the processing method
Comprise the following steps:
1) judge whether orientation curve data had 0 point of Orientation differences feature;
2) the adjacent depth point data with " crossing 0 point of Orientation differences feature " is carried out numerical transformation with 360 degree as the cycle,
Obtain with periodically variable curve data;
3) according to orientation curvilinear characteristic, by step 2) in the range of the curve data numerical value that obtains standardization [0,360].
Further, the step 1) in cross 0 point of Orientation differences feature refer to:Azimuth number when adjacent two depth point
Value α1、α2Meet | α1-α2|>During 360 °-Δ α, that is, think that the two depth point bearing datas have and " cross 0 point of Orientation differences special
Levy ", Δ α be a given angle numerical value change amount, 0 °<Δα<180°.
Further, described step 2) it is that the larger azimuth of orientation angles in two adjacent depth datas is deducted 360
Degree realizes numerical transformation.
Further, the step 2) also include to obtaining entering row interpolation and filter with periodically variable curve data
Ripple process.
Further, described step 3) be by step 2) in orientation curve data carry out 360 degree of complementation computings and obtain.
Present invention also offers a kind of orientation log numeric value display method, the display methods is comprised the following steps:
1) judge whether orientation curve data had 0 point of Orientation differences feature;
2) " crossing at 0 point " position will be found, calculated with the adjacent depth point Data Extrapolation of " crossing 0 point of Orientation differences feature "
" crossing at 0 point " depth coordinate;
3) the adjacent depth point with " crossing 0 point of Orientation differences feature " is connected with " crossing at 0 point " is obtained, realizes that curve extends
To 0 degree and/or 360 degree of orientation border.
Further, the step 1) in cross 0 point of Orientation differences feature refer to:Azimuth number when adjacent two depth point
Value α1、α2Meet | α1-α2|>During 360 °-Δ α, that is, think that the two depth point bearing datas have and " cross 0 point of Orientation differences special
Levy ", Δ α be a given angle numerical value change amount, 0 °<Δα<180°.
Further, the adjacent depth point Data Extrapolation by with " crossing 0 point of Orientation differences feature " is referred to phase
In adjacent two depth points, the less point in azimuth is extrapolated to next cycle, and point larger for azimuth in adjacent two depth point is extrapolated
Arrive the previous cycle.
Further, the step 2) in " cross 0 point " depth coordinate y0 be:
Y0=y1+x1* (y2-y1)/(x1+360-x2)
Wherein (x1, y1) and (x2, y2) is respectively the coordinate of adjacent two depth point, and y1 and y2 is respectively adjacent two depth point
Depth coordinate, x1 and x2 is respectively the azimuth angle value of adjacent two depth point, and x1<x2.
Further, described " crossing at 0 point " depth coordinate one is located at 0 degree of line, and one is located at 360 degree of lines, by adjacent two
In depth point, the less point in azimuth is connected with the depth coordinate for being located at 0 degree of line, will be larger for azimuth in adjacent two depth point
Point is connected with the depth coordinate for being located at 360 degree of lines.
The invention has the beneficial effects as follows:The present invention first determines whether whether orientation curve data has and " crosses 0 point of Orientation differences special
Levy ";Then the adjacent depth point data with " crossing 0 point of Orientation differences feature " is carried out numerical transformation with 360 degree as the cycle, is obtained
To with periodically variable curve data;Finally according to orientation curvilinear characteristic, by curve numerical value standardization [0,360] scope
Interior.Adjacent depth point data with " crossing 0 point of Orientation differences feature " is carried out numerical transformation with 360 degree as the cycle by the present invention,
Solve orientation curve to swing near 0 (360) degree and numerical value process when periodically spinning and application abnormal problem, overcome
Existing disposal methods result problem of dtmf distortion DTMF, makes process and orientation logging data application result more reasonable, accurately.
The present invention will also find " crossing at 0 point " position with the adjacent depth point Data Extrapolation of " crossing 0 point of Orientation differences feature "
Put;Adjacent depth point with " cross 0 point of Orientation differences feature " be connected with " crossing at 0 point " is obtained, realize curve extend to 0 degree with
Or 360 degree of orientation border.Overcome the problem on " orientation rotation is less than side ", it is ensured that the continuity of orientation curvilinear motion and complete
Property, reduce the truth of instrument swing.
Specific embodiment
Below in conjunction with the accompanying drawings the specific embodiment of the present invention is described further.
A kind of embodiment of orientation log Numerical Methods of the present invention
The orientation log Numerical Methods of the present invention are directed to orientation curve, and in 0 (360) degree, nearby left and right is put
Problem that is dynamic and periodically spinning, the flow process of the method is as shown in figure 1, comprise the following steps that.
1. pair pending orientation curve data carries out " crossing 0 point of Orientation differences feature " judgement.
The definition of " cross 0 point Orientation differences feature " is given first, if α1、α2For azimuthal variations (α1∈[0,Δα)、α2∈
(360- Δ α, 360]), Δ α is a less angular values variable quantity (0<Δα<180) (such as 30), when adjacent two depth point
Azimuth value is from α1It is mutated α2, or from α2It is mutated α1When, the two data points are defined as with " 0 point excessively Orientation differences are special
Levy ", the intersection point O of two depth point azimuthal coordinates lines and 0 (360) rhumb line is " crossing at 0 point ", the two adjacent depth of A, B in such as Fig. 2
Degree point and two adjacent depth point of C, D are respectively provided with " cross 0 point of Orientation differences feature ", and O points are their " crossing at 0 points ", wherein d1And d2
It is the depth coordinate of adjacent two depth point, the two numerical value differs a sampling interval.
The basis for estimation of " crossing 0 point of Orientation differences feature " is as follows:(0 ° of given less angular values variation delta α<
Δα<180 °) (such as 30), when azimuth value α of adjacent two depth point1、α2Meet | α1-α2|>During 360 °-Δ α, that is, think this
Two depth point bearing datas have " crossing 0 point of Orientation differences feature ".
2. pair adjacent depth point data with " crossing 0 point of Orientation differences feature " carries out numerical transformation with 360 degree as the cycle,
Obtain with periodically variable curve data.
Well logging orientation curve possesses following characteristics:
(1) α ∈ [0 °, 360 °]:Azimuth changes between 0 to 360 degree;
(2) 0=360 °:0 degree and 360 degree is same position of orientation;
(3) α=mod (α, 360 °), α<0 ° or α>360°:Azimuth has periodically, and mod is complementation operation function;
(4)-α=360 °-α, α=α -360 °.
When adjacent depth point data α1、α2When there is " crossing 0 point of Orientation differences feature ", numerical value can be entered according to features described above
Row such as down conversion:
If α1>α2, α1=α1-360°;
Otherwise α2=α2-360°.
Can be by the adjacent depth point data with " crossing 0 point of Orientation differences feature " with 360 degree as cycle by said process
Numerical transformation is carried out, is obtained with periodically variable curve data.
Above-mentioned ordinary curve processing mode carried out with periodically variable curve data data are inserted to obtained
The numerical value such as value, filtering are processed, and ordinary curve interpolation, filter processing method are a lot, such as linear interpolation, spline interpolation, Newton interpolation,
The rather filtering etc. of gaussian filtering, the Chinese, the present invention are no longer described in detail.If adjacent depth point data does not have " crosses 0 point of Orientation differences special
Levy ", then data keep constant, directly enter row interpolation, filtering process.
3., according to orientation curvilinear characteristic, the result data numerical value standardization that step 2 is processed is arrived in the range of [0 °, 360 °],
Application for orientation curve provides performance data.
According to orientation curvilinear characteristic (3):Val=mod (Val, 360 °), by result numerical value specification to [0 °, 360 °]
In the range of.
By said process, the present invention can overcome the defect that existing disposal methods result problem of dtmf distortion DTMF, make process and
Orientation logging data application result is more reasonable, accurate.
A kind of embodiment of orientation log numeric value display method of the present invention
1. detect between data point to be drawn with the presence or absence of " crossing 0 point of Orientation differences feature ".
The detection of " crossing 0 point of Orientation differences feature " is according to as follows:(0 ° of given less angular values variation delta α<
Δα<180 °) (such as 30 °), when azimuth value α of adjacent two depth point1、α2Meet | α1-α2|>During 360 °-Δ α, that is, think
The two depth point bearing datas have " crossing 0 point of Orientation differences feature ".
2. pair there is the adjacent depth point Data Extrapolation of " crossing 0 point of Orientation differences feature ", determine " crossing at 0 point " position.
Below to illustrate as a example by adjacent two depth points P1 (x1, y1), P2 (x2, y2), as shown in figure 3, P1 and P2 " mistakes
0 point " determination process of position is as follows:
A. the less point in azimuth in adjacent two depth point is extrapolated to next cycle, by orientation in adjacent two depth point
The larger point in angle is extrapolated to the previous cycle, with the corresponding extrapolation point for obtaining.
Azimuth of the azimuth of P1 less than P2 in adjacent two depth points P1 (x1, y1) and P2 (x2, y2) in the present embodiment,
P1 and P2 is extrapolated to next cycle and previous cycle respectively, point A1 (x1+360 °, y1), A2 (x2-360 °, y2) is obtained.
B. according to determine extrapolation point obtain adjacent two depth point respectively with 0 degree of line and the intersection point of 360 degree of lines.
Tie point A2P1, A1P2, meet at point B1 (0, y01), B2 (360 °, y02) with 0 °, 360 ° of lines respectively;
Because:Line segment A2D1=P2C2, C1P1=D2A1, C1A1//A2C2
△ B1A2D1 ≌ B2P2C2, △ C1B1P1 ≌ D2B2A1 can be released
So y01=y02, is calculated as y0, (360 °, y02) of point (0 °, y01) and point are same position of orientation, as (x1,
Y1), " crossing at the 0 point " position of (x2, y2).
C. " crossing at 0 point " depth coordinate is calculated.
Y0 computing formula can be released according to Similar Principle of Triangle:
Y0=y1+x1* (y2-y1)/(x1+360 ° of-x2)
When x1=0 °, y0=y1;When x2=360 °, y0=y2.
3. Drawing of Curve
" crossing at the 0 point " depth coordinate for obtaining is attached with adjacent two depth point, " crossing at 0 point " depth coordinate one is located at
0 degree of line, one is located at 360 degree of lines, by the less point in azimuth in adjacent two depth point and the depth coordinate phase for being located at 0 degree of line
Even, point larger for azimuth in adjacent two depth point is connected with the depth coordinate for being located at 360 degree of lines.
The zero crossing depth coordinate obtained in the present embodiment is respectively B1 and B2, B1 are located at 0 degree of line, and B2 is located at 360 degree of lines,
Azimuth of the azimuth of P1 less than P2, therefore, draws line P1B1, B2P2, curve is extended to 0 ° of (360 °) orientation border,
Ensure the continuity and integrality of orientation curvilinear motion, as shown in Figure 4.
By said process, the display methods of the present invention is calculated and curve extrapolation by " crossing at 0 point " depth coordinate, gram
Take the problem on " orientation rotation is less than side ", it is ensured that the continuity and integrality of orientation curvilinear motion, reduce instrument swing
Truth.
The effect of the present invention is described below by specific application example.
Example 1:Orientation curve numerical filtering processing method checking
In order to verify the correctness of orientation curve filtering method, the obvious example of special use bearing data numerical value change.
Accompanying drawing 9 is the pole plate orientation curve that EMI electric imaging loggings are measured, as downhole probe spins than more serious,
When routinely curve mode is to curve filtering process, cross about 0 point of orientation curve numerical value serious distortion, by the present invention in side
Method filter result is closer to true bearing plots changes, and overcomes the problem on " orientation rotation is less than side ".
Example 2:Orientation curve interpolation method application.
When EMI electric imaging loggings image is processed, need with a pole plate orientation curve to still image and dynamic image
Data carry out directional correction, and in the test data file, (10 points/rice) of orientation curve sample rate is adopted less than pole plate conductance profile
Sample rate (400 points/rice), needs to carry out resampling interpolation processing to orientation curve.Accompanying drawing 10 is that the EMI Electrical imagings to being tested are surveyed
A number pole plate orientation curve resampling Comparative result of well record, it is evident that orientation resampling mode result is closer in fact
Border orientation plots changes, and eliminate quickly spinned due to downhole probe, relatively low " less than side " effect for causing of sample rate
Should.
Example 3:Orientation curve interpolation result application
Application of the orientation curve interpolation resampling result in image procossing.
In accompanying drawing 11, left and right sides image road is respectively adopted ordinary curve interpolation method and orientation curve interpolation resampling side
Formula, carries out interpolation processing to a pole plate orientation curve, carries out orientation with the orientation curve after resampling interpolation to imaging data
Correction, it is evident that do directional correction with orientation curve interpolation method for resampling result, the strip-form eliminated in image are wrong
Position phenomenon.
Accompanying drawing 12 is orientation curve method for resampling application achievements and similar software comparison diagram both at home and abroad, and first is No. one
Pole plate orientation curve, second are the orientation curve resamplings in the sinopec Well Data Processing software Logik application present invention
The dynamic image of method process, third and fourth, five roads be that Schlumberger GeoFrame softwares and domestic other software are processed respectively
Dynamic image, by contrast it can be seen that Logik directional corrections result is consistent with GeoFrame results, and other are soft for the country
There is obvious stepped ramp type and herringbone wave mode image offset phenomenon when orientation is spinned in the image of part process.
Example 4:Orientation curves drawing display methods application
Accompanying drawing 13 is by Drawing of Curve mode in orientation of the present invention and existing Drawing of Curve mode, to pole plate orientation curves drawing
Show Comparative result, first is depth track in figure, and second curve is to show result by existing Drawing of Curve mode, the 3rd, 4
Road curve is to show that result, wherein the 3rd road curve are drawn mode by conventional curvature and drawn by Drawing of Curve mode in orientation of the present invention
0 point of closed line is crossed, the 4th road curve does not draw closed line.
Accompanying drawing 14 is the special case that orientation curve swings near 0 (360) degree, shows result pair by different drafting modes
Than, it can be seen that more can intuitively reflect the change of orientation curve by orientation Drawing of Curve mode, the mapping result for not drawing closed line
Trend, and the dotted box portion of the 4th, 5 road joints is only the true reflection of curvilinear motion.