A method and system for fine characterization of deep formation interfaces and water-rich areas

Abstract

The present invention discloses a method and system for finely depicting deep formation interfaces and water-rich areas, comprising the following steps: eliminating the static effect from the electromagnetic frequency sounding processing result data to obtain the electromagnetic frequency sounding processing result data after eliminating the static effect; using the electromagnetic frequency sounding processing result data after eliminating the static effect to improve the visualization and layered interpretation effect of the frequency sounding results, and obtaining a result pseudo-seismic profile; establishing an excitation model in the frequency domain, and obtaining the amplitude-frequency characteristic curve and phase-frequency characteristic curve of the induced polarization effect according to the excitation model; extracting polarization parameters according to the amplitude-frequency characteristic curve and phase-frequency characteristic curve, and using the polarization parameters for auxiliary detection to obtain water-rich area data; obtaining a double-frequency amplitude-frequency curve graph according to the water-rich area data, and interpreting the water-rich area by using the double-frequency amplitude-frequency curve graph. The present invention achieves the purpose of finely depicting deep formation interfaces and water-rich areas.

Description

Technical Field

[0001] The present invention relates to the technical field of geological exploration, and particularly relates to a method and system for finely depicting deep formation interfaces and water-rich areas. Background Art

[0002] Among geophysical exploration methods and techniques such as gravity, magnetics, electromagnetics, and seismic exploration, seismic exploration has advantages that cannot be compared with other methods. Especially in sedimentary rock areas, the undulating morphology and structures interpreted from reflection waves are intuitive and clear, the data format, the standardization of processing and interpretation software, and the visualization degree of results are relatively high. In electromagnetic exploration, in order to provide results similar to those of reflection seismic exploration, a lot of work has been done in pseudo-seismic processing, such as pseudo-seismic processing of magnetotelluric sounding, transient electromagnetic sounding and other data, but the actual effect is not ideal, and there is no mature technology put into production.

[0003] In view of the above problems, after formatting the processed results data of electromagnetic frequency sounding into seismic data (SEG-Y format), people use the mature seismic visualization interpretation platform to improve the visualization and stratigraphic interpretation of frequency sounding results, so as to improve the accuracy and intuitiveness of the results in some aspects. Here, the processed results data of electromagnetic frequency sounding is mainly apparent resistivity. However, due to the fact that the apparent resistivity increases with depth and is affected by static displacement, the stratification display effect is not obvious and the auxiliary interpretation effect is small, resulting in unclear interpretation of formation interfaces and water-rich areas. Summary of the Invention

[0004] In order to overcome the deficiencies of the prior art, the present invention provides a method and system for finely depicting deep formation interfaces and water-rich areas, which are used to solve the technical problem that the prior art has unclear interpretation of formation interfaces and water-rich areas as the depth increases, so as to achieve the purpose of finely depicting deep formation interfaces and water-rich areas.

[0005] To solve the above problems, the technical solutions adopted by the present invention are as follows:

[0006] A method for finely depicting deep formation interfaces and water-rich areas, comprising the following steps:

[0007] Eliminate the static effect from the processed results data of electromagnetic frequency sounding to obtain the processed results data of electromagnetic frequency sounding after eliminating the static effect;

[0008] Use the processed results data of electromagnetic frequency sounding after eliminating the static effect to improve the visualization and stratigraphic interpretation effect of frequency sounding results, and obtain a pseudo-seismic profile of the results;

[0009] Establish an excitation model in the frequency domain, and obtain the amplitude-frequency characteristic curve and phase-frequency characteristic curve of the induced polarization effect according to the excitation model;

[0010] Extract polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and use the polarization parameters for auxiliary detection to obtain water-rich area data;

[0011] Obtain a dual-frequency amplitude-frequency curve graph according to the water-rich area data, and use the dual-frequency amplitude-frequency curve graph to interpret the water-rich area.

[0012] As a preferred embodiment of the present invention, when eliminating the static effect of the electromagnetic frequency sounding processing result data, it includes:

[0013] Use the phase data to eliminate the static effect of the electromagnetic frequency sounding processing result data.

[0014] As a preferred embodiment of the present invention, when using the phase data to eliminate the static effect of the electromagnetic frequency sounding processing result data, it includes:

[0015] Obtain corresponding phase data by using the apparent resistivity, and obtain the zero position of the phase and the maximum value position of the phase according to the phase data;

