CN112347421A - Method and system for highlighting and enhancing broken-line-shaped gravity anomaly information - Google Patents
Method and system for highlighting and enhancing broken-line-shaped gravity anomaly information Download PDFInfo
- Publication number
- CN112347421A CN112347421A CN202011110088.4A CN202011110088A CN112347421A CN 112347421 A CN112347421 A CN 112347421A CN 202011110088 A CN202011110088 A CN 202011110088A CN 112347421 A CN112347421 A CN 112347421A
- Authority
- CN
- China
- Prior art keywords
- gravity
- window
- sub
- gravity anomaly
- anomaly
- 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.)
- Pending
Links
- 230000005484 gravity Effects 0.000 title claims abstract description 209
- 238000000034 method Methods 0.000 title claims abstract description 58
- 230000002708 enhancing effect Effects 0.000 title claims abstract description 36
- 238000001914 filtration Methods 0.000 claims abstract description 43
- 230000002159 abnormal effect Effects 0.000 claims abstract description 37
- JTJMJGYZQZDUJJ-UHFFFAOYSA-N phencyclidine Chemical compound C1CCCCN1C1(C=2C=CC=CC=2)CCCCC1 JTJMJGYZQZDUJJ-UHFFFAOYSA-N 0.000 claims description 9
- 238000001514 detection method Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- 230000007547 defect Effects 0.000 abstract description 10
- 238000000605 extraction Methods 0.000 abstract description 5
- 230000007812 deficiency Effects 0.000 abstract 1
- 230000000694 effects Effects 0.000 description 15
- 238000010586 diagram Methods 0.000 description 12
- 230000005856 abnormality Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 238000005755 formation reaction Methods 0.000 description 3
- 239000000463 material Substances 0.000 description 2
- 238000003672 processing method Methods 0.000 description 2
- 239000000758 substrate Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
- G01V7/02—Details
- G01V7/06—Analysis or interpretation of gravimetric records
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/20—Drawing from basic elements, e.g. lines or circles
- G06T11/206—Drawing of charts or graphs
-
- G06T5/70—
Abstract
The invention discloses a method and a system for highlighting and enhancing fracture linear gravity anomaly information, and relates to the technical field of gravity exploration in geophysical exploration. The method comprises the following steps: setting a sliding window, and dividing the sliding window into a plurality of sub-windows; taking the gravity abnormal mean value of the sub-window with the minimum mean square error as the filtering output of the center point of the sliding window; scanning the gravity anomaly data area by using a sliding window according to a preset sliding rule; and repeating the steps until the preset iteration times are reached. The invention not only realizes the reliable extraction of the fracture abnormity information, but also can identify and depict the boundary of the non-fractured geologic body, enriches the application field of the invention, effectively overcomes the defect and the deficiency that the fracture gravity abnormity information is highlighted by various related methods based on the direction derivative, adds a new processing means for gravity data processing, and has very wide application prospect.
Description
Technical Field
The invention relates to the technical field of gravity exploration in geophysical exploration, in particular to a method and a system for highlighting and enhancing broken linear gravity anomaly information.
Background
Fracturing is one of the important geologic formation phenomena, and most other geologic phenomena are related to fractured formations. Fractures often appear as linear anomaly bands extending in a certain direction on gravity anomalies, and gravity linear anomalies reflecting large-scale fractures are often represented as gradient bands of gravity anomalies. In order to better apply gravity to research fracture, the characteristic of a gradient zone of gravity anomaly is not enough, purposefully data processing of the characteristic of the gravity anomaly is required to be performed, the characteristic of the fracture is more obvious and easier to identify through processing of the characteristic expressed on the gravity of a region by specific reinforced fracture, the determined fracture position is more reliable, and fracture division is more concise.
At present, in the process of extracting and identifying fracture information of gravity, in order to extract and highlight information of a fracture in a certain direction in a gravity anomaly, the gravity is usually processed by a first-order horizontal derivative in a vertical fracture direction, and in the gravity data processing, the processing method is called as direction derivative processing of the gravity anomaly.
However, the processing of the directional derivative enhances the fracture anomaly information and also adds some unnecessary false fracture information, and the gravity anomaly generated by some other geologic bodies is processed into the fracture anomaly information in an error manner, so that the reliability of applying gravity to explain the fracture is reduced, and the effectiveness of the directional derivative on fracture information extraction is also reduced.
Disclosure of Invention
The present invention provides a method and system for highlighting and enhancing information of broken linear gravity anomaly, aiming at the defects of the prior art.
The technical scheme for solving the technical problems is as follows:
a method for highlighting and enhancing broken line-shaped gravity anomaly information comprises the following steps:
step 1, setting a sliding window for performing sub-data window filtering on a gravity anomaly data area, and dividing the sliding window into a plurality of sub-windows, wherein the shapes of the sub-windows are different from each other, but the number of included gravity anomaly points is the same;
and 4, repeating the steps 1 to 3 until the preset iteration times are reached.
Another technical solution of the present invention for solving the above technical problems is as follows:
a system for highlighting and enhancing broken line-like gravity anomaly information, comprising:
the device comprises a setting unit, a data processing unit and a processing unit, wherein the setting unit is used for setting a sliding window used for carrying out sub-data window filtering on a gravity anomaly data area, and dividing the sliding window into a plurality of sub-windows, and the shapes of the sub-windows are different from each other, but the number of included gravity anomaly points is the same;
the computing unit is used for respectively taking the mean value and the mean square error of all the gravity anomaly points of each sub-window and taking the gravity anomaly mean value of the sub-window with the minimum mean square error as the filtering output of the center point of the sliding window;
the filtering unit is used for scanning the gravity abnormal data area by using the sliding window according to a preset sliding rule and carrying out sub-data window filtering processing on the data of the gravity abnormal data area;
and the detection unit is used for detecting whether the preset iteration times are reached.
The invention has the beneficial effects that: the invention divides the sliding window in the sub-data window filtering method into a plurality of sub-windows, then the gravity anomaly mean value of the sub-window with the minimum mean square error in the sub-windows is taken as the filtering output of the central point of the sliding window, and the gravity anomaly along the fracture linear anomaly direction is determined and highlighted and enhanced by the variance of the gravity anomaly in a certain size window, thereby achieving the purpose of enhancing the fracture gravity anomaly information, not only realizing the reliable extraction of the fracture anomaly information, but also identifying and depicting the boundary of the non-fractured geologic body, enriching the application field of the invention patent, effectively overcoming the defects and the defects of highlighting the fracture gravity anomaly information based on various related methods of direction derivatives, adding a new processing means for gravity data processing, and having very wide application prospect.
Advantages of additional aspects of the invention will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention.
Drawings
FIG. 1 is a schematic flow chart diagram illustrating a method for highlighting and enhancing information on broken linear gravity anomaly according to an embodiment of the present invention;
FIG. 2 is a schematic diagram illustrating sub-window division provided by an embodiment of a method for highlighting and enhancing broken-line-shaped gravity anomaly information according to the present invention;
FIG. 3 is a schematic diagram illustrating data expansion provided by an embodiment of a method for highlighting and enhancing information about broken-line-shaped gravity anomaly according to the present invention;
FIG. 4 is a schematic diagram of a deep source region gravity anomaly provided by an embodiment of a method for highlighting and enhancing broken line-shaped gravity anomaly information according to the present invention;
FIG. 5 is a schematic diagram illustrating the effect of processing gravity anomaly in a deep source region according to an embodiment of the method for highlighting and enhancing information on broken linear gravity anomaly according to the present invention;
FIG. 6 is a schematic view of a shell source region gravity anomaly provided by an embodiment of a method for highlighting and enhancing fracture line-shaped gravity anomaly information according to the present invention;
FIG. 7 is a schematic diagram illustrating the effect of handling gravity anomaly in a shell source region according to an embodiment of the method for highlighting and enhancing information on broken-line-shaped gravity anomaly according to the present invention;
FIG. 8 is a schematic view of a Duobaoshan gravity anomaly information method according to an embodiment of the present invention for highlighting and enhancing broken line-shaped gravity anomaly information;
FIG. 9 is a schematic diagram illustrating the processing effect of the method for highlighting and enhancing the information about the broken linear gravity anomaly according to the embodiment of the present invention on the gravity anomaly of Duobaoshan;
fig. 10 is a structural frame diagram provided by an embodiment of a system for highlighting and enhancing broken-line-shaped gravity anomaly information according to the present invention.
Detailed Description
The principles and features of this invention are described below in conjunction with the following drawings, which are set forth to illustrate, but are not to be construed to limit the scope of the invention.
In oil and gas exploration, fracture not only controls the structure of a basin and the distribution of sediments, but also plays a dual role in controlling and destroying the formation of an oil and gas reservoir, and the fracture is also an important geological structure for communicating a deep part with a shallow part, and has important geological significance for the research of the fracture structure in the aspects of geological structure, regional structure, deep structure, earth kinematics and mineralogy.
Due to the geophysical characteristics of the geologic body, when the geologic body is displaced and faulted in three-dimensional space after fracture occurs, the density between stratums changes, the fracture reflects a steep zone of a stratum density interface, the larger the scale of the fracture is, the larger the scale of the steep zone of the two physical interfaces is, the larger the density difference is, and the more obvious the form of a gravity abnormal step zone (also a linear abnormal zone) is. The substrate and higher levels of fracture below it can be effectively divided by the use of regional gravity. The large-scale ultra-shell deep fracture, shell fracture, regional fracture and substrate fracture can be clearly reflected on the heavy anomaly graph.
Fractures often appear as linear anomaly bands extending in a certain direction on gravity anomalies, and gravity linear anomalies reflecting large-scale fractures are often represented as gradient bands of gravity anomalies. In order to better apply gravity to research fracture, the characteristic of a gradient zone of gravity anomaly is not enough, purposefully data processing of the characteristic of the gravity anomaly is required to be performed, the characteristic of the fracture is more obvious and easier to identify through processing of the characteristic expressed on the gravity of a region by specific reinforced fracture, the determined fracture position is more reliable, and fracture division is more concise.
At present, in the process of extracting and identifying fracture information of gravity, in order to extract and highlight information of a fracture in a certain direction in a gravity anomaly, the gravity is usually processed by a first-order horizontal derivative in a vertical fracture direction, and in the gravity data processing, the processing method is called as direction derivative processing of the gravity anomaly.
Mathematically, the directional derivative of the gravitational field is defined as follows:
wherein, g is the gravity anomaly,
is the directional derivative of the gravity anomaly,
in the case of a gradient of gravity anomaly,
is a unit vector in the direction of t,
gradient and unit vector of gravity anomaly
The inner product of (d).
Although it is not limited to
The information of the fracture in the gravity field can be highlighted to a certain extent, but the result is related to the derivation direction, the horizontal derivative of the gravity anomaly reflecting the same fracture in one direction is represented as a positive anomaly, the horizontal derivative in the other direction is represented as a negative anomaly, and the interpretation of the fracture information is ambiguous after the horizontal derivative processing is adopted for applying the gravity anomaly.
In addition, while the direction derivative is used for highlighting the gravity anomaly fracture anomaly information, some gravity anomalies which are not generated by fracture are distorted and are represented as fracture-like gravity anomaly information. For example: the directional derivative enables equiaxed anomalies created by the spheres to collapse into fracture-like linear anomaly information. The processing of the directional derivative enhances the fracture anomaly information, and simultaneously adds some unnecessary false fracture information, and incorrectly processes the gravity anomaly generated by other geologic bodies into the fracture anomaly information, thereby reducing the reliability of applying gravity to explain fracture and reducing the effectiveness of the directional derivative on fracture information extraction.
The invention provides a practical method for enhancing the position of the gravity linear abnormal information based on the defects and defects of the direction derivative prominent fracture in the gravity abnormal information, and can be seen from the following embodiment of extracting and highlighting the information of the fracture in the gravity abnormal of the invention. The method not only realizes the reliable extraction of the fracture abnormity information, but also can identify and depict the boundary of the non-fractured geologic body, thereby enriching the application field of the invention.
The scheme disclosed by the invention is realized by sub-data window filtering, the sub-data window filtering is a nonlinear filtering method, filtering is carried out by selecting a proper filtering window and iteration times, the method is applied to gravity data processing, the effect of highlighting and enhancing a gravity gradient zone can be achieved, or linear gravity anomaly can be highlighted and enhanced, and the calculation implementation process and content of the gravity sub-data window filtering are explained by combining a specific example.
As shown in fig. 1, a schematic flow chart is provided for an embodiment of a method for highlighting and enhancing broken-line-shaped gravity anomaly information according to the present invention, and the method includes:
and S1, setting a sliding window for sub-data window filtering of the gravity anomaly data area, and dividing the sliding window into a plurality of sub-windows, wherein the shapes of the sub-windows are different, but the number of the included gravity anomaly points is the same.
Before the sub-data window filtering calculation is used, the size of the sliding window, the number of rows and columns of the sliding window, the number of iterations, and the like need to be preset, and these can be set according to actual requirements. The number and shape of the divided sub-windows and the number of the gravity anomaly points contained in each sub-window can also be set according to actual requirements.
For example, as shown in fig. 2, the size of the sliding window is 5 × 5, and the number of the sliding window is denoted by
The sliding window can be used as a mother window, 8 sub-windows are divided, the number of the sub-windows can be (i) to (b), the mother window comprises 25 gravity anomaly points, each sub-window can comprise 15 gravity anomaly points, as shown in fig. 2, the shape of each sub-window can be set according to actual requirements, and it is sufficient to ensure that each sub-window comprises 15 gravity anomaly points.
Preferably, to ensure that the data is fully covered, all gravity anomaly points of the parent window should be covered when all child windows are overlapped together, for example, as shown in fig. 2, all gravity anomaly points of the parent window are covered when the child window (i) and the child window (ii) are overlapped together, and other child windows are similar and will not be described again.
Preferably, to ensure the data coverage is complete, the sub-windows may be set along the axial symmetry line, for example, as shown in fig. 2, the sub-windows (i) and (ii) are axially symmetric along the 3 rd row, and the other sub-windows are similar and will not be described again.
Preferably, in order to ensure that the data coverage is complete, a plurality of sets of child window pairs arranged along different axisymmetric lines may be provided, for example, as shown in fig. 2, a child window (i) and a child window (ii) are axisymmetric along the 3 rd row, and the two child windows cover all gravity anomaly points of the parent window and are one child window pair; the child window (c) and the child window (c) are axisymmetric along the axis between the data point at the 1 st row and the data point at the 5 th column at the 1 st row, the two child windows cover all gravity anomaly points of the parent window, and are also a child window pair, and other child window pairs are the same and are not repeated, so that the child window is selected, the data coverage can be ensured to be comprehensive, the relation among different combinations of the data can be fully excavated, and the effect of highlighting and enhancing can be further improved.
And S2, respectively taking the mean value and the mean square error of all gravity anomaly points of each sub-window, and taking the gravity anomaly mean value of the sub-window with the minimum mean square error as the filtering output of the center point of the sliding window.
It should be noted that each gravity abnormal point corresponds to a gravity abnormal value, and the average value of these gravity abnormal values can be obtained in the existing manner, and then the mean square error is obtained according to the average value.
For example, taking the sub-window shown in fig. 2 as an example, assuming that the mean square error of the gravity abnormal value of the second sub-window is found to be the smallest through calculation, the mean of the gravity abnormal values of the second sub-window may be used as the filtered output of the center point of the parent window.
And S3, scanning the gravity abnormal data area by using a sliding window according to a preset sliding rule, and performing sub-data window filtering processing on the data of the gravity abnormal data area.
It should be noted that the sliding rule may be set according to actual requirements, for example, the sliding window may scan the entire gravity anomaly data area according to a certain number of sliding rows and columns.
S4, repeating S1 to S3 until a preset number of iterations is reached.
It should be understood that since the sub-data window filtering is performed in the spatial domain of gravity anomaly, edge loss is necessarily generated. In order to obtain linear abnormal information which better highlights and enhances the broken gravity anomaly, one sub-data window filtering is not enough, and a process of repeating the sub-data window filtering is needed for multiple times, wherein the process is called iteration, each iteration is performed on the basis of the result of the last sub-data window filtering, and each iteration further enhances the linear abnormal information of the broken gravity anomaly.
Specifically, the optimal sub-data window filtering iteration number can be selected through experiments.
In the embodiment, the sliding window in the sub-data window filtering method is divided into a plurality of sub-windows, then the gravity anomaly mean value of the sub-window with the minimum mean square error in the sub-windows is used as the filtering output of the central point of the sliding window, and the gravity anomaly along the fracture linear anomaly direction is determined and highlighted and enhanced according to the variance of the gravity anomaly in the window with a certain size, so that the purpose of enhancing the fracture gravity anomaly information is achieved, not only is the fracture anomaly information reliably extracted, but also the boundary identification and drawing of the non-fractured geologic body can be carried out, the application field of the invention is enriched, the defects and the defects that the fracture gravity anomaly information is highlighted by various related methods based on the direction derivative are effectively overcome, a new processing means is added for gravity data processing, and the invention has a very wide application prospect.
Optionally, in some possible embodiments, the method further includes:
and expanding the gravity anomaly data area according to the preset iteration times and the size of the sliding window, so that the effective data area is consistent with the original data area.
It should be understood that, in order to overcome the disadvantage of the edge loss of the spatial domain filtering, a function of expanding the original gravity data area according to the number of iterations and the size of the parent window may be added. Therefore, after the sub-data window filtering is carried out for multiple times, the effective data area is consistent with the original data area, and the edge loss is effectively overcome.
The specific expansion mode may be set according to actual requirements, for example, taking fig. 3a as an example, assuming that the size of the sliding window is 2 × 2, and the right side of the 5 th column is an edge, when the sliding window slides to the 5 th column, the data in the 5 th column is only filtered once, and the output data has edge loss, at this time, the right side of the 5 th column may be subjected to data expansion, for example, as shown in fig. 3b, the data in the 5 th column may be entirely copied to the right side, and the 6 th column may be expanded, so that the effective data area may be kept consistent with the original data area. It should be understood that the mean value of the data in the 4 th column and the 5 th column in the same row may also be taken as the data in the 6 th column of the row, and other expansion manners may also be adopted, which may be set according to actual requirements.
It should be understood that the above is only an exemplary expansion manner, and specifically, what manner to expand and how much data to expand, and those skilled in the art may set the expansion manner according to the preset iteration number and the size of the sliding window, which is not described herein again.
Optionally, in some possible embodiments, the method further includes:
and storing the obtained processing result as a grid format file of SURFER drawing software, and drawing a sub-data window filter graph based on the reinforced gravity linear abnormal information protruding fracture position by using the SURFER drawing software.
The processing result is displayed through SURFER drawing software, and a sub data window filter graph based on the reinforced gravity linear abnormal information protruding fracture position can be better displayed.
Optionally, in some possible embodiments, the mean of all gravity anomaly points for each sub-window is calculated according to the following formula:
wherein，
Is the mean value of gravity anomaly points of the ith sub-window, gi(j) I is the gravity anomaly value of the jth gravity anomaly point in the ith sub-window, I is 1, 2.
Optionally, in some possible embodiments, the mean square error of all gravity anomaly points for each sub-window is calculated according to the following formula:
wherein σiIs the mean square error of the gravity anomaly point in the ith sub-window.
To further illustrate the effects of the present invention, the following description will be given with reference to examples.
Taking the effect of gravity processing of the salient and enhanced fracture linear abnormal information in the ultra-deep source region of the Songliaopelvic region as an example, fig. 4 is a schematic diagram of gravity abnormality of the deep source region, and fig. 5 is a schematic diagram of the processing effect of the gravity abnormality of the deep source region by the method disclosed by the invention, and it can be obviously seen from the diagram that the method has better effect of salient and enhanced fracture gravity linear abnormal information on the gravity abnormality of the deep source region.
Taking the effect of highlighting and enhancing the fracture linear anomaly information in the gravity of the shell source region as an example, fig. 6 is a schematic view of the gravity anomaly of the shell source region, and fig. 7 is a schematic view of the processing effect of the gravity anomaly of the shell source region by the method disclosed by the invention.
Taking the effect of highlighting and enhancing broken-line-shaped abnormal information in the gravity processing of the Duobaoshan as an example, fig. 8 is a schematic view of gravity anomaly of the Duobaoshan, and fig. 9 is a schematic view of the effect of processing gravity anomaly of the Duobaoshan by the method disclosed by the invention, and it is obvious from the figure that the method also has better effect of highlighting and enhancing broken-line-shaped abnormal information of gravity anomaly of a mountain body.
By applying the method disclosed by the invention to carry out subdata window filtering processing on the gravity data of the Songliao basin and the multiple mountains, linear abnormal information in gravity abnormity is better highlighted, important information is provided for dividing and breaking the application area gravity, and meanwhile, the method is also seen in the processing of the multiple mountains gravity data, not only is the linear gravity abnormity information reflecting the breaking highlighted and enhanced, but also the projection boundary of the three-dimensional geological body on the ground surface is better reflected, the good application effect of the patent is shown, a new processing means is added for gravity data processing, and the method has a very wide application prospect.
It is to be understood that some or all of the various embodiments described above may be included in some embodiments.
As shown in fig. 10, a structural framework diagram is provided for an embodiment of the system for highlighting and enhancing broken-line-shaped gravity anomaly information according to the present invention, the system comprising:
the device comprises a setting unit 1, a data processing unit and a processing unit, wherein the setting unit 1 is used for setting a sliding window for performing sub-data window filtering on a gravity anomaly data area, and dividing the sliding window into a plurality of sub-windows, and the shapes of the sub-windows are different from each other, but the number of included gravity anomaly points is the same;
the computing unit 2 is used for respectively taking the mean value and the mean square error of all the gravity anomaly points of each sub-window, and taking the gravity anomaly mean value of the sub-window with the minimum mean square error as the filtering output of the center point of the sliding window;
the filtering unit 3 is used for scanning the gravity abnormal data area by using a sliding window according to a preset sliding rule and performing sub-data window filtering processing on the data of the gravity abnormal data area;
and the detection unit 4 is used for detecting whether the preset iteration times are reached.
In the embodiment, the sliding window in the sub-data window filtering method is divided into a plurality of sub-windows, then the gravity anomaly mean value of the sub-window with the minimum mean square error in the sub-windows is used as the filtering output of the central point of the sliding window, and the gravity anomaly along the fracture linear anomaly direction is determined and highlighted and enhanced according to the variance of the gravity anomaly in the window with a certain size, so that the purpose of enhancing the fracture gravity anomaly information is achieved, not only is the fracture anomaly information reliably extracted, but also the boundary identification and drawing of the non-fractured geologic body can be carried out, the application field of the invention is enriched, the defects and the defects that the fracture gravity anomaly information is highlighted by various related methods based on the direction derivative are effectively overcome, a new processing means is added for gravity data processing, and the invention has a very wide application prospect.
Optionally, in some possible embodiments, the method further includes:
and the data expansion unit is used for expanding the gravity anomaly data area according to the preset iteration times and the size of the sliding window so as to keep the effective data area consistent with the original data area.
Optionally, in some possible embodiments, the method further includes:
and the export unit is used for saving the obtained processing result as a grid format file of SURFER drawing software, and drawing a sub-data window filter graph based on the reinforced gravity linear abnormal information protruding fracture position by using the SURFER drawing software.
Optionally, in some possible embodiments, the calculating unit 2 is specifically configured to calculate an average value of all gravity anomaly points of each sub-window according to the following formula:
wherein the content of the first and second substances,
is the mean value of gravity anomaly points of the ith sub-window, gi(j) I is the gravity anomaly value of the jth gravity anomaly point in the ith sub-window, I is 1, 2.
Optionally, in some possible embodiments, the calculating unit 2 is specifically configured to calculate the mean square error of all gravity anomaly points of each sub-window according to the following formula:
wherein σiIs the mean square error of the gravity anomaly point in the ith sub-window.
It is to be understood that some or all of the various embodiments described above may be included in some embodiments.
It should be noted that the above embodiments are product embodiments corresponding to previous method embodiments, and for the description of the product embodiments, reference may be made to corresponding descriptions in the above method embodiments, and details are not repeated here.
The reader should understand that in the description of this specification, reference to the description of the terms "one embodiment," "some embodiments," "an example," "a specific example," or "some examples," etc., means that a particular feature, structure, material, or characteristic described in connection with the embodiment or example is included in at least one embodiment or example of the invention. In this specification, the schematic representations of the terms used above are not necessarily intended to refer to the same embodiment or example. Furthermore, the particular features, structures, materials, or characteristics described may be combined in any suitable manner in any one or more embodiments or examples. Furthermore, various embodiments or examples and features of different embodiments or examples described in this specification can be combined and combined by one skilled in the art without contradiction.
In the several embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the above-described method embodiments are merely illustrative, and for example, the division of steps into only one logical functional division may be implemented in practice in another way, for example, multiple steps may be combined or integrated into another step, or some features may be omitted, or not implemented.
The above method, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention essentially or partially contributes to the prior art, or all or part of the technical solution can be embodied in the form of a software product stored in a storage medium and including instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: various media capable of storing program codes, such as a usb disk, a removable hard disk, a Read-only memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk.
While the invention has been described with reference to specific embodiments, it will be understood by those skilled in the art that various changes in form and details may be made therein without departing from the spirit and scope of the invention as defined by the appended claims. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (10)
1. A method for highlighting and enhancing broken linear gravity anomaly information is characterized by comprising the following steps:
step 1, setting a sliding window for performing sub-data window filtering on a gravity anomaly data area, and dividing the sliding window into a plurality of sub-windows, wherein the shapes of the sub-windows are different from each other, but the number of included gravity anomaly points is the same;
step 2, respectively taking the mean value and the mean square error of all gravity anomaly points of each sub-window, and taking the gravity anomaly mean value of the sub-window with the minimum mean square error as the filtering output of the center point of the sliding window;
step 3, scanning the gravity abnormal data area by using the sliding window according to a preset sliding rule, and performing sub-data window filtering processing on the data of the gravity abnormal data area;
and 4, repeating the steps 1 to 3 until the preset iteration times are reached.
2. The method for highlighting and enhancing fractured linear gravity anomaly information according to claim 1, further comprising:
and 5, expanding the gravity abnormal data area according to the preset iteration times and the size of the sliding window, so that the effective data area is consistent with the original data area.
3. The method for highlighting and enhancing fractured linear gravity anomaly information according to claim 1, further comprising:
and 6, storing the obtained processing result as a grid format file of SURFER drawing software, and drawing a sub data window filter graph based on the reinforced gravity linear abnormal information protruding fracture position by using the SURFER drawing software.
4. The method for highlighting and enhancing fractured linear gravity anomaly information according to any one of claims 1 to 3, wherein the mean value of all gravity anomaly points of each sub-window is calculated according to the following formula:
wherein the content of the first and second substances,
is the mean value of gravity anomaly points of the ith sub-window, gi(j) I is the gravity anomaly value of the jth gravity anomaly point in the ith sub-window, I is 1, 2.
5. The method for highlighting and enhancing fractured linear gravity anomaly information according to claim 4, wherein the mean square error of all gravity anomaly points of each sub-window is calculated according to the following formula:
wherein σiIs the mean square error of the gravity anomaly point in the ith sub-window.
6. A system for highlighting and enhancing broken thread-like gravity anomaly information, comprising:
the device comprises a setting unit, a data processing unit and a processing unit, wherein the setting unit is used for setting a sliding window used for carrying out sub-data window filtering on a gravity anomaly data area, and dividing the sliding window into a plurality of sub-windows, and the shapes of the sub-windows are different from each other, but the number of included gravity anomaly points is the same;
the computing unit is used for respectively taking the mean value and the mean square error of all the gravity anomaly points of each sub-window and taking the gravity anomaly mean value of the sub-window with the minimum mean square error as the filtering output of the center point of the sliding window;
the filtering unit is used for scanning the gravity abnormal data area by using the sliding window according to a preset sliding rule and carrying out sub-data window filtering processing on the data of the gravity abnormal data area;
and the detection unit is used for detecting whether the preset iteration times are reached.
7. The system for highlighting and enhancing broken string-like gravity anomaly information according to claim 6, further comprising:
and the data expansion unit is used for expanding the gravity abnormal data area according to the preset iteration times and the size of the sliding window, so that the effective data area is consistent with the original data area.
8. The system for highlighting and enhancing broken string-like gravity anomaly information according to claim 7, further comprising:
and the exporting unit is used for saving the obtained processing result as a grid format file of SURFER drawing software, and drawing a sub-data window filter graph based on the reinforced gravity linear abnormal information protruding fracture position by using the SURFER drawing software.
9. The system according to any one of claims 6 to 8, wherein the computing unit is specifically configured to compute a mean value of all gravity anomaly points of each of the sub-windows according to the following formula:
wherein the content of the first and second substances,
is the mean value of gravity anomaly points of the ith sub-window, gi(j) I is the gravity anomaly value of the jth gravity anomaly point in the ith sub-window, I is 1, 2.
10. The system according to claim 9, wherein the computing unit is specifically configured to compute the mean square error of all gravity anomaly points of each of the sub-windows according to the following formula:
wherein σiIs the mean square error of the gravity anomaly point in the ith sub-window.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011110088.4A CN112347421A (en) | 2020-10-16 | 2020-10-16 | Method and system for highlighting and enhancing broken-line-shaped gravity anomaly information |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011110088.4A CN112347421A (en) | 2020-10-16 | 2020-10-16 | Method and system for highlighting and enhancing broken-line-shaped gravity anomaly information |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112347421A true CN112347421A (en) | 2021-02-09 |
Family
ID=74361004
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011110088.4A Pending CN112347421A (en) | 2020-10-16 | 2020-10-16 | Method and system for highlighting and enhancing broken-line-shaped gravity anomaly information |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112347421A (en) |
Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060133206A1 (en) * | 2004-10-11 | 2006-06-22 | Landmark Graphics Corporation | Fault filter for seismic discontinuity data |
| CN101976441A (en) * | 2010-11-09 | 2011-02-16 | 东华大学 | Method for extracting Sobel operator filtering profile for representing fabric texture and fractal detail mixed characteristic vector |
| CN104122598A (en) * | 2014-07-10 | 2014-10-29 | 山西潞安环保能源开发股份有限公司 | Structure function method extracting fault anomaly from geophysical prospecting gravity anomaly |
| US20140363080A1 (en) * | 2013-06-08 | 2014-12-11 | Fujitsu Limited | Filtering method and device in image processing |
| CN105807323A (en) * | 2016-03-09 | 2016-07-27 | 哈尔滨工程大学 | Method for determining magnetic target position by using small sub-region identification algorithm |
| CN106778466A (en) * | 2016-11-11 | 2017-05-31 | 哈尔滨工程大学 | A kind of undistorted boundary extension method for the online denoising of small echo |
| CN110888144A (en) * | 2019-12-04 | 2020-03-17 | 吉林大学 | Laser radar data synthesis method based on sliding window |
| CN111045089A (en) * | 2019-12-20 | 2020-04-21 | 核工业北京地质研究院 | Method for identifying enhanced heavy magnetic field boundary |
-
2020
- 2020-10-16 CN CN202011110088.4A patent/CN112347421A/en active Pending
Patent Citations (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20060133206A1 (en) * | 2004-10-11 | 2006-06-22 | Landmark Graphics Corporation | Fault filter for seismic discontinuity data |
| CN101976441A (en) * | 2010-11-09 | 2011-02-16 | 东华大学 | Method for extracting Sobel operator filtering profile for representing fabric texture and fractal detail mixed characteristic vector |
| US20140363080A1 (en) * | 2013-06-08 | 2014-12-11 | Fujitsu Limited | Filtering method and device in image processing |
| CN104122598A (en) * | 2014-07-10 | 2014-10-29 | 山西潞安环保能源开发股份有限公司 | Structure function method extracting fault anomaly from geophysical prospecting gravity anomaly |
| CN105807323A (en) * | 2016-03-09 | 2016-07-27 | 哈尔滨工程大学 | Method for determining magnetic target position by using small sub-region identification algorithm |
| CN106778466A (en) * | 2016-11-11 | 2017-05-31 | 哈尔滨工程大学 | A kind of undistorted boundary extension method for the online denoising of small echo |
| CN110888144A (en) * | 2019-12-04 | 2020-03-17 | 吉林大学 | Laser radar data synthesis method based on sliding window |
| CN111045089A (en) * | 2019-12-20 | 2020-04-21 | 核工业北京地质研究院 | Method for identifying enhanced heavy magnetic field boundary |
Non-Patent Citations (4)
| Title |
|---|
| 刘彩云;姚长利;郑元满;: "重力异常分离的小波域优化位变滤波方法", 地球物理学报, no. 12 * |
| 周文纳;杜晓娟;李吉焱;: "改进的角度边界识别方法及其在延吉盆地的应用", 地球科学(中国地质大学学报), no. 06 * |
| 王彦国;黄远生;张瑾;: "基于重力异常增强的改进小子域滤波方法", 石油地球物理勘探, no. 01 * |
| 蔡钟;倪小东;: "重力异常处理中小子域滤波法的改进", 河北工程大学学报(自然科学版), no. 04, pages 1 - 4 * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US8350851B2 (en) | Right sizing reservoir models | |
| CN110529077B (en) | Modeling cross-faults and complex wellbores in reservoir simulation | |
| Zhang et al. | Characterization of coal porosity and permeability evolution by demineralisation using image processing techniques: A micro-computed tomography study | |
| US20160124113A1 (en) | Methods To Handle Discontinuity In Constructing Design Space Using Moving Least Squares | |
| EP3362640B1 (en) | History matching of hydrocarbon production from heterogenous reservoirs | |
| CN111796324A (en) | Seismic all-horizon tracking method and device | |
| González-Nicolás et al. | Pressure management via brine extraction in geological CO2 storage: Adaptive optimization strategies under poorly characterized reservoir conditions | |
| US20180031720A1 (en) | Method and System for Generating a Subsurface Model | |
| US20180348400A1 (en) | Method and System for Creating and Using a Subsurface Model in Hydrocarbon Operations | |
| NO20180871A1 (en) | Hybrid 3d geocellular representation of selected natural fracture network subsets | |
| CN112861890B (en) | Reservoir evaluation model construction method and reservoir identification method | |
| NO345620B1 (en) | Roughening and splitting techniques for generating grids for simulation of hydrocarbon reservoirs | |
| US9542507B2 (en) | Feature detection in seismic volumes | |
| CN112347421A (en) | Method and system for highlighting and enhancing broken-line-shaped gravity anomaly information | |
| Chen et al. | Interaction of reactive fronts during transport in a homogeneous porous medium with initial small non-uniformity | |
| Li et al. | Fracture extraction from FMI based on multiscale mathematical morphology | |
| CN108693559B (en) | Seismic wave combined imaging method and system | |
| US20160054472A1 (en) | Method of modelling a subsurface volume | |
| EP4113117A1 (en) | System and method for correlating oil distribution during drainage and imbibition using machine learning | |
| Kalisch et al. | Rayleigh–Taylor instability of immiscible fluids in porous media | |
| Kumar et al. | Qualitative Dissolution Modeling of Etch‐Pit Formation on the K‐Feldspar Surface Through Phase‐Field Approach | |
| US11353622B2 (en) | Systems and methods for hydrocarbon reservoir three dimensional unstructured grid generation and development | |
| Johnson et al. | Impact of Injection Rate for CO2 Storage Within Sedimentary Basins, a Multidisciplinary Analysis of Focused Fluid Flow | |
| US8660824B2 (en) | Discretely gridding complex fracturing for simulation | |
| Xu et al. | Effects of roughness on density-driven convection of dissolved CO2 in discrete fracture-matrix system |
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 |