CN112859181B - Local anomaly separation method of wide-area electromagnetic method - Google Patents
Local anomaly separation method of wide-area electromagnetic method Download PDFInfo
- Publication number
- CN112859181B CN112859181B CN202110023985.XA CN202110023985A CN112859181B CN 112859181 B CN112859181 B CN 112859181B CN 202110023985 A CN202110023985 A CN 202110023985A CN 112859181 B CN112859181 B CN 112859181B
- Authority
- CN
- China
- Prior art keywords
- resistivity
- local
- abnormal
- wide
- wide area
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 59
- 238000000926 separation method Methods 0.000 title claims abstract description 28
- 230000002159 abnormal effect Effects 0.000 claims abstract description 52
- 238000004364 calculation method Methods 0.000 claims description 20
- 208000010392 Bone Fractures Diseases 0.000 claims description 15
- 206010017076 Fracture Diseases 0.000 claims description 15
- 230000005856 abnormality Effects 0.000 claims description 13
- 238000011161 development Methods 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 6
- 101000872084 Danio rerio Delta-like protein B Proteins 0.000 claims description 3
- 230000008569 process Effects 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 abstract description 9
- 238000012545 processing Methods 0.000 description 10
- 238000004458 analytical method Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000004888 barrier function Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000010606 normalization Methods 0.000 description 3
- 230000000704 physical effect Effects 0.000 description 3
- 239000011435 rock Substances 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 240000007171 Imperata cylindrica Species 0.000 description 1
- 235000019738 Limestone Nutrition 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000009933 burial Methods 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000006073 displacement reaction Methods 0.000 description 1
- 239000010459 dolomite Substances 0.000 description 1
- 229910000514 dolomite Inorganic materials 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 238000002955 isolation Methods 0.000 description 1
- 239000006028 limestone Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the field of geophysical exploration, in particular to a local anomaly separation method of a wide-area electromagnetic method, which comprises the following steps: s1, acquiring inverted wide area apparent resistivity, and regarding the wide area apparent resistivity as superposition of regional abnormal resistivity and local abnormal resistivity; s2, performing repeated iterative cutting on the wide area apparent resistivity by adopting a dynamic improved interpolation cutting operator, and separating out the regional abnormal resistivity; and S3, calculating the difference between the wide area apparent resistivity and the regional abnormal resistivity, and outputting the difference as the local abnormal resistivity. Compared with the original wide-area electromagnetic apparent resistivity profile, the method can weaken the fuzzy boundary between strata, so that the boundary between adjacent electrical strata is more obvious, and for a geological complex region, the final interpretation result is more accurate and reliable due to abnormal separation, and meanwhile, the inversion imaging precision is supported and improved.
Description
Technical Field
The invention relates to the field of geophysical exploration, in particular to a local anomaly separation method of a wide-area electromagnetic method.
Background
The wide-area electromagnetic method has the advantages of strong anti-interference capability, high working efficiency, high measurement precision and large exploration depth, expands an artificial source observation area, and is very suitable for multi-layer, large-area and large-depth shale gas fine exploration under complex construction and fracture development conditions. However, for complex structures and fracture development zones, the resistivity abnormality often changes irregularly, which brings trouble to the traditional electric method depending on resistivity interpretation and severely restricts the fine interpretation of the electromagnetic method. The actual observed resistance parameters are often due to topography fluctuations, non-uniformity of the underground medium, overlapping of various rocks, crisscross fault fissures, or comprehensive reflections of ore bodies filling them. The field observation actual parameters are usually superposition of local abnormality, regional background and related noise, so that how to effectively separate the local abnormality field, regional background field and noise, and the accuracy and precision of result interpretation are improved, and the method has important research value. Various methods have been proposed for directly suppressing noise in electromagnetic signals or suppressing noise using inversion techniques to improve accuracy (Goubau et al, 1978;Gamble et al, 1979;Egbert,1997;Garcia and Jones,2008;Neukirch and Garcia,2014;Weckmann etal, 2005; gao Jinghuai, 2003,2004; alumbiugh, 1996; jing Jianen, 2012, song Shougen, 1995; weng Ai, 2012; wu Xin, 2015), which all belong to improving electromagnetic imaging accuracy using indirect methods, there is room for improvement in both accuracy and reliability of result interpretation.
Disclosure of Invention
Aiming at the defects of the prior art scheme, the invention provides a local anomaly separation method for a wide-area electromagnetic method, which is used for analyzing inversion imaging data of the wide-area electromagnetic method and directly separating local anomalies of the wide-area electromagnetic from apparent resistivity values of the wide-area electromagnetic method.
In order to achieve the above purpose, the technical scheme adopted by the invention is as follows:
a local anomaly separation method of a wide-area electromagnetic method comprises the following steps:
s1, acquiring inverted wide area apparent resistivity, and regarding the wide area apparent resistivity as superposition of regional abnormal resistivity and local abnormal resistivity;
s2, performing repeated iterative cutting on the wide area apparent resistivity by adopting a dynamic improved interpolation cutting operator, and separating out the regional abnormal resistivity;
and S3, calculating the difference between the wide area apparent resistivity and the regional abnormal resistivity, and outputting the difference as the local abnormal resistivity.
Further, the dynamic improved interpolation operator is:
wherein R (i, j) is the abnormal resistivity value of the region at the point (i, j) after iteration, Z (i, j) is the initial wide area apparent resistivity value at the point (i, j), and R (i, j) calculated by the current iteration is substituted into Z (i, j) calculated by the next iteration;
Nx and Nz the calculation windows are along the x-axis and along the x-axis respectively z The radius of the axis, i is the number of the calculation points in the calculation window, n= (2nx+1) (2nz+1) is the total number of calculation points in the calculation window,to calculate the apparent resistivity average within the window, [ F x (i,j)]Is the x-direction nonlinearity; [ F z (i,j)]Is the z-direction nonlinearity.
Generally, the larger the calculation window is, the larger the range of the regional abnormality reflected by the sliding average value is, the more the iteration times are, the smoother the obtained regional abnormal resistivity is, the more the local abnormal resistivity is obvious, and if the iteration is repeated, the false abnormality caused by fixed interference and random noise can be effectively suppressed, and the local abnormality is highlighted. The local abnormal apparent resistivity is obtained through calculation in the step S3, reflected abnormal information is extracted and analyzed, and then a low-speed abnormal zone which is difficult to find is extracted, so that weak abnormal caused by local fracture, structure and thin layer can be effectively identified.
Further, when b (i, j) is calculated: when DeltaB x When (i, j) =0, ifB (i, j) =1; if it isB (i, j) =0;
when c (i, j) is calculated, when ΔB z When (i, j) =0, ifC (i, j) =1; if->C (i, j) =0.
Further, edge expansion processing is performed on the boundary area, so that each point (M, N) within the boundary within the inner radius range can be subjected to interpolation cutting at the positions (i= -Nx, -nx+1, …, M+nx; j= -Nz, -nz+1, …, N+nz). The area edge-enlarging method includes area field edge-enlarging (regional field method/RFM), cosine edge-enlarging (cosine method/CM), polynomial fitting (polynomial fitting) and minimum curvature edge-enlarging (minimum curvature method/MCM).
Further, the iteration exit condition of the multiple iteration cutting is as follows: lim n→∞ max|R n (i,j)- R n-1 (i, j) ε, wherein R n-1 (i, j) is the region anomaly resistivity value of the previous iteration.
Further, the method is used for eliminating the influence of the underground area field on the local resistivity abnormality, and highlighting the structure and fracture development condition of the local stratum structure.
Furthermore, the method is used for further tracking and describing the thin layer extension based on the construction and fracture of the local stratum structure on the basis of the local abnormal resistivity profile output in the step S3.
A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the method described above when executing the computer program.
A computer readable storage medium storing a computer program for executing the above method.
The invention has the beneficial effects that:
according to the local anomaly separation method of the wide-area electromagnetic method, dynamic improved interpolation cutting operators are adopted to conduct anomaly separation of wide-area electromagnetic inversion resistivity, compared with an original wide-area electromagnetic apparent resistivity profile, fuzzy boundaries between strata can be weakened, boundaries between adjacent electrical layers are obvious, and interpretation results are closer to a ground model; the fault display is clearer, and the fault property and the position are more reliable and accurate; the edge effect on two sides of the end point of the measuring line is suppressed, so that the interpretation result can more truly reflect the actual geophysical model; the deep stratum relief structure form can be more clearly highlighted.
The local anomaly separation method of the wide-area electromagnetic method is based on effective inversion result data, and the local anomaly field information is further extracted and analyzed by effectively compressing the regional background field. For a geological structure complex region, the final interpretation result is more accurate and reliable due to abnormal separation, and meanwhile, the inversion imaging precision is supported and improved.
The local anomaly separation method of the wide-area electromagnetic method can be used as a supplementary means for interpretation of the conventional wide-area electromagnetic method data, can more fully mine implicit geological information in the information, can effectively improve interpretation accuracy and reliability of construction and fracture, not only gives consideration to the advantages of large depth and high precision of the conventional wide-area electromagnetic method, but also can highlight electrical anomalies caused by local geological structures or thin layers on the basis, realizes accurate positioning of fracture structures and thin layers, and forms a set of method and technical system suitable for identifying shale layers in wide-area electromagnetic exploration under complex geological conditions. The technique has important supporting function for identifying shale layers in regions with complicated geological conditions in Chuan nan area and even nationwide, and has wide application prospect.
Drawings
FIG. 1 is a flow chart of a method for local anomaly isolation for wide area electromagnetic method of example 1;
FIG. 2 is a flowchart of the wide area electromagnetic method data processing according to embodiment 1;
FIG. 3 is a conventional resistivity inversion section of a line of example 1;
FIG. 4 is a local anomaly resistivity inversion profile for example 1;
FIG. 5 is a cross-section of a known seismic result interpretation of example 1.
Detailed Description
The present invention will be described in detail with reference to the accompanying drawings.
The present invention will be described in further detail with reference to the drawings and examples, in order to make the objects, technical solutions and advantages of the present invention more apparent. It should be understood that the specific embodiments described herein are for purposes of illustration only and are not intended to limit the scope of the invention.
Example 1
Embodiment 1 provides a local anomaly separation method of a wide area electromagnetic method, as shown in fig. 1, including:
s1, acquiring inverted wide area apparent resistivity, and regarding the wide area apparent resistivity as superposition of regional abnormal resistivity and local abnormal resistivity;
s2, performing repeated iterative cutting on the wide area apparent resistivity by adopting a dynamic improved interpolation cutting operator, and separating out the regional abnormal resistivity;
and S3, calculating the difference between the wide area apparent resistivity and the regional abnormal resistivity, and outputting the difference as the local abnormal resistivity.
The dynamic improved interpolation operator is as follows:
wherein R (i, j) is the abnormal resistivity value of the region at the point (i, j) after iteration, Z (i, j) is the initial wide area apparent resistivity value at the point (i, j), and R (i, j) calculated by the current iteration is substituted into Z (i, j) calculated by the next iteration;
Nx and Nz the calculation windows are along the x-axis and along the x-axis respectively z The radius of the axis, i is the number of the calculation points in the calculation window, n= (2nx+1) (2nz+1) is the total number of calculation points in the calculation window,to calculate the apparent resistivity average within the window, [ F x (i,j)]Is the x-direction nonlinearity; [ F z (i,j)]Is the z-direction nonlinearity.
Generally, the larger the calculation window is, the larger the range of the regional abnormality reflected by the sliding average value is, the more the iteration times are, the smoother the obtained regional abnormal resistivity is, the more the local abnormal resistivity is obvious, and if the iteration is repeated, the false abnormality caused by fixed interference and random noise can be effectively suppressed, and the local abnormality is highlighted. The local abnormal apparent resistivity is obtained through calculation in the step S3, reflected abnormal information is extracted and analyzed, and then a low-speed abnormal zone which is difficult to find is extracted, so that weak abnormal caused by local fracture, structure and thin layer can be effectively identified.
In practical application, in the process of carrying out wide-area electromagnetic inversion resistivity anomaly separation by adopting a dynamic improved interpolation cutting operator, proper calculation window parameters and iteration times are required to be selected by fully combining known geological data and inversion resistivity data so as to achieve relatively accurate separation of regional anomaly resistivity and local anomaly resistivity.
When calculating b (i, j): when DeltaB x When (i, j) =0, ifB (i, j) =1; if->B (i, j) =0;
when c (i, j) is calculated, when ΔB z When (i, j) =0, ifC (i, j) =1; if->C (i, j) =0.
In order to calculate the region anomalies of all points (i, j), the edge-enlarging process is performed on the boundary area, so that each point (M, N) within the boundary within the inner radius range can be subjected to interpolation cutting at the positions (i= -Nx, -nx+1, …, M+nx; j= -Nz, -nz+1, …, N+nz). The area edge-enlarging method includes area field edge-enlarging (regional field method/RFM), cosine edge-enlarging (cosine method/CM), polynomial fitting (polynomial fitting) and minimum curvature edge-enlarging (minimum curvature method/MCM).
The iteration exit condition of the multiple iteration cutting is as follows: lim n→∞ max|R n (i,j)-R n-1 (i, j) ε, wherein R n-1 (i, j) is the region anomaly resistivity value of the previous iteration.
Step S3 calculates the local abnormal resistivity ρ Office (bureau) =Z(i,j)-R n (i,j),R n (i, j) is the region anomaly resistivity value at point (i, j) after the final iteration.
The method is used for eliminating the influence of the underground area field on the local resistivity abnormality, and highlighting the structure and fracture development condition of the local stratum structure.
The method is used for further tracking and describing the thin layer extension based on the construction and fracture of the local stratum structure by taking the local abnormal resistivity profile output in the step S3 as a support.
The specific research application cases of the local anomaly separation method adopting the wide-area electromagnetic method are described below. The formation electrical characteristics of the electrical logging data of the investigation region are shown in the following table.
Statistical table for resistivity of different stratum in research area
Shale has typical low resistance characteristics relative to limestone, siltstone and dolomite from the physical properties of the formation rock. Wherein, the upper stratum cogongrass mouth group-sauvignon group and the lower stratum pagoda group of the target layer dragon one-five peak group show obvious high resistance characteristics and are regarded as high resistance mark layers. According to the physical property analysis, the target layer and the surrounding rock have obvious electrical property difference, and a good physical property basis is provided for carrying out electromagnetic method exploration work.
Data acquisition
A certain line is arranged, the length of which is 35.4km, and the point distance is 50m. In order to achieve the desired detection depth and the requirements of the actual operating schedule, n=7 is adopted, i.e. 2 with 7 frequencies is emitted at one time n A series of pseudo-random signals are received simultaneously with 7 different depth electrical signals. The measurement frequency is 0.011-8192 Hz, 54 frequency points are used, and the method is suitable for detecting a target body ranging from the earth surface to the underground 5 km. The emission current is 100A, the emission voltage is 600-900V, the receiving and transmitting distance is about 20km, and 2 emission sources (A1-B1 and A2-B2) are distributed to cover the L4 line. The measuring device of the horizontal electric dipole source E_Ex wide area electromagnetic method is adopted, wherein the length of the horizontal electric dipole source is changed by 1.5 km-3 km.
Data processing
The data processing flow chart is shown in fig. 2, and mainly comprises the steps of field original data processing, indoor data preprocessing, two-dimensional inversion imaging, local abnormal resistivity separation and the like. The traditional apparent resistivity inversion result is processed by adopting a wide area electromagnetic data processing module in a heavy-magnetic-electricity three-dimensional inversion imaging interpretation integrated system (GME 3 DI-V6.2) which is taught by a wide area electromagnetic subject group Dai Shikun of the university of south China.
(1) And (3) processing the original data in the field, namely performing format conversion, normalization processing (current normalization and device coefficient normalization) on the original data, calculating the wide area apparent resistivity, solving amplitude anomalies of each point and the like, and removing obvious interference frequency point and distortion point data.
(2) Indoor data preprocessing, namely qualitatively analyzing a wide area apparent resistivity curve and a frequency-apparent resistivity section by combining a known electrical structure, and performing smooth filtering, flying spot or flying segment editing, static displacement correction and other processing on single-frequency point or single-channel resistivity data.
(3) And performing two-dimensional inversion imaging, namely performing one-dimensional continuous medium inversion imaging on the processed data, taking the obtained depth-apparent resistivity inversion profile as an initial model, correcting the initial model by combining known geological data or other earth physical data such as well-earthquake and the like, and finally performing two-dimensional wide area electromagnetic two-dimensional inversion imaging processing based on the model to obtain a final inversion resistivity profile.
(4) And (3) separating local abnormal resistivity, namely performing abnormal separation treatment on the reverse resistivity according to the wide-area electromagnetic abnormal separation method provided by the embodiment, eliminating the influence of the underground area field on the local resistivity abnormality, and highlighting the structure and fracture development condition of the local stratum structure.
(5) And comprehensively explaining, namely further tracking and describing the thin layer extension based on the construction and fracture of the local stratum structure by taking the local abnormal separation section as a support. And finally, combining abnormal separation results and a conventional resistivity section to make fine explanation on the combined geology.
Effect analysis
In order to further verify the effectiveness and reliability of the wide-area electromagnetic method anomaly separation method provided by the embodiment, taking a certain measuring line as an example, the distribution characteristics of a geological structure, faults and thin layers revealed by a traditional apparent resistivity section and a local anomaly section after anomaly separation are compared and analyzed, and compared with a known seismic interpretation section for research. Fig. 3 is a conventional apparent resistivity inversion section of the test line, fig. 4 is a local anomaly section after the anomaly separation, fig. 5 is an interpretation section of known seismic achievements, and the comparison analysis results are as follows:
as can be seen from FIG. 3, the electrical layer has a syncline structure characteristic with shallow ends and deep middle, and the electrical layer within 5km has a low-high-low-secondary high-low-middle-resistance-high-resistance change rule from top to bottom. The morphology of the general visual resistivity reflection is basically similar to that of the known seismic achievement interpretation section stratum.
Because of the overlapping of the regional field and the local anomaly field, it is difficult to extract and finely analyze the local information only by interpretation according to the conventional apparent resistivity. Especially fault properties and thin layers are not fine enough, and deep hidden local geological structure information is difficult to be highlighted. As can be seen from FIG. 4, after the influence of the region apparent resistivity field is removed by adopting the anomaly separation technology, the local electrical anomaly characteristics are more obvious, and particularly, the electrical horizons displayed in the region 1 and the region 3 are obviously discontinuous and have obvious electrical horizon dislocation. Due to the influence of volume effect, the local electrical abnormal characteristics of the conventional apparent resistivity inversion section (figure 3) are distributed in a bulk form, and the fault properties are difficult to accurately describe. According to the local abnormal characteristic accumulation deduction reverse faults 9 of the figure 4, the stratum of the Dragon-five peak group is tracked and described, and 1 low-resistance target thin layer is effectively identified. Analysis of the result graph of fig. 4 gave the following recognition:
(1) And (3) extracting a low-resistance thin layer, namely taking a Dragon-five peak group stratum as a target layer of the exploration, wherein the Dragon-five peak group stratum has the characteristics of deep burial and thin thickness, and meanwhile, the resistivity distribution around the fault is disordered and irregular due to the influence of the fault, so that larger misleading is brought to the thin layer division, and fine interpretation is difficult to achieve. The inverse abnormal separation local achievement map 4, whether the structure is a complex region or a simple region, can track the deep low-resistance thin layer according to the local abnormal characteristics, is more accurate in identifying and extracting the low-resistance thin layer information, and finally enables the interpretation achievement to be more accurate and reliable.
(2) The structural morphology reflected in fig. 4 is substantially identical to that revealed by the three-dimensional earthquake (fig. 5).
(3) Fault distribution three-dimensional earthquake (figure 5) reveals that region 1, region 2 and region 3 are distinct complex structural zones, namely fracture development zones. And (3) reversely observing abnormal separation technical results, wherein the stratum area 3 near 600 points is influenced by the faults of the Row 7 and the Row 31, so that local barrier characteristics are formed, and the local barrier characteristics are basically consistent with three-dimensional earthquake results. The stratum area 1 near the 60 points is affected by faults of F1, lamps 2 and lamps 1, obvious pushing and covering features are formed on the earth surface, obvious barrier features are formed in the deep part, the earthquake achievement phase axis data of the anti-observation figure 5 are poor, and the faults of lamps 2 and 1 are weak. Region 2 is shown not in fig. 3, but clearly in fig. 4, presumably as a hidden fault within the layer of interest, substantially coincident with fig. 5.
The foregoing description of the preferred embodiments of the invention is not intended to be limiting, but rather to cover all modifications, equivalents, and alternatives falling within the spirit and principles of the invention.
Claims (8)
1. The local anomaly separation method of the wide-area electromagnetic method is characterized by comprising the following steps of:
s1, acquiring inverted wide area apparent resistivity, and regarding the wide area apparent resistivity as superposition of regional abnormal resistivity and local abnormal resistivity;
s2, performing repeated iterative cutting on the wide area apparent resistivity by adopting a dynamic improved interpolation cutting operator, and separating out the regional abnormal resistivity;
s3, calculating the difference between the wide area apparent resistivity and the regional abnormal resistivity to be used as a local abnormal resistivity output;
the dynamic improved interpolation cutting operator is as follows:
a(i,j)=b(i,j)+c(i,j),
ΔB x (i,j)=Z(i+Nx,j)-Z(i-Nx,j),
ΔB z (i,j)=Z(i,j+Nz)-Z(i,j-Nz);
wherein R (i, j) is the abnormal resistivity value of the region at the point (i, j) after iteration, Z (i, j) is the original wide area apparent resistivity value at the point (i, j), and R (i, j) calculated by the current iteration is substituted into Z (i, j) calculated by the next iteration;
nx and Nz are the radii of the calculation window along the x-axis and z-axis, respectively, l is the calculation point number within the calculation window, n= (2nx+1) (2nz+1) is the total number of calculation points of the calculation window,to calculate the apparent resistivity average within the window, [ F x (i,j)]Is the x-direction nonlinearity; [ F z (i,j)]Is the z-direction nonlinearity.
2. The method according to claim 1, characterized in that: when calculating b (i, j): when DeltaB x When (i, j) =0, ifB (i, j) =1; if->B (i, j) =0;
when c (i, j) is calculated, when ΔB z When (i, j) =0, ifC (i, j) =1; if->C (i, j) =0.
3. The method according to claim 1, characterized in that: the boundary region is subjected to an edge-enlarging process, so that points within the boundary within the range of the radius (M, N) at (i= -Nx, -Nx +1, the M is equal to Nx; j= -Nz, -nz+1, ·, n+nz) can be interpolated cut.
4. A method according to claim 3, characterized in that: the iteration exit condition of the multiple iteration cutting is as follows:wherein R is n-1 (i, j) is the region anomaly resistivity value of the previous iteration.
5. The method according to any one of claims 1-4, wherein: the method is used for eliminating the influence of the underground area field on the local resistivity abnormality and highlighting the structure and fracture development condition of the local stratum structure.
6. The method according to any one of claims 1-4, wherein: the method is used for further tracking and describing the thin layer extension based on the construction and fracture of the local stratum structure by taking the local abnormal resistivity profile output in the step S3 as a support.
7. A computer device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the method of any of claims 1 to 6 when executing the computer program.
8. 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 6.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110023985.XA CN112859181B (en) | 2021-01-08 | 2021-01-08 | Local anomaly separation method of wide-area electromagnetic method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110023985.XA CN112859181B (en) | 2021-01-08 | 2021-01-08 | Local anomaly separation method of wide-area electromagnetic method |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112859181A CN112859181A (en) | 2021-05-28 |
CN112859181B true CN112859181B (en) | 2024-04-05 |
Family
ID=76005504
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110023985.XA Active CN112859181B (en) | 2021-01-08 | 2021-01-08 | Local anomaly separation method of wide-area electromagnetic method |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112859181B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113627051B (en) * | 2021-07-23 | 2024-01-30 | 中国地质科学院地球物理地球化学勘查研究所 | Gravity anomaly field separation method, system, storage medium and electronic equipment |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104375196A (en) * | 2014-11-10 | 2015-02-25 | 中国石油天然气集团公司 | Method and device for extracting three-dimensional MT weak resistivity abnormity |
JP2016085114A (en) * | 2014-10-26 | 2016-05-19 | 滋樹 水谷 | Gravity and magnetic data overall analysis method useful for geothermal resource evaluation |
CN108873085A (en) * | 2018-08-02 | 2018-11-23 | 中国石油天然气集团有限公司 | A kind of sedimentary rock thin-layer electric magnetic data processing method and device |
CN111708094A (en) * | 2020-07-20 | 2020-09-25 | 核工业北京地质研究院 | Sandstone-type uranium ore sand body identification method and system based on wide-area electromagnetic method |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20150260874A1 (en) * | 2014-03-13 | 2015-09-17 | Schlumberger Technology Corporation | System and Method for Imaging Subterranean Formations |
-
2021
- 2021-01-08 CN CN202110023985.XA patent/CN112859181B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2016085114A (en) * | 2014-10-26 | 2016-05-19 | 滋樹 水谷 | Gravity and magnetic data overall analysis method useful for geothermal resource evaluation |
CN104375196A (en) * | 2014-11-10 | 2015-02-25 | 中国石油天然气集团公司 | Method and device for extracting three-dimensional MT weak resistivity abnormity |
CN108873085A (en) * | 2018-08-02 | 2018-11-23 | 中国石油天然气集团有限公司 | A kind of sedimentary rock thin-layer electric magnetic data processing method and device |
CN111708094A (en) * | 2020-07-20 | 2020-09-25 | 核工业北京地质研究院 | Sandstone-type uranium ore sand body identification method and system based on wide-area electromagnetic method |
Non-Patent Citations (2)
Title |
---|
广域电磁法中垂直磁场分量的分析与应用;陈卫营;薛国强;;物探与化探;第39卷(第02期);147-150 * |
相对视电阻率异常的应用;孙中任;赵雪娟;滕寿仁;赵维俊;;地球物理学进展(第06期);334-339 * |
Also Published As
Publication number | Publication date |
---|---|
CN112859181A (en) | 2021-05-28 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
AU742887B2 (en) | Hydrocarbon edge detection using seismic amplitude | |
Liu et al. | Three-dimensional seismic ahead-prospecting method and application in TBM tunneling | |
CN104267442B (en) | A kind of transient electromagnetic for underground coal mine intends seismic detection method | |
Juhlin et al. | Reflection seismic studies over the end-glacial Burträsk fault, Skellefteå, Sweden | |
CN110320576A (en) | Method for forecasting advanced geology, device and storage medium | |
CN112859181B (en) | Local anomaly separation method of wide-area electromagnetic method | |
Gammaldi et al. | Fluid accumulation zone by seismic attributes and amplitude versus offset analysis at Solfatara Volcano, Campi Flegrei, Italy | |
Lv et al. | Application of GPR reverse time migration in tunnel lining cavity imaging | |
Martorana et al. | Assessment of a geological model by surface wave analyses | |
CN112433253A (en) | Three-dimensional detection method for surrounding rock tunnel with weak broken flowing water | |
Ross et al. | An algorithm for automated identification of fault zone trapped waves | |
Gordeev et al. | Detection of geodynamic activity areas based on the Earth's electromagnetic noise parameters | |
Sloan et al. | Depth estimation of voids using backscattered surface waves | |
Faleide et al. | Exploring seismic detection and resolution thresholds of fault zones and gas seeps in the shallow subsurface using seismic modelling | |
Rapetsoa et al. | Cost‐effective in‐mine seismic experiments to image platinum deposits and associated geological structures at Maseve platinum mine, South Africa | |
Schwenk* et al. | Correlation of the backscatter analysis of surface waves method (BASW) for anomaly detection | |
CN103984013A (en) | Wavelet domain pre-stack seismic trace set absorption attenuation parameter estimation algorithm | |
CN111983721A (en) | Sandstone uranium ore mud-sand-mud geological structure identification method and system | |
CN113835123B (en) | Seismic acquisition parameter analysis method based on geological target pre-stack migration imaging | |
Ma et al. | Multi-level ultra-deep fault-controlled karst reservoirs characterization methods for the Shunbei field | |
Polychronopoulou et al. | Integrating earthquake-based passive seismic methods in mineral exploration: Case study from the Gerolekas bauxite mining area, Greece | |
Amjadi et al. | An application of magnetotelluric data inversion in a stratovolcano region | |
CN113064204B (en) | Borehole seismic absorption parameter extraction and stratum lithology detection method | |
Meili et al. | Shallow-well microseismic monitoring and data processing | |
Al-Enezi et al. | 3D seismic reservoir characterization of Marrat Formation in Central Kuwait |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Country or region after: Zhong Guo Address after: No.39 Qinghua Road, Qingyang District, Chengdu, Sichuan 610072 Applicant after: Sichuan Natural Resources Investment Group Geophysical Exploration Institute Co.,Ltd. Address before: No.39 Qinghua Road, Qingyang District, Chengdu, Sichuan 610072 Applicant before: SICHUAN ZHONGCHENG COALFIELD GEOPHYSICAL ENGINEERING INSTITUTE CO.,LTD. Country or region before: Zhong Guo |
|
CB02 | Change of applicant information | ||
GR01 | Patent grant | ||
GR01 | Patent grant |