[0016] Obtain an accurate electrical interface according to the zero position of the phase and the maximum value position of the phase.

[0017] As a preferred embodiment of the present invention, when obtaining corresponding phase data by using the apparent resistivity, it includes:

[0018] Obtain the corresponding apparent resistivity according to the electromagnetic wave impedance at a frequency, and use the apparent resistivity to obtain the phase data at the frequency, as specifically shown in Formula 1:

[0019]

[0020] In the formula, ρ s (f) is the apparent resistivity calculated from the electromagnetic wave impedance at frequency f, is the phase data at frequency f.

[0021] As a preferred embodiment of the present invention, when establishing an excitation model in the frequency domain, it includes:

[0022] Establish the excitation model according to the excitation model impedance, charging rate, time constant of the excitation polarization process, and frequency correlation coefficient when the angular frequency is 0.

[0023] As a preferred embodiment of the present invention, the excitation model is specifically shown in Formula 2:

[0024]

[0025] Wherein, Z(iw) is the excitation model impedance at angular frequency w, Z(0) is the excitation model impedance at angular frequency 0, m is the chargeability, τ is the time constant characterizing the induced polarization process, and c is the frequency-dependent coefficient.

[0026] As a preferred embodiment of the present invention, when extracting polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, it includes:

[0027] Preset a polarization frequency band, select two frequencies in the polarization frequency band according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and obtain the impedances at the two frequencies;

[0028] Extract the polarization parameters according to the impedances at the two frequencies.

[0029] As a preferred embodiment of the present invention, when extracting the polarization parameters according to the impedances at the two frequencies, it includes:

[0030] By obtaining the difference between the first impedance and the second impedance, the polarization parameter is extracted, specifically as shown in Equation 3:

[0031] As shown in Equation 3:

[0032]

[0033] Wherein, Z(f L ) is the first impedance, Z(f H ) is the second impedance, and F1 is the polarization parameter;

[0034] Among them, the impedances at the two frequencies are the first impedance and the second impedance respectively, and the frequency corresponding to the second impedance is higher than the frequency corresponding to the first impedance.

[0035] As a preferred embodiment of the present invention, when extracting polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, it includes:

[0036] Preset a polarization frequency band, select two frequencies in the polarization frequency band according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and obtain the current-normalized potentials at the two frequencies;

[0037] Extract the polarization parameters according to the current-normalized potentials at the two frequencies, specifically as shown in Equation 4:

[0038]

[0039] Wherein, U(f L ) is the first current-normalized potential, U(f H ) is the second current-normalized potential, and F2 is the polarization parameter;

[0040] Among them, the current normalized potentials of the two frequencies are the first current normalized potential and the second current normalized potential respectively, and the frequency corresponding to the second current normalized potential is higher than the frequency corresponding to the first current normalized potential.

[0041] A fine characterization system for deep formation interfaces and water-rich areas, comprising:

[0042] An elimination unit: used to eliminate the static effect of the electromagnetic frequency sounding processing result data to obtain the electromagnetic frequency sounding processing result data after eliminating the static effect;

[0043] A profile construction unit: used to improve the visualization and layered interpretation effect of the frequency sounding result by using the electromagnetic frequency sounding processing result data after eliminating the static effect, and obtain a result pseudo-seismic profile;

[0044] An excitation unit: used to establish an excitation model in the frequency domain, and obtain the amplitude-frequency characteristic curve and phase-frequency characteristic curve of the induced polarization effect according to the excitation model;

[0045] An auxiliary detection unit: used to extract polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and perform auxiliary detection by using the polarization parameters to obtain water-rich area data;

[0046] A characterization unit: used to obtain a double-frequency amplitude-frequency curve graph according to the water-rich area data, and interpret the water-rich area by using the double-frequency amplitude-frequency curve graph.

[0047] Compared with the prior art, the beneficial effects of the present invention are as follows:

[0048] (1) The result pseudo-seismic profile obtained by using the method provided by the present invention can clearly show the fault and the Ordovician limestone roof interface compared with the inversion resistivity profile used in traditional interpretation, and the result is intuitive;

[0049] (2) The present invention uses polarization parameters for auxiliary detection, so that high-polarization minerals such as underground (thermal) water, metal sulfide deposits, and graphite deposits can be accurately detected, and it plays a role especially in the exploration of green energy underground hot water.

[0050] The present invention will be further described in detail below with reference to the accompanying drawings and specific embodiments. Description of the Drawings

