Data processing method and device for CO2 injection acid-rock reaction

Abstract

This specification provides a data processing method and apparatus for CO2-acid-rock reaction. Based on this method, geological parameters of a target area can be obtained first; and according to the geological parameters of the target area, a target acid-rock reaction multi-physical field model and a target chemical field model that couple a seepage field, a stress field, a temperature field, and a solute transport field for the target area can be constructed; then the initial state data and initial boundary conditions of the target acid-rock reaction multi-physical field model are determined; according to the initial state data and initial boundary conditions, by alternately using the target acid-rock reaction multi-physical field model and the target chemical field model for sequential non-iterative processing, a target acid-rock reaction simulation result for the target area is obtained. Thus, it is possible to accurately and comprehensively implement acid-rock reaction simulation under complex conditions for a long time period, and obtain a target acid-rock reaction simulation result with high reference value to guide subsequent CO2 injection for oil displacement in the target area and determination of CO2 storage.

Description

Technical Field

[0001] The present invention relates to the technical field of oil and gas development, and in particular to a data processing method and device for CO2 injection acid-rock reaction. Background Art

[0002] In the process of oil and gas development, CO2 oil flooding is often required. When CO2 is injected into the formation, it will dissolve in the formation water to generate carbonic acid and react chemically with rock minerals. The long-term CO2-water-mineral reaction will cause changes in reservoir properties, pore structure, and mineral surface properties, which in turn will affect the CO2 oil recovery efficiency and the stability of geological storage.

[0003] Based on the existing methods, most of them simply simulate the above CO2-water-mineral reaction in a short period of time through experiments. It is often impossible to accurately and comprehensively simulate the CO2 injection acid rock reaction under complex conditions for a long period of time, which will cause errors in subsequent oil and gas development, reservoir prediction, etc.

[0004] To address the above problems, no effective solution has been proposed yet. Summary of the invention

[0005] This specification provides a data processing method and device for CO2 injection acid-rock reaction, which can accurately and comprehensively realize acid-rock reaction simulation over a long period of time under complex conditions and obtain target acid-rock reaction simulation results with high reference value.

[0006] This specification provides a data processing method for CO2 injection acid rock reaction, including:

[0007] Obtain geological parameters of the target area;

[0008] According to the geological parameters of the target area, a target acid-rock reaction multi-physics field model and a target chemical field model for the target area are constructed; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, stress field, temperature field, and solute migration field;

[0009] Determine the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics model;

[0010] According to the initial state data and initial boundary conditions, the target acid-rock reaction multi-physics field model and the target chemical field model are used to perform sequential non-iterative processing to obtain target acid-rock reaction simulation results for the target area.

[0011] In one embodiment, after obtaining the target acid-rock reaction simulation results for the target area, the method further includes:

[0012] According to the simulation results of the target acid-rock reaction, predict the target change trend of reservoir physical properties, pore structure and mineral surface properties in the target area;

[0013] According to the target change trend, the oil recovery effect of CO2 injection in the target area is determined.

[0014] In one embodiment, the geological parameters of the target area include at least one of the following: reservoir parameters, formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, original formation pressure, elastic modulus of formation rock, elastic modulus of formation rock skeleton, Poisson's ratio of formation rock, density of formation rock skeleton, thermal expansion coefficient of formation rock, original formation temperature, constant-pressure specific heat capacity of formation rock, thermal conductivity of formation rock, constant-pressure specific heat capacity of formation water, thermal conductivity of formation water, diffusion coefficient of solute in formation aqueous solution, elemental composition of formation aqueous solution, pH value and PE value of formation water, and mineral composition of formation rock.

[0015] In one embodiment, according to the geological parameters of the target area, a target acid-rock reaction multi-physics field model and a target chemical field model for the target area are constructed, including:

[0016] Using the first software to construct a target acid-rock reaction multi-physics field model for the target area according to geological parameters of the target area;

[0017] The second software is used to construct a target chemical field model for the target area according to the geological parameters of the target area.

[0018] In one embodiment, a target acid-rock reaction multi-physics field model for the target area is constructed according to geological parameters of the target area, including:

[0019] Construct a geological model of the target area based on the reservoir parameters of the target area;

[0020] A seepage field model is established in the geological model of the target area according to the formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, and original formation pressure of the target area; a stress field model is established in the geological model of the target area according to the elastic modulus of the formation rock, the elastic modulus of the formation rock skeleton, the Poisson's ratio of the formation rock, the density of the formation rock skeleton, and the thermal expansion coefficient of the formation rock; a temperature field model is established in the geological model of the target area according to the original formation temperature of the target area, the constant-pressure specific heat capacity of the formation rock, the thermal conductivity of the formation rock, the constant-pressure specific heat capacity of the formation water, and the thermal conductivity of the formation water;

[0021] According to the mutual influence relationship of key parameters in the model, the seepage field model, stress field model, temperature field model, and solute migration field model are coupled to obtain the target acid-rock reaction multi-physics field model.

[0022] In one embodiment, after the target acid-rock reaction multi-physics field model is obtained by coupling the seepage field model, the stress field model, the temperature field model, and the solute migration field model according to the mutual influence relationship of the key parameters in the model, the method further includes:

[0023] According to the geological parameters of the target area, characteristic information of reactants at each grid node in the target acid-rock reaction multi-physics field model is set.

[0024] In one embodiment, according to the initial state data and the initial boundary conditions, sequential non-iterative processing is performed by using the target acid-rock reaction multi-physics field model and the target chemical field model, including:

[0025] The sequential non-iterative processing of the current round is performed as follows:

[0026] Using the first software, based on the initial state data of the current round and the initial boundary conditions of the current round, a multi-physics field coupling simulation operation of the target acid-rock reaction multi-physics field model with a specified time step of the current round is performed to obtain the physical field simulation operation result of the current round of the target acid-rock reaction multi-physics field model;

[0027] Using a connection program generated based on the third software, the first software is used to obtain the current round of physical field simulation calculation results of the target acid-rock reaction multi-physical field model; and the current round of physical field simulation calculation results of the target acid-rock reaction multi-physical field model are imported into the target chemical field model in the second software;

[0028] Using the second software, based on the current round of physical field simulation operation results of the target acid-rock reaction multi-physical field model, the target chemical field model is subjected to the current round of chemical reaction kinetics operation of the specified time step to obtain the current round of chemical field simulation operation results of the target chemical field model.

[0029] In one embodiment, after obtaining the chemical field simulation calculation result of the current round of the target chemical field model, the method further includes:

[0030] Check whether the cut-off condition is met;

[0031] When it is determined that the cutoff condition is not met, the initial state data and initial boundary conditions of the next round are updated according to the chemical field simulation operation results of the current round of the target chemical field model to perform the next round of sequential non-iterative processing.

[0032] This specification also provides a data processing device for CO2 injection acid rock reaction, including:

[0033] An acquisition module, used to acquire geological parameters of a target area;

[0034] A construction module is used to construct a target acid-rock reaction multi-physics field model and a target chemical field model for the target area according to the geological parameters of the target area; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, the stress field, the temperature field, and the solute migration field;

[0035] A determination module, used to determine the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model;

[0036] The processing module is used to obtain the target acid-rock reaction simulation result about the target area by performing sequential non-iterative processing based on the initial state data and initial boundary conditions by using the target acid-rock reaction multi-physics field model and the target chemical field model.

[0037] The present specification also provides a computer-readable storage medium having computer instructions stored thereon, which, when executed by a processor, implement the following steps: obtaining geological parameters of a target area; constructing a target acid-rock reaction multi-physics field model and a target chemical field model for the target area based on the geological parameters of the target area; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples a seepage field, a stress field, a temperature field, and a solute migration field; determining initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model; and obtaining target acid-rock reaction simulation results for the target area by performing sequential non-iterative processing using the target acid-rock reaction multi-physics field model and the target chemical field model based on the initial state data and the initial boundary conditions.

