EGS MAGNETIC NANOPARTICLE TRACER AGENT TECHNIQUE AND INTERPRETATION METHOD

Abstract

The disclosure provides an Enhanced Geothermal System (EGS) magnetic nanoparticle tracer agent technique and interpretation method. The method comprises the steps of: through a magnetic nanoparticle surface modification technique and thermal stability analysis of a high-temperature high-pressure reactor, firstly accomplishing the screening of magnetic nanoparticles, so as to prepare magnetic nanoparticles having suitable diffusivity and controllable thermal stability; upon this basis, performing a core penetration test, characterizing EGS connectivity by sampling and analyzing the change in concentration of magnetic nanoparticles, and calculating a heat exchange area between rock and injected water; and meanwhile obtaining electromagnetic signal distribution of magnetic nanoparticles entering a reservoir by utilizing an electrical measurement technology, inverting reservoir connectivity by using resistivity and calculating the heat exchange area, and calibrating the resulting reservoir connectivity and heat exchange area with the connectivity.

Claims

1. An Enhanced Geothermal System (EGS) magnetic nanoparticle tracer agent technique and interpretation method, comprising the following steps: Step 100, accomplishing selection and preparation of a magnetic nanoparticle tracer agent; Step 200, performing an indoor core penetration test by using three tracer agents namely a magnetic nanoparticle tracer agent prepared in Step 100, a conserved tracer agent NaCl and a reactive tracer agent Safraine T, detecting an electromagnetic signal using an exciting electrode, performing inversion calculation on a real component, an imaginary component and polarizability of complex resistivity, and then calculating the porosity of the core; Step 300, characterizing EGS connectivity by sampling and analyzing the change in concentration of magnetic nanoparticles, obtaining penetration curves of different peaks and trailers through a tracer agent test, respectively fitting the penetration curves using a mathematical model, and constructing a fracture solute transport model; Step 400, obtaining electromagnetic signal distribution of magnetic nanoparticles entering into a reservoir by utilizing an electrical measurement technology, and inverting the reservoir connectivity by using resisitivity; and Step 500, comparing resisitivity distribution detectiond outside the core with the penetration curve observed by sampling, comprehensively inverting the reservoir connectivity and calculating the heat exchange area.

2. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 1, wherein the Step 100 specifically comprises the following steps: the magnetic nanoparticles modified by a surface modifying agent are placed in a high-temperature high-pressure reactor, field stable temperature-pressure conditions of an EGS are given, concentration change and experience change of a magnetic nanoparticle tracer agent solution are measured so as to obtain a change relationship depending on temperatures and pressures, thereby screening an optimal surface modifying agent.

3. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 2, wherein the screening specifically comprises the following steps: Step 101, surface modification of magnetic nanoparticles, namely, preparing a certain concentration of a copolymer solution of sulfonated polystyrene and malonic acid, SiO.sub.2 modified magnetic nanoparticles and magnetic ferritin nanoparticles; Step 102, stable pressure sensitivity analysis of magnetic nanoparticles, namely, designing a high-temperature high-pressure reactor test, and analyzing a change relationship of a particle size depending on temperatures and pressures to initially select magnetic nanoparticles meeting performances; and Step 103, selection of high-temperature high-pressure diffusivity, namely, simulating reservoir conditions, displacing a tracer agent through high pressure, and determining influences of different surface modifying agents on adsorptivity and diffusivity of magnetic nanoparticles in pores through a high-pressure displacement tracer agent, thereby preferably selecting high-diffusivity magnetic nanoparticles as an ideal tracer agent.

4. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 1, wherein the Step 200 specifically comprises the following steps: the tracer agents, such as a tracer agent NaCl, a tracer agent Safraine T and a magnetic nanoparticle tracer agent, are monitored in real time respectively using an induced polarization imaging technology, the change in an imaginary part of complex resistivity of the core over time is calculated, a core penetration test result is analyzed by comparing the change with a penetration curve result, penetration time and fracture volume are calculated, the fracture aperture, diffusivity and core porosity parameters are given in a fracture solute transport model, penetration curves of different peaks and trailers are fit, and heat exchange areas are calculated.

5. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 1, wherein, the Step 400 specifically comprises the following steps: real-time distribution detection of magnetic nanoparticles is realized through an electromagnetic imaging technology, which reduces monitoring cost; and a monitoring means is innovated to replace sampling observation with physical geography detection, and the heat exchange area of the reservoir is calculated through mathematic inversion.

6. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 1, wherein the Step 500 specifically comprises the following steps: medium connectivity is inverted by utilizing an electromagnetic signal, the resistivity distribution detected outside the core is compared with the penetration curve obtained by sampling observation to integrate the inverted reservoir connectivity and calculate an effective heat exchange area between injected water and rock.

7. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 1, wherein, in the Step 400, inverting reservoir connectivity by using resisitivity is specifically as follows: the resistivity and polarizability of the core are calculated in real time utilizing a Cole-Cole parameter inversion method of a Marquette algorithm, and the calculated resistivity and polarizability are compared with the initial resistivity and polarizability of the core, so that the change value of resistivity is calculated, and the porosity and permeability parameters of the core are obtained.

8. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 1, wherein, in the Step 200, the penetration time and the fracture volume are calculated according to the test results obtained from an indoor core penetration test, the fracture aperture, diffusivity, core porosity and other parameters are given in the fracture solute transport model, the penetration curves of different peaks and trailers are fit, and the heat exchange areas are calculated.

9. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 1, wherein, in the Step 200, the inversion calculation is specifically as follows: the medium connectivity is inverted by utilizing an electromagnetic signal, the resistivity distribution detected outside the core is compared with the penetration curve obtained by sampling observation to integrate the inverted reservoir connectivity and calculate the heat exchange area between injected water and rock.

10. The EGS magnetic nanoparticle tracer agent technique and interpretation method according to claim 1, wherein in the Step 500, the penetration time and the fracture volume are calculated based on core penetration test results, the fracture aperture, diffusivity, core porosity and other parameters are given in the fracture solute transport model, the penetration curves of different peaks and trailers are fit, and the heat exchange areas are calculated; the medium connectivity is inverted by utilizing the electromagnetic signal, and the resistivity distribution detected outside the core is compared with the penetration curve obtained by sampling observation to integrate the inverted reservoir connectivity and calculate the heat exchange area between injected water and rock.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0045] FIG. 1 is a flowchart of an EGS magnetic nanoparticle tracer agent technique and interpretation method provided according to the embodiment of the disclosure;

[0046] FIG. 2 is a diagram of a core penetration test device provided according to the embodiment of the disclosure;

[0047] FIG. 3 is a structural diagram of a core pressure container provided according to the disclosure;

[0048] FIG. 4 is a diagram of a magnetic nanoparticle modified by a copolymer of sulfonated polystyrene and propandioic acid, SiO.sub.2 modified magnetic nanoparticles and magnetic ferritin;

[0049] In the drawing, reference numbers are as follows:

[0050] 1—tracer agent storage tank; 2—high-pressure spray pump; 3—tracer agent pipe regulating valve; 4—injected water pipe regulating valve; 5—high-pressure constant-current pump; 6—injected water box; 7—core pressure container; 8—collector; 9—temperature sensor; 10—pressure sensor; 11 and 13—air pump pipe stop valves; 12—air pump; 14—back pressure regulating valve; 15—solution recovery pipe stop valve; 16—tracer agent solution recovery container; 17—nitrogen bottle; 18—nitigen pipe; 19—signal transmission wire; 20—tracer agent solution pipe. 701—core body; 702—pressure container; 703—constant temperature heating ring; 704—preservation layer; 705—anti-magnetic ring.

DESCRIPTION OF THE EMBODIMENTS

[0051] Next, the disclosure will be described in detail in combination with specific examples. It should be noted that the following examples are only for further illustrating the disclosure, but the protective scope of the disclosure is not limited to the following examples.

EXAMPLES

[0052] An EGS magnetic nanoparticle tracer agent technique and interpretation method provided by according to this example comprises the following steps: see in FIG. 1:

[0053] Step 100, accomplishing selection and preparation of a magnetic nanoparticle tracer agent using a magnetic nanoparticle surface modifying technology and high-temperature high-pressure thermal stability analysis;

[0054] The Step 100 specifically comprises the following steps: the surface modification of magnetic nanoparticles adopted different magnetic nanoparticle coating materials (copolymers of sulfurized polystyrene and malonic acid, SiO.sub.2 and heat-resistant ferritin, etc.), the particle sizes of the coating materials were adjusted so that the diffusivity of the coating materials can be controlled artificially. The particle size distribution was measured to prepare magnetic nanoparticles with different diffusivities. The prepared magnetic nano tracer agent solutions were placed in the high-temperature high-pressure reactor, the change relationship of concentration of magnetic nanoparticles over temperatures and pressures was obtained, and the most proper surface modifying agent was screened.

