CN114428359A - Frequency domain induced polarization forward modeling method, device, electronic equipment and medium - Google Patents

Frequency domain induced polarization forward modeling method, device, electronic equipment and medium Download PDF

Info

Publication number
CN114428359A
CN114428359A CN202011085061.4A CN202011085061A CN114428359A CN 114428359 A CN114428359 A CN 114428359A CN 202011085061 A CN202011085061 A CN 202011085061A CN 114428359 A CN114428359 A CN 114428359A
Authority
CN
China
Prior art keywords
induced polarization
frequency domain
heat storage
deep geothermal
geothermal heat
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202011085061.4A
Other languages
Chinese (zh)
Inventor
徐凤姣
郭全仕
肖鹏飞
李弘�
齐文杰
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN202011085061.4A priority Critical patent/CN114428359A/en
Publication of CN114428359A publication Critical patent/CN114428359A/en
Pending legal-status Critical Current

Links

Images

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
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N27/00—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means
    • G01N27/02—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance
    • G01N27/04—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance
    • G01N27/041—Investigating or analysing materials by the use of electric, electrochemical, or magnetic means by investigating impedance by investigating resistance of a solid body
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/30—Circuit design

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Computer Hardware Design (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Chemical & Material Sciences (AREA)
  • General Engineering & Computer Science (AREA)
  • Analytical Chemistry (AREA)
  • Pathology (AREA)
  • Electrochemistry (AREA)
  • Chemical Kinetics & Catalysis (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Health & Medical Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention provides a frequency domain induced polarization forward modeling method, a device, electronic equipment and a medium. The forward modeling method comprises the following steps: obtaining the high-temperature high-pressure complex resistivity of a rock sample; obtaining an excitation true parameter of a research area through inversion of the excitation true parameter of the complex resistivity, and establishing an excitation model of the research area; carrying out forward calculation of frequency domain induced polarization; and determining the electromagnetic response characteristics of deep geothermal heat storage by applying the forward modeling result of the frequency domain induced polarization and combining geological and geophysical data of a research area. The method is based on the test result of the high-temperature high-pressure complex resistivity of the deep geothermal heat storage rock, establishes a deep geothermal heat storage induced polarization model through damping least square induced polarization true parameter inversion, further develops quadtree finite difference frequency domain electromagnetic forward calculation, combines geological and geophysical data of a research area, determines the response characteristic of deep geothermal heat storage induced polarization, and provides theoretical basis for fine inversion interpretation.

Description

Frequency domain induced polarization forward modeling method, device, electronic equipment and medium
Technical Field
The invention belongs to the field of deep geothermal geophysical exploration, and relates to an electromagnetic exploration method and device for providing theoretical basis for fine inversion explanation in order to clarify deep geothermal heat storage electromagnetic response characteristics.
Background
Frequency domain induced polarization has become one of the most widely used techniques in electromagnetic prospecting, and the method can better describe the morphological characteristics of the underground abnormal body from the multi-parameter perspective. The frequency domain excitation forward modeling method mainly takes a magnetotelluric Method (MT) and a controllable source audio magnetotelluric method (CSAMT) as main components, and takes a natural electromagnetic field, an artificial galvanic couple source or a magnetic couple source electromagnetic field as a field source to carry out excitation parameter response characteristic research, and determine the influence rule of excitation parameters on the apparent resistivity, thereby improving the inversion interpretation precision of the electromagnetic method.
In the prior art, there are many researches on a frequency domain induced polarization method. For example, a CSAMT one-dimensional forward study (geophysical report, 2009.52(7)) containing an excitation effect is based on a Dias model, a complex resistivity is used for replacing a direct-current resistivity without considering the earth polarization effect, a CS AMT field source one-dimensional layered model is subjected to forward simulation, a theoretical basis is provided for extracting excitation information contained in a CSAMT signal, and results show that the excitation abnormality can be visually embodied by an amplitude ratio abnormal peak value and a phase difference abnormal peak value before and after polarization, wherein the apparent abnormality (including a far-field transition field near field) occurs in a phase response curve after the excitation parameters are considered; the abnormal peak value has continuous corresponding relation with the thickness buried depth of the polarizing layer and the resistivity change, and the extraction of the induced polarization information from the frequency domain electromagnetic method signal is considered to have optimistic prospect.
A one-dimensional forward study (oil geophysical exploration, 2009.44(3)) on the MT (stimulated emission) effect of an oil and gas reservoir is used for exploring a method for extracting stimulated emission parameters from high-power artificial fixed-source low-frequency electromagnetic sounding data aiming at the stimulated emission effect of a deep oil and gas reservoir. Numerical simulation experiments show that the Dias model is suitable for the requirements of extracting the induced polarization parameters of the oil and gas reservoir in the frequency range of 0.1-300 Hz frequency and 1-10 km depth range in the frequency range of high-power artificial fixed source low-frequency electromagnetic frequency depth measurement. MT one-dimensional forward modeling is carried out based on a Dias model, and the apparent resistivity ratio and the apparent phase ratio can visually display the abnormity caused by the oil-bearing stratum induced electrical stimulation effect.
An excitation parameter forward numerical simulation (geophysical prospecting technology, 2009.31(1)) based on a finite element method is characterized in that a hexahedron and tetrahedron combined subdivision unit is used in the finite element method to disperse a ground model, a Cole-Cole model complex resistivity response under different frequencies is calculated to replace a volume polarization model in the model, and finally main parameters in a frequency domain excitation method are obtained through simulation. Important parameters of amplitude frequency Fs and phase frequency phi s in a frequency domain induced polarization method are simulated by using a finite element method, and induced polarization response characteristics of the Fs and the gamma s under a three-dimensional geoelectrical condition are contrastively researched, so that the phase frequency parameter can be used as a method for effectively distinguishing polarizers with different time constants.
Finite element numerical simulation of frequency domain induced polarization (geophysical evolution, 2008.23(4)) by substituting Cole-Cole model parameters for the dispersion effect of the earth, a finite element method was used to perform numerical simulation studies on frequency domain induced polarization. Deducing a wave equation of the complex potential of the harmonic transformation magnetic field based on Maxwell equation, obtaining the complex potential wave equation and the component problem thereof by utilizing the generalized variational principle, then dividing and discretizing the region, interpolating in the equivalent variation of the value problem under the cell boundary condition to obtain a group of linear equation sets, obtaining the complex potential value on each node by reasonably storing the rigidity matrix and solving the equation sets, and finally obtaining parameters such as amplitude frequency and the like representing the frequency domain induced polarization response. Simulation results show that frequency domain induced polarization and time domain induced polarization are equivalent in the capability of reflecting induced polarization abnormity.
However, none of the above prior arts has studied a method for establishing a rock induced electrical model for deep geothermal heat storage. With the increase of production demand, the current-stage geothermal exploration is changed from shallow surface exploration to deep geothermal exploration, so that the exploration difficulty is increased, the requirements on the extraction of electromagnetic signals, the effective acquisition of data resolution, fidelity, signal-to-noise ratio and the like are improved, and the deep geothermal exploration by the conventional frequency domain electromagnetic method faces challenges.
Therefore, the art continues with a frequency domain induced polarization forward method for deep geothermal heat storage.
Disclosure of Invention
The invention provides a deep geothermal heat storage oriented frequency domain induced polarization forward modeling method, aiming at the problems that the response of a conventional frequency domain electromagnetic method to deep geothermal heat storage is unclear and the current high-precision inversion interpretation requirements of deep geothermal heat storage cannot be effectively met. The invention provides a frequency domain induced polarization forward modeling method based on a rock physical experiment, which is characterized in that a real induced polarization model of a research area is established on the basis of an induced polarization true parameter inversion result, frequency domain electromagnetic forward response characteristic research is further carried out through a quadtree finite difference method, and a theoretical basis is provided for fine inversion interpretation by combining geological and geophysical data of the research area.
According to one aspect of the invention, a frequency domain induced polarization forward modeling method facing deep geothermal heat storage is provided, and comprises the following steps:
obtaining the high-temperature high-pressure complex resistivity of a rock sample;
obtaining an excitation true parameter of a research area through inversion of the excitation true parameter of the complex resistivity, and establishing an excitation model of the research area;
carrying out forward calculation of frequency domain induced polarization;
and determining the electromagnetic response characteristics of deep geothermal heat storage by applying the forward modeling result of the frequency domain induced polarization and combining geological and geophysical data of a research area.
Further, the rock sample in the deep geothermal research area is subjected to high-temperature and high-pressure complex resistivity tests under the conditions of different solution saturation degrees, and complex resistivity and phase data obtained through the experimental tests are collated to obtain the high-temperature and high-pressure complex resistivity of the rock sample.
Further, by a damping least square inversion method, a target function of complex resistivity induced polarization true parameter inversion and minimum damping two-multiplication inversion is carried out
Figure BDA0002720107120000031
Comprises the following steps:
Figure BDA0002720107120000032
wherein:
Figure BDA0002720107120000033
in order to set the parameters of the model,
Figure BDA0002720107120000034
is the fitted relative deviation between theoretical and measured spectra, jikIs an element of a jacobian matrix J, wherein
Figure BDA0002720107120000035
For model modifiers, the subscript k denotes the kth model parameter and i denotes the ith frequency bin.
Further, by utilizing a finite difference method of a quadtree, forward calculation of the frequency domain induced polarization is carried out, and the resistivity value is replaced by the complex resistivity in the forward calculation process.
Further, a Cole-Cole complex resistivity model is used for carrying out forward calculation of frequency domain induced polarization:
Figure BDA0002720107120000041
wherein: ρ (i ω) is the complex resistivity of the rock in ohms per meter (Ω · m) at different frequencies; omega is angular frequency in hertz (Hz); i is an imaginary unit; rho0The true resistivity of the formation in the absence of polarization, in ohms-meters (Ω · m); m is polarizability; τ is a time constant in units of microseconds(s); c is a frequency correlation coefficient.
According to another aspect of the present invention, there is provided a frequency domain induced polarization forward device for deep geothermal heat storage, comprising:
the acquisition unit is used for acquiring the high-temperature high-pressure complex resistivity of the rock sample;
the modeling unit is used for obtaining the true induced polarization parameters of the research area through the inversion of the true induced polarization parameters of the complex resistivity and establishing an induced polarization model of the research area;
the forward unit is used for carrying out forward calculation on the frequency domain induced polarization;
and the determining unit is used for determining the deep geothermal heat storage electromagnetic response characteristics by applying the frequency domain induced polarization forward modeling result and combining geological and geophysical data of a research area.
Further, by a damping least square inversion method, a target function of complex resistivity induced polarization true parameter inversion and minimum damping two-multiplication inversion is carried out
Figure BDA0002720107120000042
Comprises the following steps:
Figure BDA0002720107120000043
wherein:
Figure BDA0002720107120000044
in order to set the parameters of the model,
Figure BDA0002720107120000045
is the fitted relative deviation between theoretical and measured spectra, jikIs an element of a jacobian matrix J, wherein
Figure BDA0002720107120000046
For model modifiers, the subscript k denotes the kth model parameter and i denotes the ith frequency bin.
Further, a Cole-Cole complex resistivity model is used for carrying out forward calculation of frequency domain induced polarization:
Figure BDA0002720107120000047
wherein: ρ (i ω) is the complex resistivity of the rock in ohms per meter (Ω · m) at different frequencies; omega is angular frequency in hertz (Hz); i is an imaginary unit; rho0The true resistivity of the formation in the absence of polarization, in ohms-meters (Ω · m); m is polarizability; τ is a time constant in units of microseconds(s); c is a frequency correlation coefficient.
According to another aspect of the present invention, there is provided an electronic apparatus including:
the experimental testing device is used for testing the high-temperature high-pressure complex resistivity of the rock sample in the deep geothermal research area to obtain the high-temperature high-pressure complex resistivity of the rock sample;
a memory storing executable instructions and receiving a complex resistivity of the experimental test device;
the processor runs the executable instructions in the memory to realize the frequency domain induced polarization forward modeling method facing the deep geothermal heat storage.
According to another aspect of the present invention, a computer-readable storage medium is provided, which stores a computer program, which when executed by a processor implements the deep geothermal heat storage oriented frequency domain induced polarization forward modeling method.
The frequency domain induced polarization forward modeling method for deep geothermal heat storage has the following characteristics:
and obtaining an excitation true parameter inversion result by utilizing a rock physical test result of the deep geothermal research area and adopting a damping least square inversion algorithm, and establishing a real excitation model of the deep geothermal research area so as to develop frequency domain excitation forward response characteristic research. The method carries out forward calculation based on a real induced polarization model in a research area, and can effectively and definitely determine the deep geothermal heat storage electromagnetic response characteristics.
Drawings
The above and other objects, features and advantages of the present disclosure will become more apparent by describing in greater detail exemplary embodiments thereof with reference to the attached drawings, in which like reference numerals generally represent like parts throughout.
Fig. 1 is a flowchart of a frequency domain induced polarization forward modeling method for deep geothermal heat storage according to the present invention.
Fig. 2 is a flow chart of a frequency domain induced polarization forward modeling method for deep geothermal heat storage according to an embodiment of the invention.
FIG. 3 shows the complex resistivity test results of high-temperature and high-pressure rock according to an embodiment of the invention.
Fig. 4 shows forward results of apparent resistivity in frequency domain without induced electrical parameters according to an embodiment of the present invention.
Fig. 5 shows forward results of apparent resistivity in frequency domain including an excitation parameter according to an embodiment of the present invention.
Detailed Description
Preferred embodiments of the present disclosure will be described in more detail below with reference to the accompanying drawings. While the preferred embodiments of the present disclosure are shown in the drawings, it should be understood that the present disclosure may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the disclosure to those skilled in the art.
The invention provides a deep geothermal heat storage oriented frequency domain induced polarization forward modeling method, aiming at the problems that the response of a conventional frequency domain electromagnetic method to deep geothermal heat storage is unclear and the current high-precision inversion interpretation requirements of deep geothermal heat storage cannot be effectively met. The method is based on the test result of the high-temperature high-pressure complex resistivity of the formation rock of the geothermal research area, establishes the induced polarization model of the research area through induced polarization true parameter inversion, and then develops the forward electromagnetic modeling research of the frequency domain based on the finite difference numerical calculation method of the quadtree, so as to further clarify the response characteristic of deep geothermal heat storage induced polarization and provide a theoretical basis for the fine inversion interpretation.
As shown in fig. 1, the present disclosure provides a frequency domain induced polarization forward modeling method for deep geothermal heat storage, including:
obtaining the high-temperature high-pressure complex resistivity of a rock sample;
obtaining an excitation true parameter of a research area through inversion of the excitation true parameter of the complex resistivity, and establishing an excitation model of the research area;
carrying out forward calculation of frequency domain induced polarization;
and determining the electromagnetic response characteristics of deep geothermal heat storage by applying the forward modeling result of the frequency domain induced polarization and combining geological and geophysical data of a research area.
Specifically, rock samples in a deep geothermal research area are collected firstly, and high-temperature and high-pressure rock complex resistivity test is carried out. For example, the collected rock sample may be subjected to high temperature and high pressure complex resistivity tests under different solution saturations (e.g., clear water saturation and 5% brine saturation), and the complex resistivity and phase data obtained from the experimental tests may be collated.
Preferably, the complex resistivity induced polarization true parameter inversion can be carried out by a damped least square inversion method, induced polarization true parameters in a research area are obtained, and an induced polarization model is established. Minimum damping two-multiplication inversion objective function
Figure BDA0002720107120000071
Comprises the following steps:
Figure BDA0002720107120000072
wherein:
Figure BDA0002720107120000073
in order to set the parameters of the model,
Figure BDA0002720107120000074
is the fitted relative deviation between theoretical and measured spectra, jikIs an element of a jacobian matrix J, wherein
Figure BDA0002720107120000075
For model modifiers, the subscript k denotes the kth model parameter and i denotes the ith frequency bin.
Preferably, starting from a maxwell equation set, a frequency domain induced polarization calculation is performed by using a quadtree finite difference method, a resistivity value is replaced by a complex resistivity in the forward process, and a Cole-Cole complex resistivity model is taken as an example:
Figure BDA0002720107120000076
wherein: ρ (i ω) is the complex resistivity of the rock in ohms per meter (Ω · m) at different frequencies; omega is angular frequency in hertz (Hz); i is an imaginary unit; rho0The true resistivity of the formation in the absence of polarization, in ohms-meters (Ω · m); m is polarizability; τ is a time constant in units of microseconds(s); c is a frequency correlation coefficient.
And finally, a frequency domain induced polarization forward modeling result is applied, and the geological and geophysical data of the research area are combined to clarify deep geothermal heat storage electromagnetic response characteristics.
The method is based on the test result of the high-temperature high-pressure complex resistivity of the deep geothermal heat storage rock, establishes a deep geothermal heat storage induced polarization model through damping least square induced polarization true parameter inversion, further develops quadtree finite difference frequency domain electromagnetic forward calculation, combines geological and geophysical data of a research area, determines the response characteristic of deep geothermal heat storage induced polarization, and provides theoretical basis for fine inversion interpretation.
To facilitate understanding of the aspects of the embodiments of the present invention and the effects thereof, specific application examples are given below. It will be understood by those skilled in the art that this example is merely for the purpose of facilitating an understanding of the present invention and that any specific details thereof are not intended to limit the invention in any way.
Example 1
In this embodiment, a deep geothermal exploration block in the west of china is taken as an example, and the method provided by the present invention is used to perform a frequency domain induced polarization response characteristic study.
Fig. 2 shows a flowchart of a technique for obtaining a frequency domain induced electrical response characteristic study according to an embodiment of the invention.
In step S1, a high-temperature and high-pressure complex resistivity experiment is performed on the rock sample, and amplitude and phase data of the complex resistivity are obtained. FIG. 3 shows the complex resistivity test results of high-temperature and high-pressure rock according to an embodiment of the invention.
In step S3, complex resistivity induced polarization true parameter inversion is performed by a damped least squares inversion method to obtain an induced polarization true parameter in the study region. Minimum damping two-multiplication inversion objective function
Figure BDA0002720107120000081
Comprises the following steps:
Figure BDA0002720107120000082
wherein:
Figure BDA0002720107120000083
in order to set the parameters of the model,
Figure BDA0002720107120000084
is the fitted relative deviation between theoretical and measured spectra, jikIs an element of a jacobian matrix J, wherein
Figure BDA0002720107120000085
For model modifiers, the subscript k denotes the kth model parameter and i denotes the ith frequency bin.
In step S3, an induced polarization model is established according to the inversion result of the complex resistivity induced polarization true parameter.
In step S4, starting from maxwell equations, forward calculation of the frequency domain induced polarization is performed by using a quadtree finite difference method, and the resistivity value is replaced by the complex resistivity in the forward calculation process, taking a Cole-Cole complex resistivity model as an example:
Figure BDA0002720107120000086
wherein: ρ (i ω) is the complex resistivity of the rock in ohms per meter (Ω · m) at different frequencies; omega is angular frequency in hertz (Hz); i is an imaginary unit; rho0True resistivity of the formation in the absence of polarization, in ohmsMeters (Ω · m); m is polarizability; τ is a time constant in units of microseconds(s); c is a frequency correlation coefficient.
In step S5, the deep geothermal reservoir electromagnetic response characteristics are determined by applying the frequency domain induced polarization forward modeling result and combining the geological and geophysical data of the research area. Fig. 4 is a frequency domain electromagnetic forward result without induced polarization parameters, and fig. 5 is a frequency domain electromagnetic forward result with induced polarization parameters, and the forward result shows that the induced polarization parameters have obvious influence on the frequency domain forward response.
Example 2
The embodiment provides a frequency domain excitation forward device for deep geothermal heat storage, which includes:
the acquisition unit is used for acquiring the high-temperature high-pressure complex resistivity of the rock sample;
the modeling unit is used for obtaining the true induced polarization parameters of the research area through the inversion of the true induced polarization parameters of the complex resistivity and establishing an induced polarization model of the research area;
the forward unit is used for carrying out forward calculation on the frequency domain induced polarization;
and the determining unit is used for determining the deep geothermal heat storage electromagnetic response characteristics by applying the frequency domain induced polarization forward modeling result and combining geological and geophysical data of a research area.
The acquisition unit, the modeling unit, the forward modeling unit and the determination unit are sequentially in communication connection. The acquisition unit can acquire the high-temperature and high-pressure complex resistivity of the rock sample from the high-temperature and high-pressure rock complex resistivity experimental testing instrument and send the high-temperature and high-pressure complex resistivity to the modeling unit, and is used for acquiring the true induced polarization parameters of the research area through inversion of the true induced polarization parameters of the complex resistivity and establishing an induced polarization model of the research area. The forward unit carries out frequency domain induced polarization calculation based on the induced polarization model, and sends a frequency domain induced polarization result to the determining unit, and the forward unit is used for determining the deep geothermal heat storage electromagnetic response characteristics by combining geological and geophysical data of a research area.
Example 3
The present embodiment provides an electronic device including: the experimental testing device is used for testing the high-temperature high-pressure complex resistivity of the rock sample in the deep geothermal research area to obtain the high-temperature high-pressure complex resistivity of the rock sample; a memory storing executable instructions; and the processor runs the executable instructions in the memory to realize the frequency domain induced polarization forward modeling method for deep geothermal heat storage.
An electronic device according to an embodiment of the present disclosure includes a memory and a processor.
The memory is to store non-transitory computer readable instructions. In particular, the memory may include one or more computer program products that may include various forms of computer-readable storage media, such as volatile memory and/or non-volatile memory. The volatile memory may include, for example, Random Access Memory (RAM), cache memory (cache), and/or the like. The non-volatile memory may include, for example, Read Only Memory (ROM), hard disk, flash memory, etc.
The processor may be a Central Processing Unit (CPU) or other form of processing unit 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 disclosure, the processor is configured to execute the computer readable instructions stored in the memory.
Those skilled in the art should 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, an interface, and the like, and these well-known structures should also be included in the protection scope of the present disclosure.
For the detailed description of the present embodiment, reference may be made to the corresponding descriptions in the foregoing embodiments, which are not repeated herein.
Example 4
The present embodiment provides a computer-readable storage medium, which stores a computer program, and when the computer program is executed by a processor, the computer program implements the deep geothermal heat storage oriented frequency domain induced polarization forward modeling method.
A computer-readable storage medium according to an embodiment of the present disclosure has non-transitory computer-readable instructions stored thereon. The non-transitory computer readable instructions, when executed by a processor, perform all or a portion of the steps of the methods of the embodiments of the disclosure previously described.
The computer-readable storage media include, but are not limited to: optical storage media (e.g., CD-ROMs and DVDs), magneto-optical storage media (e.g., MOs), magnetic storage media (e.g., magnetic tapes or removable disks), media with built-in rewritable non-volatile memory (e.g., memory cards), and media with built-in ROMs (e.g., ROM cartridges).
It will be appreciated by persons skilled in the art that the above description of embodiments of the invention is intended only to illustrate the benefits of embodiments of the invention and is not intended to limit embodiments of the invention to any examples given.
Having described embodiments of the present invention, the foregoing description is intended to be exemplary, not exhaustive, and not limited to the embodiments disclosed. Many modifications and variations will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments.

Claims (10)

1. A deep geothermal heat storage oriented frequency domain induced polarization forward modeling method is characterized by comprising the following steps:
obtaining the high-temperature high-pressure complex resistivity of a rock sample;
obtaining an excitation true parameter of a research area through inversion of the excitation true parameter of the complex resistivity, and establishing an excitation model of the research area;
carrying out forward calculation of frequency domain induced polarization;
and determining the electromagnetic response characteristics of deep geothermal heat storage by applying the forward modeling result of the frequency domain induced polarization and combining geological and geophysical data of a research area.
2. The deep geothermal heat storage oriented frequency domain excitation forward modeling method according to claim 1, characterized in that the high-temperature high-pressure complex resistivity of the rock sample is obtained by performing high-temperature high-pressure complex resistivity tests on the rock sample in the deep geothermal research area under different solution saturation conditions and by collating the complex resistivity and phase data obtained by the experimental tests.
3. The deep geothermal heat storage oriented frequency domain induced polarization forward modeling method according to claim 1, characterized in that a damping least square inversion method is adopted to carry out complex resistivity induced polarization true parameter inversion and a minimum damping two-multiplication inversion target function
Figure FDA0002720107110000011
Comprises the following steps:
Figure FDA0002720107110000012
wherein:
Figure FDA0002720107110000013
in order to set the parameters of the model,
Figure FDA0002720107110000014
is the fitted relative deviation between theoretical and measured spectra, jikIs an element of a jacobian matrix J, wherein
Figure FDA0002720107110000015
For model modifiers, the subscript k denotes the kth model parameter and i denotes the ith frequency bin.
4. The deep geothermal heat storage oriented frequency domain induced polarization forward modeling method according to claim 1, characterized in that a finite difference quadtree method is used for carrying out frequency domain induced polarization forward modeling calculation, and resistivity values are replaced by complex resistivity values in the forward modeling process.
5. The deep geothermal heat storage oriented frequency domain induced polarization forward modeling method according to claim 4, characterized in that a Cole-Cole complex resistivity model is used for frequency domain induced polarization forward modeling calculation:
Figure FDA0002720107110000021
wherein: ρ (i ω) is the complex resistivity of the rock in ohms per meter (Ω · m) at different frequencies; omega is angular frequency in hertz (Hz); i is an imaginary unit; rho0The true resistivity of the formation in the absence of polarization, in ohms-meters (Ω · m); m is polarizability; τ is a time constant in units of microseconds(s); c is a frequency correlation coefficient.
6. A frequency domain induced polarization forward device facing deep geothermal heat storage is characterized by comprising:
the acquisition unit is used for acquiring the high-temperature high-pressure complex resistivity of the rock sample;
the modeling unit is used for obtaining the true induced polarization parameters of the research area through the inversion of the true induced polarization parameters of the complex resistivity and establishing an induced polarization model of the research area;
the forward unit is used for carrying out forward calculation on the frequency domain induced polarization;
and the determining unit is used for determining the deep geothermal heat storage electromagnetic response characteristics by applying the frequency domain induced polarization forward modeling result and combining geological and geophysical data of a research area.
7. The deep geothermal heat storage oriented frequency domain excitation forward modeling device according to claim 6, characterized in that a damping least square inversion method is adopted to perform complex resistivity excitation true parameter inversion and a minimum damping two-multiplication inversion target function
Figure FDA0002720107110000022
Comprises the following steps:
Figure FDA0002720107110000023
wherein:
Figure FDA0002720107110000024
in order to set the parameters of the model,
Figure FDA0002720107110000025
is the fitted relative deviation between theoretical and measured spectra, jikIs an element of a jacobian matrix J, wherein
Figure FDA0002720107110000031
For model modifiers, the subscript k denotes the kth model parameter and i denotes the ith frequency bin.
8. The deep geothermal heat storage oriented frequency domain induced polarization forward device according to claim 6, wherein a Cole-Cole complex resistivity model is used for frequency domain induced polarization forward calculation:
Figure FDA0002720107110000032
wherein: ρ (i ω) is the complex resistivity of the rock in ohms per meter (Ω · m) at different frequencies; omega is angular frequency in hertz (Hz); i is an imaginary unit; rho0The true resistivity of the formation in the absence of polarization, in ohms-meters (Ω · m); m is polarizability; τ is a time constant in units of microseconds(s); c is a frequency correlation coefficient.
9. An electronic device, characterized in that the electronic device comprises:
the experimental testing device is used for testing the high-temperature high-pressure complex resistivity of the rock sample in the deep geothermal research area to obtain the high-temperature high-pressure complex resistivity of the rock sample;
a memory storing executable instructions and receiving a complex resistivity of the experimental test device;
a processor executing the executable instructions in the memory to implement the frequency domain induced polarization forward method for deep geothermal heat storage of any one of claims 1-5.
10. A computer-readable storage medium, wherein the computer-readable storage medium stores a computer program, and the computer program is executed by a processor to implement the deep geothermal heat storage oriented frequency domain induced polarization forward modeling method of any one of claims 1 to 5.
CN202011085061.4A 2020-10-12 2020-10-12 Frequency domain induced polarization forward modeling method, device, electronic equipment and medium Pending CN114428359A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011085061.4A CN114428359A (en) 2020-10-12 2020-10-12 Frequency domain induced polarization forward modeling method, device, electronic equipment and medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011085061.4A CN114428359A (en) 2020-10-12 2020-10-12 Frequency domain induced polarization forward modeling method, device, electronic equipment and medium

Publications (1)

Publication Number Publication Date
CN114428359A true CN114428359A (en) 2022-05-03

Family

ID=81310345

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011085061.4A Pending CN114428359A (en) 2020-10-12 2020-10-12 Frequency domain induced polarization forward modeling method, device, electronic equipment and medium

Country Status (1)

Country Link
CN (1) CN114428359A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114791633A (en) * 2022-06-21 2022-07-26 西安石油大学 Method, system and medium for monitoring shale gas fracturing
CN117520730A (en) * 2024-01-04 2024-02-06 成都理工大学 Transient electromagnetic data Cole-Cole model parameter jacobian matrix calculation method

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114791633A (en) * 2022-06-21 2022-07-26 西安石油大学 Method, system and medium for monitoring shale gas fracturing
CN114791633B (en) * 2022-06-21 2022-09-30 西安石油大学 Method, system and medium for monitoring shale gas fracturing
CN117520730A (en) * 2024-01-04 2024-02-06 成都理工大学 Transient electromagnetic data Cole-Cole model parameter jacobian matrix calculation method

Similar Documents

Publication Publication Date Title
US11098565B2 (en) Method for estimating permeability of fractured rock formations from induced slow fluid pressure waves
Mueller-Petke et al. QT inversion—Comprehensive use of the complete surface NMR data set
Behroozmand et al. Improvement in MRS parameter estimation by joint and laterally constrained inversion of MRS and TEM data
Jardani et al. Stochastic joint inversion of 2D seismic and seismoelectric signals in linear poroelastic materials: A numerical investigation
WO2009067330A2 (en) Methods and systems for evaluating fluid movement related reservoir properties via correlation of low-frequency part of seismic data with borehole measurements
CN114428359A (en) Frequency domain induced polarization forward modeling method, device, electronic equipment and medium
Liang et al. Estimating petrophysical parameters and average mud-filtrate invasion rates using joint inversion of induction logging and pressure transient data
Commer et al. MPiTOUGH2: A parallel parameter estimation framework for hydrological and hydrogeophysical applications
Kavian et al. A new empirical complex electrical resistivity model
CN113775326B (en) Method and device for evaluating movable water saturation, electronic equipment and medium
EP2929375B1 (en) Method of using a data weighted electromagnetic source for time-lapse monitoring of a reservoir production or hydraulic fracturing
Nadri et al. Estimation of stress-dependent anisotropy from P-wave measurements on a spherical sample
Jänicke et al. Viscoelastic substitute models for seismic attenuation caused by squirt flow and fracture leak off
Buddo et al. Joint interpretation of seismic and TEM data from the Kovykta gas-condensate field, East Siberia
MacLennan Monitoring CO2 storage using well-casing source electromagnetics
Ashlock et al. Critical depths for higher modes by minimally-invasive multimodal surface wave (MMSW) method: Simulations and field test
Proestakis et al. Petrophysical characterization of Lower Cretaceous marly chalk through 1H-NMR and complex conductivity
Calandro et al. Numerical simulation of the response of laterolog tool using the integral equation technique in three dimensional reservoirs
CN114427436A (en) Method and device for evaluating organic carbon content in reservoir well, electronic equipment and medium
Zhdanov et al. Carbonate reservoir rocks show induced polarization effects, based on generalized effective medium theory
Zhao et al. Joint Inversion of Saturation and Q in Low-Permeability Sandstones Using v Spontaneous Potential and Resistivity Logs
Xu et al. An extraction optimization algorithm design for IP parameters in CSAMT
CN116774300A (en) Polarization anomaly information extraction technology based on resistivity inversion result
Rugari Electrokinetic Methods and Applications in Australian Aquifer Settings: High-Dimension Electrical Tomography Imaging and Neural Network Filtration Techniques
Colombo et al. Quantificação de incertezas do acúmulo de sobrepressão utilizando modelagem inversa dos processos de compactação em bacias sedimentares

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