[0038] Based on the data processing method and device for the CO2 injection acid-rock reaction provided in this specification, the geological parameters of the target area can be obtained first; and according to the geological parameters of the target area, a target acid-rock reaction multi-physics field model and a target chemical field model for the target area can be constructed; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, stress field, temperature field, and solute migration field; and then the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model are determined; according to the initial state data and initial boundary conditions, the target acid-rock reaction multi-physics field model and the target chemical field model are alternately used for sequential non-iterative processing to obtain the target acid-rock reaction simulation results for the target area. Thus, the acid-rock reaction simulation for a long period of time under complex conditions can be accurately and comprehensively realized, and the target acid-rock reaction simulation results with high reference value and good effect can be obtained, and then the above-mentioned target acid-rock reaction simulation results can be used to better guide the subsequent determination of CO2 injection for oil recovery and CO2 storage in the target area. Furthermore, COMSOL is used as the first software to construct and run the target acid-rock reaction multi-physics field model, and PHREEQC is used as the second software to construct and run the target chemical field model. At the same time, MATLAB is used as the third software to connect the first software and the second software, so as to make full use of the characteristics and advantages of the above-mentioned different software and realize the above-mentioned sequential non-iterative processing more efficiently and accurately. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] In order to more clearly illustrate the embodiments of this specification, the drawings required for use in the embodiments will be briefly introduced below. The drawings described below are only some embodiments recorded in this specification. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying creative work.

[0040] Figure 1 It is a flow chart of a data processing method for CO2 injection acid rock reaction provided by an embodiment of this specification;

[0041] Figure 2 It is a schematic diagram of an embodiment of a data processing method for CO2 acid rock reaction provided by an embodiment of this specification in a scenario example;

[0042] Figure 3 It is a schematic diagram of the structure of a server provided by an embodiment of this specification;

[0043] Figure 4 This is a schematic diagram of the structure of a data processing device for CO2 injection acid rock reaction provided by an embodiment of this specification;

[0044] Figure 5It is a schematic diagram of a scenario example, in which the data processing method for CO2 acid rock injection provided in the embodiment of this specification is applied to perform a multi-field and multi-scale coupled simulation of heat-fluid-solid-chemical CO2 acid rock injection reaction;

[0045] Figure 6 It is a schematic diagram of the verification result of the molar amount of calcite in a one-dimensional reaction migration model obtained by applying the data processing method for the CO2 injection acid rock reaction provided in the embodiment of this specification in a scenario example;

[0046] Figure 7 It is a schematic diagram of the verification result of the molar amount of dolomite in a one-dimensional reaction migration model obtained by applying the data processing method for the CO2 injection acid rock reaction provided in the embodiment of this specification in a scenario example;

[0047] Figure 8 It is a schematic diagram of pH verification results in an aqueous solution of a one-dimensional reaction migration model obtained by applying the data processing method for CO2 injection acid rock reaction provided in an embodiment of this specification in a scenario example;

[0048] Fig. 9 This is a schematic diagram of the verification result of the molar amount of Ca element in a one-dimensional reaction migration model obtained by applying the data processing method for the CO2 injection acid rock reaction provided in the embodiment of this specification in a scenario example;

[0049] Fig.10 This is a schematic diagram of the verification result of the molar amount of Mg element in a one-dimensional reaction migration model obtained by applying the data processing method for CO2 injection acid rock reaction provided in an embodiment of this specification in a scenario example;

[0050] Fig.11 This is a schematic diagram of the verification results of the molar amount of Cl element in a one-dimensional reaction migration model obtained by applying the data processing method for CO2 injection acid rock reaction provided in an embodiment of this specification in a scenario example. DETAILED DESCRIPTION

[0051] In order to enable those skilled in the art to better understand the technical solutions in this specification, the technical solutions in the embodiments of this specification will be clearly and completely described below in conjunction with the drawings in the embodiments of this specification. Obviously, the described embodiments are only part of the embodiments of this specification, not all of the embodiments. Based on the embodiments in this specification, all other embodiments obtained by ordinary technicians in this field without creative work should fall within the scope of protection of this specification.

[0052] See also Figure 1 As shown, the embodiment of this specification provides a data processing method for CO2 injection acid rock reaction. Among them, when the method is implemented specifically, it can include the following contents:

[0053] S101: Acquire geological parameters of the target area;

[0054] S102: constructing a target acid-rock reaction multi-physics field model and a target chemical field model for the target area according to geological parameters of the target area; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples a seepage field, a stress field, a temperature field, and a solute migration field;

[0055] S103: Determine initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model;

[0056] S104: According to the initial state data and initial boundary conditions, sequential non-iterative processing is performed by using the target acid-rock reaction multi-physics field model and the target chemical field model to obtain a target acid-rock reaction simulation result for the target area.

[0057] Based on the above embodiments, a target acid-rock reaction multi-physical field model and a target chemical field model for the target area that are simultaneously coupled with multiple physical fields including seepage field, stress field, temperature field, and solute migration scenarios can be constructed first; then, the above-mentioned standard acid-rock reaction model and target chemical field model can be used to perform sequential non-iterative processing to conduct a more comprehensive acid-rock reaction simulation for a long period of time under complex conditions in a simulation rather than an experimental manner, so as to accurately and quickly obtain target acid-rock reaction simulation results with high reference value, thereby better guiding subsequent oil and gas development work such as CO2 oil drive in the target area.

[0058] In some embodiments, the above-mentioned target area can be specifically understood as a reservoir area to be predicted.

[0059] In some embodiments, the geological parameters of the above-mentioned target area may specifically include at least one of the following: reservoir parameters, formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, original formation pressure, elastic modulus of formation rock, elastic modulus of formation rock skeleton, Poisson's ratio of formation rock, density of formation rock skeleton, thermal expansion coefficient of formation rock, original formation temperature, constant-pressure specific heat capacity of formation rock, thermal conductivity of formation rock, constant-pressure specific heat capacity of formation water, thermal conductivity of formation water, diffusion coefficient of solute in formation aqueous solution, elemental composition of formation aqueous solution, pH value and PE value of formation water, mineral composition of formation rock, etc.

[0060] The pH value may specifically refer to the negative logarithm of the hydrogen ion concentration (activity) in the solution, which is used to characterize the concentration of hydrogen ions in the solution, that is, the acidity or alkalinity of the solution. The pe value may specifically refer to the negative logarithm of the electron activity in the solution, which is used to characterize the electron activity in the solution, that is, the oxidation / reduction ability of the system.

[0061] Of course, it should be noted that the geological parameters listed above are only a schematic illustration, and in specific implementation, according to specific circumstances and processing requirements, the geological parameters may also include other types of parameter data, which is not limited in this specification.

[0062] In some embodiments, the target acid-rock reaction multi-physics field model and the target chemical field model for the target area are constructed according to the geological parameters of the target area. When implemented specifically, the following contents may be included:

[0063] S1: using the first software to construct a target acid-rock reaction multi-physics field model for the target area according to the geological parameters of the target area;

[0064] S2: Using the second software to construct a target chemical field model for the target area according to the geological parameters of the target area.

[0065] The first software may specifically include: COMSOL. The second software may specifically include: PHREEQC.

[0066] The above-mentioned COMSOL specifically refers to a scientific research and engineering auxiliary software that mainly uses finite element numerical analysis and is supplemented by mixed numerical calculation and equation mechanism research. The above-mentioned PHREEQC specifically refers to a computer software used to calculate various low-temperature hydrogeochemical reactions.

[0067] In some embodiments, the target acid-rock reaction multi-physics field model for the target area is constructed according to the geological parameters of the target area. When implemented specifically, the following contents may be included:

[0068] S1: Construct a geological model of the target area based on the reservoir parameters of the target area;

[0069] S2: Establish a seepage field model in the geological model of the target area according to the formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, and original formation pressure of the target area; establish a stress field model in the geological model of the target area according to the elastic modulus of the formation rock, the elastic modulus of the formation rock skeleton, the Poisson's ratio of the formation rock, the density of the formation rock skeleton, and the thermal expansion coefficient of the formation rock; establish a temperature field model in the geological model of the target area according to the original formation temperature of the target area, the constant-pressure specific heat capacity of the formation rock, the thermal conductivity of the formation rock, the constant-pressure specific heat capacity of the formation water, and the thermal conductivity of the formation water;

[0070] S3: According to the mutual influence relationship of key parameters in the model, the seepage field model, stress field model, temperature field model, and solute migration field model are coupled to obtain the target acid-rock reaction multi-physics field model.

[0071] In specific implementation, a geological model of the target area may be created in the first software based on reservoir parameters such as the length, width and height of the reservoir in the target area.

[0072] Furthermore, a seepage field model can be established in the geological model of the above-mentioned target area according to geological parameters related to seepage, such as formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, and original formation pressure in the target area; a stress field model can be established in the geological model of the above-mentioned target area according to geological parameters related to stress, such as elastic modulus of formation rock in the target area, elastic modulus of formation rock skeleton, Poisson's ratio of formation rock, density of formation rock skeleton, thermal expansion coefficient of formation rock; a temperature field model can be established in the geological model of the target area according to geological parameters related to temperature, such as original formation temperature in the target area, constant pressure specific heat capacity of formation rock, thermal conductivity of formation rock, constant pressure specific heat capacity of formation water, thermal conductivity of formation water. In addition, a solute migration field model can be further established in the geological model of the target area according to geological parameters related to solute migration.

[0073] Specifically, the seepage process of the fluid in the porous medium can be described based on Darcy's law. By using the first software, the seepage field model is established in the geological model of the target area by constructing the following seepage field control equation:

[0074]

[0075] Where φ is the porosity, dimensionless; ρ f is the fluid density, kg / m 3 ;Q m is the source and sink term of the fluid, kg / (m 3 ·s); t is time; u is the velocity vector.

[0076] Furthermore, the velocity vector u can be calculated according to the following formula:

[0077]

[0078] Where p is the fluid pressure, Pa; k is the permeability, m 2 ; μ is the dynamic viscosity of the fluid, Pa·s.

[0079] Specifically, the first software can be used to construct the following temperature control equation to establish a temperature field model in the geological model of the target area:

[0080]

[0081] Where T is temperature, K; Cp,f is the constant pressure heat capacity of the fluid, J / (kg·K); u is the velocity vector, which can be obtained from the seepage field model; Q is the heat source term, W / m 3 ;k eff is the effective thermal conductivity; (ρC p ) eff is the effective volume heat capacity.

[0082] Furthermore, the above effective volume heat capacity can be calculated according to the following formula:

[0083] (ρC p ) eff =θ s ρ s C p,s +θ f ρ f C p,f

[0084] Among them, θ s is the volume fraction of porous media, dimensionless; θ f is the fluid volume fraction, dimensionless; ρ s is the density of porous media, kg / m 3 ; C p,s is the constant pressure heat capacity of the porous medium, J / (kg·K); k eff is the effective thermal conductivity.

[0085] The above effective thermal conductivity can be calculated by calculating the volume average according to the following formula:

[0086] k eff =θ s k s +θ f k f

[0087] Among them, k s is the thermal conductivity of porous media, W / (m·K); k f is the thermal conductivity of the fluid, W / (m·K).

[0088] Specifically, the first software can be used to construct the following solid mechanics control equation to establish a stress field model in the geological model of the target area:

[0089]

[0090] Where s is the displacement vector; σ is the stress, Pa; f v is the body force, which includes gravity, body load, and body force caused by thermal stress, fluid pressure, and adsorption stress.

[0091] Among them, the relevant porous media porosity and permeability evolution formula can refer to the following:

[0092]

[0093]

[0094] Where φ0 is the initial porosity, dimensionless; k0 is the initial permeability, m 2 ; α is the Biot coefficient, dimensionless; ε v is the volume strain of porous media, dimensionless; ε s is the adsorption strain of porous media, dimensionless; α s is the thermal expansion coefficient of porous media, K -1 ; K s is the bulk modulus of the porous medium, Pa.

[0095] Specifically, the first software can be used to construct the following solute migration field control equation to establish a solute migration field model in the geological model of the target area:

[0096]

[0097] Where u is the velocity vector, obtained from the seepage field; C i is the concentration of component i in the solution, mol / m 3 ;D i is the diffusion coefficient of component i in the solution, m 2 / s; Rate i is the reaction rate of component i in the solution, mol / (m 3 ·s).

[0098] Based on the control equations of the above-mentioned multiple physical field models, the following mutual influence relationship between the key parameters of different physical field models can be obtained: the density and viscosity of the fluid in the seepage field model will change with temperature, and the porosity and permeability will also be affected by temperature and volume strain; the convective heat transfer part in the temperature field model (or thermal field) will be affected by the fluid flow rate; and factors such as effective heat capacity are related to the stress field model (or solid mechanics field), and the stress field model contains thermal stress and volume force generated by fluid pressure. Based on the mutual influence relationship of the above-mentioned key parameters, the seepage field model, stress field model, and temperature field model can be effectively coupled to realize the coupling of the three physical fields of heat, fluid and solid, and obtain the corresponding THM coupling model. On the basis of the above-mentioned THM coupling model, the solute transport field model can be coupled again, so as to obtain a target acid-rock reaction multi-physical field model that can better couple a variety of different physical fields, can be applied to complex reaction conditions, is more comprehensive, and has higher accuracy.

[0099] In some embodiments, after the target acid-rock reaction multi-physics field model is obtained by coupling the seepage field model, stress field model, temperature field model, and solute migration field model according to the mutual influence relationship of the key parameters in the model, the method may also include the following contents when it is implemented: according to the geological parameters of the target area, the characteristic information of the reactants at each grid node in the target acid-rock reaction multi-physics field model is set.

[0100] Specifically, the characteristic information of the reactants at each grid node in the target acid-rock reaction multi-physics field model can be set according to the mineral composition of the formation rock and its substance amount, the pH value and PE value of the formation water in the target area, and the elemental composition and concentration of the formation water solution.

[0101] For example, the mineral composition of the formation rock and the amount of its substance can be: 0.1mol, calcite: 0.05mol, albite: 0.05mol, illite 0.1mol. The pH value of the formation water is 7; the pe value is 4. The element composition and concentration of the formation water solution are: Na 0.1mol / L; Cl 0.3mol / L; Mg 0.1mol / L.

[0102] In addition, during specific implementation, the amount of CO2 to be injected into the target area can be converted into the corresponding hydrogen ion concentration and carbonate ion concentration in the formation aqueous solution based on the geological parameters of the target area, so as to obtain the elemental composition and concentration of the formation aqueous solution after the CO2 injection; and then the elemental composition and concentration of the formation aqueous solution after the CO2 injection can be used to simulate the acid-rock reaction of CO2 injection in the target area, that is, the CO2-water-mineral reaction over a long period of time.

[0103] In specific implementation, a related program written based on the third software can be used to obtain and set the node information of each grid node in the target acid-rock reaction multi-physics field model in the first software, and the corresponding characteristic information of the reactants such as the molar amount of minerals involved in the acid-rock reaction, the ion composition and concentration in the aqueous solution, etc. can be set at each grid node.

[0104] The third software may be specifically MATLAB. Further, MATLAB is specifically a data processing software used in fields such as data analysis, wireless communication, deep learning, image processing and computer vision, signal processing, quantitative finance and risk management, robotics, and control systems.

[0105] Here, it should be noted that since MATLAB can access both COMSOL and PHREEQC software at the same time, using MATLAB as the third software can better realize the data transmission between the first software (COMSOL) and the second software (PHREEQC), and has relatively more powerful post-processing capabilities and scalability.

