CN111309295A - Inversion visualization method, device and readable storage medium for three-dimensional variable density constraint underground interface - Google Patents
Inversion visualization method, device and readable storage medium for three-dimensional variable density constraint underground interface Download PDFInfo
- Publication number
- CN111309295A CN111309295A CN202010172103.1A CN202010172103A CN111309295A CN 111309295 A CN111309295 A CN 111309295A CN 202010172103 A CN202010172103 A CN 202010172103A CN 111309295 A CN111309295 A CN 111309295A
- Authority
- CN
- China
- Prior art keywords
- variable density
- inversion
- visualization
- operation port
- interface
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F8/00—Arrangements for software engineering
- G06F8/20—Software design
- G06F8/24—Object-oriented
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
Abstract
The invention provides a three-dimensional variable density constraint underground interface inversion visualization method, equipment and a readable storage medium, wherein a calculation formula for obtaining variable density parameters is configured; configuring a formula into a Matlab numerical simulation calculation module; realizing visual packaging of variable density parameters by using a Matlab GUI component; and configuring an inversion visualization operation port. And a reusable calculation mode can be realized by the three-dimensional variable density constraint underground interface inversion visualization method. And packaging the interfaces of the function codes by using a Matlab GUI component to form an interface, so as to form a mode convenient for operation, and facilitate the use of research and development personnel. The problems that the cost for carrying out large-area exploration by utilizing seismic reflection waves is too high, effective exploration cannot be carried out, and the exploration process has limitation are solved.
Description
Technical Field
The invention relates to the field of three-dimensional data modeling visualization, in particular to a three-dimensional variable density constraint underground interface inversion visualization method, equipment and a readable storage medium.
Background
The inversion research of the underground density interface has important significance for disclosing underground geological structures, for example, the underground important density interface-Mohuo surface which is the interface between the crust and the mantle can be obtained through the research on the depth of the Mohuo surface, the thickness change and the distribution condition of the crust in the research area range can be obtained, and the inversion research of the underground density interface has important significance for the structure of the mantle and the deep dynamics.
At present, the three-dimensional fluctuation form of a mujojoba surface is obtained by utilizing the reflection of an earthquake, but the cost for large-area exploration by utilizing the reflection of the earthquake is too high, effective exploration cannot be carried out, and the exploration process has limitation.
Disclosure of Invention
In order to overcome the defects in the prior art, the invention provides a three-dimensional variable density constraint underground interface inversion visualization method, which comprises the following steps:
step one, configuring the following calculation formula for obtaining the variable density parameter:
step two, the formula is configured in a Matlab numerical simulation calculation module;
thirdly, realizing the visual packaging of the variable density parameters by using a Matlab GUI component;
and step four, configuring an inversion visualization operation port.
It should be further noted that the second step further includes:
step 11: configuring an operating environment for realizing a Matlab numerical simulation calculation module;
step 12: finishing the functional code compiling of the variable density parameter calculation formula based on a Matlab numerical simulation calculation module, and forming an operation file of the variable density parameter calculation formula;
step 13: configuring a Matlab GUI component under the current operating environment;
step 14: compiling a visual packaging interface by using a Matlab GUI component, and forming a visual file;
step 15: configuring a variable density parameter calculation formula operation file and a visual file into a preset folder;
step 16: and configuring an inversion visualization operation port, so that a user can calculate the variable density parameters through the inversion visualization operation port.
The invention also provides equipment for realizing the inversion visualization method of the three-dimensional variable density constraint underground interface, which comprises the following steps: the memory is used for storing a computer program and a three-dimensional variable density constraint underground interface inversion visualization method; and the processor is used for executing the computer program and the three-dimensional variable density constraint underground interface inversion visualization method so as to realize the steps of the three-dimensional variable density constraint underground interface inversion visualization method.
The invention also provides a readable storage medium with a three-dimensional variable density constraint underground interface inversion visualization method, wherein a computer program is stored on the readable storage medium, and the computer program is executed by a processor to realize the steps of the three-dimensional variable density constraint underground interface inversion visualization method.
According to the technical scheme, the invention has the following advantages:
the invention improves Parker-Oldenburg, and perfects the theory; the invention also completes the code realization of the Parker-Oldenburg theory and designs a visualization system for the Parker-Oldenburg theory.
And a reusable calculation mode can be realized through the three-dimensional variable density constraint underground interface inversion visualization method. And packaging the interfaces of the function codes by using a Matlab GUI component to form an interface, so as to form a mode convenient for operation, and facilitate the use of research and development personnel. The problems that the cost for carrying out large-area exploration by utilizing seismic reflection waves is too high, effective exploration cannot be carried out, and the exploration process has limitation are solved.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings used in the description will be briefly introduced, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without creative efforts.
FIG. 1 is a flow chart of a three-dimensional variable density constraint subsurface interface inversion visualization method;
fig. 2 is a schematic diagram of an embodiment of Matlab GUI component in the present invention.
Detailed Description
Those of ordinary skill in the art will appreciate that the elements and algorithm steps of the examples described in connection with the embodiments disclosed herein may be embodied in electronic hardware, computer software, or combinations of both, and that the components and steps of the examples have been described in a functional general in the foregoing description for the purpose of illustrating clearly the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
The block diagrams shown in the figures are functional entities only and do not necessarily correspond to physically separate entities. I.e. these functional entities may be implemented in the form of software, or in one or more hardware modules or integrated circuits, or in different networks and/or processor means and/or microcontroller means.
In the several embodiments provided in the present application, it should be understood that the disclosed system, apparatus and method may be implemented in other ways. For example, the above-described apparatus embodiments are merely illustrative, and for example, the division of the units is only one logical division, and other divisions may be realized in practice, for example, a plurality of units or components may be combined or integrated into another system, or some features may be omitted, or not executed. In addition, the shown or discussed mutual coupling or direct coupling or communication connection may be an indirect coupling or communication connection through some interfaces, devices or units, and may also be an electric, mechanical or other form of connection.
The invention provides a three-dimensional variable density constraint underground interface inversion visualization method, as shown in figure 1, the method comprises the following steps:
s11, configuring the following calculation formula for obtaining the variable density parameter:
abnormal F [ 2 ] by gravity]For Fourier transform, σ0(ξ) is the density difference of the transverse two-dimensional change, ξ and η are constant values, z is0Is the average depth of the interface, Δ h is the average depth of the interface relative to z0Delta G is residual gravity anomaly, G is a universal gravitation constant, k is a frequency domain circular wave number, and mu is an attenuation coefficient of residual density along with depth;
s12, configuring the formula into a Matlab numerical simulation calculation module;
the meaning of the Matlab numerical simulation computation module related by the invention is the combination of two words of matrix & laboratory, which means a matrix factory. Is a high-tech computing environment facing scientific computing, visualization, and interactive programming. The method integrates a plurality of powerful functions of data analysis, matrix calculation, scientific data visualization, modeling and simulation of a nonlinear dynamic system and the like into an easy-to-use window environment, provides a comprehensive solution for scientific research, engineering design and a plurality of scientific fields which need to carry out effective numerical calculation, and gets rid of the editing mode of the traditional non-interactive programming language to a great extent.
The MATLAB can perform matrix operation, draw functions and data, realize algorithms, create user interfaces, connect programs of other programming languages and the like, and is mainly applied to the fields of engineering calculation, control design, signal processing and communication, image processing, signal detection, financial modeling design and analysis and the like.
The basic data unit of MATLAB is matrix, its instruction expression is very similar to the form commonly used in mathematics and engineering, so that it is much simpler and more convenient to use MATLAB to solve the problem than to use C, FORTRAN and other languages to implement the same things, and the MATLAB also absorbs the advantages of software like Maple and so on, and makes MATLAB become a powerful mathematical software. Support for C, FORTRAN, C + +, JAVA was also added to the new version.
Namely, the variable density parameter calculation formula is configured in the Matlab numerical simulation calculation module to realize automatic calculation and obtain the relevant data.
The specific arrangement process is as follows:
step 11: configuring an operating environment for realizing a Matlab numerical simulation calculation module;
step 12: finishing the functional code compiling of the variable density parameter calculation formula based on a Matlab numerical simulation calculation module, and forming an operation file of the variable density parameter calculation formula;
the formed variable density parameter calculation formula operation file is an m file. The m-file is a file format, the Chinese name is an executable code file, and the file extension name is ". m". Belonging to a package in the Wolfram language.
Step 13: configuring a Matlab GUI component under the current operating environment;
step 14: compiling a visual packaging interface by using a Matlab GUI component, and forming a visual file;
step 15: configuring a variable density parameter calculation formula operation file and a visual file into a preset folder;
step 16: and configuring an inversion visualization operation port, so that a user can calculate the variable density parameters through the inversion visualization operation port.
S13, realizing visual packaging of variable density parameters by using a Matlab GUI component;
and S14, configuring an inversion visualization operation port.
The inversion visualization operation port provided by the invention comprises: a parameter operation port and a calculation operation port. As shown in particular in figure 2 of the drawings,
the parameter operation port includes: the device comprises an average depth input port, an iteration number input port, a filtering frequency low value input port, a filtering frequency high value input port, an X-direction wavelength input port, a Y-direction wavelength input port, a residual density constant input port and a residual density e exponential change input port. The calculation operation port comprises an opening operation port, a closing operation port and an output operation port.
The user can input corresponding data into the port based on the current requirement, and the starting, closing and output functions of the method are realized by calculating the starting operation port, the closing operation port and the output operation port of the operation port.
Based on the method, the invention also provides equipment for realizing the inversion visualization method of the three-dimensional variable density constraint underground interface, which comprises the following steps: the memory is used for storing a computer program and a three-dimensional variable density constraint underground interface inversion visualization method; and the processor is used for executing the computer program and the three-dimensional variable density constraint underground interface inversion visualization method so as to realize the steps of the three-dimensional variable density constraint underground interface inversion visualization method.
Based on the method, the invention also provides a readable storage medium with the three-dimensional variable density constraint underground interface inversion visualization method, and the readable storage medium stores a computer program which is executed by a processor to realize the steps of the three-dimensional variable density constraint underground interface inversion visualization method.
The apparatus implementing the method for three-dimensional variable density constrained subsurface inversion visualization is the elements and algorithm steps of the examples described in connection with the embodiments disclosed herein, which can be implemented in electronic hardware, computer software, or combinations thereof, the components and steps of the examples having been generally described in the foregoing description by function for the purpose of clearly illustrating the interchangeability of hardware and software. Whether such functionality is implemented as hardware or software depends upon the particular application and design constraints imposed on the implementation. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present invention.
Through the above description of the embodiments, those skilled in the art will readily understand that the apparatus for implementing the three-dimensional variable density constraint subsurface interface inversion visualization method described herein can be implemented by software, or by software in combination with necessary hardware. Therefore, the technical solution according to the disclosed embodiments of the method for implementing inversion visualization of a three-dimensional variable density constrained subsurface interface may be embodied in the form of a software product, which may be stored in a non-volatile storage medium (which may be a CD-ROM, a usb disk, a mobile hard disk, etc.) or on a network, and includes several instructions to enable a computing device (which may be a personal computer, a server, a mobile terminal, or a network device, etc.) to execute the indexing method according to the disclosed embodiments.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (8)
1. A three-dimensional variable density constraint underground interface inversion visualization method is characterized by comprising the following steps:
step one, configuring the following calculation formula for obtaining the variable density parameter:
abnormal F [ 2 ] by gravity]For Fourier transform, σ0(ξ) is the density difference of the transverse two-dimensional change, ξ and η are constant values, z is0Is the average depth of the interface, Δ h is the average depth of the interface relative to z0Delta G is residual gravity anomaly, G is a universal gravitation constant, k is a frequency domain circular wave number, and mu is an attenuation coefficient of residual density along with depth;
step two, the formula is configured in a Matlab numerical simulation calculation module;
thirdly, realizing the visual packaging of the variable density parameters by using a Matlab GUI component;
and step four, configuring an inversion visualization operation port.
2. The method of claim 1,
the second step further comprises:
step 11: configuring an operating environment for realizing a Matlab numerical simulation calculation module;
step 12: finishing the functional code compiling of the variable density parameter calculation formula based on a Matlab numerical simulation calculation module, and forming an operation file of the variable density parameter calculation formula;
step 13: configuring a Matlab GUI component under the current operating environment;
step 14: compiling a visual packaging interface by using a Matlab GUI component, and forming a visual file;
step 15: configuring a variable density parameter calculation formula operation file and a visual file into a preset folder;
step 16: and configuring an inversion visualization operation port, so that a user can calculate the variable density parameters through the inversion visualization operation port.
3. The method of claim 2,
in Step12, the operating file of the variable density parameter calculation formula is formed as an m file.
4. The method according to claim 1 or 2,
the inversion visualization operation port comprises: a parameter operation port and a calculation operation port.
5. The method of claim 4,
the parameter operation port includes: the device comprises an average depth input port, an iteration number input port, a filtering frequency low value input port, a filtering frequency high value input port, an X-direction wavelength input port, a Y-direction wavelength input port, a residual density constant input port and a residual density e exponential change input port.
6. The method of claim 4,
the calculation operation port comprises an opening operation port, a closing operation port and an output operation port.
7. An apparatus for implementing a visualization method for inversion of a three-dimensional variable density constrained subsurface interface, comprising:
the memory is used for storing a computer program and a three-dimensional variable density constraint underground interface inversion visualization method;
a processor for executing the computer program and the three-dimensional variable density constrained subsurface inversion visualization method to realize the steps of the three-dimensional variable density constrained subsurface inversion visualization method according to any one of claims 1 to 6.
8. A readable storage medium having a three-dimensional variable density constrained subsurface inversion visualization method, wherein the readable storage medium has stored thereon a computer program, which is executed by a processor to implement the steps of the three-dimensional variable density constrained subsurface inversion visualization method according to any one of claims 1 to 6.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010172103.1A CN111309295B (en) | 2020-03-12 | 2020-03-12 | Three-dimensional variable density constraint underground interface inversion visualization method, equipment and readable storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202010172103.1A CN111309295B (en) | 2020-03-12 | 2020-03-12 | Three-dimensional variable density constraint underground interface inversion visualization method, equipment and readable storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111309295A true CN111309295A (en) | 2020-06-19 |
| CN111309295B CN111309295B (en) | 2023-04-18 |
Family
ID=71158568
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202010172103.1A Active CN111309295B (en) | 2020-03-12 | 2020-03-12 | Three-dimensional variable density constraint underground interface inversion visualization method, equipment and readable storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111309295B (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101051394A (en) * | 2007-04-11 | 2007-10-10 | 中国科学院地质与地球物理研究所 | Three dimension visual system based on geological body of geophysical field data |
| US20130265851A1 (en) * | 2012-02-06 | 2013-10-10 | Kees Faber | Sensor System of Buried Seismic Array |
| CN103576212A (en) * | 2012-07-19 | 2014-02-12 | 中国石油天然气集团公司 | Complex structural well constrained three dimensional density sequence inversion method |
| CN104297814A (en) * | 2014-10-24 | 2015-01-21 | 张维加 | Machine perspective vision implementation method |
| US20170132339A1 (en) * | 2015-08-20 | 2017-05-11 | Fracgeo Llc | System for predicting induced seismicity potential resulting from injection of fluids in naturally fractured reservoirs |
| CN107748399A (en) * | 2017-09-12 | 2018-03-02 | 中国石油化工股份有限公司 | Mountain front infrastructure layer method is identified using Interface Inversion of Gravity Field |
| CN109543353A (en) * | 2018-12-29 | 2019-03-29 | 广东电网有限责任公司 | Three-dimensional water vapor retrieval method, apparatus, equipment and computer readable storage medium |
| CN109976729A (en) * | 2019-05-05 | 2019-07-05 | 东北大学 | One kind depositing calculation and shows globally configurable Data Analysis Software architecture design method |
-
2020
- 2020-03-12 CN CN202010172103.1A patent/CN111309295B/en active Active
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101051394A (en) * | 2007-04-11 | 2007-10-10 | 中国科学院地质与地球物理研究所 | Three dimension visual system based on geological body of geophysical field data |
| US20130265851A1 (en) * | 2012-02-06 | 2013-10-10 | Kees Faber | Sensor System of Buried Seismic Array |
| CN103576212A (en) * | 2012-07-19 | 2014-02-12 | 中国石油天然气集团公司 | Complex structural well constrained three dimensional density sequence inversion method |
| CN104297814A (en) * | 2014-10-24 | 2015-01-21 | 张维加 | Machine perspective vision implementation method |
| US20170132339A1 (en) * | 2015-08-20 | 2017-05-11 | Fracgeo Llc | System for predicting induced seismicity potential resulting from injection of fluids in naturally fractured reservoirs |
| CN107748399A (en) * | 2017-09-12 | 2018-03-02 | 中国石油化工股份有限公司 | Mountain front infrastructure layer method is identified using Interface Inversion of Gravity Field |
| CN109543353A (en) * | 2018-12-29 | 2019-03-29 | 广东电网有限责任公司 | Three-dimensional water vapor retrieval method, apparatus, equipment and computer readable storage medium |
| CN109976729A (en) * | 2019-05-05 | 2019-07-05 | 东北大学 | One kind depositing calculation and shows globally configurable Data Analysis Software architecture design method |
Non-Patent Citations (3)
| Title |
|---|
| 张盛等: "约束变密度界面反演方法", 《地球物理学进展》 * |
| 袁东;雷玉山;张涛;黄磊;: "基于Voxler的高密度电法数据三维可视化的实现" * |
| 陈康;李世峰;张政杰;张晋磊;李晓波;: "高密度电法勘探数据三维处理技术及城市勘察中的应用" * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111309295B (en) | 2023-04-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Magoules et al. | Asynchronous optimized Schwarz methods with and without overlap | |
| Momani et al. | Numerical investigation for Caputo-Fabrizio fractional Riccati and Bernoulli equations using iterative reproducing kernel method | |
| Tanaka | Verification of Hasselmann's energy transfer among surface gravity waves by direct numerical simulations of primitive equations | |
| Aagaard et al. | A domain decomposition approach to implementing fault slip in finite‐element models of quasi‐static and dynamic crustal deformation | |
| US10691410B2 (en) | Neural network computing | |
| Petersson et al. | Model reduction using a frequency-limited H2-cost | |
| Ganis et al. | Stochastic collocation and mixed finite elements for flow in porous media | |
| El Bouajaji et al. | A quasi-optimal domain decomposition algorithm for the time-harmonic Maxwell's equations | |
| Xu et al. | Conditional simulation of nonstationary fluctuating wind speeds for long-span bridges | |
| Lemarié et al. | Toward an optimized global-in-time schwarz algorithm for diffusion equations with discontinuous and spatially variable coefficients, part 2: the variable coefficients case | |
| Park | Guaranteed-stable sliding DFT algorithm with minimal computational requirements | |
| Antoine et al. | An analysis of Schwarz waveform relaxation domain decomposition methods for the imaginary-time linear Schrödinger and Gross-Pitaevskii equations | |
| Krishnaiah | Developments in statistics | |
| Du et al. | Isotropic variogram matrix functions on spheres | |
| CN115685334A (en) | Anisotropic elastic wave field decomposition method and device and computer equipment | |
| Wang et al. | Passivity enforcement for descriptor systems via matrix pencil perturbation | |
| Yahia et al. | Description of turbulent dynamics in the interstellar medium: multifractal-microcanonical analysis-I. Application to Herschel observations of the musca filament | |
| Régaldo-Saint Blancard et al. | Generative Models of Multichannel Data from a Single Example—Application to Dust Emission | |
| Li et al. | Numerical solutions of coupled systems of fractional order partial differential equations | |
| CN111309295B (en) | Three-dimensional variable density constraint underground interface inversion visualization method, equipment and readable storage medium | |
| He et al. | A numerical dispersion-dissipation analysis of discontinuous Galerkin methods based on quadrilateral and triangular elements | |
| Han et al. | Classification of all single traveling wave solutions of fractional coupled Boussinesq equations via the complete discrimination system method | |
| US11914671B2 (en) | Performing uncertainty quantification analysis with efficient two dimensional random fields | |
| CN112764105B (en) | HTI medium quasi-longitudinal wave forward simulation method and device, storage medium and processor | |
| Garcia-Donoro et al. | Higher-order finite element electromagnetics code for HPC environments |
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 | ||
| CB02 | Change of applicant information |
Address after: 250101 no.2877 Kehang Road, Suncun Town, high tech Zone, Jinan City, Shandong Province Applicant after: Chaoyue Technology Co.,Ltd. Address before: 250101 no.2877 Kehang Road, Suncun Town, high tech Zone, Jinan City, Shandong Province Applicant before: SHANDONG CHAOYUE DATA CONTROL ELECTRONICS Co.,Ltd. |
|
| CB02 | Change of applicant information | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |