Rock induced polarization parameter inversion method and device, electronic equipment and medium
Through the inversion method of rock excitation polarization parameters based on the MGEMTIP model, the problem that traditional models cannot effectively characterize the characteristics of rock excitation polarization, and more accurate interpretation of excitation reservoirs and electrical data processing are achieved.
Patent Information
- Application Number
- CN202311586043.8
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2023-11-24
- Publication Date
- 2025-05-27
AI Technical Summary
The traditional Cole-Cole model cannot effectively characterize the excitation polarization characteristics of underground rocks, and cannot describe the multipolarization mechanism, cannot adapt to high-dimensional problems, and lacks quantitative relationships. между excitation parameters and structural properties.
Based on the MGEMTIP model, by establishing the relationship between resistivity dispersion, zero-frequency resistivity, polarization rate, and time constant, combining experimental data to judge the excitation frequency band, calculate the polarization spectrum, and iteratively adjust the initial value of the true conductivity until specific conditions are met, the true conductivity, IP polarization rate and IP equivalent time constant are then calculated.
Effectively remove the low-frequency polarization caused by electrode polarization, accurately extract the characteristics of rock excitation polarization, and provide a theoretical basis and electrical data processing method for the interpretation of excitation reservoirs.
Smart Images
Figure CN120044632A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to the field of oil and natural gas exploration and development, and more specifically, to a rock induced polarization parameter inversion method, device, electronic equipment and medium. Background Art
[0002] The MGEMTIP model is a unified and effective representation of the multipolarization mechanism, and a complex resistivity model that constructs induced polarization parameters and structural parameters. It can effectively represent the induced polarization characteristics of underground rocks. The IP parameters obtained by inversion based on the MGEMTIP model are beneficial to the interpretation of electromagnetic exploration reservoirs. The traditional Cole-Cole model has the defects of lack of quantitative relationship between induced polarization parameters and structural properties, inability to adapt to high-dimensional problems, and inability to describe the multipolarization mechanism.
[0003] Therefore, it is necessary to develop a rock induced polarization parameter inversion method, device, electronic equipment and medium based on the MGEMTIP model.
[0004] The information disclosed in the background technology section of the present invention is only intended to deepen the understanding of the general background technology of the present invention, and should not be regarded as acknowledging or suggesting in any form that the information constitutes the prior art already known to those skilled in the art. Summary of the invention
[0005] The present invention proposes a method, device, electronic device and medium for inversion of rock induced polarization parameters, and its MGEMTIP model is more suitable for the complex polarization characterization after the superposition of multiple polarization mechanisms of actual rocks. The corresponding inversion algorithm focuses on the influence of electrode effect on rock testing during the experiment, which will effectively suppress or remove the low-frequency polarization caused by electrode polarization, and more accurately and effectively extract the characteristics of rock induced polarization, providing an effective theoretical basis and electrical data processing method for induced polarization reservoir interpretation.
[0006] In a first aspect, an embodiment of the present disclosure provides a method for inverting rock induced polarization parameters, comprising:
[0007] The MGEMTIP model is established based on resistivity dispersion, zero-frequency resistivity, polarizability, and time constant;
[0008] Determine the IP frequency band based on experimental data, determine the high frequency band, IP frequency band and low frequency band, and input the initial value of true conductivity;
[0009] Calculate polarizability spectra;
[0010] Determine whether the low frequency band of the polarizability spectrum satisfies η H =0; if not, adjust the initial value of true conductivity and repeatedly calculate the polarizability spectrum until η is satisfied. H =0;
[0011] The effective polarizability spectrum is obtained, and then the true conductivity σ is calculated according to the MGEMTIP model. 0 M 0 , IP polarizability and IP equivalent time constant.
[0012] As a specific implementation of the embodiment of the present disclosure, the MGEMTIP model is:
[0013]
[0014]
[0015]
[0016]
[0017] Among them, ρ 0 , l Represent the resistivity of the background medium and the first disturbance medium, a l is the spherical radius of the lth disturbance medium, f l is the volume content of the first disturbance medium, α l is the surface polarization coefficient, ω=2πf is the circular frequency, C l is the frequency correlation coefficient.
[0018] As a specific implementation of the embodiment of the present disclosure, the initial value of true conductivity is:
[0019] σ 0 M 0 =bσ e (min(ω m )).
[0020] As a specific implementation of the embodiment of the present disclosure, the polarizability spectrum is:
[0021]
[0022] Among them, η l is the polarizability spectrum,
[0023] As a specific implementation of the embodiment of the present disclosure, the initial value of the true conductivity is adjusted as follows:
[0024]
[0025] in, ΓH is the high time constant range set corresponding to low frequency, η l For different time constants τ l The corresponding polarizability.
[0026] As a specific implementation of the embodiment of the present disclosure, the IP polarization rate is:
[0027]
[0028] As a specific implementation of the embodiment of the present disclosure, the IP equivalent time constant is:
[0029]
[0030] In a second aspect, the embodiments of the present disclosure further provide a rock induced polarization parameter inversion device, comprising:
[0031] Modeling module, which establishes MGEMTIP model based on resistivity dispersion, zero-frequency resistivity, polarizability, and time constant;
[0032] The assignment module determines the IP frequency band based on experimental data, determines the high frequency band, IP frequency band and low frequency band, and inputs the initial value of true conductivity;
[0033] A first calculation module calculates a polarizability spectrum;
[0034] Iteration module, determining whether the low frequency band of the polarizability spectrum satisfies η H =0; if not, adjust the initial value of true conductivity and repeatedly calculate the polarizability spectrum until η is satisfied. H =0;
[0035] The second calculation module obtains the effective polarizability spectrum and then calculates the true conductivity σ according to the MGEMTIP model 0 M 0 , IP polarizability and IP equivalent time constant.
[0036] As a specific implementation of the embodiment of the present disclosure, the MGEMTIP model is:
[0037]
[0038]
[0039]
[0040]
[0041] Among them, ρ 0 , l Represent the resistivity of the background medium and the first disturbance medium, a l is the spherical radius of the lth disturbance medium, f l is the volume content of the first disturbance medium, α l is the surface polarization coefficient, ω=2πf is the circular frequency, C l is the frequency correlation coefficient.
[0042] As a specific implementation of the embodiment of the present disclosure, the initial value of true conductivity is:
[0043] σ 0 M 0 =bσ e (min(ω m )).
[0044] As a specific implementation of the embodiment of the present disclosure, the polarizability spectrum is:
[0045]
[0046] Among them, η l is the polarizability spectrum,
[0047] As a specific implementation of the embodiment of the present disclosure, the initial value of the true conductivity is adjusted as follows:
[0048]
[0049] in, Γ H is the high time constant range set corresponding to low frequency, η l For different time constants τ l The corresponding polarizability.
[0050] As a specific implementation of the embodiment of the present disclosure, the IP polarization rate is:
[0051]
[0052] As a specific implementation of the embodiment of the present disclosure, the IP equivalent time constant is:
[0053]
[0054] In a third aspect, an embodiment of the present disclosure further provides an electronic device, the electronic device comprising:
[0055] A memory storing executable instructions;
[0056] A processor runs the executable instructions in the memory to implement the rock induced polarization parameter inversion method.
[0057] In a fourth aspect, an embodiment of the present disclosure further provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the rock induced polarization parameter inversion method is implemented.
[0058] Its beneficial effects are:
[0059] The present invention takes into account the interference of macroscopic factors such as rock structure and test electrodes on the complex resistivity of rock, effectively distinguishes the polarization characteristics under different induced polarization mechanisms, is conducive to the quantitative analysis of rock induced polarization and structural parameters, and provides a theoretical basis and analysis tools for further induced polarization reservoir interpretation.
[0060] The methods and apparatus of the present invention have other features and advantages that will be apparent from, or will be described in detail in, the accompanying drawings and subsequent detailed descriptions incorporated herein, which together serve to explain the specific principles of the present invention. BRIEF DESCRIPTION OF THE DRAWINGS
[0061] The above and other objects, features and advantages of the present invention will become more apparent through a more detailed description of exemplary embodiments of the present invention in conjunction with the accompanying drawings, wherein like reference numerals generally represent like components throughout the exemplary embodiments of the present invention.
[0062] Figure 1 A flow chart showing the steps of a rock induced polarization parameter inversion method according to an embodiment of the present invention.
[0063] Figure 2 A schematic diagram showing rock phases under different temperature and mineralogy conditions according to an embodiment of the present invention is shown.
[0064] Figure 3 A schematic diagram showing the comparison between the complex resistivity test curve of a shaly sandstone sample and suppressed low-frequency inversion data according to an embodiment of the present invention is shown.
[0065] Figure 4 A schematic diagram showing the comparison between the complex resistivity test curve of a fractured shale sample according to an embodiment of the present invention and the unsuppressed low-frequency inversion data is shown.
[0066] Figure 5 A block diagram of a rock induced polarization parameter inversion device according to an embodiment of the present invention is shown.
[0067] Description of reference numerals:
[0068] 201, modeling module; 202, assignment module; 203, first calculation module; 204, iteration module; 205, second calculation module. DETAILED DESCRIPTION
[0069] The preferred embodiments of the present invention will be described in more detail below. Although the preferred embodiments of the present invention are described below, it should be understood that the present invention can be implemented in various forms and should not be limited to the embodiments set forth herein.
[0070] To facilitate understanding of the solutions and effects of the embodiments of the present invention, six specific application examples are given below. Those skilled in the art should understand that the examples are only for facilitating understanding of the present invention, and any specific details thereof are not intended to limit the present invention in any way.
[0071] Example 1
[0072] Figure 1 A flow chart showing the steps of a rock induced polarization parameter inversion method according to an embodiment of the present invention.
[0073] like Figure 1 As shown in FIG. 1 , the rock induced polarization parameter inversion method includes: step 101, establishing an MGEMTIP model according to resistivity dispersion, zero-frequency resistivity, polarizability, and time constant; step 102, judging the induced polarization frequency band based on experimental data, determining the high-frequency band, the induced polarization frequency band, and the low-frequency band, and inputting the initial value of true conductivity; step 103, calculating the polarizability spectrum; step 104, judging whether the low-frequency band of the polarizability spectrum satisfies η H =0; if not, adjust the initial value of true conductivity and repeat the calculation of polarizability spectrum until η is satisfied H = 0; Step 105, obtain the effective polarizability spectrum, and then calculate the true conductivity σ according to the MGEMTIP model 0 M 0 , IP polarizability and IP equivalent time constant.
[0074] In one example, the MGEMTIP model is:
[0075]
[0076]
[0077]
[0078]
[0079] Among them, ρ 0 , l Represent the resistivity of the background medium and the first disturbance medium, a l is the spherical radius of the lth disturbance medium, f l is the volume content of the first disturbance medium, α l is the surface polarization coefficient, ω=2πf is the circular frequency, C l is the frequency correlation coefficient.
[0080] In one example, the initial true conductivity value is:
[0081] σ 0 M 0 =bσ e (min(ω m )).
[0082] In one example, the polarizability spectrum is:
[0083]
[0084] Among them, η l is the polarizability spectrum,
[0085] In one example, the initial value of true conductivity is adjusted as follows:
[0086]
[0087] in, Γ H is the high time constant range set corresponding to low frequency, η l For different time constants τ l The corresponding polarizability.
[0088] In one example, the IP polarization ratio is:
[0089]
[0090] In one example, the IP equivalent time constant is:
[0091]
[0092] Figure 2 A schematic diagram showing rock phases under different temperature and mineralogy conditions according to an embodiment of the present invention is shown.
[0093] Specifically, the low-frequency polarization analysis in experimental measurement. The electrode effect will produce low-frequency polarization during the test, interfering with the accurate measurement of rock complex resistivity. The electrode material, temperature, pressure, rock solution mineralogy and other factors will affect the low-frequency induced polarization characteristics of the rock. Especially in high temperature and low salinity environments, the polarization caused by the electrode will be significantly increased, such as Figure 2 shown.
[0094] Table 1 shows the MGEMTIP model.
[0095] Table 1
[0096]
[0097] Among them, ρ 0 , l Represent the resistivity of the background medium and the first disturbance medium, a l is the spherical radius of the lth disturbance medium, f l is the volume content of the first disturbance medium, α l is the surface polarization coefficient, ω=2πf is the circular frequency, C l is the frequency correlation coefficient, C l =1 corresponds to the Debye model, C l The case of <1 can be approximated by a series of linear combinations of Debye models.
[0098] C l =1 The corresponding conductivity form is:
[0099]
[0100] in, A series of τ l Covering from low frequency to high frequency, effectively covering the experimental frequency band of data measurement (0.01~10 4 Hz). The model parameters can be obtained by Debye decomposition (DD) inversion method.
[0101] remember The linear relationship can be obtained:
[0102]
[0103] According to the measured frequency band data, a series of data corresponding to y(ω m ), we can get l The linear equations with unknowns are:
[0104] A = Y
[0105] in, a l ′(ω m )=re(a l (ω m ));a l ″(ω m )=Im(a l (ω m ));y′(ω m )=re(y(ω m ));y″(ω m )=Im(y(ω m )). By solving the system of equations A TAη=A T Y, constraint η l >0,l=1,2…N can obtain a corresponding series of τ l and η l .
[0106] In the actual calculation process, the initial value setting is analyzed in combination with the test environment and rock structure analysis. If suppression is not required, the DC measurement result can be selected. If low-frequency polarization needs to be removed or suppressed, the amplitude σ 0 M 0 The initial value is a lower value, or bσ is assigned iteratively e (min(ω m )), where b < 1, the effective DC conductivity is obtained by removing the low-frequency component (high time constant)
[0107]
[0108] in, The calculated time constant spectrum τ l and its corresponding polarizability η l , the independent parameters equivalent to the Cole-Cole model can be extracted, but it needs to be within a fixed frequency range, that is, the time constant spectrum interval Γ corresponding to the excitation polarization frequency considered in the experiment is defined IP =(τ L ,τ H ). Thus, the corresponding induced polarization equivalent polarizability and equivalent time constant are obtained:
[0109]
[0110] Combined with the DC conductivity σ after removing the low-frequency component 0 M 0 =σ e (ω m ) / (1-η L ) to obtain the rock IP parameters.
[0111] The specific process of the present invention is:
[0112] Combine the experimental data to determine the IP frequency band, determine the high frequency band, IP frequency band and low frequency band, and input the initial value of true conductivity bσ according to the data e (min(ω m ));
[0113] Combination Calculate the polarizability spectrum η l ;
[0114] Investigate whether the low frequency band of the polarizability spectrum (corresponding to the high time constant) satisfies η H=0; if not satisfied, adjust the initial value coefficient to Repeat the calculation of the polarizability spectrum until η is satisfied H =0;
[0115] Obtain the effective polarizability spectrum η l , and obtain the corresponding true conductivity σ 0 M 0 、IP polarizability and IP equivalent time constant
[0116] Example 2
[0117] The present invention also provides a rock induced polarization parameter inversion device, comprising:
[0118] Modeling module, which establishes MGEMTIP model based on resistivity dispersion, zero-frequency resistivity, polarizability, and time constant;
[0119] The assignment module determines the IP frequency band based on experimental data, determines the high frequency band, IP frequency band and low frequency band, and inputs the initial value of true conductivity;
[0120] A first calculation module calculates a polarizability spectrum;
[0121] Iteration module, to determine whether the low frequency band of the polarizability spectrum satisfies η H =0; if not, adjust the initial value of true conductivity and repeat the calculation of polarizability spectrum until η is satisfied H =0;
[0122] The second calculation module obtains the effective polarizability spectrum and then calculates the true conductivity σ according to the MGEMTIP model 0 M 0 , IP polarizability and IP equivalent time constant.
[0123] In one example, the MGEMTIP model is:
[0124]
[0125]
[0126]
[0127]
[0128] Among them, ρ 0 , l Represent the resistivity of the background medium and the first disturbance medium, a l is the spherical radius of the lth disturbance medium, f l is the volume content of the first disturbance medium, α l is the surface polarization coefficient, ω=2πf is the circular frequency, C l is the frequency correlation coefficient.
[0129] In one example, the initial true conductivity value is:
[0130] σ 0 M 0 =bσ e (min(ω m )).
[0131] In one example, the polarizability spectrum is:
[0132]
[0133] Among them, η l is the polarizability spectrum,
[0134] In one example, the initial value of true conductivity is adjusted as follows:
[0135]
[0136] in, Γ H is the high time constant range set corresponding to low frequency, η l For different time constants τ l The corresponding polarizability.
[0137] In one example, the IP polarization ratio is:
[0138]
[0139] In one example, the IP equivalent time constant is:
[0140]
[0141] Specifically, the low-frequency polarization analysis in experimental measurement. The electrode effect will produce low-frequency polarization during the test, interfering with the accurate measurement of rock complex resistivity. The electrode material, temperature, pressure, rock solution mineralogy and other factors will affect the low-frequency induced polarization characteristics of the rock. Especially in high temperature and low salinity environments, the polarization caused by the electrode will be significantly increased, such as Figure 2 shown.
[0142] Table 1 shows the MGEMTIP model. l =1 corresponds to the Debye model, C l The case of <1 can be approximated by a series of linear combinations of Debye models.
[0143] C l=1 The corresponding conductivity form is:
[0144]
[0145] in, A series of τ l Covering from low frequency to high frequency, effectively covering the experimental frequency band of data measurement (0.01~10 4 Hz). The model parameters can be obtained by Debye decomposition (DD) inversion method.
[0146] remember The linear relationship can be obtained:
[0147]
[0148] According to the measured frequency band data, a series of data corresponding to y(ω m ), we can get l The linear equations with unknowns are:
[0149] A = Y
[0150] in, a l ′(ω m )=re(a l (ω m ));a l ″(ω m )=Im(a l (ω m ));y′(ω m )=re(y(ω m ));y″(ω m )=Im(y(ω m )). By solving the system of equations A T Aη=A T Y, constraint η l >0,l=1,2…N can obtain a corresponding series of τ l and η l .
[0151] In the actual calculation process, the initial value setting is analyzed in combination with the test environment and rock structure analysis. If suppression is not required, the DC measurement result can be selected. If low-frequency polarization needs to be removed or suppressed, the amplitude σ 0 M 0 The initial value is a lower value, or bσ is assigned iteratively e (min(ω m )), where b < 1, the effective DC conductivity is obtained by removing the low-frequency component (high time constant)
[0152]
[0153] in, The calculated time constant spectrum τ l and its corresponding polarizability η l , the independent parameters equivalent to the Cole-Cole model can be extracted, but it needs to be within a fixed frequency range, that is, the time constant spectrum interval Γ corresponding to the excitation polarization frequency considered in the experiment is defined IP =(τ L ,τ H ). Thus, the corresponding induced polarization equivalent polarizability and equivalent time constant are obtained:
[0154]
[0155] Combined with the DC conductivity σ after removing the low-frequency component 0 M 0 =σ e (ω m ) / (1-η L ) to obtain the rock IP parameters.
[0156] The specific process of the present invention is:
[0157] Combine the experimental data to determine the IP frequency band, determine the high frequency band, IP frequency band and low frequency band, and input the initial value of true conductivity bσ according to the data e (min(ω m ));
[0158] Combination Calculate the polarizability spectrum η l ;
[0159] Investigate whether the low frequency band of the polarizability spectrum (corresponding to the high time constant) satisfies η H =0; if not satisfied, adjust the initial value coefficient to Repeat the calculation of the polarizability spectrum until η is satisfied H =0;
[0160] Obtain the effective polarizability spectrum η l , and obtain the corresponding true conductivity σ 0 M 0 、IP polarizability and IP equivalent time constant
[0161] Example 3
[0162] The rock complex resistivity curve is characterized based on the MGEMTIP model. The influence of rock temperature, pressure, electrode, rock macro cracks and other factors on low-frequency polarization are comprehensively analyzed, and the factors affecting low-frequency polarization of rock complex resistivity are comprehensively analyzed to remove or suppress low-frequency polarization caused by electrode effect (not related to the rock itself). Comprehensively analyze the rock test environment and structure, perform rapid inversion of MGEMTIP model induced polarization parameters, and obtain effective induced polarization parameters such as true resistivity, polarizability, and time constant.
[0163] The specific plan is explained in combination with rock experiments and analysis.
[0164] A shale sample with a diameter of 2.5 cm was obtained by sampling the core of a well in a certain area, and the rock complex resistivity curve was obtained through experimental measurement.
[0165] Figure 3 A schematic diagram showing the comparison between the complex resistivity test curve of a shaly sandstone sample and suppressed low-frequency inversion data according to an embodiment of the present invention is shown.
[0166] Test sample 1 is argillaceous sandstone, and the test environment is: 5% NaCl solution, temperature is 120℃, and pressure is 40MPa. According to the test environment analysis, the rock particle diameter is less than 0.1mm, and the polarization time constant is not more than 0.1s. The analysis shows that the low-frequency phase rise is an electrode effect, which can be suppressed by the algorithm during the inversion process. Figure 3 As shown, the low-frequency polarization of the sample is mainly caused by the external test electrode, and the inversion results are: true resistivity 112.57Ω·m, polarization rate 56.47%, and time constant 36.57s.
[0167] Figure 4 A schematic diagram showing the comparison between the complex resistivity test curve of a fractured shale sample according to an embodiment of the present invention and the unsuppressed low-frequency inversion data is shown.
[0168] Test sample 2 is fractured shale, and the test environment is: 5% NaCl solution, temperature 40℃, pressure 10Mpa. According to the test environment analysis, the fracture scale can reach the sample scale, the polarization time constant can be higher than 0.01s, the test temperature is low, the electrode effect is not obvious, and there is no need to suppress low-frequency polarization, such as Figure 4 As shown in the figure, it is believed that low-frequency polarization is an inherent feature of rocks formed by cracks, and the inversion data basically retains low-frequency polarization. The inversion results are: true resistivity 29.71Ω·m, polarizability 5.80%, and time constant 0.01s.
[0169] Example 4
[0170] Figure 5 A block diagram of a rock induced polarization parameter inversion device according to an embodiment of the present invention is shown.
[0171] like Figure 5 As shown, the rock induced polarization parameter inversion device comprises:
[0172] Modeling module 201, establishing MGEMTIP model according to resistivity dispersion, zero frequency resistivity, polarizability and time constant;
[0173] The assignment module 202 determines the IP frequency band based on the experimental data, determines the high frequency band, IP frequency band and low frequency band, and inputs the initial value of true conductivity;
[0174] A first calculation module 203 calculates a polarizability spectrum;
[0175] Iteration module 204, determining whether the low frequency band of the polarizability spectrum satisfies η H =0; if not, adjust the initial value of true conductivity and repeat the calculation of polarizability spectrum until η is satisfied H =0;
[0176] The second calculation module 205 obtains the effective polarizability spectrum and then calculates the true conductivity σ according to the MGEMTIP model 0 M 0 , IP polarizability and IP equivalent time constant.
[0177] As an alternative, the MGEMTIP model is:
[0178]
[0179]
[0180]
[0181]
[0182] Among them, ρ 0 , l Represent the resistivity of the background medium and the first disturbance medium, a l is the spherical radius of the lth disturbance medium, f l is the volume content of the first disturbance medium, α l is the surface polarization coefficient, ω=2πf is the circular frequency, C l is the frequency correlation coefficient.
[0183] As an option, the initial value of true conductivity is:
[0184] σ 0 M 0 =bσ e (min(ω m )).
[0185] As an alternative, the polarizability spectrum is:
[0186]
[0187] Among them, η l is the polarizability spectrum,
[0188] As an option, adjust the initial true conductivity value:
[0189]
[0190] in, Γ H is the high time constant range set corresponding to low frequency, η l For different time constants τ l The corresponding polarizability.
[0191] As an alternative, the IP polarization ratio is:
[0192]
[0193] As an alternative, the IP equivalent time constant is:
[0194]
[0195] Example 5
[0196] The present disclosure provides an electronic device, which includes: a memory storing executable instructions; a processor, which runs the executable instructions in the memory to implement the above-mentioned rock induced polarization parameter inversion method.
[0197] An electronic device according to an embodiment of the present disclosure includes a memory and a processor.
[0198] The memory is used to store non-temporary computer-readable instructions. Specifically, the memory may include one or more computer program products, which may include various forms of computer-readable storage media, such as volatile memory and / or non-volatile memory. The volatile memory may, for example, include random access memory (RAM) and / or cache memory (cache), etc. The non-volatile memory may, for example, include read-only memory (ROM), hard disk, flash memory, etc.
[0199] The processor may be a central processing unit (CPU) or other forms of processing units having data processing capabilities and / or instruction execution capabilities, and may control other components in the electronic device to perform desired functions. In one embodiment of the present disclosure, the processor is used to run the computer-readable instructions stored in the memory.
[0200] Those skilled in the art should be able to understand that in order to solve the technical problem of how to obtain a good user experience, the present embodiment may also include well-known structures such as a communication bus and an interface, and these well-known structures should also be included in the protection scope of the present disclosure.
[0201] For detailed description of this embodiment, reference may be made to the corresponding descriptions in the aforementioned embodiments, which will not be repeated here.
[0202] Example 6
[0203] An embodiment of the present disclosure provides a computer-readable storage medium storing a computer program, which implements the rock induced polarization parameter inversion method when executed by a processor.
[0204] According to the computer-readable storage medium of the embodiment of the present disclosure, non-transitory computer-readable instructions are stored thereon. When the non-transitory computer-readable instructions are executed by a processor, all or part of the steps of the above-mentioned methods of each embodiment of the present disclosure are executed.
[0205] The above-mentioned computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROM and DVD), magneto-optical storage media (e.g., MO), magnetic storage media (e.g., magnetic tape or mobile hard disk), media with built-in rewritable non-volatile memory (e.g., memory card) and media with built-in ROM (e.g., ROM box).
[0206] Those skilled in the art should understand that the purpose of the above description of the embodiments of the present invention is only to exemplarily illustrate the beneficial effects of the embodiments of the present invention, and is not intended to limit the embodiments of the present invention to any given examples.
[0207] The embodiments of the present invention have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those skilled in the art without departing from the scope and spirit of the described embodiments.
Claims
1. A method for inverting rock induced polarization parameters, characterized in that, it includes: establishing an MGEMTIP model based on resistivity dispersion, zero-frequency resistivity, polarizability, and time constant; judging the induced polarization frequency band based on experimental data, determining the high-frequency band, the induced polarization band, and the low-frequency band, and inputting the initial value of true conductivity; calculating the polarizability spectrum; Determine whether the low-frequency band of the polarizability spectrum satisfies η H = 0; if not, adjust the initial value of the true conductivity and repeat the calculation of the polarizability spectrum until η H = 0; Obtain an effective polarizability spectrum, and then calculate the true conductivity σ according to the MGEMTIP model 0 M 0 , IP polarizability and IP equivalent time constant.
2. The method for inverting rock induced polarization parameters according to claim 1, wherein, the MGEMTIP model is: where ρ 0 and ρ l represent the resistivity of the background medium and the resistivity of the l-th disturbed medium, respectively, a l is the sphere radius of the l-th disturbed medium, f l is the volume fraction of the l-th disturbed medium, α l is the surface polarization coefficient, ω = 2πf is the angular frequency, and C l is the frequency-dependent coefficient.
3. The method for inverting rock induced polarization parameters according to claim 1, wherein, the initial value of true conductivity is: σ 0 M 0 = bσ e (min(ω m ))。 4. The method for inverting rock induced polarization parameters according to claim 3, wherein, the polarizability spectrum is: Among them, η l is the polarizability spectrum, 5. The method for inverting rock induced polarization parameters according to claim 1, wherein, adjusting the initial value of true conductivity to: Among them, Γ H is a set of high time constant ranges corresponding to low frequencies, and η l is the polarizability corresponding to different time constants τ l respectively.
6. The method for inverting rock induced polarization parameters according to claim 1, wherein, the IP polarizability is:
7. The method for inverting rock induced polarization parameters according to claim 1, wherein, the IP equivalent time constant is:
8. A device for inverting rock induced polarization parameters, characterized in that, it includes: a modeling module that establishes an MGEMTIP model based on resistivity dispersion, zero-frequency resistivity, polarizability, and time constant; an assignment module that judges the induced polarization frequency band based on experimental data, determines the high-frequency band, the induced polarization band, and the low-frequency band, and inputs the initial value of true conductivity; a first calculation module that calculates the polarizability spectrum; An iterative module determines whether the low-frequency band of the polarizability spectrum satisfies η H = 0; if not, adjust the initial value of the true conductivity and repeat the calculation of the polarizability spectrum until η H = 0; A second calculation module, which obtains an effective polarizability spectrum and then calculates the true conductivity σ according to the MGEMTIP model 0 M 0 , IP polarizability and IP equivalent time constant.
9. An electronic device, characterized in that, the electronic device includes: a memory that stores executable instructions; a processor that runs the executable instructions in the memory to implement the method for inverting rock induced polarization parameters according to any one of claims 1-7.
10. A computer-readable storage medium, characterized in that, the computer-readable storage medium stores a computer program, and when the computer program is executed by a processor, it implements the method for inverting rock induced polarization parameters according to any one of claims 1-7.