[0106] In some embodiments, the above-mentioned target chemical field model for the target area is constructed based on the geological parameters of the target area. When implemented specifically, it may include: constructing a target chemical field model for the target area based on the diffusion coefficient of the solute in the formation water solution, the ion composition of the formation water (or the elemental composition of the formation water solution), the pH value and PE value of the formation water, the mineral composition of the formation rock, the physical properties of each mineral (for example, density, molar mass, specific surface area, etc.), and the acid-rock reaction parameters of each mineral (for example, the reaction rate constant at 25°C, the reaction equilibrium constant, the activation energy, etc.).

[0107] In some embodiments, the above-mentioned initial state data may also be referred to as initial conditions, which may specifically be the initial values ​​of each dependent variable under the initial state of each physical field in the target acid-rock reaction multi-physics field model. Specifically, for example, it may include the initial pressure of the seepage field model, the initial temperature of the temperature field model, the initial displacement of the stress field model, and the initial concentration of the solute in the solute transport field model in the target acid-rock reaction multi-physics field model.

[0108] Among them, the initial pressure of the seepage field model, the initial temperature of the temperature field model, the initial displacement of the stress field model, and the initial concentration of the solute in the solute transport field model can be determined based on the initial parameters of the seepage field model, the initial parameters of the temperature field model, and the initial parameters of the stress field model, combined with the geological parameters of the target area. The initial concentration of the solute in the solute transport field model can be obtained by performing thermodynamic equilibrium calculations based on the initial parameters of the above models.

[0109] In specific implementation, the initial concentration of the solute in the above-mentioned solute transport field model can be determined in the following manner: use a connection program written based on a third software (e.g., MATLAB) to extract the initial parameters of the relevant model in the target acid-rock reaction multi-physics field model and input them into the second software (e.g., PHREEQC) for thermodynamic equilibrium calculation; then input the calculation results from the second software into the target acid-rock reaction multi-physics field model in the first software (e.g., COMSOL) through the connection program.

[0110] The above-mentioned thermodynamic equilibrium calculation may specifically refer to taking the law of mass action as a condition for judging the equilibrium of the entire chemical system, and then combining the element equilibrium equation and the charge equilibrium equation to solve the equilibrium state of the entire chemical system through an iterative method.

[0111] Specifically, the following mass action equation for aqueous solution can be constructed according to the law of mass action:

[0112]

[0113] Among them, K i is the equilibrium constant of the chemical reaction, which is related to temperature; a i is the activity of substance i in aqueous solution, mol / m 3 ; a m is the activity of the main substance m in aqueous solution, mol / m 3 ;M aq is the amount of the main substance in the aqueous solution during the reaction, dimensionless; c m,i It is the stoichiometric coefficient of the main substance m in substance i. The specific value can be positive or negative. In this model, the term on the right side of the association reaction is designated as a negative value, and the term on the left side is a positive value and dimensionless.

[0114] The element balance equation is constructed according to the following formula to describe the balance of the total amount of substances in the system, specifically, the total number of moles of the elements in the system is the sum of the number of moles initially present in the pure phase and the aqueous phase:

[0115]

[0116] Among them, T m is the total number of moles of elements in the system, mol; N p is the number of phases in the pure phase combination, dimensionless; n p is the molar number of each phase in the pure phase combination, mol; b m,p is the number of moles of element m per mole of pure phase, mol; N aq is the total amount of substance in the aqueous solution, dimensionless; n i is the number of moles of substance i in aqueous solution, mol; b m,i is the number of moles of element m in the substance per mole of aqueous solution, mol; where function f m When the value of is zero, the system reaches equilibrium in terms of total mass.

[0117] In specific implementation, based on the above equation, the charge balance equation can be used to adjust the pH value of the aqueous solution to make the aqueous solution electrically neutral, or to calculate the activity of the main substance in the aqueous solution, thereby obtaining the total concentration of each element. The charge balance equation in the entire chemical system can be expressed as follows:

[0118]

[0119] Among them, T z is the charge imbalance in the system, mol; zi is the ionic charge number of substance i in aqueous solution, dimensionless; where, when the function f z When the value of is zero, the system reaches charge equilibrium.

[0120] In the above processing, the parameters output from the first software and input to the second software may specifically include: temperature, porosity, water content, mineral components, elements in the aqueous solution, and ion composition and concentration, etc. The parameters output from the second software and input to the first software may specifically include: porosity, mineral components, ion composition and concentration, etc.

[0121] In some embodiments, the above boundary conditions may specifically refer to the values ​​of various field dependent variables at the boundary of the target acid-rock reaction multi-physics model. Similar to the initial state data, the pressure, temperature, displacement, etc. at the boundary may be determined according to the boundary parameters of the target acid-rock reaction multi-physics model. The concentration of the solute at the boundary may be calculated by thermodynamic equilibrium based on the above data values.

[0122] In specific implementation, a connection program written based on a third software (e.g., MATLAB) can be used to input boundary parameters such as pressure, temperature, displacement, etc. in the target acid-rock reaction multi-physics field model into a second software (e.g., PHREEQC) for thermodynamic equilibrium calculation, and the calculation results can be input into the target acid-rock reaction multi-physics field model of the first software (e.g., COMSOL), thereby determining the boundary conditions of the target acid-rock reaction multi-physics field model.

[0123] In some embodiments, the above-mentioned sequential non-iterative (SNIA) processing may specifically refer to implementing multiple physical fields such as the heat-fluid-solid-solute transport field in the target acid-rock reaction multi-physics field model through a fully coupled iterative approach, that is, using an iterative method to solve all physical field equations at the same time. As for the coupling of multiple physical fields such as the chemical field and the heat-fluid-solid-solute transport field, the fully coupled model of the heat-fluid-solid-solute transport field is first solved for calculation, and then the physical field calculation result data (for example, solute data) is input into the chemical field for reaction kinetics calculation at the same time step; it should be noted that the solution here is not simultaneous, but sequential, and naturally no iterative method is used, so it is sequential non-iterative.

[0124] In some embodiments, see Figure 2 As shown, the above-mentioned sequential non-iterative processing is performed by using the target acid-rock reaction multi-physics field model and the target chemical field model according to the initial state data and the initial boundary conditions. Specifically, the implementation may include: performing the current round of sequential non-iterative processing in the following manner:

[0125] S1: using the first software, based on the initial state data of the current round and the initial boundary conditions of the current round, performing a multi-physics field coupling simulation operation of the target acid-rock reaction multi-physics field model at a specified time step of the current round, and obtaining a physical field simulation operation result of the target acid-rock reaction multi-physics field model of the current round;

[0126] S2: using a connection program generated based on the third software, through the first software, obtaining the current round of physical field simulation calculation results of the target acid-rock reaction multi-physical field model; and importing the current round of physical field simulation calculation results of the target acid-rock reaction multi-physical field model into the target chemical field model in the second software;

[0127] S3: Using the second software, based on the current round of physical field simulation operation results of the target acid-rock reaction multi-physical field model, perform chemical reaction kinetics operation of the current round of the specified time step on the target chemical field model to obtain the current round of chemical field simulation operation results of the target chemical field model.

[0128] In specific implementation, the above steps can be repeated for multiple rounds of simulation until the cutoff condition is reached.

[0129] Among them, in the same round of simulation, the specified time step in the target chemical long model is the same time step as the specified time step in the target acid-rock reaction multi-physics field model.

[0130] During specific implementation, the above-specified time step may be set to a smaller time step according to specific business needs and the total simulation duration.

[0131] Specifically, for example, if the total simulation duration is 1000 years, the specified time step may be set to 1 year. For another example, if the total simulation duration is 30 days, the specified time step may be set to 1 hour.

[0132] Alternatively, you can set a shorter time step as the specified time step in the early stages of the simulation, and a longer time step as the specified time step in the middle and late stages of the simulation.

[0133] The above physical field simulation calculation results may specifically include the temperature, porosity, water content, concentration of each ion in the aqueous solution, etc. at each grid node. In specific implementation, the above physical field simulation calculation results and the corresponding mineral composition parameters of each grid node may be input into the target chemical field model of the second software using a connection program generated based on the third software.