[0055] The screening steps are as follows:

[0056] Step 101, surface modification of magnetic nanoparticles, namely, preparing a certain concentration of a copolymer (PSS-co-MA) solution of sulfonated polystyrene and propane diacid, self-assembling the copolymer solution into a cage-like structure, simulating biomimetic minerlization conditions (pH 8.5, 65° C.) of ferritin, synthesizing monodisperse magnetite (Fe.sub.3O.sub.4) nanoparticles in the cage of the copolymer, and removing aggregated magnetic nanoparticles through centrifugation and concentration; obtaining SiO.sub.2 modified magnetic nanoparticles by using the same method; and by controlling the number of iron atoms entering the genetically engineered recombinant heat-resisting ferritin shell, synthesizing the magnetic ferritin nanoparticles with magnetite core through biomimetic mineralization, and controlling the particle size of the magnetic nanoparticle to be between 10 nm and 12 nm (see FIG. 4);

[0057] Step 102, stable-pressure sensitivity analysis, namely, designing a high-temperature high-pressure reactor test, setting temperatures of 90° C. 150° C. and 200° C. and the pressures of 1 Mpa, 10 MPa and 30 Mpa one by one, respectively adding the modified magnetic nanoparticle solutions modified by the different modification technique into a reactor and a displacement devices, wherein the particle size distribution of magnetic nanoparticles was measured respectively using high-resolution transmission electron microscope before and after addition, analyzing a change relationship of particle sizes depending on temperatures and pressures, wherein the surface modifying agent having the most stable particle size distribution and the least coagulation behavior was used as a magnetic nanoparticle for core penetration test; and

[0058] Step 103, selection of high-temperature high-pressure diffusivity, simulating reservoir conditions, displacing the tracer agent through high pressure, determining influences of different surface modifying agents on adsorptivity and diffusivity of magnetic nanoparticles in pores, thereby selecting high-diffusivity magnetic nanoparticles as an ideal tracer agent.

[0059] Step 200, performing an indoor core penetration test by using three tracer agents namely the magnetic nanoparticle tracer agent and the conserved tracer agent NaCl as well as the reactive tracer agent Safraine T, detecting an electromagnetic signal with an exciting electrode, performing inversion calculation on a real component, an imaginary component and polarizability of complex resistivity and then calculating the porosity of the core.

[0060] In the Step 200, the indoor core penetration test was performed by using three tracer agents namely the magnetic nanoparticle tracer agent, the conserved tracer agent NaCl and the reactive tracer agent Safraine T. The core penetration test device is as shown in FIG. 2, and the structure of the core pressure container is as shown in FIG. 3.

[0061] Step 201, before test, pressing the fractures of the core 701 of the selected reservoir sample;

[0062] Step 202, rinsing the core 701 using deionized water to remove any mineral fine particles which may subsequently cause plugging, and then placing the core into a pressure container 702, wherein a constant-temperature heating ring 703 with a thermal insulation material 704 was arranged outside the container, and an anti-magnetic ring 705 was arranged outside the thermal insulation material 704, and three tracer agent solutions were placed in a tracer agent storage box 1;

[0063] Step 203, starting the test, opening an air pump 12 so that the pressure container 7 is kept negative pressure, closing the air pump 12 and stop valves 11 and 13 after the tracer agent solution entered the core 701; setting the pressure parameters of the backpressure regulating valve 14, opening a high-pressure constant-current pump 5 to pressurize the pressure container 7, heating the pressure container 7 using the constant-temperature heating device 703, and monitoring the temperature and pressure values in the pressure container 7 using a temperature sensor 9 and a pressure sensor 10. After temperature and pressure values were stable, setting the pressure parameters of the backpressure regulating valve 14 so that the set pressure parameters were smaller than pressure parameters of the pressure container 7, allowing the tracer agent to flow back to the partial pressure end of the backpressure regulating valve 14, and collecting effluent in a tracer agent solution recovery container 16 at regular intervals.

[0064] Step 204, in the process of performing the core penetration test with three tracer agents, measurement frequency was 1-1000 Hz, 20 frequencies was selected. According to a fluid flowing-in rate, measurement time interval was selected as per min for once.