[0051] Figure 1 - is the inversion resistivity profile of using the wide-area electromagnetic method to detect geothermal energy in a certain place in Shandong in the embodiment of the present invention;

[0052] Figure 2 - is the result pseudo-seismic profile of using the wide-area electromagnetic method to detect geothermal energy in a certain place in Shandong in the embodiment of the present invention;

[0053] Figure 3 - is the apparent resistivity - frequency pseudo - section map of a certain place in the embodiment of the present invention;

[0054] Figure 4 - is the apparent resistivity - frequency pseudo - section map of a certain place in the embodiment of the present invention after static displacement correction;

[0055] Figure 5 - is the inverted resistivity section map of a certain place in the embodiment of the present invention;

[0056] Figure 6 - is the result pseudo - seismic profile map of a certain place in the embodiment of the present invention;

[0057] Figure 7 - is the amplitude - frequency characteristic curve and phase - frequency characteristic curve graph in the embodiment of the present invention;

[0058] Figure 8 - is the amplitude - frequency curve graph extracted by using two frequencies in a certain place in the embodiment of the present invention;

[0059] Figure 9 - is the method step diagram for fine characterization of deep - formation interfaces and water - rich areas in the embodiment of the present invention. Detailed implementation manners

[0060] The method for fine characterization of deep - formation interfaces and water - rich areas provided by the present invention, as Figure 9 shown, includes the following steps:

[0061] Step S1: Eliminate the static effect from the electromagnetic frequency sounding processing result data to obtain the electromagnetic frequency sounding processing result data after eliminating the static effect;

[0062] Step S2: Use the electromagnetic frequency sounding processing result data after eliminating the static effect to improve the visualization and layered interpretation effect of the frequency sounding result, and obtain a result pseudo - seismic profile;

[0063] Step S3: Establish an excitation model in the frequency domain, and obtain the amplitude - frequency characteristic curve and phase - frequency characteristic curve of the induced polarization effect according to the excitation model;

[0064] Step S4: Extract polarization parameters according to the amplitude - frequency characteristic curve and phase - frequency characteristic curve, and use the polarization parameters for auxiliary detection to obtain water - rich area data;

[0065] Step S5: Obtain a dual - frequency amplitude - frequency curve graph according to the water - rich area data, and interpret the water - rich area by using the dual - frequency amplitude - frequency curve graph.

[0066] In the above step S1, when eliminating the static effect from the electromagnetic frequency sounding processing result data, it includes:

[0067] Eliminating the static effect of the electromagnetic frequency sounding processing result data by using phase data.

[0068] Furthermore, when eliminating the static effect of the electromagnetic frequency sounding processing result data by using phase data, it includes:

[0069] Obtaining the corresponding phase data by using apparent resistivity, and obtaining the zero position of the phase and the maximum position of the phase according to the phase data;

[0070] Obtaining an accurate electrical interface according to the zero position of the phase and the maximum position of the phase.

[0071] Even further, when obtaining the corresponding phase data by using apparent resistivity, it includes:

[0072] Obtaining the corresponding apparent resistivity according to the electromagnetic wave impedance at a frequency, and obtaining the phase data at the frequency by using the apparent resistivity, as specifically shown in Formula 1:

[0073]

[0074] In the formula, ρ s (f) is the apparent resistivity calculated from the electromagnetic wave impedance at frequency f, is the phase data at frequency f.

[0075] In the above step S3, when establishing an excitation model in the frequency domain, it includes:

[0076] Establishing an excitation model according to the excitation model impedance, chargeability, time constant of the excitation polarization process, and frequency correlation coefficient when the angular frequency is 0.

[0077] Furthermore, the excitation model is specifically shown in Formula 2:

[0078]

[0079] In the formula, Z(iw) is the excitation model impedance at angular frequency w, Z(0) is the excitation model impedance at angular frequency 0, m is the chargeability, τ is the time constant characterizing the excitation polarization process, and c is the frequency correlation coefficient.

[0080] Even further, the value range of the frequency correlation coefficient is 0 to 1, generally 0.1 to 0.6, and the typical value is 0.25.

[0081] In the above step S4, when extracting polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, it includes:

[0082] Presetting a polarization frequency band, selecting two frequencies in the polarization frequency band according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and obtaining the impedances at the two frequencies;

[0083] Extract polarization parameters based on the impedances at two frequencies.