[0134] In specific implementation, the above-mentioned specified time step can be set to a relatively short time step to improve the simulation accuracy, and then the chemical field and the heat-fluid-solid-solute transport field can be more accurately and effectively coupled through sequential non-iterative (SNIA) processing to obtain a relatively good simulation effect.

[0135] The above chemical field simulation calculation results may specifically include: the porosity at each grid node, the element composition and concentration of the aqueous solution (or the ion composition and concentration of the aqueous solution), etc.

[0136] In specific implementation, the third software can also be used to call the second software to perform parallel operations on reactants at multiple grid nodes in the target chemical field model to shorten data processing time and more efficiently obtain the chemical field simulation operation results of the current round of the target chemical field model.

[0137] In some embodiments, when using the second software to perform chemical reaction kinetics calculations for the specified time step of the current round on the target chemical field model based on the physical field simulation calculation results of the current round of the target acid-rock reaction multi-physics field model, it is also necessary to calculate the dissolution / precipitation reaction rate of each mineral, which can be calculated specifically according to the following formula:

[0138]

[0139] Where r is the mineral dissolution / precipitation reaction rate, mol / (kgw·s); S A is the surface area of ​​minerals per unit water, m 2 / kgw; A is the prefactor, mol / (m 2 ·s); E a is the reaction activation energy, J / mol; Ω is the mineral saturation quotient, dimensionless; n is the reaction order, dimensionless; p and q are empirical exponents, dimensionless; H, w, and OH represent acidic environment, neutral environment, and alkaline environment, respectively.

[0140] In some embodiments, after obtaining the current round of chemical field simulation operation results of the target chemical field model, the method may also include the following contents when implemented: detecting whether the cutoff condition is met; if it is determined that the cutoff condition is not met, updating the initial state data and initial boundary conditions of the next round according to the current round of chemical field simulation operation results of the target chemical field model to perform the next round of sequential non-iterative processing.

[0141] Among them, detecting whether the cut-off condition is met may specifically include: detecting whether the sum of the accumulated specified time steps is greater than or equal to the total simulation time; and determining that the cut-off condition is met when it is determined that the sum of the accumulated specified time steps is greater than or equal to the total simulation time. It may also include: detecting whether the number of rounds of the current round simulation reaches the specified number of rounds; and determining that the cut-off condition is met when it is determined that the number of rounds of the current round simulation reaches the specified number of rounds.

[0142] When it is determined that the cutoff conditions are met, the simulation is terminated and the next round of sequential non-iterative processing is no longer triggered; the final target acid-rock reaction simulation results for the target area are determined based on the current round of physical field simulation results of the target acid-rock reaction multi-physical field model and the current round of chemical field simulation results of the target chemical field model.

[0143] The above-mentioned target acid-rock reaction simulation results can be specifically understood as an acid-rock reaction result under multi-field coupling, which can specifically include the reaction process and the final temperature field data, pressure field data, displacement field data, stress field data, as well as porosity distribution, permeability distribution, rock mineral component changes, solute component changes in aqueous solution and other data.

[0144] In some embodiments, after obtaining the target acid-rock reaction simulation results for the target area, the method may further include the following when implemented:

[0145] S1: Based on the target acid-rock reaction simulation results, predict the target change trends of reservoir physical properties, pore structure, and mineral surface properties in the target area;

[0146] S2: Determine the oil recovery effect of CO2 injection in the target area according to the target change trend.

[0147] Subsequently, the above-mentioned oil recovery effects can be used to more effectively guide CO2 injection for oil recovery in the target area.

[0148] In addition, the target change trends of reservoir physical properties, pore structure, and mineral surface properties in the target area can be used to predict the stability of geological storage in the target area to guide CCUS (carbon capture, utilization and storage) in the target area.

[0149] As can be seen from the above, based on the data processing method of the CO2 injection acid-rock reaction provided in the embodiment of this specification, the geological parameters of the target area can be obtained first; and according to the geological parameters of the target area, a target acid-rock reaction multi-physics field model and a target chemical field model for the target area are constructed; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, stress field, temperature field, and solute migration field; and then the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model are determined; according to the initial state data and initial boundary conditions, the target acid-rock reaction multi-physics field model and the target chemical field model are used for sequential non-iterative processing to obtain the target acid-rock reaction simulation results for the target area. Thus, the acid-rock reaction simulation for a long period of time under complex conditions can be accurately and comprehensively realized, and the target acid-rock reaction simulation results with high reference value can be obtained, which can better guide the subsequent determination of CO2 injection for oil recovery and CO2 storage in the target area. Furthermore, the stability of the geological storage of CO2 saltwater layers in the target area can also be accurately evaluated based on CO2 storage.

[0150] The embodiment of the present specification also provides a server, including a processor and a memory for storing processor executable instructions. When the processor is implemented, it can perform the following steps according to the instructions: obtain geological parameters of the target area; construct a target acid-rock reaction multi-physics field model and a target chemical field model for the target area according to the geological parameters of the target area; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, stress field, temperature field, and solute migration field; determine the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model; according to the initial state data and initial boundary conditions, obtain the target acid-rock reaction simulation results for the target area by using the target acid-rock reaction multi-physics field model and the target chemical field model for sequential non-iterative processing.

[0151] In order to complete the above instructions more accurately, refer to Figure 3 As shown, the embodiment of this specification also provides another specific server, wherein the server includes a network communication port 301, a processor 302 and a memory 303, and the above structures are connected through internal cables so that each structure can perform specific data interaction.

[0152] The network communication port 301 can be used to obtain geological parameters of the target area.

[0153] The processor 302 can be specifically used to construct a target acid-rock reaction multi-physics field model and a target chemical field model for the target area according to the geological parameters of the target area; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, stress field, temperature field, and solute migration field; determine the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model; and obtain the target acid-rock reaction simulation results for the target area by performing sequential non-iterative processing using the target acid-rock reaction multi-physics field model and the target chemical field model according to the initial state data and the initial boundary conditions.

[0154] The memory 303 may be specifically used to store corresponding instruction programs.

[0155] In this embodiment, the network communication port 301 can be a virtual port that is bound to different communication protocols so that different data can be sent or received. For example, the network communication port can be a port responsible for web data communication, a port responsible for FTP data communication, or a port responsible for email data communication. In addition, the network communication port can also be a physical communication interface or communication chip. For example, it can be a wireless mobile network communication chip, such as GSM, CDMA, etc.; it can also be a Wifi chip; it can also be a Bluetooth chip.

[0156] In this embodiment, the processor 302 may be implemented in any appropriate manner. For example, the processor may take the form of a microprocessor or processor and a computer-readable medium storing a computer-readable program code (such as software or firmware) executable by the (micro)processor, a logic gate, a switch, an application specific integrated circuit (ASIC), a programmable logic controller, and an embedded microcontroller, etc. This specification does not limit this.

[0157] In this embodiment, the memory 303 may include multiple levels. In a digital system, anything that can store binary data can be a memory; in an integrated circuit, a circuit with a storage function that has no physical form is also called a memory, such as RAM, FIFO, etc.; a storage device with a physical form in the system is also called a memory, such as a memory stick, TF card, etc.

[0158] The embodiments of this specification also provide a computer-readable storage medium based on the data processing method of the above-mentioned CO2 injection acid-rock reaction, wherein the computer-readable storage medium stores computer program instructions, and when the computer program instructions are executed, the following are achieved: obtaining geological parameters of the target area; constructing a target acid-rock reaction multi-physical field model and a target chemical field model for the target area according to the geological parameters of the target area; wherein the target acid-rock reaction multi-physical field model is a physical field model that couples the seepage field, stress field, temperature field, and solute migration field; determining the initial state data and initial boundary conditions of the target acid-rock reaction multi-physical field model; and obtaining the target acid-rock reaction simulation results for the target area by sequentially non-iteratively processing using the target acid-rock reaction multi-physical field model and the target chemical field model according to the initial state data and the initial boundary conditions.