[0065] Step 205, in the process of monitoring the electromagnetic signal, phases with different frequencies and amplitude data were measured. The real component, imaginary component and polarizability of complex resistivity were calculated by using Marquette inversion method. The calculation process is as follows:

[0066] The complex resistivity of the core can be represented as ρ*=ρ′(ω)+iρ*(ω)

[0067] The complex resistivity spectrum caused by IP effect satisfies a Cole-Cole model:

[00001]$\rho \left(i\omega \right)={\rho }_0\left\{1-m\left[1-\frac{1}{1+\left(i\mathrm{\omega \tau }\right)}\right]\right\}$

[0068] In the formula, ρ-resistivity (excluding IP effect), ρ0-zero frequency resistivity (including IP effect), m-chargeability (polarizability), τ-time constant (unit s), c-frequency-associated coefficient.

[0069] Expressions of imaginary component, real component, phase and amplitude of complex resistivity:

[00002]$\begin{array}{c}\begin{array}{c}\mathrm{Imaginary}\\ \mathrm{component}\end{array}:\mathrm{Im}\rho \left(i\omega \right)=\frac{-\rho \left(0\right){m\left(\mathrm{\omega \tau }\right)}^c\sin \frac{c\pi }{2}}{1+2{\left(\mathrm{\omega \tau }\right)}^c\cos \frac{c\pi }{2}+{\left(\mathrm{\omega \tau }\right)}^{2c}}\\ \begin{array}{c}\mathrm{Real}\\ \mathrm{component}\end{array}:\mathrm{Re}\rho \left(i\omega \right)=\rho \left(0\right)\frac{1+\left(2-m\right){\left(\mathrm{\omega \tau }\right)}^c\cos \frac{c\pi }{2}+\left(1-m\right){\left(\mathrm{\omega \tau }\right)}^{2c}}{1+2{\left(\mathrm{\omega \tau }\right)}^c\cos \frac{c\pi }{2}+{\left(\mathrm{\omega \tau }\right)}^{2c}}\\ \mathrm{Phase}:\varphi \left(\omega \right)=\mathrm{arc}\mathrm{tg}\frac{1+\left(2-m\right){\left(\mathrm{\omega \tau }\right)}^c\cos -{m\left(\mathrm{\omega \tau }\right)}^c\sin \frac{c\pi }{2}}{\frac{c\pi }{2}+\left(1-m\right){\left(\mathrm{\omega \tau }\right)}^{2c}}\\ \mathrm{Amplitude}:\sqrt[\rho \left(0\right)]{\frac{1+\left(2-m\right){\left(\mathrm{\omega \tau }\right)}^c\cos \frac{c\pi }{2}+\left(1-m\right){\left(\mathrm{\omega \tau }\right)}^{2c}}{1+2{\left(\mathrm{\omega \tau }\right)}^c\cos \frac{c\pi }{2}+{\left(\mathrm{\omega \tau }\right)}^{2c}}}\end{array}$

[0070] Step 206, calculating the porosity of the core according to amended Archie equation based on the resisitivity value obtained by inversion.

[0071] The effective resistivity and porosity of fluid-containing statured rock and fluid resistor meet the following relationship formula:

[00003]$\begin{array}{c}{\sigma }_{\mathrm{eff}}={{\sigma }_1\left(1-{\chi }_2\right)}^p+{\sigma }_2{\chi }_2^m\end{array}\mathrm{Where},p=\frac{\log \left(1-{\chi }_2^m\right)}{\log \left(1-{\chi }_2\right)}$

[0072] The porosity and permeability of the rock meet the following formula according to RGPZ model:

[00004]$K_{\mathrm{RGPZ}}=\frac{d^2{\phi }^{3m}}{4{\mathrm{am}}^2}$

[0073] Where, KPGPZ is permeability, unit: m.sup.2, d is a geometrical mean of particle diameter, ϕ is porosity, m is cementation index, usually empirical constant, a is parameter constant, for a tree-dimensional geological body composed of quasi spherical particles, a=8/3.