[0084] Further, when extracting polarization parameters based on the impedances at two frequencies, it includes:

[0085] Extract the polarization parameter by obtaining the difference between the first impedance and the second impedance, as specifically shown in Equation 3:

[0086]

[0087] In the formula, Z(f L ) is the first impedance, Z(f H ) is the second impedance, and F1 is the polarization parameter;

[0088] Among them, the impedances at the two frequencies are the first impedance and the second impedance respectively, and the frequency corresponding to the second impedance is higher than the frequency corresponding to the first impedance.

[0089] In the above step S4, when extracting polarization parameters based on the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, it includes:

[0090] Preset a polarization frequency band, select two frequencies in the polarization frequency band according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and obtain the current normalized potentials at the two frequencies;

[0091] Extract the polarization parameter according to the current normalized potentials at the two frequencies, as specifically shown in Equation 4:

[0092]

[0093] In the formula, U(f L ) is the first current normalized potential, U(f H ) is the second current normalized potential, and F2 is the polarization parameter;

[0094] Among them, the current normalized potentials at the two frequencies are the first current normalized potential and the second current normalized potential respectively, and the frequency corresponding to the second current normalized potential is higher than the frequency corresponding to the first current normalized potential.

[0095] The fine characterization system for deep formation interfaces and water-rich areas provided by the present invention includes: an elimination unit, a profile construction unit, an excitation unit, an auxiliary detection unit, and a characterization unit.

[0096] The elimination unit is used to eliminate the static effect from the electromagnetic frequency sounding processing result data to obtain the electromagnetic frequency sounding processing result data after eliminating the static effect.

[0097] The profile construction unit is used to improve the visualization and stratification interpretation effect of the frequency sounding result by using the electromagnetic frequency sounding processing result data after eliminating the static effect, and obtain a result pseudo-seismic profile.

[0098] The excitation unit is used to establish an excitation model in the frequency domain, and obtain the amplitude-frequency characteristic curve and phase-frequency characteristic curve of the induced polarization effect according to the excitation model.

[0099] The auxiliary detection unit is used to extract polarization parameters according to the amplitude-frequency characteristic curve and phase-frequency characteristic curve, and perform auxiliary detection using the polarization parameters to obtain water-rich area data.

[0100] The characterization unit is used to obtain a dual-frequency amplitude-frequency curve graph according to the water-rich area data, and interpret the water-rich area using the dual-frequency amplitude-frequency curve graph.

[0101] The following embodiments are further descriptions of the present invention, but the scope of the present invention is not limited thereto. Embodiment 1 (Quasi-seismic profile display and interpretation)

[0102] In electromagnetic frequency sounding, as long as quantities related to the horizontal electric field are measured for sounding exploration, such as CSAMT, wide-area electromagnetic, etc., they will be affected by static displacement. In the apparent resistivity profile, the resistivity isolines are vertical, and the sounding curves are difficult to interpret and layer. There are three existing static displacement correction methods: one is to perform theoretical calculations on the effect; the second is to use various processing methods such as filtering and phase integration; the third is to use independent measurement methods without static effects. However, no matter which method is used, it is impossible to correct the static displacement well.

[0103] The basis for the method of using phase data to eliminate static effects in the present invention is that static effects generally do not affect the phase curve, and the phase is related to the derivative of the apparent resistivity with respect to frequency under certain conditions. Therefore, using phase data for static correction can achieve good results. The phase of CSAMT and MT is the phase of the wave impedance, and the phase can be expressed by the above formula 1. In the present invention, the phase is not measured, and there is no wave impedance. The corresponding phase is obtained using the apparent resistivity calculated by the horizontal electric field. The frequency sounding method is a volume exploration method. The obtained apparent resistivity does not represent the true resistivity of a specific depth (geological body) below the measuring point, but is a comprehensive reflection of the electrical properties of each layer of the underground medium near the measuring point and is related to the electromagnetic wave frequency f. It can be seen from the above formula 1 that the phase is obtained from the differential of the apparent resistivity, has higher resolution than the apparent resistivity data, and can obtain a relatively accurate electrical interface. In horizontally layered media, the zero position and maximum value position of the phase roughly correspond to the extreme points of the apparent resistivity curve, which is equivalent to the interface of the electrical layer.