[0159] In this embodiment, the storage medium includes, but is not limited to, a random access memory (RAM), a read-only memory (ROM), a cache, a hard disk (HDD), or a memory card. The memory may be used to store computer program instructions. The network communication unit may be an interface for network connection communication set in accordance with the standard specified by the communication protocol.

[0160] In this embodiment, the functions and effects specifically implemented by the program instructions stored in the computer-readable storage medium can be explained in comparison with other implementations and will not be described in detail here.

[0161] See also Figure 4 As shown, at the software level, the embodiment of this specification also provides a data processing device for CO2 injection acid rock reaction, which may specifically include the following structural modules:

[0162] An acquisition module 401 may be used to acquire geological parameters of a target area;

[0163] The construction module 402 can be specifically used to construct a target acid-rock reaction multi-physics field model and a target chemical field model for the target area according to the geological parameters of the target area; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, the stress field, the temperature field, and the solute migration field;

[0164] The determination module 403 may be specifically used to determine the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model;

[0165] The processing module 404 can be specifically used to obtain the target acid-rock reaction simulation results for the target area by performing sequential non-iterative processing based on the initial state data and initial boundary conditions by using the target acid-rock reaction multi-physics field model and the target chemical field model.

[0166] In some embodiments, after obtaining the target acid-rock reaction simulation results for the target area, the device, when implemented, can also be used to predict the target change trends of reservoir physical properties, pore structure, and mineral surface properties in the target area based on the target acid-rock reaction simulation results; and determine the oil recovery effect of CO2 injection in the target area based on the target change trends.

[0167] In some embodiments, the geological parameters of the target area may specifically include at least one of the following: reservoir parameters, formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, original formation pressure, elastic modulus of formation rock, elastic modulus of formation rock skeleton, Poisson's ratio of formation rock, density of formation rock skeleton, thermal expansion coefficient of formation rock, original formation temperature, constant-pressure specific heat capacity of formation rock, thermal conductivity of formation rock, constant-pressure specific heat capacity of formation water, thermal conductivity of formation water, diffusion coefficient of solute in formation aqueous solution, elemental composition of formation aqueous solution, pH value and PE value of formation water, mineral composition of formation rock, etc.

[0168] In some embodiments, when the above-mentioned construction module 402 is implemented, a target acid-rock reaction multi-physical field model and a target chemical field model for the target area can be constructed according to the geological parameters of the target area in the following manner: using the first software to construct the target acid-rock reaction multi-physical field model for the target area according to the geological parameters of the target area; using the second software to construct the target chemical field model for the target area according to the geological parameters of the target area.

[0169] In some embodiments, when the above-mentioned construction module 402 is specifically implemented, a target acid-rock reaction multi-physics field model for the target area can be constructed according to the geological parameters of the target area in the following manner: construct a geological model of the target area according to the reservoir parameters of the target area; establish a seepage field model in the geological model of the target area according to the formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, and original formation pressure of the target area; establish a stress field model in the geological model of the target area according to the elastic modulus of the formation rock, the elastic modulus of the formation rock skeleton, the Poisson's ratio of the formation rock, the density of the formation rock skeleton, and the thermal expansion coefficient of the formation rock in the target area; establish a temperature field model in the geological model of the target area according to the original formation temperature of the target area, the constant-pressure specific heat capacity of the formation rock, the thermal conductivity of the formation rock, the constant-pressure specific heat capacity of the formation water, and the thermal conductivity of the formation water; and couple the seepage field model, the stress field model, the temperature field model, and the solute migration field model according to the mutual influence relationship of the key parameters in the model to obtain the target acid-rock reaction multi-physics field model.

[0170] In some embodiments, after the target acid-rock reaction multi-physics field model is obtained by coupling the seepage field model, stress field model, temperature field model, and solute migration field model according to the mutual influence relationship of the key parameters in the model, the construction module 402, when implemented specifically, can also be used to set the characteristic information of the reactants at each grid node in the target acid-rock reaction multi-physics field model according to the geological parameters of the target area.

[0171] In some embodiments, when the above-mentioned processing module 404 is specifically implemented, the current round of sequential non-iterative processing can be performed in the following manner according to the initial state data and initial boundary conditions by using the target acid-rock reaction multi-physics field model and the target chemical field model: using the first software, based on the initial state data of the current round and the initial boundary conditions of the current round, the target acid-rock reaction multi-physics field model is subjected to multi-physics field coupling simulation operation of the specified time step of the current round to obtain the physical field simulation operation result of the current round of the target acid-rock reaction multi-physics field model; using the connection program generated based on the third software, through the first software, the physical field simulation operation result of the current round of the target acid-rock reaction multi-physics field model is obtained; and the physical field simulation operation result of the current round of the target acid-rock reaction multi-physics field model is imported into the target chemical field model in the second software; using the second software, based on the physical field simulation operation result of the current round of the target acid-rock reaction multi-physics field model, the target chemical field model is subjected to chemical reaction kinetic operation of the specified time step of the current round to obtain the chemical field simulation operation result of the current round of the target chemical field model.

[0172] In some embodiments, after obtaining the current round of chemical field simulation operation results of the target chemical field model, the device can also be used to detect whether the cutoff condition is met during its specific implementation; when it is determined that the cutoff condition is not met, the initial state data and initial boundary conditions of the next round are updated according to the current round of chemical field simulation operation results of the target chemical field model to perform the next round of sequential non-iterative processing.

[0173] It should be noted that the units, devices or modules described in the above embodiments can be implemented by computer chips or entities, or by products with certain functions. For the convenience of description, the above devices are described separately by functions divided into various modules. Of course, when implementing this specification, the functions of each module can be implemented in the same or more software and / or hardware, or the modules that implement the same function can be implemented by a combination of multiple sub-modules or sub-units. The device embodiments described above are only schematic. For example, the division of the units is only a logical function division. There may be other division methods in actual implementation, such as multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interfaces, indirect coupling or communication connection of devices or units, which can be electrical, mechanical or other forms.

[0174] As can be seen from the above, the data processing device for the CO2 injection acid-rock reaction provided in the embodiments of this specification can accurately and comprehensively realize the acid-rock reaction simulation for a long period of time under complex conditions, and obtain the target acid-rock reaction simulation results with high reference value. And thus, it can better guide the subsequent CO2 injection oil flooding in the target area.

[0175] In a specific scenario example, the data processing method for CO2 injection acid rock reaction provided in this specification can be used to perform multi-field and multi-scale coupled simulation of CO2 injection acid rock reaction heat-fluid-solid-chemical (THMC). Figure 5 , refer to the following content for implementation.

[0176] In this scenario example, see Figure 5As shown, a fully coupled numerical model of the temperature field, seepage field, solid mechanics field and solute transport field of the CO2 injection acid-rock reaction can be first established in COMSOL software (corresponding to the first software) (corresponding to the target acid-rock reaction multi-physics field model); then the initial and boundary parameters (corresponding to the initial state data and initial boundary conditions) of the COMSOL model are obtained through MATLAB software (corresponding to the third software) and input into PHREEQC software (corresponding to the third software) to perform geochemical reaction thermodynamic equilibrium calculation to obtain the initial conditions and boundary conditions of the model; then the initial conditions and boundary conditions are input into the COMSOL model through MATLAB software and a simulation of a specified time step is performed; the COMSOL model calculation results are input into PHREEQC software to perform geochemical reaction kinetics calculation of the same time step to obtain the chemical field results of the time step; the chemical field results are input into the COMSOL model as the initial conditions for the next time step model calculation to calculate, and the process is repeated for sequential non-iterative processing until the simulation is performed according to all time steps to realize the THMC coupling simulation of the CO2-water-rock reaction and obtain the final THMC coupling simulation results of the CO2 injection acid-rock reaction (corresponding to the target acid-rock reaction simulation results). Specifically, you can follow the steps below.

