CN116611267B - All-constraint configuration method for priori boundary structure in geophysical prospecting data regularized inversion - Google Patents
All-constraint configuration method for priori boundary structure in geophysical prospecting data regularized inversion Download PDFInfo
- Publication number
- CN116611267B CN116611267B CN202310884362.0A CN202310884362A CN116611267B CN 116611267 B CN116611267 B CN 116611267B CN 202310884362 A CN202310884362 A CN 202310884362A CN 116611267 B CN116611267 B CN 116611267B
- Authority
- CN
- China
- Prior art keywords
- boundary
- structural
- configuration
- constraint
- units
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 31
- 238000011156 evaluation Methods 0.000 claims abstract description 87
- 238000012545 processing Methods 0.000 claims description 19
- 239000011159 matrix material Substances 0.000 claims description 12
- 238000004364 calculation method Methods 0.000 claims description 8
- 239000002131 composite material Substances 0.000 claims description 7
- 238000012937 correction Methods 0.000 claims description 4
- 230000008521 reorganization Effects 0.000 claims description 3
- 230000003313 weakening effect Effects 0.000 claims description 3
- 238000009499 grossing Methods 0.000 claims description 2
- 238000001514 detection method Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 6
- 210000004027 cell Anatomy 0.000 description 4
- 230000007547 defect Effects 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 239000000546 pharmaceutical excipient Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 238000000053 physical method Methods 0.000 description 1
- 238000013139 quantization Methods 0.000 description 1
- 210000003537 structural cell Anatomy 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000002123 temporal effect Effects 0.000 description 1
Classifications
-
- 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
- G06F18/00—Pattern recognition
- G06F18/10—Pre-processing; Data cleansing
- G06F18/15—Statistical pre-processing, e.g. techniques for normalisation or restoring missing data
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/23—Clustering techniques
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/29—Graphical models, e.g. Bayesian networks
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F18/00—Pattern recognition
- G06F18/20—Analysing
- G06F18/29—Graphical models, e.g. Bayesian networks
- G06F18/295—Markov models or related models, e.g. semi-Markov models; Markov random fields; Networks embedding Markov models
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
Landscapes
- Engineering & Computer Science (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- General Engineering & Computer Science (AREA)
- Evolutionary Computation (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Evolutionary Biology (AREA)
- Life Sciences & Earth Sciences (AREA)
- Artificial Intelligence (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Geometry (AREA)
- Probability & Statistics with Applications (AREA)
- Computer Hardware Design (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application discloses a priori boundary structure full-constraint configuration method in geophysical prospecting data regularization inversion, which comprises the steps of firstly quantifying evaluation coefficients of regularization items according to medium information on two sides of a boundary; then, judging the type of the structure boundary by relying on prior structure information; formulating a configuration combination of the structure full constraint evaluation coefficients according to different structure types; calibrating the evaluation coefficient configuration of the boundary of the special structure; finally, the structure full constraint assembly and format standardization are carried out to form regularization items containing prior structure full constraint configuration.
Description
Technical Field
The application relates to the technical field of data processing, in particular to a priori boundary structure full-constraint configuration method in geophysical prospecting data regularization inversion.
Background
Geophysical prospecting techniques mainly develop geospatial medium detection by physical methods. The detection result has continuity, so that the working efficiency can be greatly improved, the cost is reduced, and the method is widely applied to the aspects of exploration of various ore body resources such as coal, petroleum, metal and the like, environmental pollution investigation, engineering stability prediction and the like. Geophysical prospecting techniques infer the dielectric profile of a target body by measuring the spatial and temporal variations of the physical field of the target and performing an inversion analysis on a large number of measured data. In the inversion process, a regularization method is often utilized, and model roughness is introduced to solve the optimal solution of the inversion model. Under the condition that no priori information exists in the detection target area, the roughness of the whole inversion model is set to be the same value by default, and the inversion result is represented as a smooth model. However, most of the detection targets are mostly non-homogeneous bodies, so that the inversion results deviate from the same model roughness. Generally, the stronger the detected target is heterogeneous, the more obvious such defects are. Therefore, how to utilize prior structure information to configure a reasonable model roughness, and further improve accuracy of geophysical prospecting data inversion is an urgent problem to be solved.
Disclosure of Invention
The application overcomes the defects of the prior art and provides a prior boundary structure full-constraint configuration method in geophysical prospecting data regularization inversion.
The technical scheme adopted by the application for achieving the purpose is as follows:
the application discloses a prior boundary structure full-constraint configuration method in geophysical prospecting data regularization inversion, which comprises the following steps:
s1: quantizing the evaluation coefficients of the regularization items according to the medium information on two sides of the boundary;
s2: judging the type of the structure boundary according to the prior structure boundary information;
s3: formulating a configuration combination of the structure full constraint evaluation coefficients according to different structure types;
s4: configuring and calibrating the evaluation coefficient of the boundary of the special structure;
s5: the structure full constraint assembly and format normalization form a regularized term comprising a priori structure full bundle configuration.
Further, in a preferred embodiment of the present application, the regularization term calculation formula is:
；
in the formula ,is a regularization term; />Regularizing global coefficients; />Is a model physical attribute matrix;
initial physical attribute matrix for the model; />Transpose the matrix; />Smoothing matrixes respectively in the horizontal direction, the vertical direction, the first diagonal direction and the second diagonal direction of the element grid;
the evaluation coefficients of regularization items of the unit grid in the horizontal direction, the vertical direction, the first diagonal direction and the second diagonal direction are respectively classified into 1 and 1 respectively>Three kinds of>The value of (2) is smaller than 1, in order to reduce the inter-cell spaceDegree of association (s)/(s)>In order to increase the degree of association between units, the evaluation coefficient is assigned a value of 1 in order to maintain the degree of association between conventional units.
Further, in a preferred embodiment of the present application,the calculation formula of (2) is as follows: />The calculation formula of (2) is as follows: />, wherein />Is the large value in the physical field characteristic values of the medium on both sides of the same physical field lower boundary, wherein +.>Is the small value in the physical field characteristic value of the medium at both sides of the same physical field lower boundary, +.>As the correction coefficient, its value range: />。
Further, in a preferred embodiment of the present application, the determining the type of the structural boundary specifically includes:
judging the type of the structural boundary according to the approximation degree of the trend of the structural boundary and the passing unit grid in the horizontal direction, the vertical direction, the first type diagonal direction and the second type diagonal direction; the structural boundary types are divided into four boundary types, namely a horizontal boundary, a vertical boundary, a first type inclined boundary and a second type inclined boundary.
Further, in a preferred embodiment of the present application, the determining the type of the structural boundary further includes:
according to the cell grids with different sizes, defining the included angle between the diagonal direction of the second class and the horizontal direction asThe included angle between the diagonal direction of the first type and the horizontal direction is +.>The angle between the first diagonal direction and the vertical direction is +.>The method comprises the steps of carrying out a first treatment on the surface of the The angle between the trend of the boundary of the prior structure and the vertical direction is defined as +.>；
When (when)When the prior structure boundary is a horizontal boundary;
when (when)When the prior structure boundary is a vertical boundary;
when (when)When the prior structure boundary is the first type of inclined boundary;
when (when)The a priori structure boundaries are then second class tilt boundaries.
Further, in a preferred embodiment of the present application, in the formulating structural fully-constrained evaluation coefficient configuration combination, the structural fully-constrained evaluation coefficient configuration principle is as follows: for the two-dimensional model of the structured grid, defining units on the boundary of the structure as structural units, carrying out regularization constraint processing between each structural unit and 8 adjacent units, and configuring special evaluation coefficients only for regularization parameters related to the structural units; excipient (C)Evaluation of coefficients for structural units on the same side of the boundaryEnhancing the association in the structural direction; structural element evaluation coefficient assigned across structural boundary +.>Weakening the correlation in the structural direction; setting the evaluation coefficient value to 1 between the structural unit and the non-structural unit, and keeping normal association.
Further, in a preferred embodiment of the present application, in the configuration combination of the fully constrained evaluation coefficients of the structure, four different configuration combinations of the evaluation coefficients are adopted for four types of boundary types of the structure:
for the horizontal boundary, as two sides of the structural boundary are respectively provided with a layer of structural unit layer, the regularized constraint processing between the structural unit and 8 adjacent units is completed by respectively applying a layer of evaluation coefficient configuration to the structural units at the upper layer and the structural units at the lower layer of the horizontal boundary according to the structural full constraint evaluation coefficient configuration principle;
for the vertical boundary, as two sides of the structural boundary are respectively provided with a layer of structural unit layer, the regularized constraint processing between the structural unit and 8 adjacent units is completed by respectively applying a layer of evaluation coefficient configuration to the structural units on the left side and the structural units on the right side of the vertical boundary according to the structural full constraint evaluation coefficient configuration principle;
for the first type inclined boundary, as two layers of structural unit layers are respectively arranged at two sides of the structural boundary, two layers of evaluation coefficient configuration are respectively applied to the upper left structural unit and the lower right structural unit of the first type inclined boundary according to the structural full-constraint evaluation coefficient configuration principle for completing regularization constraint processing between the structural unit and 8 adjacent structural units;
for the second type of inclined boundary, two layers of structural unit layers are respectively arranged on two sides of the structural boundary, so that two layers of evaluation coefficient configuration are respectively applied to the structural unit at the upper right side and the structural unit at the lower left side of the first type of inclined boundary according to the structural full-constraint evaluation coefficient configuration principle for completing regularization constraint processing between the structural unit and 8 adjacent structural units.
Further, in a preferred embodiment of the present application, the configuration calibration of the evaluation coefficient of the boundary of the special structure is specifically: for the composite boundary region which does not belong to four types of conventional structure boundary types, the configuration principle is required to be configured according to the structure full constraint evaluation coefficient, and special processing is required to be carried out on special structural units according to specific composite boundary combination conditions.
Further, in a preferred embodiment of the present application, the structure full constraint assembly and format normalization form a regularized item comprising a priori structure full bundle configuration, specifically:
according to the evaluation coefficient configuration completed by the prior structure boundary, the unit grid dimension of the inversion model is calculated according to the formula:
format standard reorganization in (3) will finally +.>And sequentially carrying the above formulas to form regularization terms containing the prior structure full-constraint configuration.
The application solves the technical defects existing in the background technology, and has the following beneficial effects: firstly, quantifying an evaluation coefficient of a regularization term according to medium information on two sides of a boundary; then, judging the type of the structure boundary by relying on prior structure information; formulating a configuration combination of the structure full constraint evaluation coefficients according to different structure types; calibrating the evaluation coefficient configuration of the boundary of the special structure; and finally, carrying out full constraint assembly and format standardization on the structure to form a regularization term containing prior full constraint configuration of the structure. The application provides a method for extracting structural full constraint by using priori structural information, and brings the structural full constraint into geophysical prospecting data inversion so as to improve accuracy of geophysical prospecting data inversion.
Drawings
In order to more clearly illustrate the embodiments of the application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the application, and that other embodiments of the drawings can be obtained according to these drawings without inventive effort for a person skilled in the art.
FIG. 1 is a flow chart of a prior boundary structure full constraint configuration method in geophysical prospecting data regularization inversion;
FIG. 2 is a schematic diagram of four structural boundary types in the present application;
FIG. 3 is a schematic diagram of a fully constrained evaluation coefficient configuration of the horizontal boundary structure type of the present application;
FIG. 4 is a schematic diagram of a fully constrained evaluation coefficient configuration of the vertical boundary structure type of the present application;
FIG. 5 is a schematic diagram of a fully constrained evaluation coefficient configuration of a first type of sloped boundary structure according to the present application;
FIG. 6 is a schematic diagram of a fully constrained evaluation coefficient configuration for a second type of inclined boundary structure of the present application;
FIG. 7 is a schematic diagram of the configuration of evaluation coefficients of a boundary of a specific structure according to the present application;
in the figure: 1. horizontal direction, 2, vertical direction, 3, first class diagonal direction, 4, second class diagonal direction, 5, large value in physical field characteristic value of medium at two sides of the same physical field lower boundary, 6, small value in physical field characteristic value of medium at two sides of the same physical field lower boundary, 7, horizontal boundary, 8, vertical boundary, 9, first class inclined boundary, 10, second class inclined boundary, 11, structural unit, 12, upper structural unit of horizontal boundary, 13, structural unit of lower layer of horizontal boundary, 14, left structural unit of vertical boundary, 15, right structural unit of vertical boundary, 16, upper right structural unit of first class inclined boundary, 17, lower left structural unit of first class inclined boundary, 18, upper left structural unit of second class inclined boundary, 19, lower right structural unit of second class inclined boundary.
Detailed Description
In order that the above-recited objects, features and advantages of the present application will be more clearly understood, a more particular description of the application will be rendered by reference to the appended drawings and appended detailed description. It should be noted that, without conflict, the embodiments of the present application and features in the embodiments may be combined with each other.
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present application, however, the present application may be practiced in other ways than those described herein, and therefore the scope of the present application is not limited to the specific embodiments disclosed below.
As shown in fig. 1 and 2, the application provides a prior boundary structure full-constraint configuration method in geophysical prospecting data regularization inversion, which comprises the following steps:
according to the evaluation coefficients of the medium information quantization regularization terms at two sides of the boundary, the regularization term formula is as follows:
wherein->For regularized item->For regularizing global coefficients, +.>For the model physical attribute matrix +.>For model initial physical property matrix +.>A smoothed matrix of the cell grid in the horizontal direction 1, in the vertical direction 2, in the first diagonal direction 3 and in the second diagonal direction 4, respectively. />Regularization of cell grid horizontal direction 1, vertical direction 2, first class diagonal direction 3 and second class diagonal direction 4, respectivelyEvaluation coefficient of the term. The four evaluation coefficient assignment types are classified as 1, < >>Three kinds of>The value of (2) is smaller than 1, in order to cut down the degree of association between units, < >>In order to increase the degree of association between units, the evaluation coefficient is assigned a value of 1 in order to maintain the degree of association between conventional units.
Optionally, the evaluating coefficient of the regularization term is quantized according to the medium information on two sides of the boundary, and is characterized in that,the calculation formula of (2) is as follows: />，/>The calculation formula of (2) is as follows: />, wherein />Is a large value 5 in the physical field characteristic values of the medium on both sides of the same physical field lower boundary, wherein +.>Is 6, which is a small value in the physical field characteristic value of the medium on both sides of the same physical field lower boundary>As the correction coefficient, its value range: />. Here we defineFor a dielectric resistivity value of 20Ω×m +.>For a medium resistivity value of 1 Ω m, correction factor +.>0.5, then->Calculated as 10>Calculated as 0.1.
Optionally, the structure boundary type is determined according to the prior structure information, and the structure boundary type is divided into four boundary types of a horizontal boundary 7, a vertical boundary 8, a first type inclined boundary 9 and a second type inclined boundary 10. And judging the type of the structural boundary according to the approximation degree of the trend of the structural boundary and the passing unit grids of the structural boundary in the horizontal direction 1, the vertical direction 2, the first type diagonal direction 3 and the second type diagonal direction 4.
Optionally, the structure boundary type is determined according to the prior structure information, and according to the unit grids with different sizes, we define the included angle between the diagonal direction of the second type and the horizontal direction asThe included angle between the diagonal direction of the first type and the horizontal direction is +.>The angle between the first diagonal direction and the vertical direction is +.>. The angle between the trend of the boundary of the prior structure and the vertical direction is defined as +.>。
When (when)When the prior structure boundary is a horizontal boundary 7;
when (when)When, the prior structure boundary is a vertical boundary 8;
when (when)When the prior structure boundary is the first type of inclined boundary 9;
when (when)The a priori structure boundaries are then the second type of tilt boundaries 10.
Here we take a square-sized grid as an example, then45 deg., then:
when (when)When the prior structure boundary is a horizontal boundary 7;
when (when)When, the prior structure boundary is a vertical boundary 8;
when (when)When the prior structure boundary is the first type of inclined boundary 9;
when (when)The a priori structure boundaries are then the second type of tilt boundaries 10.
Optionally, the configuration combination of the structural full-constraint evaluation coefficients is formulated according to different structural types, and the configuration principle of the structural full-constraint evaluation coefficients is as follows: for a structured grid two-dimensional model, we define cells on the structural boundary as structural cells 11, eachThe structural unit 11 should perform regularization constraint processing between 8 units (units with serial numbers of 1-8 in fig. 3-6) adjacent to the structural unit 11, and we only configure special evaluation coefficients for regularization parameters related to the structural unit 11. Structural element evaluation coefficient assigned to boundary ipsilateralEnhancing the association in the structural direction; structural element evaluation coefficient assigned across structural boundary +.>Weakening the correlation in the structural direction; setting the evaluation coefficient value to 1 between the structural unit and the non-structural unit, and keeping normal association.
Optionally, the configuration combination of the structural full constraint evaluation coefficients is formulated according to different structural types, and is characterized in that for four types of structural boundary types, four different configuration combinations of the evaluation coefficients are adopted.
For the horizontal boundary 7, since there is one layer of the structural unit layer 11 on each side of the structural boundary, the regularization constraint processing between the structural unit 11 and 8 units (units numbered 1 to 8 in fig. 3) adjacent to the structural unit needs to apply one layer of the configuration of the evaluation coefficients to each of the structural units 12 and 13 on the upper layer of the horizontal boundary according to the above evaluation coefficient application principle.
For the vertical boundary 8, since there is one layer of structural unit layer on each side of the structural boundary, the task of performing regularization constraint processing between the structural unit 11 and 8 units (units numbered 1-8 in fig. 4) adjacent to the structural unit needs to apply one layer of evaluation coefficient configuration to the structural unit 14 on the left side and the structural unit 15 on the right side of the vertical boundary according to the above evaluation coefficient application principle.
For the first type inclined boundary 9, since two layers of structural unit layers are respectively arranged on two sides of the structural boundary, two layers of evaluation coefficient configuration are required to be respectively applied to the structural unit 16 at the right upper side and the structural unit 17 at the left lower side of the first type inclined boundary 9 according to the above evaluation coefficient application principle to complete the regularization constraint processing between the structural unit 11 and 8 units (the units with the serial numbers of 1-8 in fig. 5) adjacent to the structural unit.
For the second type of inclined boundary 10, since two layers of structural unit layers are provided on both sides of the structural boundary, two layers of evaluation coefficient arrangements are required to be applied to the upper left structural unit 18 and the lower right structural unit 19 of the second type of inclined boundary 10 according to the above evaluation coefficient application principle to complete the regularization constraint processing between the structural unit 11 and 8 units (the units with the numbers 1 to 8 in fig. 6) adjacent to the structural unit.
Optionally, the configuration and calibration of the evaluation coefficients of the boundaries of the special structure require special treatment of the special structural units according to the specific composite boundary combination condition according to the application principle of the evaluation coefficients as above for the composite boundary regions not belonging to four types of conventional structural boundary types. The composite boundary shown in fig. 7, optionally, the structure full constraint assembly and format normalization is characterized in that, for the configuration of the evaluation coefficients completed according to the prior structure boundary, the dimension of the unit grid of the inversion model is required according to the formula:
format standard reorganization in (3) will finally +.>And sequentially carrying the above formulas to form regularization terms containing the prior structure full-constraint configuration.
In addition, the full constraint configuration method further comprises the following steps: acquiring inverted geophysical prospecting data, constructing a Markov random field model, importing the inverted geophysical prospecting data into the Markov random field model, and calculating the Markov distance of each piece of inverted geophysical prospecting data based on the Markov distance metric;
presetting a clustering center, and calculating the Euclidean distance between each piece of inverted geophysical prospecting data and the clustering center based on the clustering center;
calculating an update coefficient according to the Markov distance and the Euclidean distance, importing the update coefficient into the Markov random field model, and carrying out data update on each piece of inverted geophysical prospecting data based on the update coefficient so as to obtain an updated cluster center;
building a membership matrix based on the updated clustering center, importing the inverted geophysical prospecting data into the membership matrix, and carrying out secondary deduction on the inverted geophysical prospecting data in the membership matrix to obtain priori probability values of the inverted geophysical prospecting data;
and removing the inverted geophysical prospecting data with the priori probability value larger than the preset probability value.
The method is used for verifying the geophysical prospecting data after modeling, comparing the inversion data with the actual data, evaluating the reliability and accuracy of the inversion data, and eliminating the inversion data which do not accord with the physical meaning and the actual situation so as to establish a more reliable inversion model.
In addition, the full constraint configuration method further comprises the following steps:
acquiring inverted geophysical prospecting data and acquiring inverted preset geophysical prospecting data;
constructing an evaluation system, determining an evaluation index based on the inverted geophysical prospecting data, and constructing an evaluation score based on the inverted preset geophysical prospecting data;
calculating the weight occupancy rate between the evaluation index and the evaluation score by using an analytic hierarchy process; comparing the weight occupancy rate with a preset weight occupancy rate;
if the weight duty ratio is larger than the preset weight duty ratio, acquiring working parameters of the geophysical prospecting equipment during detection, and generating a relevance text based on the weight duty ratio and the working parameters;
and importing the relevance text into a Bayesian network for fault deduction, deducting to obtain corresponding fault sub-equipment, generating fault information based on the fault sub-equipment, and outputting the fault information.
When geophysical prospecting data are collected through geophysical prospecting equipment such as a geomagnetic instrument and an electrode group, the phenomenon of drift and loss of the data obtained through initial collection can be caused due to equipment faults, inversion data are abnormal synchronously, and therefore whether the equipment breaks down in the collecting process is judged through deduction of the inversion data.
Compared with the prior art, the application has the beneficial effects that:
according to the application, through the evaluation coefficient of the regularization term in the geophysical regularization inversion process of the system configuration, the known structure constraint information is completely imported into the regularization formula, the model roughness conforming to the prior structure information is constructed to constrain inversion, the inversion accuracy of geophysical measurement data is improved, and the interpretation of geophysical detection results and the identification work of detection targets can be effectively optimized.
In the several embodiments provided by the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. The above described device embodiments are only illustrative, e.g. the division of the units is only one logical function division, and there may be other divisions in practice, such as: multiple units or components may be combined or may be integrated into another system, or some features may be omitted, or not performed. In addition, the various components shown or discussed may be coupled or directly coupled or communicatively coupled to each other via some interface, whether indirectly coupled or communicatively coupled to devices or units, whether electrically, mechanically, or otherwise. The units described above as separate components may or may not be physically separate, and components shown as units may or may not be physical units; can be located in one place or distributed to a plurality of network units; some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional unit in each embodiment of the present application may be integrated in one processing unit, or each unit may be separately used as one unit, or two or more units may be integrated in one unit; the integrated units may be implemented in hardware or in hardware plus software functional units.
Those of ordinary skill in the art will appreciate that: all or part of the steps for implementing the above method embodiments may be implemented by hardware related to program instructions, and the foregoing program may be stored in a computer readable storage medium, where the program, when executed, performs steps including the above method embodiments; and the aforementioned storage medium includes: a mobile storage device, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk or an optical disk, or the like, which can store program codes. Alternatively, the above-described integrated units of the present application may be stored in a computer-readable storage medium if implemented in the form of software functional modules and sold or used as separate products. Based on such understanding, the technical solutions of the embodiments of the present application may be embodied in essence or a part contributing to the prior art in the form of a software product stored in a storage medium, including several instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to execute all or part of the methods of the embodiments of the present application. And the aforementioned storage medium includes: a removable storage device, ROM, RAM, magnetic or optical disk, or other medium capable of storing program code.
The foregoing is merely illustrative embodiments of the present application, but the scope of the present application is not limited thereto, and any person skilled in the art can easily think about variations or substitutions within the technical scope of the present application, and the application should be covered. Therefore, the protection scope of the application is subject to the protection scope of the claims.
Claims (3)
1. A prior boundary structure full-constraint configuration method in geophysical prospecting data regularization inversion is characterized by comprising the following steps:
s1: quantizing the evaluation coefficients of the regularization items according to the medium information on two sides of the boundary;
s2: judging the type of the structure boundary according to the prior structure boundary information;
s3: formulating a configuration combination of the structure full constraint evaluation coefficients according to different structure types;
s4: configuring and calibrating the evaluation coefficient of the boundary of the special structure;
s5: the method comprises the steps of full constraint assembly and format standardization of a structure to form a regularization term containing prior structure full beam configuration;
the regularization term calculation formula is as follows:
；
in the formula ,is a regularization term; />Regularizing global coefficients; />Is a model physical attribute matrix; />Initial physical attribute matrix for the model; />Transpose the matrix; />，/>，/> and />Smoothing matrixes respectively in the horizontal direction, the vertical direction, the first diagonal direction and the second diagonal direction of the element grid; />，/>，/>，/>The evaluation coefficients of regularization items of the unit grid in the horizontal direction, the vertical direction, the first diagonal direction and the second diagonal direction are respectively classified into 1 and 1 respectively> and />Three kinds of>The value of (2) is smaller than 1, in order to cut down the degree of association between units, < >>In order to increase the degree of association between units, the evaluation coefficient is assigned 1 in order to maintain the degree of association between conventional units;
wherein ,the calculation formula of (2) is as follows: />，/>The calculation formula of (2) is as follows: />, wherein Is the large value in the physical field characteristic values of the medium on both sides of the same physical field lower boundary, wherein +.>Is the small value in the physical field characteristic value of the medium at both sides of the same physical field lower boundary, +.>As the correction coefficient, its value range: />；
In the configuration combination of the formulated structure full constraint evaluation coefficients, the configuration principle of the structure full constraint evaluation coefficients is as follows: for the two-dimensional model of the structured grid, defining units on the boundary of the structure as structural units, carrying out regularization constraint processing between each structural unit and 8 adjacent units, and configuring special evaluation coefficients only for regularization parameters related to the structural units; structural element evaluation coefficient assigned to boundary ipsilateralEnhancing the association in the structural direction; structural element evaluation coefficient assigned across structural boundary +.>Weakening the correlation in the structural direction; setting an evaluation coefficient value to be 1 between the structural unit and the non-structural unit, and keeping normal association;
in the configuration combination of the fully constrained evaluation coefficients of the prepared structure, four different configuration combinations of the evaluation coefficients are adopted for four types of structure boundary types:
for the horizontal boundary, as two sides of the structural boundary are respectively provided with a layer of structural unit layer, the regularized constraint processing between the structural unit and 8 adjacent units is completed by respectively applying a layer of evaluation coefficient configuration to the structural units at the upper layer and the structural units at the lower layer of the horizontal boundary according to the structural full constraint evaluation coefficient configuration principle;
for the vertical boundary, as two sides of the structural boundary are respectively provided with a layer of structural unit layer, the regularized constraint processing between the structural unit and 8 adjacent units is completed by respectively applying a layer of evaluation coefficient configuration to the structural units on the left side and the structural units on the right side of the vertical boundary according to the structural full constraint evaluation coefficient configuration principle;
for the first type inclined boundary, as two layers of structural unit layers are respectively arranged at two sides of the structural boundary, two layers of evaluation coefficient configuration are respectively applied to the upper left structural unit and the lower right structural unit of the first type inclined boundary according to the structural full-constraint evaluation coefficient configuration principle for completing regularization constraint processing between the structural unit and 8 adjacent structural units;
for the second type of inclined boundary, as two layers of structural unit layers are respectively arranged on two sides of the structural boundary, two layers of evaluation coefficient configuration are respectively applied to the upper right structural unit and the lower left structural unit of the first type of inclined boundary according to the structural full-constraint evaluation coefficient configuration principle for completing regularization constraint processing between the structural unit and 8 adjacent structural units;
the configuration calibration of the evaluation coefficient of the boundary of the special structure is specifically as follows: for a composite boundary region which does not belong to four types of conventional structure boundary types, special processing is required to be carried out on special structural units according to a specific composite boundary combination condition according to the structure full constraint evaluation coefficient configuration principle;
the method comprises the steps of fully constrained assembly and format standardization of a structure to form a regularization item containing prior structural fully beam configuration, and specifically comprises the following steps:
according to the evaluation coefficient configuration completed by the prior structure boundary, the unit grid dimension of the inversion model is calculated according to the formula:
format standard reorganization in (3) will finally +.>，/>，/>，/>And sequentially carrying the above formulas to form regularization terms containing the prior structure full-constraint configuration.
2. The method for configuring the prior boundary structure in the geophysical prospecting data regularization inversion according to claim 1, wherein the determining the structure boundary type specifically includes:
judging the type of the structural boundary according to the approximation degree of the trend of the structural boundary and the passing unit grid in the horizontal direction, the vertical direction, the first type diagonal direction and the second type diagonal direction; the structural boundary types are divided into four boundary types, namely a horizontal boundary, a vertical boundary, a first type inclined boundary and a second type inclined boundary.
3. The method for fully constrained configuration of a priori boundary structure in geophysical prospecting data regularized inversion according to claim 2, wherein the determining the type of the boundary of the structure further comprises:
according to the cell grids with different sizes, defining the included angle between the diagonal direction of the second class and the horizontal direction asThe included angle between the diagonal direction of the first type and the horizontal direction is +.>The angle between the first diagonal direction and the vertical direction is +.>The method comprises the steps of carrying out a first treatment on the surface of the The angle between the trend of the boundary of the prior structure and the vertical direction is defined as +.>；
When (when)When the prior structure boundary is a horizontal boundary;
when (when)Or->When the prior structure boundary is a vertical boundary;
when (when)When the prior structure boundary is the first type of inclined boundary;
when (when)The a priori structure boundaries are then second class tilt boundaries.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310884362.0A CN116611267B (en) | 2023-07-19 | 2023-07-19 | All-constraint configuration method for priori boundary structure in geophysical prospecting data regularized inversion |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310884362.0A CN116611267B (en) | 2023-07-19 | 2023-07-19 | All-constraint configuration method for priori boundary structure in geophysical prospecting data regularized inversion |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116611267A CN116611267A (en) | 2023-08-18 |
| CN116611267B true CN116611267B (en) | 2023-09-19 |
Family
ID=87685690
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310884362.0A Active CN116611267B (en) | 2023-07-19 | 2023-07-19 | All-constraint configuration method for priori boundary structure in geophysical prospecting data regularized inversion |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116611267B (en) |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112434439A (en) * | 2020-12-01 | 2021-03-02 | 中国自然资源航空物探遥感中心 | Multi-segment inversion method and device for constraint physical property multi-interval distribution |
| CN115469376A (en) * | 2022-11-14 | 2022-12-13 | 四川省冶勘设计集团有限公司 | Tunnel advanced geological forecast comprehensive detection method and system based on non-homologous data |
| CN115586573A (en) * | 2022-09-15 | 2023-01-10 | 河海大学 | Dynamic constraint physical property parameter seismic inversion method for tight sandstone reservoir |
| CN116050046A (en) * | 2022-09-21 | 2023-05-02 | 海南经贸职业技术学院 | Magnetotelluric fuzzy constraint inversion method based on cluster analysis |
| CN116088044A (en) * | 2023-01-17 | 2023-05-09 | 中海石油(中国)有限公司海南分公司 | Deep water shallow layer seismic data structure constraint attenuation compensation speed modeling method |
| CN116165722A (en) * | 2023-03-29 | 2023-05-26 | 中国地质大学（武汉） | Loop source transient electromagnetic three-dimensional rapid inversion method adopting Gaussian Newton method |
Family Cites Families (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2539744A4 (en) * | 2010-03-19 | 2017-11-22 | Schlumberger Technology B.V. | Uncertainty estimation for large-scale nonlinear inverse problems using geometric sampling and covariance-free model compression |
| US10591638B2 (en) * | 2013-03-06 | 2020-03-17 | Exxonmobil Upstream Research Company | Inversion of geophysical data on computer system having parallel processors |
-
2023
- 2023-07-19 CN CN202310884362.0A patent/CN116611267B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112434439A (en) * | 2020-12-01 | 2021-03-02 | 中国自然资源航空物探遥感中心 | Multi-segment inversion method and device for constraint physical property multi-interval distribution |
| CN115586573A (en) * | 2022-09-15 | 2023-01-10 | 河海大学 | Dynamic constraint physical property parameter seismic inversion method for tight sandstone reservoir |
| CN116050046A (en) * | 2022-09-21 | 2023-05-02 | 海南经贸职业技术学院 | Magnetotelluric fuzzy constraint inversion method based on cluster analysis |
| CN115469376A (en) * | 2022-11-14 | 2022-12-13 | 四川省冶勘设计集团有限公司 | Tunnel advanced geological forecast comprehensive detection method and system based on non-homologous data |
| CN116088044A (en) * | 2023-01-17 | 2023-05-09 | 中海石油(中国)有限公司海南分公司 | Deep water shallow layer seismic data structure constraint attenuation compensation speed modeling method |
| CN116165722A (en) * | 2023-03-29 | 2023-05-26 | 中国地质大学（武汉） | Loop source transient electromagnetic three-dimensional rapid inversion method adopting Gaussian Newton method |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116611267A (en) | 2023-08-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112990330B (en) | User energy abnormal data detection method and device | |
| CN112884320B (en) | Foundation pit risk assessment method, device, equipment and medium based on entropy model | |
| CN116150676A (en) | Equipment fault diagnosis and identification method and device based on artificial intelligence | |
| CN115033591B (en) | Intelligent detection method, system, storage medium and computer equipment for electric charge data abnormality | |
| CN112766537A (en) | Short-term electric load prediction method | |
| CN112949735A (en) | Liquid hazardous chemical substance volatile concentration abnormity discovery method based on outlier data mining | |
| CN116611267B (en) | All-constraint configuration method for priori boundary structure in geophysical prospecting data regularized inversion | |
| CN114639102A (en) | Cell segmentation method and device based on key point and size regression | |
| CN117335409B (en) | Power consumer load prediction system based on artificial intelligence | |
| CN116365519B (en) | Power load prediction method, system, storage medium and equipment | |
| CN117710824A (en) | Core image RQD intelligent analysis method based on improved Cascade Mask R-CNN model | |
| CN117077513A (en) | Noise analysis method for solving neutron diffusion eigenvalue problem based on DEPINN | |
| CN117422195A (en) | Water quality evaluation method, device, computer equipment and storage medium | |
| CN111863124B (en) | Copy number variation detection method, system, storage medium and computer equipment | |
| CN115616408A (en) | Battery thermal management data processing method and system | |
| CN115953578A (en) | Method and device for evaluating inference result of fracture semantic segmentation model | |
| CN109840479B (en) | Health state matching method and device | |
| CN118432952B (en) | Abnormality detection method under zero trust environment, electronic equipment and storage medium | |
| CN111459838A (en) | Software defect prediction method and system based on manifold alignment | |
| CN112257586B (en) | Truth box selection method, device, storage medium and equipment in target detection | |
| CN118583347B (en) | Stress detection method of organic electroluminescent device | |
| CN111626586B (en) | Data quality detection method, device, computer equipment and storage medium | |
| CN114519168A (en) | Spectrum resolving method of multi-sphere neutron spectrometer based on parameterized energy spectrum | |
| CN114677000A (en) | Sandstone thermal reservoir development and utilization suitability quantification method and device and storage medium | |
| CN117455893A (en) | GIS defect classification and assessment method and system based on ultraviolet image |
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 |