CN113312695A - Ancient landform restoration method based on trend surface analysis - Google Patents
Ancient landform restoration method based on trend surface analysis Download PDFInfo
- Publication number
- CN113312695A CN113312695A CN202110575415.1A CN202110575415A CN113312695A CN 113312695 A CN113312695 A CN 113312695A CN 202110575415 A CN202110575415 A CN 202110575415A CN 113312695 A CN113312695 A CN 113312695A
- Authority
- CN
- China
- Prior art keywords
- trend surface
- trend
- equation
- thickness
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 53
- 238000005211 surface analysis Methods 0.000 title claims abstract description 14
- 238000007689 inspection Methods 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims description 21
- 238000011084 recovery Methods 0.000 claims description 7
- 230000000052 comparative effect Effects 0.000 claims description 3
- 238000001514 detection method Methods 0.000 claims description 2
- 230000015572 biosynthetic process Effects 0.000 abstract description 9
- 238000004364 calculation method Methods 0.000 abstract description 3
- 238000004458 analytical method Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 230000015556 catabolic process Effects 0.000 description 4
- 238000006731 degradation reaction Methods 0.000 description 4
- 230000001419 dependent effect Effects 0.000 description 4
- 230000000694 effects Effects 0.000 description 4
- 230000008021 deposition Effects 0.000 description 3
- 239000006185 dispersion Substances 0.000 description 3
- 238000005056 compaction Methods 0.000 description 2
- 238000000611 regression analysis Methods 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000012937 correction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005429 filling process Methods 0.000 description 1
- 229910052500 inorganic mineral Inorganic materials 0.000 description 1
- 238000012067 mathematical method Methods 0.000 description 1
- 239000011707 mineral Substances 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000004576 sand Substances 0.000 description 1
- 238000007493 shaping process Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/10—Geometric CAD
- G06F30/13—Architectural design, e.g. computer-aided architectural design [CAAD] related to design of buildings, bridges, landscapes, production plants or roads
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/02—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness
- G01B11/06—Measuring arrangements characterised by the use of optical techniques for measuring length, width or thickness for measuring thickness ; e.g. of sheet material
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01B—MEASURING LENGTH, THICKNESS OR SIMILAR LINEAR DIMENSIONS; MEASURING ANGLES; MEASURING AREAS; MEASURING IRREGULARITIES OF SURFACES OR CONTOURS
- G01B11/00—Measuring arrangements characterised by the use of optical techniques
- G01B11/24—Measuring arrangements characterised by the use of optical techniques for measuring contours or curvatures
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
- G06F2119/02—Reliability analysis or reliability optimisation; Failure analysis, e.g. worst case scenario performance, failure mode and effects analysis [FMEA]
Landscapes
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Geometry (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- General Engineering & Computer Science (AREA)
- Structural Engineering (AREA)
- Computational Mathematics (AREA)
- Civil Engineering (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Architecture (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses an ancient landform restoration method based on trend surface analysis, which collects enough single-well target stratum thickness data h for stratum ancient landforms needing to be restored1Single well target overburden thickness data h2And single well geographic coordinates x, y; and then using the calculated superposed stratum thickness data z-h1+h2Constructing a trend surface equation by using a least square method principle, programming a calculation program by using MATLAB software to obtain different-order trend surface model equations, and substituting geographic coordinates x and y into the equations to obtain a series of trend surface valuesThereafter using R2Carrying out moderate inspection on the model by three different methods of inspection, F inspection and moderate successive inspection of the trend surface to obtain the trend surface of the optimal model, and finally carrying out the obtained trend surface and the target formation thickness h1And overlapping to obtain the final restored ancient landform. Compared with the traditional residual thickness method, the method combines the residual thickness method and the trend surface method, the technical means are properly fused, and the reliability of the comprehensive calculation result is high.
Description
Technical Field
The invention relates to the field of ancient landform restoration in mineral exploration technology, in particular to a method for restoring weathering crust type ancient landform.
Background
Ancient landforms are units of landforms whose basic form is completely decoupled from the modern landform shaping process and which are not compatible with modern tectonic conditions. The ancient landscape is often not easily left intact, especially in denuded zones where elevations are constructed. The ancient landform can be exposed on the earth surface in a lifting area, such as a quasi-plain of a mountain top, an ancient valley of a wind gap and the like; in subsidence areas, they are buried under the surface of the existing landform, such as ancient river channels, ancient buried mountains, ancient sand dunes, etc. The ancient landform is an important part of ancient geographic environment research and has important significance for finding oil gas.
The general idea of recovering the ancient landform is to find out the ancient structural surface of the underlying stratum by various methods and then accumulate the ancient structural surface and the thickness of the current stratum to obtain the ancient landform. The method mainly used for recovering the ancient landform at present is a residual thickness method, and the principle of the method is that stratum deposition is a filling and filling process, the deposition thickness at the high places of the terrain is small, the deposition thickness at the low places of the terrain is large, therefore, the top interface of the stratum is leveled, and then the ancient structural surface of the underlying stratum can be obtained according to the thickness of the stratum, as shown in figure 1.
However, this method has its limitations, and works best when the formation is well preserved, and once the formation has degraded, it can be less accurate. In particular, when ancient landforms are restored in a weathering crust type of ground layer, the residual thickness method is not applicable at all. The reason is that the weathering crust is strongly degraded, and although the degradation intensity is influenced by the terrain, the influence factors of the degradation intensity on the top boundary of the stratum are various, so that the weathering crust often shows non-uniformity at the same height. The unconsolidated formation thus recovered has more localized undulations, which is not consistent with the actual situation. In fact, the paleo-tectonic surface should be a relatively smooth surface reflecting the spatially varying tendency of the formation, which may fluctuate over a large range, but is usually locally stable, without significant fluctuations. As shown in fig. 2, for the weathering crust type formation, the ancient structural surface recovered from the direct flat top interface is greatly different from the actual structural surface, and therefore, the ancient geomorphology recovery of the weathering crust formation needs to be performed by improving the original method or using a new method.
Disclosure of Invention
In order to solve the problems in the prior art, the invention aims to provide an improved method for recovering ancient landforms based on a residual thickness method.
The specific technical scheme is as follows:
an ancient landform restoration method based on trend surface analysis comprises the following steps:
step 1, collecting thickness h of weathering crust of target stratum1Thickness h of overburden2And geographic coordinates x and y, calculating the thickness z of the superposed stratum according to the obtained data, constructing a trend surface equation by using a least square method principle, and programming a calculation program by using MATLAB software to obtain a primary simple trend surface equation z11Once bilinear trend surface z12Equation z of the quadratic trend surface model2;
Step 2, substituting the independent variables x and y into the three trend surface model formulas to obtain a series of trend surface values
Step 3, the result is processedLine R2Checking; will zi、Andsubstituting into correlation formula to obtain SSD、SSRAnd SSTThen, R is found2Comparison ofThe higher the numerical value is, the higher the fitting degree is;
step 4, carrying out F inspection on the result; substituting the correlation data into a formula to obtain an F value of a corresponding model, and looking up a table to obtain the F of the single linear trend surface when the confidence level alpha is 0.050.05(p11,n-p11-1), a one-time bilinear trend surface F0.05(p12,n-p12-1), a secondary trend surface F0.05(p11,n-p11-1); comparison F0.05And F, if F is larger, the trend surface equation is obvious, otherwise, the trend surface equation is not obvious;
step 6, integrating the three model test results of the step 3-5, and selecting an optimal model;
step 7, utilizing a trend surface model formula z2A series of resultsValue and weathering crust thickness h1And (5) overlapping, and importing related mapping software to finally obtain the final ancient landform recovery map.
The trend surface equation in the step (1) can also use a higher-order trend surface model equation of a cubic trend surface and a quartic trend surface according to the complexity of the ancient landform to be restored.
Z in the above steps (2) and (3)i、Andrespectively indicating the stratum stacking thickness, the trend face value and the stratum stacking thickness average value,in the form of a remaining sum of squares,to return the sum of squares, SST=SSD+SSRIs the sum of the squares of the total deviations.
The formula in the step (4) isp is the number of coefficients except for a constant term in the trend surface equation; n is the number of observation points.
The method for successive detection in the step (5) comprises the following steps: first, the regression sum of squares of the higher order polynomial equation is foundRegression sum of squares with lower degree polynomial equationThe difference between the two; this difference is then divided by the difference in degrees of freedom p-q of the sum of squares of the regression to yield the regression mean square error due to the increase in the polynomial degreeFinally, the mean square error is divided by the residual mean square error of the higher-order polynomialObtaining the comparative appropriateness check value F of the two successive order trend surface modelsK+1→K(ii) a Wherein, MSR=SSR/p,MSD=SSD/(n-p-1)。
The method provided by the invention is based on the residual thickness method, combines the overburden stratum thickness and trend surface analysis method, performs more refined recovery on the ancient landform, can perform relatively accurate recovery under the condition that the stratum is degraded, and has wider application range compared with the residual thickness method.
Drawings
FIG. 1 is a diagram of a recovery pattern of ancient landform by a residual thickness method in the prior art;
FIG. 2 is a schematic diagram of the restoration of the weathering crust ancient landform by using the residual thickness method according to the embodiment;
FIG. 3 is a diagram of an example weathered shell and overlying strata;
FIG. 4 is a graph illustrating differential compaction of an overburden of an example weathering crust;
FIG. 5 is a stratigraphic diagram after an example leveling of the overburden top boundary;
FIG. 6 is a graph of paleotopographic structures calculated by the example trend surface analysis method;
FIG. 7 is a schematic diagram of the final restoration of ancient geomorphic geometry of the embodiment;
FIG. 8 is a diagram of the final ancient landform restoration results in the example.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention will be described in further detail with reference to the following embodiments and accompanying drawings. The exemplary embodiments and descriptions of the present invention are provided to explain the present invention, but not to limit the present invention.
The method is mainly used for improving the residual thickness method from two aspects, so that the method is suitable for restoring the weathering crust type ancient landform. First, the target formation thickness plus the overburden thickness is used as a recovery basis to reduce the effects of degradation without weathering degradation of the overburden, as shown in fig. 3, where the top interface is stable and flat compared to the weathering crust top interface and is suitable for leveling. Second, although the overlying strata is stable and flat relative to the weathering crust, due to differential compaction, local undulations may still be present at the same terrain elevation, as shown in FIG. 4, and the topographical surface obtained after flattening it along the top interface may still be inaccurate, as shown in FIG. 5. Therefore, further correction was performed using trend surface analysis.
The trend surface analysis is a mathematical method for simulating the distribution and variation trend of geographic system elements in space by using a mathematical curved surface. The method is characterized in that a two-dimensional nonlinear function is fitted by using a least square method through a regression analysis principle, the distribution rule of the geographic elements on the space is simulated, and the change trend of the geographic elements on the region space is displayed. The ancient earth structural surface is exactly the change trend of the stratum in space, so that more accurate ancient structural surface can be obtained by using trend surface analysis.
The observation surface of the trend surface analysis is composed of a trend surface part and a residual part. The trend surface part reflects the change condition in a regional large range, the residual error part is the difference between the measured value and the corresponding value of the trend function, the local change condition is reflected, and the combination of the two parts is helpful for deep analysis. One basic requirement of trend surface analysis is that the selected trend surface model should have the smallest residual value and the largest trend value, so that fitting accuracy can be sufficient. The spatial trend surface analysis is to analyze trend values and residual values from actual data of geographic element distribution, so as to reveal the trend and law of the geographic element spatial distribution. For data with mild changes, analysis can be performed by using a low-frequency trend surface; for complex and fluctuating data, a polynomial order higher trend surface can be used.
First, a model of the trend surface is built. Thickness z of superposed stratai(xi,yi) (i ═ 1, 2.. times.n) as a dependent variable, geographic position coordinates (x) are determinedi,yi) As an independent variable, the trend fit value isεiAs residual, then:
the core of the trend surface analysis method is to calculate the trend surface from the actual observed value, and a regression analysis method is generally adopted, so that the sum of squares of residual errors Q tends to be minimum, namely:
this is a trend surface fit in the least squares sense.
Fitting values to trendsFormulas of different orders can be set according to needs, and generally polynomial functions of different orders are used for calculation, such as:
in the above formulas, formula (3) corresponds to a first simple linear trend surface, formula (4) corresponds to a first bilinear trend surface, formula (5) corresponds to a second trend surface, formula (6) corresponds to a third trend surface, and a higher-order standard preset formula is analogized, and the formula can also be customized.
After the models are established by using different formulas, the fitting degree of the models needs to be checked, and the models with the highest fitting degree are selected as the basis for the final ancient landform restoration. The model test mainly consists of the following three parts.
1)R2Examination of
Fitting degree coefficient R of trend surface and practical surface2Is an important index for determining the goodness of fit of the regression model. The goodness of fit of the regression model is generally expressed in terms of the proportion of the regression sum of squares in the total dispersion sum of the variable z.
The total sum of squared deviations equals the sum of the regression sum of squares and the remaining sum of squares, i.e.:
wherein,is the residual sum of squares, which represents the effect of random factors on dispersion,
it is the regression sum of squares, which represents the total effect of the dispersion of the independent variables on the dependent variables.
SSRThe larger (or SS)D) Smaller values indicate closer relationship between the dependent variable and the independent variable, and stronger regression regularity and better effect are achieved. This gives equation (8):
R2the larger the fit of the trend surface.
2) Significance F test
The trend surface moderated F test is a significance test for the whole trend surface regression model. The specific method is to determine whether the regression relationship between the variable z and the independent variables x and y is significant or not by utilizing the ratio of the residual square sum to the regression square sum in the total deviation square sum of the variable z. As shown in formula (9):
in the formula, p is the number of coefficients except for a constant term in the trend surface; n is the number of observation points.
At significance level α, look up F of the F distribution Tableα(p, n-p-1) if the calculated F value is greater than the threshold value FαThen, the trend surface equation is considered to be significant; otherwise, it is not significant.
3) Successive examination
First, the regression sum of squares of the higher order polynomial equation is foundRegression sum of squares with lower degree polynomial equationThe difference between the two; this difference is then divided by the difference in degrees of freedom p-q of the sum of squares of the regression to yield the regression mean square error due to the increase in the polynomial degreeFinally, the mean square error is divided by the residual mean square error of the higher-order polynomialObtaining the comparative appropriateness check value F of the two successive order trend surface modelsK+1→K. The specific formula can be seen in table 1.
TABLE 1 regression significance test with K-fold increase to (K + 1-fold) fold polynomial trend surface
Note: in the table, p and q are the number of coefficients except for constant terms in the trend surface equation; n is the number of observation points.
Looking up the table to obtain Fα(p-q, n-p-1) and F obtainedK+1→KIn comparison, if the resulting F value is significant, the higher order polynomial contributes to the regression, and if the F value is not significant, the higher order polynomial does not contribute new to the regression.
After the trend surface model is tested by the 3 methods, a proper model is selected, and then the earth structure as shown by a red dotted line in fig. 6 can be obtained. And then, the weathering crust stratum is superposed on the earth structure, and the ancient landform restoration map (shown in the attached figure 7) can be obtained.
The method according to the invention will be described by taking the example of a weathering crust in a certain area. As shown in Table 2, the data of the geographic coordinate positions of 100 wells and the thickness of the earth strata are combined, the distribution situation of the stratum thickness on the geographic space position is used as a trend surface model to be constructed for restoring the ancient landform, wherein the thickness h of the weathering crust1And thickness h of overburden2The thickness z of the stack of (a) as a dependent variable and the geographic coordinates as independent variables x, y.
TABLE 2 thickness and geographical coordinates of the strata in a certain area
Step 1, constructing a trend surface equation by using a least square method principle according to data in table 2, and obtaining a primary simple trend surface, a primary bilinear trend surface and a secondary trend surface model equation by using an MATLAB software to compile a calculation program as follows:
z11=-5111+0.0004049x-0.0006157y
z12=0.5797+0.0001457x-0.001306y+(3.464e-11)xy
z2=(-8.264e+06)+0.6645x+0.8433y+(-1.44e-08)x2+(-2.451e-08)xy+(-4.272e-08)y2
step 2, substituting the independent variables x and y into the three trend surface model formulas (9), (10) and (11) to obtain a series of trend surface values
Step 3, performing R on the result2And (6) checking. Will zi、Andsubstituting into correlation formula to obtain SSD、SSRAnd SSTThen, R is found2The correlation results are shown in table 3.
TABLE 3 analysis of the trend surface of the formation in a certain area R2Test sum F test results
Watch with watch3 according to R2As a result of the examination, the determination coefficients of the primary linear trend surface, the primary bilinear trend surface, and the secondary trend surface are 0.3047, 0.3760, and 0.5852, respectively, and the fitness of the secondary trend surface is higher.
And 4, carrying out F test on the result. The correlation data was substituted into equation (9) to find F values of the three models. In the case of confidence level α being 0.05, F of the one-time single-linear trend surface is obtained by table lookup0.05(2,97)=3.09<Calculated value F1121.25; one time bilinear trend surface F0.05(3,96)=2.70<Calculated value F1219.28; second trend surface F0.05(9,90)=1.98<Calculated value F214.11. Therefore, the three equations are all significant from the F test, and the fitting result is reasonable.
And 5, carrying out successive testing on the trend surface of the result, and respectively comparing the primary single-linear trend surface with the primary bilinear trend surface, the primary single-linear trend surface with the secondary trend surface and the primary bilinear trend surface with the secondary trend surface.
The results of the successive-test analysis are shown in table 4 for the one-time single-linear trend surface and the one-time bilinear trend surface.
TABLE 4 successive testing analysis Table of Once-Mono-Linear and Once-bilinear trend surface regression model
Increase from a one-time single-linear trend surface to a one-time double-linear trend surface, F12→1110.97. In the case of the confidence level alpha being 0.05, F is obtained by looking up the F distribution table0.05(1,96)=3.96<F12→1110.97. Therefore, the increase of the primary linear trend surface to the primary bilinear trend surface makes a new contribution to the regression equation, namely, the primary bilinear trend surface has higher fitting degree than the primary single linear trend surface.
The results of the sequential test analysis for the primary single linear trend surface and the secondary trend surface are shown in table 5.
TABLE 5 successive testing analysis of the Primary Mono-Linear and Secondary Trend surface regression models
Increase from a primary single linear trend surface to a secondary trend surface, F2→118.70. In the case of the confidence level alpha being 0.05, F is obtained by looking up the F distribution table0.05(7,90)=2.11<F2→118.70. The rising of the primary, single-linear trend surface to the secondary trend surface makes a new contribution to the regression equation, i.e., the secondary trend surface has a higher fitness than the primary, single-linear trend surface.
The results of the sequential testing analysis for the primary bilinear trend surface and the secondary trend surface are shown in table 6.
TABLE 6 successive testing analysis of the Primary bilinear and Secondary Trend surface regression models
Increase from a primary single linear trend surface to a secondary trend surface, F2→127.57. In the case of the confidence level alpha being 0.05, F is obtained by looking up the F distribution table0.05(6,90)=2.2<F2→127.57. Therefore, the rising of the primary bilinear trend surface to the secondary trend surface makes a new contribution to the regression equation, i.e., the secondary trend surface has a higher fitness than the primary bilinear trend surface.
By combining the three analyses, the model with the highest fitness of the secondary trend surface can be obtained.
Step 6, synthesizing the three model test results of steps 3-5, and obtaining the result model of the secondary trend surface with the highest fitting degree, so that the trend surface model formula is used:
z2=(-8.264e+06)+0.6645x+0.8433y+(-1.44e-08)x2+(-2.451e-08)xy+(-4.272e-08)y2
step 7, utilizing a trend surface model formula z2A series of resultsValue and weathering crust thickness h1And (4) superposing, and importing related mapping software to finally obtain the final ancient landform restoration map (figure 8).
Note: in this example, all the thickness values used for convenience in the calculation part are absolute values, and actually, if the ground is taken as a zero interface and the upward direction is positive according to a general situation, the trend face value isShould be negative, so proceedAnd h1Coping before superpositionAnd taking the negative value.
The above description is only an embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that are not thought of through the inventive work should be included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope defined by the claims.
Claims (5)
1. An ancient landform restoration method based on trend surface analysis is characterized by comprising the following steps:
step 1, collecting thickness h of weathering crust of target stratum1Thickness h of overburden2And geographic coordinates x and y, calculating the thickness z of the superposed stratum according to the obtained data, constructing a trend surface equation by using a least square method principle, and programming a calculation program by using MATLAB software to obtain a primary simple trend surface equation z11Once bilinear trend surface z12Equation z of the quadratic trend surface model2;
Step 2Substituting the independent variables x and y into the three trend surface model formulas to obtain a series of trend surface values
Step 3, performing R on the result2Checking; will zi、Andsubstituting into correlation formula to obtain SSD、SSRAnd SSTThen, R is found2Comparison ofThe higher the numerical value is, the higher the fitting degree is;
step 4, carrying out F inspection on the result; substituting the correlation data into a formula to obtain an F value of a corresponding model, and looking up a table to obtain the F of the single linear trend surface when the confidence level alpha is 0.050.05(p11,n-p11-1), a one-time bilinear trend surface F0.05(p12,n-p12-1), a secondary trend surface F0.05(p11,n-p11-1); comparison F0.05And F, if F is larger, the trend surface equation is obvious, otherwise, the trend surface equation is not obvious;
step 5, carrying out successive testing on the trend surface of the result, and respectively comparing the primary single-linear trend surface with the primary bilinear trend surface, the primary single-linear trend surface with the secondary trend surface and the primary bilinear trend surface with the secondary trend surface;
step 6, integrating the three model test results of the step 3-5, and selecting an optimal model;
2. The method according to claim 1, wherein the method comprises: the trend surface equation in the step (1) can also use a higher-order trend surface model equation of a cubic trend surface and a quartic trend surface according to the complexity of the ancient landform to be restored.
3. The method according to claim 1, wherein the method comprises: z in the above steps (2) and (3)i、Andrespectively indicating the stratum stacking thickness, the trend face value and the stratum stacking thickness average value,in the form of a remaining sum of squares,to return the sum of squares, SST=SSD+SSRIs the sum of the squares of the total deviations.
5. The method according to claim 1, wherein the ancient landform restoration method is based on trend surface analysisThe method is characterized in that: the method for successive detection in the step (5) comprises the following steps: first, the regression sum of squares of the higher order polynomial equation is foundRegression sum of squares with lower degree polynomial equationThe difference between the two; this difference is then divided by the difference in degrees of freedom p-q of the sum of squares of the regression to yield the regression mean square error due to the increase in the polynomial degreeFinally, the mean square error is divided by the residual mean square error of the higher-order polynomialObtaining the comparative appropriateness check value F of the two successive order trend surface modelsK+1→K(ii) a Wherein, MSR=SSR/p,MSD=SSD/(n-p-1)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110575415.1A CN113312695B (en) | 2021-05-26 | 2021-05-26 | Ancient landform restoration method based on trend surface analysis |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110575415.1A CN113312695B (en) | 2021-05-26 | 2021-05-26 | Ancient landform restoration method based on trend surface analysis |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113312695A true CN113312695A (en) | 2021-08-27 |
CN113312695B CN113312695B (en) | 2022-03-11 |
Family
ID=77374877
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110575415.1A Expired - Fee Related CN113312695B (en) | 2021-05-26 | 2021-05-26 | Ancient landform restoration method based on trend surface analysis |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113312695B (en) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114415233A (en) * | 2022-01-20 | 2022-04-29 | 北京师范大学 | Method for recovering landform form of exploited fossil energy stratum |
CN118037135A (en) * | 2024-04-11 | 2024-05-14 | 生态环境部卫星环境应用中心 | Qualitative analysis method for influence of vegetation coverage change on ecosystem service |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425315A (en) * | 2015-12-21 | 2016-03-23 | 中国石油大学(北京) | Inversion method of microscale deposition stage palaeogeomorphology in weak construction deformation range |
CN105740574A (en) * | 2016-03-03 | 2016-07-06 | 中国石油大学(北京) | Palaeo-geomorphic cross section recovery method |
CN105956238A (en) * | 2016-04-22 | 2016-09-21 | 中国石油天然气股份有限公司 | Ancient landform restoration method and device |
CN107315195A (en) * | 2017-05-31 | 2017-11-03 | 中国石油天然气股份有限公司 | Method and device for determining micro landform |
CN108196302A (en) * | 2017-11-28 | 2018-06-22 | 中国石油天然气股份有限公司 | Forecasting method and device for dessert region of dolomite fracture-cave development reservoir |
CN110361781A (en) * | 2019-08-14 | 2019-10-22 | 西南石油大学 | A kind of deposition phase palaeogeomorphology based on the analysis of seismic reflection amplitude |
AU2020102025A4 (en) * | 2020-07-02 | 2020-10-01 | China University Of Geosciences, Beijing | Method for Reconstructing Basin Paleogeomorphology |
-
2021
- 2021-05-26 CN CN202110575415.1A patent/CN113312695B/en not_active Expired - Fee Related
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105425315A (en) * | 2015-12-21 | 2016-03-23 | 中国石油大学(北京) | Inversion method of microscale deposition stage palaeogeomorphology in weak construction deformation range |
CN105740574A (en) * | 2016-03-03 | 2016-07-06 | 中国石油大学(北京) | Palaeo-geomorphic cross section recovery method |
CN105956238A (en) * | 2016-04-22 | 2016-09-21 | 中国石油天然气股份有限公司 | Ancient landform restoration method and device |
CN107315195A (en) * | 2017-05-31 | 2017-11-03 | 中国石油天然气股份有限公司 | Method and device for determining micro landform |
CN108196302A (en) * | 2017-11-28 | 2018-06-22 | 中国石油天然气股份有限公司 | Forecasting method and device for dessert region of dolomite fracture-cave development reservoir |
CN110361781A (en) * | 2019-08-14 | 2019-10-22 | 西南石油大学 | A kind of deposition phase palaeogeomorphology based on the analysis of seismic reflection amplitude |
AU2020102025A4 (en) * | 2020-07-02 | 2020-10-01 | China University Of Geosciences, Beijing | Method for Reconstructing Basin Paleogeomorphology |
Non-Patent Citations (3)
Title |
---|
OLEG SIZOV: ""Predictive Mapping of Glacial and Fluvioglacial Landforms in the Nadym River Basin (North of West Siberia) With TanDEM-X DEM"", 《IEEE JOURNAL OF SELECTED TOPICS IN APPLIED EARTH OBSERVATIONS AND REMOTE SENSING ( VOLUME: 14)》 * |
左丽群: ""古地貌恢复方法综述"", 《石油地质与工程》 * |
贾会冲;丁晓琪: ""鄂尔多斯盆地大牛地气田马五1+2段白云岩喀斯特储层特征"", 《成都理工大学学报(自然科学版)》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN114415233A (en) * | 2022-01-20 | 2022-04-29 | 北京师范大学 | Method for recovering landform form of exploited fossil energy stratum |
CN118037135A (en) * | 2024-04-11 | 2024-05-14 | 生态环境部卫星环境应用中心 | Qualitative analysis method for influence of vegetation coverage change on ecosystem service |
CN118037135B (en) * | 2024-04-11 | 2024-07-09 | 生态环境部卫星环境应用中心 | Qualitative analysis method for influence of vegetation coverage change on ecosystem service |
Also Published As
Publication number | Publication date |
---|---|
CN113312695B (en) | 2022-03-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN113312695B (en) | Ancient landform restoration method based on trend surface analysis | |
Kang et al. | Topographic mapping on large-scale tidal flats with an iterative approach on the waterline method | |
Shahid et al. | GIS integration of remote sensing and electrical sounding data for hydrogeological exploration | |
CN103234920B (en) | Based on the underground water enriching appraisal procedure of sensor information | |
CN110163965B (en) | Solid three-dimensional model construction method for monoclinic rock stratum structure landform | |
CN102707332A (en) | Interpretation and evaluation method for reservoir area engineering geological survey | |
CN111161415B (en) | Three-dimensional space simulation method for karst development channel of underground water in karst area | |
CN103412333B (en) | A kind of static correction base level defining method | |
Zhang et al. | A case study on integrated modeling of spatial information of a complex geological body | |
Adcroft | Representation of topography by porous barriers and objective interpolation of topographic data | |
CN112231818B (en) | Mountain city three-dimensional geological modeling method based on BIM | |
Chitsazan et al. | Land subsidence susceptibility mapping using PWRSTFAL framework and analytic hierarchy process: fuzzy method (case study: Damaneh-Daran Plain in the west of Isfahan Province, Iran) | |
Fregoso et al. | Sediment deposition, erosion, and bathymetric change in Central San Francisco Bay: 1855-1979 | |
CN112800518B (en) | Stratum surface model correction method utilizing adjacent stratum cross-correlation constraint | |
CN113985490B (en) | Method and device for carrying out surface gravity simulation by utilizing terrain and crust density data | |
Gerhart et al. | Evaluation of the ground-water resources of the lower Susquehanna River basin, Pennsylvania and Maryland | |
Simo et al. | Geostatistical investigations for suitable mapping of the water table: the Bordeaux case (France) | |
Cafaro et al. | Landslide hazard assessment and judgment of reliability: a geomechanical approach | |
Gao et al. | Mathematical Interpolation and Correction of Three-Dimensional Modelling of High-Speed Railway. | |
CN110687591B (en) | Method for determining physical property parameters of coal bed and surrounding rock based on waveform matching of prior data | |
CN106125149A (en) | The Point-mass Model middle-shallow layer optimal buried depth of high-resolution point mass determines method | |
Saleh et al. | Geometric investigation of Al-Wind Dam reservoir Northeastern Iraq, using digital elevation models and spatial analyses system | |
CN115469361B (en) | Clastic rock stratum three-dimensional geological modeling method | |
Desjarlais | Analysis of piezometer measurements and stratigraphy at Sugar Lake Dam | |
Perton et al. | A New Method to Assess the Volume of Volcanic Landforms and Associated Deposits |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20220311 |