[0177] Step 1: Use COMSOL software to establish a fully coupled mathematical model of the acid-rock reaction seepage field, stress field, temperature field, and solute migration field (for example, construct a target acid-rock reaction multi-physics field model).

[0178] Step 2: Set the main reaction species on each grid node of the acid-rock reaction model.

[0179] Specifically, a program can be written through MATLAB software to obtain the grid information of the COMSOL model and set the molar amount of minerals involved in the acid-rock reaction, the ion composition and concentration in the aqueous solution, etc. at each grid node.

[0180] Step 3: Set the initial conditions for the acid-rock reaction model.

[0181] Specifically, a program was written using MATLAB software to input the initial parameters in the acid-rock reaction model into the PHREEQC software for thermodynamic equilibrium calculation, and the calculation results were input into the COMSOL model as its initial conditions. Among them, the thermodynamic equilibrium in the PHREEQC software uses the law of mass action as the condition for judging the equilibrium of the entire chemical system, and then combines the element equilibrium equation and the charge balance equation to solve the equilibrium state of the entire chemical system through an iterative method.

[0182] Specifically, the parameters input into the PHREEQC software include temperature, porosity, water content, mineral composition, elements in the aqueous solution, and ion composition and concentration; the parameters input into the COMSOL model as initial conditions after the PHREEQC thermodynamic equilibrium calculation include porosity, mineral composition, ion composition and concentration in the aqueous solution.

[0183] Step 4: Set the boundary conditions for the acid-rock reaction model.

[0184] Specifically, a program was written using MATLAB software to input the boundary parameters in the acid-rock reaction model into PHREEQC software for thermodynamic equilibrium calculation, and the calculation results were input into the COMSOL model as its boundary conditions.

[0185] Step 5: Use the COMSOL model to perform a fully coupled simulation of THM and solute transport at a specified time step, and input the calculated parameters of each grid node into the PHREEQC software.

[0186] Specifically, a program was written using MATLAB software to obtain the temperature, porosity, water content, and concentration of each ion in the aqueous solution of each grid node in the COMSOL model calculation results, and input them into the PHREEQC software together with the mineral component parameters of each grid node.

[0187] Step 6: Use PHREEQC software to perform chemical reaction kinetics calculations at the same time step and input the calculation results into the COMSOL model as the initial conditions for the next time step simulation.

[0188] The same time step means that the time step used in the PHREEQC chemical reaction kinetics simulation is consistent with the time step used in the COMSOL model calculation in step 3.

[0189] After the PHREEQC reaction kinetics calculation, the parameters input into the COMSOL model as the initial conditions for the next time step simulation include the porosity of each node, the ion composition and concentration in the aqueous solution.

[0190] Specifically, MATLAB software was used to write a program to establish the PHREEQC model and perform chemical reaction calculations. After the calculations, the calculation results were obtained and input into the COMSOL model. In order to save calculation time, MATLAB programming was used to implement parallel calculations of the reaction kinetics of each node on the model grid of PHREEQC.

[0191] Step 7: Repeat steps 5 and 6 until the THMC coupling simulation results of CO2 injection acid rock reaction are obtained according to all time steps.

[0192] Since this simulation method uses the sequential non-iterative method (SNIA) to couple the chemical field and the heat-fluid-solid-solute transport field, the time step should be set as small as possible. Too large a time step will lead to inaccurate simulation results.

[0193] In this scenario example, the one-dimensional reaction transport model is used for verification results. Figures 6 to 11 As shown. Among them, Figure 6 This is a schematic diagram of the verification results of the molar amount of calcite; Figure 7 This is a schematic diagram of the verification results of the molar amount of dolomite; Figure 8 It is a schematic diagram of pH verification results in aqueous solution; Fig. 9 This is a schematic diagram of the verification results of the molar amount of Ca element; Fig.10 This is a schematic diagram of the verification results of the molar amount of Mg element; Fig.11 This is a schematic diagram of the verification results of the molar amount of Cl element.

[0194] Through the above scenario examples, combined with specific verification results, it can be proved that based on the data processing method of CO2 injection acid-rock reaction provided in this specification, on the one hand, it is possible to establish the THMC four-field coupling model of CO2-water-rock reaction, which can more accurately describe the CO2 injection acid-rock reaction process and its impact on reservoir properties, and provide theoretical support for CCUS projects; on the other hand, it is possible to connect COMSOL and PHREEQC software through MATLAB software programming, so that the present invention has good applicability and greater scalability, and can realize the simulation of water-rock reaction under various conditions; on the other hand, since CO2-water-rock reaction is usually slow, the method provided in this specification can obtain long-term water-rock reaction results that cannot be obtained by physical simulation experiments.

[0195] Although the present specification provides method operation steps as described in the embodiments or flow charts, more or less operation steps may be included based on conventional or non-creative means. The order of steps listed in the embodiments is only one way of executing the order of many steps, and does not represent a unique execution order. When the device or client product in practice is executed, it can be executed in sequence or in parallel according to the method shown in the embodiments or the drawings (for example, a parallel processor or a multi-threaded processing environment, or even a distributed data processing environment). The term "include", "include" or any other variant thereof is intended to cover non-exclusive inclusion, so that the process, method, product or device including a series of elements includes not only those elements, but also includes other elements that are not explicitly listed, or also includes elements inherent to such a process, method, product or device. In the absence of more restrictions, it is not excluded that there are other identical or equivalent elements in the process, method, product or device including the elements. The first, second, etc. words are used to represent the name, and do not represent any particular order.

[0196] Those skilled in the art also know that, in addition to implementing the controller in a purely computer-readable program code, the controller can be made to implement the same function in the form of logic gates, switches, application-specific integrated circuits, programmable logic controllers, and embedded microcontrollers by logically programming the method steps. Therefore, such a controller can be considered as a hardware component, and the devices for implementing various functions included therein can also be considered as structures within the hardware component. Or even, the devices for implementing various functions can be considered as both software modules for implementing the method and structures within the hardware component.

[0197] This specification may be described in the general context of computer-executable instructions executed by a computer, such as program modules. Generally, program modules include routines, programs, objects, components, data structures, classes, etc. that perform specific tasks or implement specific abstract data types. This specification may also be practiced in distributed computing environments where tasks are performed by remote processing devices connected through a communication network. In a distributed computing environment, program modules may be located in local and remote computer-readable storage media including storage devices.

[0198] Through the description of the above embodiments, it can be known that those skilled in the art can clearly understand that the present specification can be implemented by means of software plus a necessary general hardware platform. Based on such an understanding, the technical solution of the present specification can essentially be embodied in the form of a software product, which can be stored in a storage medium such as ROM / RAM, a disk, an optical disk, etc., and includes a number of instructions for enabling a computer device (which can be a personal computer, a mobile terminal, a server, or a network device, etc.) to execute the methods described in each embodiment of the present specification or some parts of the embodiments.

[0199] The various embodiments in this specification are described in a progressive manner, and the same or similar parts between the various embodiments can be referred to each other, and each embodiment focuses on the differences from other embodiments. This specification can be used in many general or special computer system environments or configurations. For example: personal computers, server computers, handheld devices or portable devices, tablet devices, multi-processor systems, microprocessor-based systems, set-top boxes, programmable electronic devices, network PCs, minicomputers, mainframe computers, distributed computing environments including any of the above systems or devices, etc.

[0200] Although the present specification is described through embodiments, those skilled in the art will appreciate that there are many modifications and changes to the present specification without departing from the spirit of the present specification, and it is intended that the appended claims include these modifications and changes without departing from the spirit of the present specification.

Claims