[0074] Step 300, characterizing EGS connectivity by sampling and analyzing the change in concentration of magnetic nanoparticles, obtaining penetration curves of different peaks and trailers through a tracer agent test, respectively fitting the penetration curves using a mathematical model, and constructing a fracture solute transport model based on a mathematical model for the transfer migration process of a tracer agent in a fracture medium, it was needed to comprehensively consider the flow of the fluid in the fracture medium, heat transfer in the fracture and migration of the tracer agent.

[0075] For underground water flow in the fracture medium, stimulation was performed using simplified Navier-Stokes equation: μ∇.sup.2v=∇P−ρ.sub.wg

[0076] Where, μ is fluid viscosity coefficient, v is fluid velocity, p is pressure, ρ.sub.w is the fluid density. Meanwhile, a mass conservation equation of fluid is combined:

[00005]$\stackrel{\text{?}}{v}.Math.\left({\rho }_{w}b\left(v\right)\right)=\underset{1}{.Math.}{\stackrel{.Math.}{m}}_1\delta \left(r-r_i\right)$$\text{?}\text{indicates text missing or illegible when filed}$

[0077] Where, b represents the aperture of the fracture, and the right item of the equation represents a source and sink term flowing through the fracture. Combined with the above two sets of equations and Poiseuille's fluid law, we can obtain a main governing equation of fluid flow in the fracture medium:

[00006]$-\nabla .Math.\left(\frac{{\rho }_w{w}^3}{12\mu }\nabla P\right)=\underset{1}{.Math.}{\stackrel{.Math.}{m}}_1\delta \left(r-r_i\right)$

[0078] This equation is solved to obtain the distribution of a pressure field of a fracture flow field, and the distribution of a velocity field can be obtained by substituting the simplified Navier-Stokes equation to solve the migration of the tracer agent below. To describe the influence of temperature change, it is necessary to consider an energy conservation equation based on the flow field so as to describe the heat transfer process in the fracture medium:

[00007]${\rho }_w{c}_{w}b\frac{\partial T}{\partial t}+{\rho }_w{c}_{w}b.Math.v.Math..Math.\nabla T-k_{w}b{\nabla }^{2}T=q_s$

[0079] Where, c.sub.w and k.sub.w represent the heat capacity and thermal conductivity of fracture fluid respectively, and represent heat exchanged between a fracture surface and a rock matrix. The heat transfer process in the rock matrix can be solved by the following energy conservation equation:

[00008]$\frac{\partial \left(H_t\right)}{\partial t}=-\nabla \left({\rho }_L{h}_L{v}_L\right)+\nabla \left(\lambda \nabla T\right)+Q_{\mathrm{heat}}$

[0080] Where, h.sub.L is enthalpy of liquid phase, λ is the heat conductivity coefficient, Q.sub.heat is the source and sink of heat, H.sub.t is the total enthalpy in the system, including the contributions of fluid and rock:

H.sub.t=ϕρ.sub.Lh.sub.L+(1−ϕ)ρ.sub.Rc.sub.pRT

[0081] Where, ρ.sub.R is the density of rock, c.sub.pR is the specific heat capacity of rock. On the basis of the flow field and heat transfer field in the fracture medium obtained above, the flow of the tracer agent in the fracture medium is further solved. Based on the mass conservation equation of the tracer agent in the fracture medium, we can establish the main flow control equation of the tracer agent:

[00009]${\frac{\partial C}{\partial t}+.Math.v.Math.\nabla C-\nabla .Math.\left(D_t.Math.\nabla c\right)-\frac{{\varphi }_m}{b}{D}_m\nabla C_m.Math.}_{Z=b}=0$

[0082] Where, C is the concentration of the tracer agent in the fracture, D is the diffusion coefficient of the tracer agent in the fracture, and v is the flow rate of solute in the fracture, which can be obtained by solving the flow field of the fracture medium described above. Cm is the concentration of solute in the rock matrix, Dm is the diffusion coefficient of solute in the rock matrix. The following control equation can be established and solved by the diffusion of the tracer agent in the rock matrix:

[00010]$\frac{\partial C_m}{\partial t}-D_m\frac{{\partial }^2{C}_m}{\partial y^2}=0$

[0083] The whole solution of the mathematical model is discretized in space using a Galerkin finite element method and in time using an Euler difference method. Based on a global implicit coupling algorithm, the nonlinearity of the coupling equation is handled using Newton Raphson method, and SparseLU is directly used to solve a sparse matrix.

[0084] Step 400, obtaining the electromagnetic signal distribution of magnetic nanoparticles entering the reservoir by using an electrical measurement technology, and inverting the reservoir connectivity of the reservoir by resistivity;

