CN112785910A - Large-dynamic-range nonlinear geophysical contour map drawing method and device - Google Patents
Large-dynamic-range nonlinear geophysical contour map drawing method and device Download PDFInfo
- Publication number
- CN112785910A CN112785910A CN201911079926.3A CN201911079926A CN112785910A CN 112785910 A CN112785910 A CN 112785910A CN 201911079926 A CN201911079926 A CN 201911079926A CN 112785910 A CN112785910 A CN 112785910A
- Authority
- CN
- China
- Prior art keywords
- geophysical data
- spatial distribution
- type
- geophysical
- gradient
- 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
- G09—EDUCATION; CRYPTOGRAPHY; DISPLAY; ADVERTISING; SEALS
- G09B—EDUCATIONAL OR DEMONSTRATION APPLIANCES; APPLIANCES FOR TEACHING, OR COMMUNICATING WITH, THE BLIND, DEAF OR MUTE; MODELS; PLANETARIA; GLOBES; MAPS; DIAGRAMS
- G09B29/00—Maps; Plans; Charts; Diagrams, e.g. route diagram
- G09B29/003—Maps
- G09B29/006—Representation of non-cartographic information on maps, e.g. population distribution, wind direction, radiation levels, air and sea routes
- G09B29/007—Representation of non-cartographic information on maps, e.g. population distribution, wind direction, radiation levels, air and sea routes using computer methods
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Remote Sensing (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geology (AREA)
- Theoretical Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Acoustics & Sound (AREA)
- Ecology (AREA)
- Educational Technology (AREA)
- Educational Administration (AREA)
- Business, Economics & Management (AREA)
- Mathematical Physics (AREA)
- General Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a method and a device for drawing a nonlinear geophysical contour map with a large dynamic range, wherein the method comprises the following steps: acquiring geophysical data; calculating a gradient non-uniformity of the geophysical data; determining a spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution; determining a contour interval sequence of the geophysical data according to the spatial distribution type; and drawing a geophysical contour map according to the contour interval sequence of the geophysical data. The method avoids the problem that the distance between adjacent contour lines in the contour map is abnormally dense or sparse when the geophysical data change very severely, and effectively retains the detail characteristics of the geophysical data on the basis of clearly expressing the overall form of the geophysical data.
Description
Technical Field
The invention relates to the technical field of geophysical exploration, in particular to a method and a device for drawing a large-dynamic-range nonlinear geophysical contour map.
Background
Geophysical prospecting (geophysical prospecting for short) refers to a technique for detecting the characteristics of geological structures, lithology, fluids and the like of the stratum by researching and observing the changes of various geophysical fields. Generally, the different formation media that make up the earth's crust differ in density, elasticity, electrical conductivity, magnetism, radioactivity, and thermal conductivity, which cause corresponding local changes in the geophysical field. The geological character can be inferred by measuring the distribution and change characteristics of the geophysical field and analyzing and researching by combining with known geological data.
The basis of geophysical exploration is the difference in physical properties, density, magnetization properties, conductivity and radioactivity between rocks in a stratum and surrounding rocks, and on the basis of the difference, the exploration and the research are carried out on the medium structures, substance compositions, formation and evolution of an earth body and a near-earth space, various related natural phenomena and change rules thereof are sought, and theories, methods and technologies are provided for detecting the internal structure and structure of the earth, seeking energy sources, resources and environment monitoring or important basis is provided for geological disaster prediction.
The geophysical exploration methods commonly used at present comprise gravity exploration, magnetic exploration (comprising a geoelectromagnetic method, a time-frequency electromagnetic method and the like), electrical exploration (comprising an electrical sounding method, an electrical section method, a high-density electrical method, a natural electric field method, a charging method, an induced polarization method, a controllable source audio magnetotelluric sounding method, a transient electromagnetic method and the like), seismic exploration (comprising a refracted wave method, a reflected wave method, a Rayleigh wave method and the like), geochemical exploration and the like.
Geophysical exploration is an indirect observation method, which utilizes the physical principle and specific instruments to obtain the physical parameters and the rules of a known rock ore specimen or geological model, and then analyzes, explains and extracts the parameter values which can reflect the space morphology of a geologic body, the rock physical properties or the mineral characteristics and are contained in geophysical observation information according to the established physical rules (mathematical physical models) to form a geophysical result and a geological exploration result. In geophysical exploration, the main expression form of a geophysical observation result is geophysical data, and the geophysical data needs to be drawn in a graphic form for observation, analysis and deep research, namely the field of geophysical graphic drawing and industrial mapping.
The geophysical data is a generic term for geophysical observation information collected by a geophysical prospecting instrument and for geological result information extracted by analyzing and interpreting the geophysical observation information. Due to the numerous geophysical prospecting methods and types, geophysical data are complex in type and diverse in form. According to the geophysical exploration method, the geophysical data can be divided into seismic, gravity, magnetic force, electromagnetic and other data; geophysical prospecting can be divided into one, two and three dimensions depending on their dimensionality. The geophysical map is a map for expressing geophysical response and geologic morphological characteristics contained in geophysical data. Corresponding to the geophysical data, the geophysical map can be divided into seismic response map, gravity response map, magnetic response map, electromagnetic response map and other types according to the geophysical response type; according to the dimensionality of geophysical exploration, the geophysical graph can be divided into one-dimensional, two-dimensional, three-dimensional and other types; according to the graph form, the geophysical graph can be divided into a section view, a plan view, a three-dimensional visualization graph and the like; the geophysical figures may be classified into variable density map, contour map, bubble map, spider-web map, and the like, according to the form of graphical representation. In geophysical exploration, the geophysical map with the most important position and the most widely applied position belongs to a geophysical contour map.
A geophysical contour map is a graph in which projections of connecting lines of equal numerical points on a horizontal plane represent continuously distributed and gradually changing data features, which is known for the fact that the adopted data representation form is a contour line. Geophysical contour maps are commonly used to represent features such as elevation of the earth's surface, formation space morphology, earth isotherms, earth barometric pressure, and the like. In the field of geophysical exploration, contour maps can be divided into terrain contour maps, isotherm maps, isobologram maps, stratigraphic structure maps, stratigraphic isobath maps, mineral composition contour maps, mineral reserve contour maps, and the like, according to the difference in the geophysical and geological characteristics that it expresses.
The geophysical contour map is suitable for representing the phenomena of geography, physics, chemistry, even humanity and the like which are continuously distributed and have quantitative characteristics (height, size, strength, speed and the like), and the spatial change rule of the geophysical data can be intuitively obtained by utilizing the contour map. For example, on an isotherm diagram, the dense isotherms indicate that the temperature difference of the region is large, and the sparse isotherms indicate that the temperature difference of the region is small; on the contour map, topographic features and slope slopes can be judged according to density degree and form of contour lines, wherein the density of the contour lines indicates high slope steepness, and the sparsity of the contour lines indicates low slope slowness. The contour is typically drawn in the form of a thin line segment plus a numerical mark. On contour maps, in addition to numerically marking the values represented by the contours, colors are often used to enhance intuitiveness and highlight differences between quantities. If the color with different shades or the density of the halo lines is drawn between the contour lines, and the shade of the color or the density of the halo lines corresponds to the numerical value of the data, the law and the regional difference of the quantity change can be reflected more obviously, for example, on the global isotherm diagram, the warm degree is represented by red, and the cold degree is represented by blue.
In the prior art, when the geophysical data change very violently, the distance between adjacent contour lines in a drawn contour map is abnormally dense or sparse, so that the contour lines are overlapped together and are difficult to distinguish, or the contour lines are too sparse and are difficult to obtain the macroscopic characteristics of the corresponding geophysical data, and the obvious detailed characteristics of the geophysical data are difficult to keep on the basis of clearly expressing the integral form of the geophysical data.
Disclosure of Invention
The embodiment of the invention provides a large-dynamic-range nonlinear geophysical contour map drawing method, which is used for drawing a geophysical contour map and reserving obvious detail characteristics of geophysical data on the basis of clearly expressing the integral form of the geophysical data, and the method comprises the following steps:
acquiring geophysical data;
calculating a gradient non-uniformity of the geophysical data;
determining a spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution;
determining a contour interval sequence of the geophysical data according to the spatial distribution type;
and drawing a geophysical contour map according to the contour interval sequence of the geophysical data.
The embodiment of the invention provides a large dynamic range nonlinear geophysical contour map drawing device, which is used for drawing a geophysical contour map and reserving obvious detail characteristics of geophysical data on the basis of clearly expressing the integral form of the geophysical data, and comprises the following components:
the acquisition module is used for acquiring geophysical data;
a calculation module for calculating gradient non-uniformity of the geophysical data;
a type determination module for determining a spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution;
a sequence determination module for determining a contour interval sequence of the geophysical data according to the spatial distribution type;
and the drawing module is used for drawing the geophysical contour map according to the contour line interval sequence of the geophysical data.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the geophysical contour map drawing method when executing the computer program.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program which is stored on the memory and can run on the processor, wherein the processor realizes the geophysical contour map drawing method when executing the computer program.
The embodiment of the invention obtains the geophysical data; calculating a gradient non-uniformity of the geophysical data; determining a spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution; determining a contour interval sequence of the geophysical data according to the spatial distribution type; and drawing a geophysical contour map according to the contour interval sequence of the geophysical data. The embodiment of the invention calculates the gradient non-uniformity of the geophysical data, determines the spatial distribution type of the geophysical data, further determines the contour line interval sequence, sets different contour line interval sequences for the geophysical data of each spatial distribution type, and adjusts the distance between adjacent contour lines in the contour map, thereby avoiding the problem that the distance between adjacent contour lines in the contour map is abnormally dense or sparse when the geophysical data are changed severely, and effectively retaining the obvious detailed characteristics of the geophysical data on the basis of clearly expressing the integral form of the geophysical data.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts. In the drawings:
FIG. 1 is a schematic diagram of a large dynamic range non-linear geophysical contour map rendering method in an embodiment of the present invention;
FIG. 2 is a block diagram of a large dynamic range non-linear geophysical contour map drawing apparatus in an embodiment of the present invention;
FIG. 3 is a diagram illustrating gravity grid anomaly data in an embodiment of the present invention;
FIG. 4 is a table with the contour spacing sequence as a constant sequence in an embodiment of the present invention;
FIG. 5 is a geophysical contour map plotted against a sequence of constants in an embodiment of the present invention;
FIG. 6 is a table in which the contour interval sequence is a linear sequence in an embodiment of the present invention;
FIG. 7 is a geophysical contour map plotted according to a linear-type sequence in accordance with an embodiment of the present invention;
FIG. 8 is a table in which the contour interval sequence is a logarithmic sequence in an embodiment of the present invention;
FIG. 9 is a geophysical contour map plotted against a logarithmic sequence in accordance with an embodiment of the present invention;
FIG. 10 is a table of contour spacing sequences as free type sequences in an embodiment of the present invention;
FIG. 11 is a geophysical contour map plotted against a freeform sequence in an embodiment of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the embodiments of the present invention are further described in detail below with reference to the 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.
In order to draw a geophysical contour map and retain obvious detailed characteristics of geophysical data on the basis of clearly expressing the overall form of the geophysical data, an embodiment of the invention provides a large-dynamic-range nonlinear geophysical contour map drawing method, as shown in fig. 1, which may include:
102, calculating gradient non-uniformity of the geophysical data;
103, determining the spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution;
104, determining a contour line interval sequence of the geophysical data according to the spatial distribution type;
and 105, drawing a geophysical contour map according to the contour line interval sequence of the geophysical data.
As can be appreciated from FIG. 1, embodiments of the present invention operate by acquiring geophysical data; calculating a gradient non-uniformity of the geophysical data; determining a spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution; determining a contour interval sequence of the geophysical data according to the spatial distribution type; and drawing a geophysical contour map according to the contour interval sequence of the geophysical data. The embodiment of the invention calculates the gradient non-uniformity of the geophysical data, determines the spatial distribution type of the geophysical data, further determines the contour line interval sequence, sets different contour line interval sequences for the geophysical data of each spatial distribution type, and adjusts the distance between adjacent contour lines in the contour map, thereby avoiding the problem that the distance between adjacent contour lines in the contour map is abnormally dense or sparse when the geophysical data are changed severely, and effectively retaining the detail characteristics of the geophysical data on the basis of clearly expressing the integral form of the geophysical data.
In specific implementation, geophysical data are acquired.
The geophysical data is a general term for geophysical observation information collected by a geophysical prospecting instrument and extracted geological result information obtained by analyzing and interpreting the geophysical observation information. Due to the wide variety of geophysical prospecting methods and types, geophysical data are complex in type and diverse in form. According to the geophysical exploration method, the geophysical data can be divided into seismic, gravity, magnetic force, electromagnetic and other data; geophysical prospecting can be divided into one, two and three dimensions depending on their dimensionality. For two-dimensional geophysical data, a two-dimensional surface in space is typically represented, consisting of X, Y, Z three components, where X and Y are absolute or relative geographic coordinates, and Z is a feature value representing the morphology of the two-dimensional surface. The geophysical data described in the embodiments of the present invention represent only characteristic values of a two-dimensional surface form.
In particular implementations, gradient non-uniformities of the geophysical data are calculated.
In an embodiment, after acquiring geophysical data, performing normalization processing on the geophysical data; the calculating of gradient inhomogeneity of the geophysical data comprises: regularization processing is carried out on the geophysical data after normalization processing; computational gridding deals with the gradient non-uniformity of geophysical data. The normalization process is a type of data normalization process, which refers to transforming a numerical range of geophysical data into the [ -1.0,1.0] interval by a linear normalization algorithm.
Normalization is performed according to the following formula:
wherein, yiTo normalize the processed geophysical data, xiFor normalizing the pre-processed geophysical data, x1Is the minimum value, x, of the pre-normalization geophysical data2Is the maximum of the geophysical data before normalization processing.
In an embodiment, the regularizing the normalized geophysical data refers to performing a gridding process on the geophysical data to form regularized geophysical data which are uniformly distributed in the horizontal direction and the vertical direction.
In an embodiment, the calculating the gradient non-uniformity of the geophysical data comprises: computing transverse and longitudinal gradient fields of the geophysical data; determining a gradient strength field of the geophysical data according to the transverse gradient field and the longitudinal gradient field of the geophysical data; calculating the non-uniformity of the geophysical data according to the gradient intensity field of the geophysical data.
In the present embodiment, in calculating the gradient nonuniformity of the geophysical data, first, the transverse gradient field and the longitudinal gradient field of the geophysical data are calculated. The gradient field is a vector field for indicating the degree of data change. For a vector y ═ f (x), the gradient is the derivative of the vector, expressed as:
wherein x is an independent variable, y is a dependent variable, and y is a gradient of the dependent variable y. Since the gradient field is a two-dimensional vector field, it can be divided into a transverse gradient field and a longitudinal gradient field according to the coordinate axis direction.
In this embodiment, the transverse gradient field and the longitudinal gradient field of the geophysical data are obtained by filtering the geophysical data by using a transverse Sobel gradient operator and a longitudinal Sobel gradient operator, respectively, and the expression of the transverse gradient field of the geophysical data is as follows:
the longitudinal gradient field expression for geophysical data is:
in this embodiment, after obtaining the transverse gradient field and the longitudinal gradient field of the geophysical data, the gradient strength field of the geophysical data is determined from the transverse gradient field and the longitudinal gradient field of the geophysical data. The gradient strength field is obtained by utilizing the operation of a transverse gradient field and a longitudinal gradient field and is used for expressing the strength degree of the gradient, and the gradient strength field of the geophysical data is determined according to the following formula:
wherein G is the gradient intensity field of the geophysical data, GxTransverse gradient field, G, being geophysical datayIs a longitudinal gradient field of geophysical data.
In this embodiment, after determining the gradient intensity field of the geophysical data, the non-uniformity of the geophysical data is calculated from the gradient intensity field of the geophysical data. Non-uniformity is a property that measures the degree of non-uniformity of a data field distribution. In general, referring to the abnormal value judgment method based on the standard deviation, it can be considered that the data field is uniform when the non-uniformity is less than 2.5; the data field is moderately non-uniform when the non-uniformity is greater than 2.5 and less than 5.0; when the non-uniformity is greater than 5.0, the data field is extremely non-uniform. Determining the non-uniformity of the geophysical data according to the following formula:
HG＝|(G2-G1)/Ga| (6)
wherein HGIs the non-uniformity of geophysical data, G2Is the maximum value of the gradient intensity field G of geophysical data, G1Is the minimum of the gradient intensity field G of the geophysical data, GaIs the average of the gradient intensity field G of the geophysical data.
In specific implementation, the spatial distribution type of the geophysical data is determined according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution.
In an embodiment, the determining the type of spatial distribution of the geophysical data from the gradient non-uniformity comprises: and comparing the gradient non-uniformity with a first preset value, and determining the spatial distribution type of the geophysical data according to a comparison result.
In this embodiment, the determining the type of the spatial distribution of the geophysical data according to the comparison result includes: when the gradient non-uniformity is smaller than a first preset value, determining that the type of the spatial distribution of the geophysical data is linear spatial distribution; and when the gradient non-uniformity is larger than or equal to a first preset value, determining that the type of the spatial distribution of the geophysical data is nonlinear spatial distribution.
In this embodiment, the first preset value is set according to actual needs, and in this embodiment, the first preset value is 2.5.
In specific implementation, the contour interval sequence of the geophysical data is determined according to the spatial distribution type.
In an embodiment, determining a sequence of contour intervals of the geophysical data according to the spatial distribution type comprises: when the spatial distribution type of the geophysical data is linear spatial distribution, determining that a contour line interval sequence of the geophysical data is a constant sequence; when the type of the spatial distribution of the geophysical data is nonlinear spatial distribution, determining that the contour interval sequence of the geophysical data is a non-constant sequence. The constant sequence means that each value in the contour interval sequence is the same, and the non-constant sequence means that each value in the contour interval sequence is different.
In this embodiment, determining the sequence of contour intervals of the geophysical data according to the spatial distribution type further includes: when the type of the spatial distribution of the geophysical data is nonlinear spatial distribution, comparing the gradient non-uniformity with a second preset value, and determining the type of the non-constant sequence according to a comparison result, wherein the type of the non-constant sequence comprises: linear type sequences and logarithmic type sequences.
In this embodiment, the determining the type of the non-constant sequence according to the alignment result includes: when the gradient non-uniformity is smaller than a second preset value, determining that the type of the non-constant sequence is a linear sequence; and when the gradient non-uniformity is greater than or equal to a second preset value, determining that the type of the non-constant sequence is a logarithmic sequence.
In this embodiment, the second preset value is set according to actual needs, and in this embodiment, the second preset value is 5.0. The linear sequence means that the median of the isoline interval sequence is in linear distribution, and the logarithmic sequence means that the median of the isoline interval sequence is in logarithmic distribution.
In an embodiment, the type of the non-constant sequence may further include: the free type sequence is that each value in the contour interval sequence is set according to the requirement so as to adapt to any complex space distribution type.
In specific implementation, a geophysical contour map is drawn according to the contour interval sequence of the geophysical data.
In an embodiment, the geophysical data is processed according to the sequence of maxima, minima, and contour intervals of the geophysical data, and a geophysical contour map is then rendered from the processed geophysical data.
In an embodiment, the processed geophysical data may also be mapped and then a geophysical contour map may be rendered from the mapped geophysical data. The mapping refers to a mathematical transformation process, or number domain transformation. Such as transforming the geophysical data from the true domain (i.e., untransformed geophysical data) to the log domain (i.e., logarithmically computing the measurements for each measurement point in one of the geophysical data and replacing the original measurements for the corresponding measurement point with them).
In an embodiment, the mapped geophysical data is subjected to planar meshing to form a mapping grid. Specifically, the whole drawing area is divided into a regular grid according to geographic coordinates, discrete geophysical data distributed in an irregular form are subjected to interpolation processing by adopting algorithms such as rectangular gridding and triangular gridding, the geophysical data of the position corresponding to each grid node in the regular grid are obtained, and drawing grid data are formed. And carrying out contour tracing on the gridded geophysical data to form an unequally-spaced contour set. Specifically, firstly, dividing equally according to a contour line interval geophysical data range to obtain a contour line tracking sequence; all contours in the sequence of contours are then traced from the drawing grid data to form a set of contours. The contour line is a smooth curve formed by connecting grid points with equal values of certain geophysical data in the drawing object and used for representing the spatial form of the characteristic value of the geophysical data, and one contour line has and only has one corresponding value of the geophysical data; the contour tracing sequence refers to a set of all contour values; the set of contours refers to the set of all traced contours. The set of contours is rendered on a computer screen or other drawing medium according to specified contour rendering parameters. Wherein the contour drawing parameters include: line type, thickness, color and inter-contour color fill one or any combination of palettes.
Based on the same inventive concept, the embodiment of the invention also provides a large-dynamic-range nonlinear geophysical contour map drawing device, which is described in the following embodiment. Because the principles for solving the problems are similar to the method for drawing the nonlinear geophysical contour map with the large dynamic range, the implementation of the device can be referred to the implementation of the method, and repeated parts are not described again.
Fig. 2 is a block diagram of a large dynamic range non-linear geophysical contour map drawing apparatus according to an embodiment of the present invention, as shown in fig. 2, the apparatus comprising:
an obtaining module 201, configured to obtain geophysical data;
a calculation module 202 for calculating gradient non-uniformity of the geophysical data;
a type determining module 203, configured to determine a spatial distribution type of the geophysical data according to the gradient non-uniformity, where the spatial distribution type includes: linear spatial distribution and non-linear spatial distribution;
a sequence determination module 204, configured to determine a contour interval sequence of the geophysical data according to the spatial distribution type;
and the drawing module 205 is configured to draw a geophysical contour map according to the contour interval sequence of the geophysical data.
In one embodiment, the calculation module 202 is specifically configured to:
computing transverse and longitudinal gradient fields of the geophysical data;
determining a gradient strength field of the geophysical data according to the transverse gradient field and the longitudinal gradient field of the geophysical data;
calculating gradient heterogeneity of the geophysical data according to the gradient intensity field of the geophysical data.
In one embodiment, the calculation module 202 is specifically configured to:
normalizing the geophysical data;
and calculating the gradient non-uniformity of the normalized geophysical data.
In an embodiment, the type determining module 203 is specifically configured to:
and comparing the gradient non-uniformity with a first preset value, and determining the spatial distribution type of the geophysical data according to a comparison result.
In an embodiment, the type determining module 203 is specifically configured to:
when the gradient non-uniformity is smaller than a first preset value, determining that the type of the spatial distribution of the geophysical data is linear spatial distribution;
and when the gradient non-uniformity is larger than or equal to a first preset value, determining that the type of the spatial distribution of the geophysical data is nonlinear spatial distribution.
In one embodiment, the sequence determination module 204 is further configured to:
when the spatial distribution type of the geophysical data is linear spatial distribution, determining that a contour line interval sequence of the geophysical data is a constant sequence;
when the type of the spatial distribution of the geophysical data is nonlinear spatial distribution, determining that the contour interval sequence of the geophysical data is a non-constant sequence.
In one embodiment, the sequence determination module 204 is further configured to:
when the type of the spatial distribution of the geophysical data is nonlinear spatial distribution, comparing the gradient non-uniformity with a second preset value, and determining the type of the non-constant sequence according to a comparison result, wherein the type of the non-constant sequence comprises: linear type sequences and logarithmic type sequences.
In one embodiment, the sequence determining module 204 is specifically configured to:
when the gradient non-uniformity is smaller than a second preset value, determining that the type of the non-constant sequence is a linear sequence;
and when the gradient non-uniformity is greater than or equal to a second preset value, determining the type of the non-constant sequence to be a logarithmic sequence.
A specific embodiment is given below to illustrate a specific application of the large dynamic range non-linear geophysical contour map drawing method in the embodiment of the present invention. In this embodiment, as shown in fig. 3, the gravity grid anomaly data collected from the gravity survey is a two-dimensional data reflecting the spatial variation characteristics of the gravity acceleration caused by the variation of the density of the underground rock. The data has 2774 measurement points which are regularly arranged, the value range of the gravity grid distribution abnormal data is distributed between 3.3877 and 842.8824, the value range of the gradient intensity is 0.00 to 701.86, the average value is 105.38, and the non-uniformity is 5.66, so the spatial distribution type belongs to non-linear spatial distribution.
Fig. 4 is a table in which the contour interval sequence is a constant sequence, and in fig. 4, the contour interval sequence is a constant sequence with an interval of 50.00. FIG. 5 is a geophysical contour map plotted according to a sequence of constants. As can be seen from fig. 5, because the distribution of the numerical range of the gravity anomaly data is not uniform, a large area of blank space exists at the lower part of the center, and the details of the distribution characteristics of the gravity anomaly space cannot be shown.
FIG. 6 is a geophysical contour plot plotted according to a linear type sequence, in FIG. 6, the sequence of contour intervals exhibits a linear variation from 5.00 to 85.00 with increasing values of contour intervals and increments of adjacent contour intervals of 5.00. Figure 7 is a geophysical contour map plotted according to a linear-type sequence. Comparing fig. 7 with fig. 5, it can be seen that the large-area blank space existing in the lower part of the center of fig. 5 is filled up by the contour lines 2, 3, and 4, and the gravity abnormal distribution characteristics are more detailed.
FIG. 8 is a geophysical contour plot plotted according to a logarithmic sequence, in FIG. 8 the sequence of contour intervals exhibits a logarithmic change, from 1.00 to 100.00, with values of contour intervals exhibiting a step-wise increase. FIG. 9 is a geophysical contour map plotted according to a logarithmic sequence. Comparing fig. 9 with fig. 5, it can be seen that the large-area blank space existing in the lower part of the center of fig. 5 is filled by the contour lines 4-11, and the details of the gravity abnormal distribution feature shown are clearer.
FIG. 10 is a geophysical contour plot plotted against a freeform sequence, in FIG. 10 the sequence of contour spacings exhibits an irregular variation from 4.00 to 50.00. FIG. 11 is a geophysical contour map plotted according to a freeform sequence. Comparing fig. 11 with fig. 5, it can be seen that two contour lines are added to fig. 11, so that the details of the displayed gravity anomaly distribution feature are clearer. The geophysical contour map drawn based on the free form sequence can control and adjust the gravity abnormal distribution characteristics according to the needs, and can further highlight the detail characteristics of the geophysical data on the basis of clearly expressing the overall form of the geophysical data.
In summary, the embodiments of the present invention obtain geophysical data; calculating a gradient non-uniformity of the geophysical data; determining a spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution; determining a contour interval sequence of the geophysical data according to the spatial distribution type; and drawing a geophysical contour map according to the contour interval sequence of the geophysical data. The embodiment of the invention calculates the gradient non-uniformity of the geophysical data, determines the spatial distribution type of the geophysical data, further determines the contour line interval sequence, sets different contour line interval sequences for the geophysical data of each spatial distribution type, and adjusts the distance between adjacent contour lines in the contour map, thereby avoiding the problem that the distance between adjacent contour lines in the contour map is abnormally dense or sparse when the geophysical data are changed severely, and effectively retaining the detail characteristics of the geophysical data on the basis of clearly expressing the integral form of the geophysical data.
As will be appreciated by one skilled in the art, embodiments of the present invention may be provided as a method, system, or computer program product. Accordingly, the present invention may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present invention may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present invention is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the invention. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Claims (18)
1. A large dynamic range non-linear geophysical contour map drawing method is characterized by comprising the following steps:
acquiring geophysical data;
calculating a gradient non-uniformity of the geophysical data;
determining a spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution;
determining a contour interval sequence of the geophysical data according to the spatial distribution type;
and drawing a geophysical contour map according to the contour interval sequence of the geophysical data.
2. The method of claim 1, wherein calculating gradient non-uniformities for the geophysical data comprises:
computing transverse and longitudinal gradient fields of the geophysical data;
determining a gradient strength field of the geophysical data according to the transverse gradient field and the longitudinal gradient field of the geophysical data;
calculating gradient heterogeneity of the geophysical data according to the gradient intensity field of the geophysical data.
3. The method of claim 1, wherein after acquiring the geophysical data, performing a normalization process on the geophysical data;
the calculating of gradient inhomogeneity of the geophysical data comprises: and calculating the gradient non-uniformity of the normalized geophysical data.
4. The method of claim 1, wherein said determining a type of spatial distribution of the geophysical data based on the gradient non-uniformity comprises:
and comparing the gradient non-uniformity with a first preset value, and determining the spatial distribution type of the geophysical data according to a comparison result.
5. The method of claim 4, wherein determining the type of spatial distribution of the geophysical data from the comparison comprises:
when the gradient non-uniformity is smaller than a first preset value, determining that the type of the spatial distribution of the geophysical data is linear spatial distribution;
and when the gradient non-uniformity is larger than or equal to a first preset value, determining that the type of the spatial distribution of the geophysical data is nonlinear spatial distribution.
6. The method of claim 1, wherein determining a sequence of contour intervals for the geophysical data based on the spatial distribution type comprises:
when the spatial distribution type of the geophysical data is linear spatial distribution, determining that a contour line interval sequence of the geophysical data is a constant sequence;
when the type of the spatial distribution of the geophysical data is nonlinear spatial distribution, determining that the contour interval sequence of the geophysical data is a non-constant sequence.
7. The method of claim 6, wherein determining a sequence of contour intervals for the geophysical data from the spatial distribution types further comprises:
when the type of the spatial distribution of the geophysical data is nonlinear spatial distribution, comparing the gradient non-uniformity with a second preset value, and determining the type of the non-constant sequence according to a comparison result, wherein the type of the non-constant sequence comprises: linear type sequences and logarithmic type sequences.
8. The method of claim 7, wherein determining the type of the non-constant sequence from the alignment comprises:
when the gradient non-uniformity is smaller than a second preset value, determining that the type of the non-constant sequence is a linear sequence;
and when the gradient non-uniformity is greater than or equal to a second preset value, determining that the type of the non-constant sequence is a logarithmic sequence.
9. A large dynamic range non-linear geophysical contour map drawing apparatus comprising:
the acquisition module is used for acquiring geophysical data;
a calculation module for calculating gradient non-uniformity of the geophysical data;
a type determination module for determining a spatial distribution type of the geophysical data according to the gradient non-uniformity, wherein the spatial distribution type comprises: linear spatial distribution and non-linear spatial distribution;
a sequence determination module for determining a contour interval sequence of the geophysical data according to the spatial distribution type;
and the drawing module is used for drawing the geophysical contour map according to the contour line interval sequence of the geophysical data.
10. The apparatus of claim 9, wherein the computing module is specifically configured to:
computing transverse and longitudinal gradient fields of the geophysical data;
determining a gradient strength field of the geophysical data according to the transverse gradient field and the longitudinal gradient field of the geophysical data;
calculating the non-uniformity of the gradient intensity field of the geophysical data according to the gradient intensity field of the geophysical data.
11. The apparatus of claim 9, wherein the computing module is specifically configured to:
normalizing the geophysical data;
and calculating the gradient non-uniformity of the normalized geophysical data.
12. The apparatus of claim 9, wherein the type determination module is specifically configured to:
and comparing the gradient non-uniformity with a first preset value, and determining the spatial distribution type of the geophysical data according to a comparison result.
13. The apparatus of claim 12, wherein the type determination module is specifically configured to:
when the gradient non-uniformity is smaller than a first preset value, determining that the type of the spatial distribution of the geophysical data is linear spatial distribution;
and when the gradient non-uniformity is larger than or equal to a first preset value, determining that the type of the spatial distribution of the geophysical data is nonlinear spatial distribution.
14. The apparatus of claim 9, wherein the sequence determination module is further to:
when the spatial distribution type of the geophysical data is linear spatial distribution, determining that a contour line interval sequence of the geophysical data is a constant sequence;
when the type of the spatial distribution of the geophysical data is nonlinear spatial distribution, determining that the contour interval sequence of the geophysical data is a non-constant sequence.
15. The apparatus of claim 14, wherein the sequence determination module is further to:
when the type of the spatial distribution of the geophysical data is nonlinear spatial distribution, comparing the gradient non-uniformity with a second preset value, and determining the type of the non-constant sequence according to a comparison result, wherein the type of the non-constant sequence comprises: linear type sequences and logarithmic type sequences.
16. The apparatus of claim 15, wherein the sequence determination module is specifically configured to:
when the gradient non-uniformity is smaller than a second preset value, determining that the type of the non-constant sequence is a linear sequence;
and when the gradient non-uniformity is greater than or equal to a second preset value, determining the type of the non-constant sequence to be a logarithmic sequence.
17. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the processor implements the method of any one of claims 1 to 8 when executing the computer program.
18. A computer-readable storage medium, characterized in that the computer-readable storage medium stores a computer program for executing the method of any one of claims 1 to 8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911079926.3A CN112785910B (en) | 2019-11-07 | 2019-11-07 | Large dynamic range nonlinear geophysical contour drawing method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201911079926.3A CN112785910B (en) | 2019-11-07 | 2019-11-07 | Large dynamic range nonlinear geophysical contour drawing method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN112785910A true CN112785910A (en) | 2021-05-11 |
| CN112785910B CN112785910B (en) | 2023-07-25 |
Family
ID=75747753
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201911079926.3A Active CN112785910B (en) | 2019-11-07 | 2019-11-07 | Large dynamic range nonlinear geophysical contour drawing method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112785910B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN119293708A (en) * | 2024-12-13 | 2025-01-10 | 应急管理部国家自然灾害防治研究院 | Spatial data local dense anomaly analysis method, device and storage medium |
| CN119625202A (en) * | 2024-11-27 | 2025-03-14 | 深层探索(北京)科技有限公司 | Geophysical data visualization method and device based on subatomic radiation characteristics |
| CN119625202B (en) * | 2024-11-27 | 2025-10-17 | 深层探索(北京)科技有限公司 | Geophysical data visualization method and device based on subatomic radiation characteristics |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090087034A1 (en) * | 2007-10-01 | 2009-04-02 | Magellan Navigation, Inc. | Static and dynamic contours |
| CN103995944A (en) * | 2014-06-10 | 2014-08-20 | 中国地质大学（武汉） | Space abnormal extreme value unbiased isogram drawing method |
| CN104166158A (en) * | 2013-05-17 | 2014-11-26 | 中国石油天然气集团公司 | Earthquake data non-uniform interval contour line interaction tracking method |
| CN106228589A (en) * | 2016-08-04 | 2016-12-14 | 北京金双狐油气技术有限公司 | Gridding method with direction constraint scatterplot data |
| CN106981084A (en) * | 2016-10-28 | 2017-07-25 | 阿里巴巴集团控股有限公司 | A kind of method and device of drawing isoline |
| CN107492303A (en) * | 2017-08-28 | 2017-12-19 | 中国南方电网有限责任公司电网技术研究中心 | Drawing method and system for equivalent salt deposit density distribution map of power transmission line in coastal region |
| CN108197347A (en) * | 2017-12-08 | 2018-06-22 | 北京农业信息技术研究中心 | The geographic space distribution uniformity detection method of continuous irregular area sampled point |
| CN108776999A (en) * | 2018-06-06 | 2018-11-09 | 天津大学 | Grid isoplethes drawing method based on ocean Internet of Things |
| CN109741446A (en) * | 2018-12-12 | 2019-05-10 | 四川华控图形科技有限公司 | A kind of method of the fine coastal landform of dynamic generation in three-dimensional digital earth |
-
2019
- 2019-11-07 CN CN201911079926.3A patent/CN112785910B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090087034A1 (en) * | 2007-10-01 | 2009-04-02 | Magellan Navigation, Inc. | Static and dynamic contours |
| CN104166158A (en) * | 2013-05-17 | 2014-11-26 | 中国石油天然气集团公司 | Earthquake data non-uniform interval contour line interaction tracking method |
| CN103995944A (en) * | 2014-06-10 | 2014-08-20 | 中国地质大学（武汉） | Space abnormal extreme value unbiased isogram drawing method |
| CN106228589A (en) * | 2016-08-04 | 2016-12-14 | 北京金双狐油气技术有限公司 | Gridding method with direction constraint scatterplot data |
| CN106981084A (en) * | 2016-10-28 | 2017-07-25 | 阿里巴巴集团控股有限公司 | A kind of method and device of drawing isoline |
| CN107492303A (en) * | 2017-08-28 | 2017-12-19 | 中国南方电网有限责任公司电网技术研究中心 | Drawing method and system for equivalent salt deposit density distribution map of power transmission line in coastal region |
| CN108197347A (en) * | 2017-12-08 | 2018-06-22 | 北京农业信息技术研究中心 | The geographic space distribution uniformity detection method of continuous irregular area sampled point |
| CN108776999A (en) * | 2018-06-06 | 2018-11-09 | 天津大学 | Grid isoplethes drawing method based on ocean Internet of Things |
| CN109741446A (en) * | 2018-12-12 | 2019-05-10 | 四川华控图形科技有限公司 | A kind of method of the fine coastal landform of dynamic generation in three-dimensional digital earth |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN119625202A (en) * | 2024-11-27 | 2025-03-14 | 深层探索(北京)科技有限公司 | Geophysical data visualization method and device based on subatomic radiation characteristics |
| CN119625202B (en) * | 2024-11-27 | 2025-10-17 | 深层探索(北京)科技有限公司 | Geophysical data visualization method and device based on subatomic radiation characteristics |
| CN119293708A (en) * | 2024-12-13 | 2025-01-10 | 应急管理部国家自然灾害防治研究院 | Spatial data local dense anomaly analysis method, device and storage medium |
| CN119293708B (en) * | 2024-12-13 | 2025-03-21 | 应急管理部国家自然灾害防治研究院 | Spatial data local dense anomaly analysis method, device and storage medium |
Also Published As
| Publication number | Publication date |
|---|---|
| CN112785910B (en) | 2023-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Dannemiller et al. | A new method for determination of magnetization direction | |
| Pilkington et al. | Potential field continuation between arbitrary surfaces—Comparing methods | |
| Paoletti et al. | Self-constrained inversion of potential fields | |
| Paoletti et al. | Inversion of gravity gradient tensor data: Does it provide better resolution? | |
| Fedi et al. | MHODE: a local-homogeneity theory for improved source-parameter estimation of potential fields | |
| Abbas et al. | Automatic DEXP imaging of potential fields independent of the structural index | |
| CN107221028B (en) | A 3D reconstruction method of closed surface of geological body based on seismic interpretation data | |
| Huang et al. | Bidimensional empirical mode decomposition (BEMD) for extraction of gravity anomalies associated with gold mineralization in the Tongshi gold field, Western Shandong Uplifted Block, Eastern China | |
| US10422900B2 (en) | Analyzing seismic data | |
| Goff et al. | Conditional simulation of Thwaites Glacier (Antarctica) bed topography for flow models: Incorporating inhomogeneous statistics and channelized morphology | |
| Dai et al. | The forward modeling of 3D gravity and magnetic potential fields in space-wavenumber domains based on an integral method | |
| CA2936156A1 (en) | System and method for two dimensional gravity modeling with variable densities | |
| Feng et al. | Gravity inversion of blocky basement relief using L0 norm constraint with exponential density contrast variation | |
| Mastellone et al. | Volume continuation of potential fields from the minimum-length solution: An optimal tool for continuation through general surfaces | |
| CN112785910B (en) | Large dynamic range nonlinear geophysical contour drawing method and device | |
| Liu et al. | Multiscale fractal characterization of hierarchical heterogeneity in sandstone reservoirs | |
| Wang et al. | A constrained scheme for high precision downward continuation of potential field data | |
| Nava-Flores et al. | 3D gravity modeling of complex salt features in the southern Gulf of Mexico | |
| Hidalgo-Gato et al. | Fast 3D magnetic inversion of a surface relief in the space domain | |
| AU2013303068A1 (en) | Methods and corresponding software module for quantifying risks or likelihoods of hydrocarbons being present in a geological basin or region | |
| Chowdhury et al. | Geometric shapes derived from airborne gravity gradiometry data: New tools for the explorationist | |
| Yañez et al. | Topographic correction of magnetic data on rugged topography with application to Río Blanco-Los Bronces and El Teniente porphyry copper districts, Southern Andes, Chile | |
| Wang et al. | Exponential Approximate Invariance: mapping the depth to basement from gravity anomalies | |
| Guo et al. | Apparent magnetization mapping in the presence of strong remanent magnetization: The space-domain inversion approach | |
| Qi* et al. | Segmentation of salt domes, mass transport complexes on 3D seismic data volumes using Kuwahara windows and multiattribute cluster analysis |
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 |