[0104] Figure 1 is the inversion resistivity profile of using the wide-area electromagnetic method to detect geothermal energy in a certain place in Shandong. It explores and interprets the positions of Yishan Branch 1, Yishan Branch 2, Liqiao Fault extending from outside the exploration area to the exploration area, two new faults, the top interface of the Ordovician limestone, and the favorable positions for geothermal energy development. The traditional interpretation uses the inversion resistivity profile ( Figure 1)It reflects the fault quite obviously, but the interpretation of the stratigraphic interface is not clear. Figure 2 It is a pseudo-seismic profile of the results of applying the wide-area electromagnetic method to detect geothermal energy in a certain place in Shandong. It can clearly show the fault and the Ordovician limestone top interface, and the results are intuitive.

[0105] Figure 3 It is the apparent resistivity-frequency pseudo-section map of a certain place, which is composed of Figure 3 It can be seen that the static displacements of measuring points such as 800, 1200, 1700, 2400, and 2600 are obvious. Figure 4 It is the apparent resistivity-frequency pseudo-section map after static displacement correction of the above-mentioned certain place, which is composed of Figure 4 It can be seen that the apparent resistivity isolines after correction basically reflect the stratigraphic undulation, and the display of the northwest shallowness and southeast depth of the high-resistance basement of the Ordovician limestone is more intuitive. Figure 5 It is the inversion resistivity section map of the above-mentioned certain place, Figure 5 The resistivity at different positions in the underground space is given in it. According to the resistivity distribution characteristics, electrical stratification interpretation is carried out. The resistivity of the strata above the elevation of -1700 to -1400m is low, which is interpreted as the bottom boundary of the Cenozoic (Neogene) strata. The elevation of the Ordovician limestone is between -2750 and -2250m, and there is an obvious low-resistance anomaly at point 1500, which is interpreted as a fault. Figure 6 It is the pseudo-seismic profile of the results of the above-mentioned certain place. The reference plane selected for the pseudo-seismic profile of the results is 450m. By comparing Figure 6 and Figure 5 it can be found that with the help of the display of the pseudo-seismic profile of the results, the electrical stratigraphic boundary, faults and other structures can be interpreted more accurately. Three faults are interpreted near the Ordovician limestone interface, including 1 large fault and 2 small faults.

[0106] Example 2 (Pseudo-dual-frequency IP water-richness interpretation)

[0107] The induced polarization effect of the polarized geological body can be represented by an excitation model, and the excitation model in the frequency domain can be represented by the above formula 2. According to formula 2, the amplitude-frequency characteristic curve and phase-frequency characteristic curve of the induced polarization effect can be obtained, as shown in Figure 7 In the electromagnetic frequency sounding of the present invention, by using the impedances Z(f H ) and Z(f L ) at two lower frequencies or the normalized potentials U(f H ) and U(f L ) between MN, polarization parameters similar to the amplitude frequency F are extracted to represent the polarization anomaly, so as to assist in detecting high-polarization minerals such as underground (thermal) water, metal sulfide deposits, and graphite deposits, and particularly play a role in the exploration of green energy underground hot water. Figure 8 It is the amplitude-frequency curve graph extracted by using two frequencies in the above-mentioned certain place, which is composed of Figure 8It can be seen that the amplitude frequency within 800 points is the largest, greater than the average value of 0.46; the amplitude frequency near 2000 points is the second, and there are 3 measurement point values greater than 0.46, which are the sections with relatively strong water-richness.

[0108] The above embodiments are only the preferred embodiments of the present invention, and cannot be used to limit the scope of protection of the present invention. Any non-substantive changes and substitutions made by those skilled in the art on the basis of the present invention fall within the scope of protection required by the present invention.

Claims

1. A fine characterization method for deep formation interfaces and water-rich areas, characterized in that Including the following steps: Eliminate the static effect from the electromagnetic frequency sounding processing result data to obtain the electromagnetic frequency sounding processing result data after eliminating the static effect; Utilize the electromagnetic frequency sounding processing result data after eliminating the static effect to improve the visualization and layered interpretation effect of the frequency sounding result, and obtain a result pseudo-seismic profile; Establish an excitation model in the frequency domain, and obtain the amplitude-frequency characteristic curve and phase-frequency characteristic curve of the induced polarization effect according to the excitation model; Extract polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and perform auxiliary detection using the polarization parameters to obtain water-rich area data; Obtain a double-frequency amplitude-frequency curve graph according to the water-rich area data, and interpret the water-rich area using the double-frequency amplitude-frequency curve graph.