[0085] The electromagnetic signal distribution of magnetic nanoparticles entering the reservoir was obtained by electrical logging technology. The fracture connectivity of the reservoir was inversed by resistivity, and the penetration curves of different peaks and trailers were obtained. The penetration curves were simulated and obtained by using the above equation, and the average residence time, recovery rate and fluid flow rate of the tracer agent were calculated, the simulation results of the fracture solute transport model were calibrated, and the flow path and permeability of fluid were evaluated.

[0086] Step 500, calibrating the simulated penetration curve with the really measured penetration curve, and fitting to obtain a concentrated parameter ø/b√{square root over (D.sub.m)}, which was represented by ΔC:

ΔC=ø/b√{square root over (D.sub.m)}

[0087] In the formula, φ is porosity, b is the half opening of the fracture, and Dm is dispersion coefficient of the tracer agent. φ is only associated with the core, Dm is only associated with the tracer agent, and b has the following relationship with the heat exchange area (SA/V) between water and rock in the fracture in unit volume (V):

SA/V=1/b

[0088] The fracture volume V can be obtained from the average residence time (τ) of the tracer agent curve:

V=Q*τ

[0089] In the formula, Q is the flow of the fluid in the core, and Q and τ can be obtained by a single tracer agent curve.

[0090] The resistivity distribution detected outside the core was compared with the penetration curve obtained by sampling observation, and the reservoir connectivity was comprehensively inversed and the heat exchange area was calculated.

[0091] As a preferred technical solution of the disclosure, in fracture connectivity evaluation and heat exchange area calculation of an EGS heat reservoir based on an electromagnetic signal, electromagnetic signal measurement was performed by using the induced polarization imaging technology to respectively monitor three tracer agents in real time, the change in the imaginary part of core complex resistivity over time was calculated, and the change was compared with the results of the penetration curve for analysis.

[0092] For evaluation of fracture connectivity, the resistivity and polarizability of the core were calculated in real time by using the Cole-Cole parameter inversion method of Marquette algorithm, the resulting resistivity and polarizability were compared with the initial resistivity and polarizability of the core, the change value of resistivity was calculated, and the parameters such as porosity and permeability of the core were obtained.

[0093] For the calculation of heat exchange area, parameters such as penetration time, porosity and permeability were calculated based on the real-time monitoring results of core geophysics real-time monitoring results, the fracture aperture and diffusion coefficient were given in the fracture solute transport model, and the change in the imaginary part of complex resistivity over time was fitted, and the heat exchange area was calculated.

[0094] Compared with the prior art, the disclosure has the following advantages:

[0095] (1) The disclosure adopts a method of combining an indoor test and a mathematical model and meanwhile is based on the established EGS project site data, magnetic nanoparticles are used as the tracer agents, the wild geological conditions are simulated through an indoor high-temperature high-pressure reactor test, different temperatures and pressures are given, surface modifying agents with significant differences in diffusivity and adsorptivity are selected to perform surface modification on magnetic nanoparticles, and the change in particle size and diffusivity is measured, thereby selecting the surface modifying agent to solve the problems of thermal stability and adsorptivity of magnetic nanoparticles.

[0096] (2) In the disclosure, the conserved tracer agent NaCl is combined with the adsorption tracer agent Safraine T to establish an indoor test platform for tracer agent test; the medium connectivity is inverted by using the electromagnetic signal with relatively low monitoring cost, the resistivity distribution detected outside the core is compared with the penetration curve obtained by sampling observation to integrate the inverted reservoir connectivity and calculate the effective heat exchange area between injected water and rock.

[0097] (3) In the disclosure, performance exhibition of the magnetic nanoparticle as the tracer agent at high temperature and high pressure is analyzed through the indoor core penetration test, a quantitative analysis method for key production parameters such as fracture aperture, connectivity and heat exchange area is established to obtain a new understanding of the magnetic nanoparticles tracer agent technology, and a new method for explaining the fracture connectivity in the reservoir by using the electromagnetic imaging technology provides a new technical means for future geothermal energy development.

[0098] The specific embodiments of the disclosure are described above. It should be understood that the disclosure is not limited to the above specific embodiments, and those skilled in the art can make various deformations or modifications within the scope of the claims, which does not affect the essence of the disclosure.