1. A data processing method for CO2 injection acid rock reaction, characterized in that: include: Obtain geological parameters of the target area; According to the geological parameters of the target area, a target acid-rock reaction multi-physics field model and a target chemical field model for the target area are constructed; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, stress field, temperature field, and solute migration field; wherein the stress field model is constructed according to the following formula: Where s is the displacement vector; σ is the stress, Pa; f v is the body force, which includes gravity, body load, and body force caused by thermal stress, fluid pressure, and adsorption stress; Determine the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics model; According to the initial state data and initial boundary conditions, by using the target acid-rock reaction multi-physics field model and the target chemical field model to perform sequential non-iterative processing, a target acid-rock reaction simulation result about the target area is obtained; according to the target acid-rock reaction simulation result, the determination of CO2 injection for oil recovery and CO2 storage in the target area is guided; Among them, according to the geological parameters of the target area, a target acid-rock reaction multi-physics field model for the target area is constructed, including: according to the geological parameters of the target area, a seepage field model, a stress field model, and a temperature field model for the target area are respectively established; the relationship between the density and viscosity of the fluid in the seepage field model and the temperature is obtained, and according to the influence of temperature and volume strain on the porosity and permeability, the influence of the convective heat transfer part in the temperature field model on the fluid flow rate, and the correlation between the effective heat capacity and the stress field model, the seepage field model, the stress field model, and the temperature field model are coupled to obtain the THM coupling model; on the basis of the THM coupling model, the solute migration model is coupled to obtain the target acid-rock reaction multi-physics field model; Furthermore, based on the amount of CO2 to be injected into the target area and the geological parameters of the target area, the hydrogen ion concentration and carbonate ion concentration in the corresponding formation aqueous solution are calculated, and the elemental composition and concentration of the formation aqueous solution after CO2 injection are obtained and utilized, and the target acid-rock reaction multi-physics field model and the target chemical field model are combined to simulate the CO2-water-mineral reaction in the target area for a long period of time.

2. The method according to claim 1, characterized in that After obtaining the target acid-rock reaction simulation results for the target area, the method further includes: According to the simulation results of the target acid-rock reaction, predict the target change trend of reservoir physical properties, pore structure and mineral surface properties in the target area; According to the target change trend, the oil recovery effect of CO2 injection in the target area is determined.

3. The method according to claim 1, characterized in that The geological parameters of the target area include at least one of the following: reservoir parameters, formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, original formation pressure, elastic modulus of formation rock, elastic modulus of formation rock skeleton, Poisson's ratio of formation rock, density of formation rock skeleton, thermal expansion coefficient of formation rock, original formation temperature, constant-pressure specific heat capacity of formation rock, thermal conductivity of formation rock, constant-pressure specific heat capacity of formation water, thermal conductivity of formation water, diffusion coefficient of solute in formation aqueous solution, elemental composition of formation aqueous solution, pH value and PE value of formation water, and mineral composition of formation rock.

4. The method according to claim 3, characterized in that According to the geological parameters of the target area, a target acid-rock reaction multi-physics field model and a target chemical field model for the target area are constructed, including: Using the first software to construct a target acid-rock reaction multi-physics field model for the target area according to geological parameters of the target area; The second software is used to construct a target chemical field model for the target area according to the geological parameters of the target area.

5. The method according to claim 3, characterized in that: According to the geological parameters of the target area, a multi-physics field model of the target acid-rock reaction is constructed for the target area, including: Construct a geological model of the target area based on the reservoir parameters of the target area; A seepage field model is established in the geological model of the target area according to the formation porosity, formation permeability, formation rock density, formation water density, formation water viscosity, and original formation pressure of the target area; a stress field model is established in the geological model of the target area according to the elastic modulus of the formation rock, the elastic modulus of the formation rock skeleton, the Poisson's ratio of the formation rock, the density of the formation rock skeleton, and the thermal expansion coefficient of the formation rock; a temperature field model is established in the geological model of the target area according to the original formation temperature of the target area, the constant-pressure specific heat capacity of the formation rock, the thermal conductivity of the formation rock, the constant-pressure specific heat capacity of the formation water, and the thermal conductivity of the formation water; According to the mutual influence relationship of key parameters in the model, the seepage field model, stress field model, temperature field model, and solute migration field model are coupled to obtain the target acid-rock reaction multi-physics field model.

6. The method according to claim 5, characterized in that After the target acid-rock reaction multi-physics field model is obtained by coupling the seepage field model, the stress field model, the temperature field model, and the solute migration field model according to the mutual influence relationship of the key parameters in the model, the method further includes: According to the geological parameters of the target area, characteristic information of reactants at each grid node in the target acid-rock reaction multi-physics field model is set.

7. The method according to claim 6, characterized in that According to the initial state data and initial boundary conditions, sequential non-iterative processing is performed by using the target acid-rock reaction multi-physics field model and the target chemical field model, including: The sequential non-iterative processing of the current round is performed as follows: Using the first software, based on the initial state data of the current round and the initial boundary conditions of the current round, a multi-physics field coupling simulation operation of the target acid-rock reaction multi-physics field model with a specified time step of the current round is performed to obtain the physical field simulation operation result of the current round of the target acid-rock reaction multi-physics field model; Using a connection program generated based on the third software, the first software is used to obtain the current round of physical field simulation calculation results of the target acid-rock reaction multi-physical field model; and the current round of physical field simulation calculation results of the target acid-rock reaction multi-physical field model are imported into the target chemical field model in the second software; Using the second software, based on the current round of physical field simulation operation results of the target acid-rock reaction multi-physical field model, the target chemical field model is subjected to the current round of chemical reaction kinetics operation of the specified time step to obtain the current round of chemical field simulation operation results of the target chemical field model.

8. The method according to claim 7, characterized in that After obtaining the chemical field simulation calculation result of the current round of the target chemical field model, the method further includes: Check whether the cut-off condition is met; When it is determined that the cutoff condition is not met, the initial state data and initial boundary conditions of the next round are updated according to the chemical field simulation operation results of the current round of the target chemical field model to perform the next round of sequential non-iterative processing.

9. A data processing device for CO2 injection acid rock reaction, characterized in that: include: An acquisition module, used to obtain geological parameters of the target area; A construction module is used to construct a target acid-rock reaction multi-physics field model and a target chemical field model for the target area according to the geological parameters of the target area; wherein the target acid-rock reaction multi-physics field model is a physical field model that couples the seepage field, the stress field, the temperature field, and the solute migration field; wherein the stress field model is constructed according to the following formula: Where s is the displacement vector; σ is the stress, Pa; f v is the body force, which includes gravity, body load, and body force caused by thermal stress, fluid pressure, and adsorption stress; A determination module, used to determine the initial state data and initial boundary conditions of the target acid-rock reaction multi-physics field model; A processing module is used to obtain a target acid-rock reaction simulation result about a target area by performing sequential non-iterative processing based on the initial state data and initial boundary conditions by using the target acid-rock reaction multi-physics field model and the target chemical field model; and guide the determination of CO2 injection for oil recovery and CO2 storage in the target area based on the target acid-rock reaction simulation result; The construction module is specifically used to: establish the seepage field model, stress field model and temperature field model of the target area according to the geological parameters of the target area; obtain and couple the seepage field model, stress field model and temperature field model according to the relationship between the density and viscosity of the fluid in the seepage field model and the temperature, the influence of temperature and volume strain on the porosity and permeability, the influence of fluid flow rate on the convective heat transfer part in the temperature field model, and the correlation between the effective heat capacity and the stress field model to obtain the THM coupling model; on the basis of the THM coupling model, couple the solute migration model to obtain the multi-physics field model of the target acid-rock reaction; In addition, the device is also used to calculate the hydrogen ion concentration and carbonate ion concentration in the corresponding formation water solution according to the amount of CO2 to be injected into the target area and the geological parameters of the target area, obtain and use the elemental composition and concentration of the formation water solution after the CO2 injection, and simulate the CO2-water-mineral reaction in the target area for a long time in combination with the target acid-rock reaction multi-physical field model and the target chemical field model.

10. A computer-readable storage medium, characterized in that: Computer instructions are stored thereon, and when the instructions are executed by a processor, the steps of the method according to any one of claims 1 to 8 are implemented.

Images