2. The fine characterization method of deep formation interfaces and water-rich areas according to claim 1, characterized in that When eliminating the static effect from the electromagnetic frequency sounding processing result data, it includes: Utilize the phase data to eliminate the static effect of the electromagnetic frequency sounding processing result data.

3. The fine characterization method of deep formation interfaces and water-rich areas according to claim 2, wherein, When utilizing the phase data to eliminate the static effect of the electromagnetic frequency sounding processing result data, it includes: Utilize the apparent resistivity to obtain the corresponding phase data, and obtain the zero position of the phase and the maximum value position of the phase according to the phase data; Obtain an accurate electrical interface according to the zero position of the phase and the maximum value position of the phase.

4. The fine characterization method for deep formation interfaces and water-rich zones according to claim 3, characterized in that When utilizing the apparent resistivity to obtain the corresponding phase data, it includes: Obtain the corresponding apparent resistivity according to the electromagnetic wave impedance at a frequency, and utilize the apparent resistivity to obtain the phase data at the frequency, as specifically shown in Formula 1: where ρ s (f) is the apparent resistivity for calculating the electromagnetic wave impedance at frequency f, is the phase data at frequency f.

5. The fine characterization method of deep formation interfaces and water-rich areas according to claim 1, characterized in that When establishing an excitation model in the frequency domain, it includes: Establish the excitation model according to the excitation model impedance at angular frequency 0, the chargeability, the time constant of the induced polarization process, and the frequency correlation coefficient.

6. The fine characterization method of deep formation interface and water-rich area according to claim 5, characterized in that The excitation model is specifically shown in Formula 2: In the formula, Z(iw) is the excitation model impedance at angular frequency w, Z(0) is the excitation model impedance at angular frequency 0, m is the chargeability, τ is the time constant representing the induced polarization process, and c is the frequency correlation coefficient.

7. The fine characterization method of deep formation interface and water-rich area according to claim 1, characterized in that When extracting polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, it includes: Preset a polarization frequency band, select two frequencies in the polarization frequency band according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and obtain the impedances of the two frequencies; Extract the polarization parameters according to the impedances of the two frequencies.

8. The fine characterization method of deep formation interface and water-rich area according to claim 7, characterized in that, When extracting the polarization parameters according to the impedances of the two frequencies, it includes: Extract the polarization parameters by obtaining the difference between the first impedance and the second impedance, as specifically shown in Formula 3: Wherein, Z(f L ) is the first impedance, Z(f H ) is the second impedance, and F1 is the polarization parameter; Wherein, the impedances of the two frequencies are respectively the first impedance and the second impedance, and the frequency corresponding to the second impedance is higher than the frequency corresponding to the first impedance.

9. The fine characterization method of deep formation interfaces and water-rich areas according to claim 1, wherein When extracting polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, it includes: Preset a polarization frequency band, select two frequencies in the polarization frequency band according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and obtain the current-normalized potentials of the two frequencies; Extract the polarization parameters according to the current-normalized potentials of the two frequencies, as specifically shown in Formula 4: where U(f L ) is the first current-normalized potential, U(f H ) is the second current-normalized potential, and F2 is the polarization parameter; Among them, the current normalized potentials of the two frequencies are the first current normalized potential and the second current normalized potential respectively, and the frequency corresponding to the second current normalized potential is higher than the frequency corresponding to the first current normalized potential.

10. A fine characterization system for deep formation interfaces and water-rich areas, characterized in that, It includes: Elimination unit: used to eliminate the static effect of the electromagnetic frequency sounding processing result data to obtain the electromagnetic frequency sounding processing result data after eliminating the static effect; Profile construction unit: used to improve the visualization and layered interpretation effect of the frequency sounding result by using the electromagnetic frequency sounding processing result data after eliminating the static effect, and obtain a result pseudo-seismic profile; Excitation unit: used to establish an excitation model in the frequency domain, and obtain the amplitude-frequency characteristic curve and phase-frequency characteristic curve of the induced polarization effect according to the excitation model; Auxiliary detection unit: used to extract polarization parameters according to the amplitude-frequency characteristic curve and the phase-frequency characteristic curve, and perform auxiliary detection by using the polarization parameters to obtain water-rich area data; Characterization unit: used to obtain a double-frequency amplitude-frequency curve graph according to the water-rich area data, and interpret the water-rich area by using the double-frequency amplitude-frequency curve graph.

Images