High-precision dynamic acquisition method for barite precipitation amount in reservoir modeling

Abstract

The invention relates to the technical field of resource exploration and development, in particular to a high-precision dynamic acquisition method for barite precipitation capacity in reservoir modeling, which comprises the following steps of: loading key reaction parameters of a fluid-rock system of a reservoir into a three-dimensional pore model; and carrying out numerical simulation of mineral dissolution and barite precipitation processes by using the three-dimensional pore model loaded with the key reaction parameters, and calculating the barite precipitation amount based on the numerical simulation of the mineral dissolution and barite precipitation processes. Through calculation of the barite precipitation amount, the reservoir damage risk caused by barite precipitation can be effectively identified, decision support is provided for reservoir protection measures, and potential loss in the reservoir development process is reduced; and an injection-production scheme can be optimized, the oil and gas recovery efficiency is improved, and the economic life of a reservoir is prolonged. Besides, by calculating the precipitation amount in the barite precipitation simulation process, development and evolution of the reservoir can be accurately predicted, and a reliable basis is provided for reservoir modeling and simulation.

Description

Technical Field

[0001] The present invention relates to the technical field of resource exploration and development, and is a high-precision dynamic acquisition method for barite precipitation in reservoir modeling. Background Art

[0002] Underground barite precipitation and crystallization are important processes in oil and gas resource development and geological reservoir research. The driving force behind barite precipitation is the amount of solute in a solution that exceeds its solubility, a phenomenon known as supersaturation. Supersaturation is typically caused by fluid mixing, pressure changes, temperature changes, or chemical reactions. The crystallization process primarily involves two mechanisms: nucleation and growth. Nucleation occurs when barite molecules in a solution aggregate to form new crystal nuclei, while growth occurs when barite molecules fuse into existing nuclei, gradually increasing the size of the crystals. Factors influencing barite precipitation include fluid composition, temperature and pressure, fluid flow, and the mineral composition of the rock.

[0003] Barite precipitation can have multiple impacts on reservoirs and equipment during oil and gas production. First, barite precipitation reduces reservoir porosity and permeability, hindering oil and gas flow and leading to reservoir damage, reducing well production and increasing production costs. Second, the precipitation can form solid particles near the wellbore, causing blockage and further hindering oil and gas flow. In severe cases, clearing operations may be necessary. Furthermore, barite precipitation can be accompanied by other chemical reactions, leading to equipment corrosion, shortening equipment life, and increasing maintenance costs. The precipitation can also form scale layers with other minerals, increasing flow resistance in pipelines and equipment and reducing system efficiency. During waterflooding, barite precipitation can alter reservoir wettability, affecting flooding efficiency and reducing oil and gas recovery. Furthermore, barite precipitation increases operational complexity, requiring additional monitoring and treatment measures, such as chemical treatment or mechanical clearing. Furthermore, improper handling and disposal of the precipitation can lead to environmental contamination and increase environmental risks. Ultimately, these factors increase economic costs, both direct operating expenses and indirect losses due to reduced production. Therefore, barite precipitation has a significant negative impact on oil and gas production, and effective monitoring and management measures are needed to mitigate its impact and improve production efficiency and economic benefits.

[0004] Numerical simulation methods can accurately calculate the precipitation process of barite in underground reservoirs. The simulation results can be applied to reservoir modeling, development plan design, and reservoir protection, thereby improving oil and gas recovery and reducing the risk of reservoir damage. This process is crucial for understanding and predicting changes in reservoir physical properties.

[0005] In the field of reservoir modeling and simulation, the mineral dissolution-precipitation process is an important factor affecting the physical properties of the reservoir (such as porosity and permeability). In the prior art, there has been some research on the calculation methods of the mineral dissolution-precipitation process, but it mainly focuses on the simulation of carbonate reservoirs. For example, the prior art combines experiments and numerical simulations to analyze the effects of parameters such as temperature, pressure, and fluid composition on the mineral dissolution-precipitation process. However, the prior art still has the following problems and shortcomings in the calculation method of the amount of barite (i.e., barium sulfate) precipitation:

[0006] The model is overly simplified: the nonlinear dynamic characteristics of barite precipitation are ignored.

[0007] High computational complexity: Traditional methods have difficulty balancing accuracy and efficiency.

[0008] Therefore, the deficiencies in the prior art regarding the method for calculating the amount of barite precipitation provide an improvement direction for the present invention. Summary of the Invention

[0009] The present invention provides a high-precision dynamic acquisition method for barite precipitation in reservoir modeling, which overcomes the shortcomings of the above-mentioned existing technologies. By calculating the precipitation amount during the barite precipitation simulation process, the development and evolution of the reservoir can be accurately predicted, providing a reliable basis for reservoir modeling and simulation.

[0010] The technical solution of the present invention is achieved by the following measures: a high-precision dynamic acquisition method of barite precipitation in reservoir modeling, comprising:

[0011] The key reaction parameters of the fluid-rock system of the reservoir are loaded into the three-dimensional pore model. The three-dimensional pore model loaded with the key reaction parameters is used to carry out numerical simulation of the mineral dissolution and barite precipitation process. Based on the numerical simulation of the mineral dissolution and barite precipitation process, the amount of barite precipitation is calculated. The barite precipitation amount is calculated as follows:

[0012]

[0013] m=r×t

[0014] Where m represents the amount of barite precipitation, mol / m 3 ; t represents time, s; r represents the chemical reaction rate of barium sulfate crystal formation, mol / m 3 / s;B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s; G represents the growth rate of barium sulfate crystals, m / s; n represents the particle number density, 1 / m 4 ρ c Indicates the density of barium sulfate crystals, kg / m 3; L0 represents the minimum crystal size of barium sulfate, m; k V The shape factor, which represents the volume, is used to relate the geometric properties of a particle to its size; ΔL i Indicates the length of the barium sulfate crystal size interval, m; L i represents the barium sulfate crystal size, m.

[0015] The following are further optimizations and / or improvements to the above technical solutions:

[0016] Preferably, the nucleation rate B of the barium sulfate particle source is N Calculate as follows:

[0017]

[0018] Where B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s;D AB represents the apparent diffusion coefficient, m 2 / s, describes the diffusion ability of solute in solution; K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; N A represents Avogadro's constant, 1 / mol, which is used to relate the amount of a substance to the number of particles; γ CL represents the interfacial energy, mJ / m 2 , describing the interfacial energy between the crystal and the solution; k B represents the Boltzmann constant, J / K, which is used to describe the thermal motion of microscopic particles; V m Represents the molecular volume, m 3 , reflecting the size of the space occupied by the molecules; T represents the temperature, K, which affects the thermodynamics and kinetics of the crystallization process.

[0019] Preferably, the supersaturation is calculated as follows:

[0020]

[0021] Where S represents supersaturation; γ represents activity coefficient; c Ba2+ Indicates the mass concentration of barium ions, kg / m 3 ; Indicates the mass concentration of sulfate ions, kg / m 3 ;K SP represents the solubility product, mol 2 / m 6 .

[0022] Preferably, the barium sulfate crystal growth rate is calculated as follows:

[0023]

[0024] G represents the growth rate of barium sulfate crystals, m / s; k a and k V are the shape factors for area and volume, respectively, used to relate the geometric properties of the particle to its size; Sh is the Sherwood number, which describes the mass transfer efficiency of solutes on the particle surface; M W,C represents the molar mass of barium sulfate crystals, kg / mol, reflecting the chemical composition of the crystals; ρ c Indicates the density of barium sulfate crystals, kg / m 3 , describing the physical properties of crystals; D AB represents the apparent diffusion coefficient, m 2 / s, describes the diffusion ability of solute in solution; K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; L represents the size of barium sulfate crystal particles, m.

[0025] Preferably, the particle number density is calculated as follows:

[0026]

[0027] Where L represents the size of the barium sulfate crystal particles, m; N represents the total number of barium sulfate crystal particles per unit volume, 1 / m 3 ; n represents the particle number density, 1 / m 4 , which can more accurately describe the distribution characteristics of the particle group.

[0028] Preferably, the minimum crystal size of barium sulfate is calculated as follows:

[0029]

[0030] Where L0 represents the minimum crystal size of barium sulfate, m; ν d represents the dissociation number, which describes the degree of dissociation of crystal molecules in solution; γ CL represents the interface energy, mJ / m 2 , describing the interfacial energy between the crystal and the solution; k B represents the Boltzmann constant, J / K, which is used to describe the thermal motion of microscopic particles; Vm represents the molecular volume, m 3 , reflecting the size of the space occupied by the molecules; T represents temperature, K, which affects the thermodynamics and kinetics of the crystallization process; S represents supersaturation.

[0031] The key reaction parameters of the fluid-rock system of the above reservoir include physical parameters of barium sulfate, equilibrium constant parameters, fluid dynamics parameters, interface and diffusion parameters, temperature parameters (°C), and concentration parameters at the inlet of the three-dimensional pore model.

[0032] The concentration parameters at the entrance of the three-dimensional pore model include the barium ion concentration at the entrance of the three-dimensional pore model (kmol / m 3 ), hydrosulfate ion concentration (kmol / m 3 ), sulfate ion concentration (kmol / m 3 ), chloride ion concentration (kmol / m 3 ), hydrogen ion concentration (kmol / m 3 ).

[0033] The physical parameters of barium sulfate include the molar mass (g / mol), density (g / cm 3 ), the molar mass of sulfate ion (g / mol), the molar mass of barium ion (g / mol).

[0034] The above equilibrium constant parameters include the equilibrium constant k1 (kmol / m 3 ), equilibrium constant k2(m 3 / kmol), dissolved products (mol 2 / l 2 ).

[0035] The above fluid dynamics parameters include the mixing tube diameter D (mm) of the three-dimensional pore model, the channel length Lc (mm) of the three-dimensional pore model, the Reynolds number Re, the solvent viscosity (Pa·s), the solvent density (kg / m 3 ), solvent molar mass (kg / mol), average outlet velocity of the three-dimensional pore model (m / s), average inlet velocity of the three-dimensional pore model (m / s), fluid kinematic viscosity (m 2 / s), Schmidt number, material diffusion coefficient (m 2 / s), particle diffusion coefficient (m 2 / s).

[0036] The above interface and diffusion parameters include interface energy (mJ / m 2 ), Sherwood number, apparent diffusion coefficient (m 2 / s).

[0037] The above three-dimensional pore model is a three-dimensional Y-shaped pore throat structure model or a three-dimensional irregular pore structure model.

[0038] The present invention applies simulation results to reservoir modeling, development plan design and production performance prediction.

[0039] Compared with the existing technology, it has the following significant beneficial effects:

[0040] When simulating and calculating, the present invention uses key reaction parameters close to the actual geological environment, so that the calculated barite precipitation amount is closer to the barite precipitation amount of the actual geological environment; through the calculation of the barite precipitation amount, the reservoir damage risk caused by barite precipitation can be effectively identified, and decision support can be provided for reservoir protection measures, reducing potential losses in the reservoir development process; and the injection and production scheme can be optimized, the oil and gas recovery rate can be increased, and the economic life of the reservoir can be extended. In addition, by calculating the precipitation amount during the barite precipitation simulation process, the development and evolution of the reservoir can be accurately predicted, providing a reliable basis for reservoir modeling and simulation. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] Attachment Figure 1 The modeling geometry of the 3D pore model is shown (the left figure is a 3D Y-shaped pore throat structure model, and the right figure is a 3D irregular pore structure model).

[0042] Attachment Figure 2 The modeling grid of the three-dimensional pore model modeling geometry is shown (the left picture is the modeling grid of the three-dimensional Y-shaped pore throat structure model, and the right picture is the modeling grid of the three-dimensional irregular pore structure model).

[0043] Attachment Figure 3 The time evolution of supersaturation is shown.

[0044] Attachment Figure 4 The time evolution of the mass concentrations of the reaction species and the crystalline products is shown.

[0045] Attachment Figure 5A The effect of saturation change on crystal size distribution in a three-dimensional Y-type pore throat structure model is demonstrated.

[0046] Attachment Figure 5B The effect of saturation variation on crystal size distribution in a three-dimensional irregular pore structure model is demonstrated.

[0047] Attachment Figure 6 The velocity field, particle field, and barium sulfate concentration field in the three-dimensional Y-shaped pore throat structure model grid are displayed.

[0048] Attachment Figure 7 The particle number density distribution at different positions in the three-dimensional Y-shaped pore throat structure model channel is shown.

[0049] Attachment Figure 8 Chemical substance concentration and streamlines in the grid of a three-dimensional irregular pore structure model.

[0050] Attachment Figure 9A is the pressure field corresponding to the three-dimensional irregular pore structure model.

[0051] Attachment Figure 9B is the particle field corresponding to the three-dimensional irregular pore structure model.

[0052] Figure 7 In the equation, population density refers to the particle number density. DETAILED DESCRIPTION

[0053] The present invention is not limited by the following embodiments, and specific implementation methods can be determined based on the technical solutions and actual conditions of the present invention. In the present invention, the crystals refer to barium sulfate crystals.

[0054] The supersaturation ratio refers to the actual solute concentration and solubility of a supersaturated solution.

[0055] The present invention will be further described below in conjunction with the embodiments:

[0056] Example 1: A high-precision dynamic acquisition method for barite precipitation in reservoir modeling includes:

[0057] The key reaction parameters of the fluid-rock system of the reservoir are loaded into the three-dimensional pore model. The three-dimensional pore model loaded with the key reaction parameters is used to carry out numerical simulation of the mineral dissolution and barite precipitation process. Based on the numerical simulation of the mineral dissolution and barite precipitation process, the amount of barite precipitation is calculated. The barite precipitation amount is calculated as follows:

[0058]

[0059] m=r×t

[0060] Where m represents the amount of barite precipitation, mol / m 3 ; t represents time, s; r represents the chemical reaction rate of barium sulfate crystal formation, mol / m 3 / s;B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s; G represents the growth rate of barium sulfate crystals, m / s; n represents the particle number density, 1 / m 4 ρ c Indicates the density of barium sulfate crystals, kg / m 3 ; L0 represents the minimum crystal size of barium sulfate, m; k V The shape factor, which represents the volume, is used to relate the geometric properties of a particle to its size; ΔL i Indicates the length of the barium sulfate crystal size interval, m; L i represents the barium sulfate crystal size, m.

[0061] The present invention utilizes COMSOL Software: Loads the key reaction parameters of the reservoir fluid-rock system into a three-dimensional pore model, and uses the three-dimensional pore model loaded with the key reaction parameters to carry out numerical simulations of mineral dissolution and barite precipitation processes.

[0062] Example 2: As an optimization of the above example, the nucleation rate B of the barium sulfate particle source is N Calculate as follows:

[0063]

[0064] Where B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s;D AB represents the apparent diffusion coefficient, m 2 / s, describes the diffusion ability of solute in solution; K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; N A represents Avogadro's constant, 1 / mol, which is used to relate the amount of a substance to the number of particles; γ CL represents the interfacial energy, mJ / m 2 , describing the interfacial energy between the crystal and the solution; k B represents the Boltzmann constant, J / K, which is used to describe the thermal motion of microscopic particles; V m Represents the molecular volume, m 3 , reflecting the size of the space occupied by the molecules; T represents the temperature, K, which affects the thermodynamics and kinetics of the crystallization process.

[0065] Example 3: As an optimization of the above example, the supersaturation is calculated as follows:

[0066]

[0067] Where S represents supersaturation; γ represents activity coefficient; c Ba2+ Indicates the mass concentration of barium ions, kg / m 3 ; Indicates the mass concentration of sulfate ions, kg / m 3 ;K SP represents the solubility product, mol 2 / m 6 .

[0068] Example 4: As an optimization of the above example, the barium sulfate crystal growth rate is calculated as follows:

[0069]

[0070] G represents the growth rate of barium sulfate crystals, m / s; ka and k V are the shape factors for area and volume, respectively, used to relate the geometric properties of the particle to its size; Sh is the Sherwood number, which describes the mass transfer efficiency of solutes on the particle surface; M W,C represents the molar mass of barium sulfate crystals, kg / mol, reflecting the chemical composition of the crystals; ρ C Indicates the density of barium sulfate crystals, kg / m 3 , describing the physical properties of crystals; D AB represents the apparent diffusion coefficient, m 2 / s, describes the diffusion ability of solute in solution; K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; L represents the size of barium sulfate crystal particles, m.

[0071] Example 5: As an optimization of the above example, the particle number density is calculated as follows:

[0072]

[0073] Where L represents the size of the barium sulfate crystal particles, m; N represents the total number of barium sulfate crystal particles per unit volume (i.e., the number of particles), 1 / m 3 ; The particle number density n can more accurately describe the distribution characteristics of the particle group.

[0074] Example 6: As an optimization of the above example, the minimum crystal size of barium sulfate is calculated as follows:

[0075]

[0076] Where L0 represents the minimum crystal size of barium sulfate, m; ν d represents the dissociation number, which describes the degree of dissociation of crystal molecules in solution; γ CL represents the interface energy, mJ / m 2 , describing the interfacial energy between the crystal and the solution; k B represents the Boltzmann constant, J / K, which is used to describe the thermal motion of microscopic particles; Vm represents the molecular volume, m 3 , reflecting the size of the space occupied by the molecules; T represents temperature, K, which affects the thermodynamics and kinetics of the crystallization process; S represents supersaturation.

[0077] Example 7: As an optimization of the above-mentioned Example 1, the key reaction parameters of the fluid-rock system of the reservoir include the physical parameters of barium sulfate, equilibrium constant parameters, fluid dynamics parameters, interface and diffusion parameters, temperature parameters (°C), and concentration parameters at the entrance of the three-dimensional pore model.

[0078] Example 8: As an optimization of the above-mentioned Example 7, the concentration parameters at the entrance of the three-dimensional pore model include the barium ion concentration (kmol / m 3 ), hydrosulfate ion concentration (kmol / m 3 ), sulfate ion concentration (kmol / m 3 ), chloride ion concentration (kmol / m 3 ), hydrogen ion concentration (kmol / m 3 ).

[0079] Example 9: As an optimization of Example 7, the physical parameters of barium sulfate include the molar mass (g / mol) of barium sulfate, the density (g / cm 3 ), the molar mass of sulfate ion (g / mol), the molar mass of barium ion (g / mol).

[0080] Example 10: As an optimization of Example 7 above, the equilibrium constant parameters include the equilibrium constant k1 (kmol / m 3 ), equilibrium constant k2(m 3 / kmol), dissolved products (mol 2 / l 2 ).

[0081] Example 11: As an optimization of Example 7, the fluid dynamics parameters include the mixing tube diameter D (mm) of the three-dimensional pore model, the channel length Lc (mm) of the three-dimensional pore model, the Reynolds number Re, the solvent viscosity (Pa·s), the solvent density (kg / m 3 ), solvent molar mass (kg / mol), average outlet velocity of the three-dimensional pore model (m / s), average inlet velocity of the three-dimensional pore model (m / s), fluid kinematic viscosity (m 2 / s), Schmidt number, material diffusion coefficient (m 2 / s), particle diffusion coefficient (m 2 / s).

[0082] Example 12: As an optimization of Example 7 above, the interface and diffusion parameters include interface energy (mJ / m 2 ), Sherwood number, apparent diffusion coefficient (m 2 / s).

[0083] Example 13: As an optimization of the above-mentioned Example 1, the three-dimensional pore model is a three-dimensional Y-type pore throat structure model or a three-dimensional irregular pore structure model.

[0084] The present invention utilizes COMSOL Software: The key reaction parameters of the fluid-rock system of the reservoir are loaded into a three-dimensional pore model, and the three-dimensional pore model loaded with the key reaction parameters is used to carry out numerical simulation of the mineral dissolution and barite precipitation process.

[0085] During the numerical simulation, the three-dimensional pore model assumed that barium ions and sulfate ions were completely mixed, focusing on the rapid evaluation of the kinetic characteristics of barite precipitation.

[0086] 3D Y-shaped pore-throat structure model: This model meticulously simulates the flow and sedimentation mechanisms of fluids in microscopic pore throats. Its structure consists of a main mixing tube (10 mm long, 0.1 mm diameter) and two symmetrical Y-shaped pore-throat inlets (0.1 mm diameter). The inlet flow rate is precisely set to a Reynolds number of 250 to ensure that the simulated flow state is highly consistent with actual reservoir conditions.

[0087] Three-dimensional irregular pore structure model: reconstructed based on the irregular pore network of the actual reservoir. The overall size of the model is 10 mm long, 5 mm wide, and 0.2 mm thick. It can finely depict the transmission and sedimentation process of fluids in complex pore structures, providing a high-precision simulation environment for studying the dynamic changes of reservoir physical properties.

[0088] Reaction mechanism and kinetic characteristics

[0089] During the simulation, the models all assumed that barium chloride dissociated completely, and the first step of sulfuric acid dissociation was also considered to have occurred completely. The second step of sulfuric acid dissociation and the formation of water-soluble barium sulfate follow the following reaction formula:

[0090]

[0091] Dissociation process of sulfuric acid:

[0092] The first step of dissociation of sulfuric acid (H2SO4→H + +HSO4 - ) is considered to occur completely in aqueous solution, so we mainly focus on its second dissociation step (HSO4 - →H + +SO4 2- ). This process is controlled by the balance between the concentration of hydrogen ions and the concentration of sulfate ions in the solution.

[0093] The reaction formula for the formation of water-soluble barium sulfate is:

[0094]

[0095] In solution, barium ions (Ba 2+ ) and sulfate ions (SO4 2-) combine to form water-soluble barium sulfate (BaSO4aq). This process is a key step before the formation of precipitation and determines the concentration change of barium ions and sulfate ions in the solution.

[0096] The chemical reaction formula of barium sulfate precipitation is:

[0097]

[0098] When the concentration of barium and sulfate ions in solution exceeds their solubility product (Ksp), barium sulfate begins to precipitate as a solid (BaSO4s). This process is one of the main causes of reservoir damage because the precipitates block pores and reduce the permeability of the reservoir.

[0099] In the model, we configured the aforementioned reactions using the software's Reaction Engineering and Chemistry interfaces to ensure accurate capture of the formation and precipitation of barium sulfate during simulation. These reaction configurations provided a solid theoretical foundation for subsequent numerical simulations, enabling detailed investigation of the impact of barite precipitation on reservoir physical properties.

[0100] Reaction parameters:

[0101] Supersaturation is the driving force behind crystallization or precipitation. When the concentration of a solute in a solution exceeds its solubility, supersaturation occurs, prompting the solute to precipitate from the solution, forming crystals or precipitates. The degree of supersaturation directly influences the crystallization rate and crystal growth behavior. Therefore, supersaturation is an essential parameter for studying and controlling crystallization processes.

[0102] Supersaturation is calculated as follows:

[0103]

[0104] Where S represents supersaturation; γ represents activity coefficient; c Ba2+ Indicates the mass concentration of barium ions, kg / m 3 ; Indicates the mass concentration of sulfate ions, kg / m 3 ;K SP represents the solubility product, mol 2 / m 6 .

[0105] In chemical reactions, supersaturation can influence reaction rates and reaction pathways. For example, in precipitation reactions, supersaturation determines whether the concentrations of reactants are high enough to overcome the activation energy barrier and initiate the reaction. By controlling supersaturation, reaction conditions can be optimized and reaction efficiency can be improved.

[0106] The extended Debye-Hückel model, proposed by Bromley, is a theoretical model for calculating activity coefficients in solutions. It is an improvement on the classic Debye-Hückel model, aiming to more accurately describe the interactions between ions in ionic solutions, particularly those at high ionic strength or in complex electrolyte solutions.

[0107] The classic Debye-Hückel theory is a classic theory describing the activity coefficients of ions in dilute electrolyte solutions. It is based on the following assumptions:

[0108] The solution is dilute and the interactions between ions are mainly dominated by electrostatic interactions.

[0109] Ions are point charges, so their size and shape are ignored.

[0110] The concentration of ions in the solution is very low and the interaction between ions can be approximated as linear.

[0111] This study adopts the extended Debye-Huckel model proposed by Bromley to consider the activity coefficient.

[0112] The activity coefficient is calculated as follows:

[0113]

[0114] Where γ represents the activity coefficient; A is taken as 0.511; Z represents the charge of the ion; F represents the correction term for the activity coefficients of barium ions and sulfate ions in the model, which is used to quantify their interactions with other ions in the solution; I m Indicates the ionic strength of the solution.

[0115] The activity coefficient γ is an important parameter that describes the actual behavior of ions in solution. It reflects the difference between the effective concentration (activity) of ions in solution and their actual concentration. The magnitude of the activity coefficient is affected by the following factors:

[0116] Charge of the ion: the higher the charge, the lower the activity coefficient.

[0117] Concentration of ions: The higher the concentration, the lower the activity coefficient.

[0118] Ionic strength of the solution: The higher the ionic strength, the lower the activity coefficient.

[0119] The population equilibrium equation is:

[0120] During barite crystallization, the dynamics of the crystal population are influenced by a variety of factors, including nucleation rate, growth rate, and interparticle interactions. To quantitatively describe these changes, we use the Population Balance Equation (PBE) to model the temporal evolution of particle number density. The PBE comprehensively accounts for processes such as particle formation, growth, aggregation, and fragmentation, providing a theoretical basis for studying changes in particle distribution during barite crystallization.

[0121] In order to simulate the distribution of barite crystal particles as the particle size changes, the particle number density n (unit: 1 / m 4 ). The particle number density n is defined as the number of particles within a specific particle size range per unit volume, and its mathematical expression is as follows:

[0122]

[0123] Where L is the particle size (in meters), N is the total number of particles per unit volume (in 1 / m 3 ). The distribution characteristics of the particle group can be described more accurately through the particle number density n.

[0124] The population equilibrium equation is as follows:

[0125]

[0126] G represents the crystal growth rate (unit: m / s), which describes the growth rate of crystal particles over time; ν represents the kinematic viscosity (unit: m 2 / s), reflecting the viscosity of the fluid; S c represents the Schmidt number, which is the ratio of momentum diffusion to mass diffusion in the fluid; B0 represents the nucleation rate (unit: 1 / m 4 / s), describes the rate of generation of new crystal particles per unit volume; L c0 represents the minimum stable crystal size (in m), which is the minimum particle size at which a crystal can exist stably; ▽(un) represents the convection top; u represents the flow rate, in m / s.

[0127] The left side of the population equilibrium equation is expressed from left to right:

[0128] Changes in population density over time;

[0129] Convective crystal transport effects;

[0130] Effect of crystal growth on particle size distribution.

[0131] The right side of the population equilibrium equation is expressed from left to right:

[0132] Diffusion crystal transport according to Fick's second law;

[0133] The generation rate of the crystal nucleation process.

[0134] Here, the nucleation rate B0 as the number density source and the nucleation rate B as the particle number source N (Unit is 1 / m 3 / s) are related, that is: B0=B N / L,B N According to the classical nucleation theory, it is defined as follows:

[0135]

[0136] Where B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s;D AB represents the apparent diffusion coefficient, m 2 / s, describes the diffusion ability of solute in solution; K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; N A represents Avogadro's constant, 1 / mol, which is used to relate the amount of a substance to the number of particles; γ CL represents the interface energy, mJ / m 2 , describing the interfacial energy between the crystal and the solution; k B represents the Boltzmann constant, J / K, which is used to describe the thermal motion of microscopic particles; V m Represents the molecular volume, m 3 , reflecting the size of the space occupied by the molecules; T represents the temperature, K, which affects the thermodynamics and kinetics of the crystallization process.

[0137] Among them, B N Classical nucleation theory describes the rate at which new crystal particles form per unit volume. This theory considers the thermodynamic driving force and energy barriers in the crystal nucleation process, providing a theoretical basis for calculating the nucleation rate.

[0138] Nucleation process:

[0139] The nucleation process determines the initial size and number of crystals, while the growth process describes how the crystals grow over time. Nucleation occurs at the minimum barium sulfate crystal size L0, which is defined as:

[0140]

[0141] Where L0 represents the minimum crystal size of barium sulfate, m; ν drepresents the dissociation number, which describes the degree of dissociation of crystal molecules in solution; γ CL represents the interfacial energy, mJ / m 2 , describing the interfacial energy between the crystal and the solution; k B represents the Boltzmann constant, J / K, which is used to describe the thermal motion of microscopic particles; V m Represents the molecular volume, m 3 , reflecting the size of the space occupied by the molecules; T represents temperature, K, which affects the thermodynamics and kinetics of the crystallization process; S represents supersaturation.

[0142] L0 is regarded as a constant, which means that only when the crystal size reaches L0 can the crystal exist stably and grow further.

[0143] Crystal growth process:

[0144] The growth rate of the crystal is controlled by solute transport, which is expressed as:

[0145]

[0146] G represents the growth rate of barium sulfate crystals, m / s; k a and k V are the shape factors for area and volume, respectively, used to relate the geometric properties of the particle to its size; Sh is the Sherwood number, which describes the mass transfer efficiency of solutes on the particle surface; M w,c represents the molar mass of barium sulfate crystals, kg / mol, reflecting the chemical composition of the crystals; ρ c Indicates the density of barium sulfate crystals, kg / m 3 , describing the physical properties of crystals; D AB represents the apparent diffusion coefficient, m 2 / s, describes the diffusion ability of solute in solution; K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; L represents the size of barium sulfate crystal particles, m.

[0147] The size range of barium sulfate crystals is divided into multiple intervals i = 1, 2, 3, ..., I, each interval corresponds to a specific crystal size range, from L i-1 to L i , the interval length is ΔL i .

[0148] The minimum crystal size of barium sulfate, L0, is defined as the minimum stable crystal size L C0 , which is the minimum size at which a crystal can exist stably.

[0149] In using COMSOL During numerical simulations, the population equilibrium equation for each barium sulfate crystal size interval is modeled in the software's zero-dimensional model using the Global Ordinary Differential and Differential Algebraic Equations interface, while in the three-dimensional pore model, it is modeled using a separate Steady Convection-Diffusion Equation interface.

[0150] A transient solver was used with a time step size of 0.001s to 0.06s to ensure that the rapid changes in the crystal growth and nucleation processes were captured.

[0151] Output results include supersaturation, crystal size distribution, and mass concentration curves. These results can describe the dynamic behavior of barite crystallization in detail and provide basic data for subsequent analysis.

[0152] In the three-dimensional pore model, a steady-state laminar flow solver is used, combined with parallel computing to optimize the grid resolution. The minimum cell size is set to 2e -5 m, the maximum element size is set to 5e -5 m, to ensure the best balance between computational efficiency and accuracy.

[0153] During the simulation process, changes in reservoir flow rate, pressure and crystal number are monitored in real time to ensure the consistency of simulation results with the actual process.

[0154] Through the dynamic feedback mechanism, the model response parameters can be adjusted in time, the simulation process can be optimized, and the accuracy of the prediction can be improved.

[0155] The simulation generates detailed spatiotemporal distribution maps, including changes in supersaturation, crystal number density, porosity, and permeability. These maps visually demonstrate the dynamic changes in the reservoir during barite crystallization, providing key data support for subsequent analysis.

[0156] Supersaturation analysis

[0157] Supersaturation is a key driving force for precipitation reactions. It not only reflects the extent to which the solute concentration in the solution exceeds its solubility, but also directly determines the amount of available reactants in the solution.

[0158] 1. Initial stage:

[0159] At the beginning of the reaction, there is a high degree of supersaturation in the solution, which is due to the concentration of barium sulfate in the solution exceeding its solubility.

[0160] High supersaturation provides sufficient driving force for nucleation and promotes the formation of a large number of tiny crystals.

[0161] 2. Nucleation stage:

[0162] Over time, the precipitation process of barium sulfate begins to consume solutes in the solution, resulting in a gradual decrease in supersaturation.

[0163] During this stage, the nucleation rate of crystals is high and the number of crystals in the solution increases rapidly.

[0164] 3. Growth stage:

[0165] When the supersaturation decreases to a certain level, the nucleation rate gradually decreases and crystal growth becomes the dominant process.

[0166] The crystals gradually increase in size by absorbing solutes from the solution, and the supersaturation continues to decrease.

[0167] 4. Balanced state:

[0168] Eventually, the supersaturation stabilizes when the ion concentration in the solution reaches equilibrium with the solubility of barium sulfate.

[0169] At this point, the crystal growth rate and dissolution rate reach a dynamic equilibrium, and the crystal size distribution tends to be stable.

[0170] Application of model results

[0171] Optimizing the crystallization process:

[0172] By monitoring and controlling supersaturation, the crystallization process can be optimized to ensure crystal growth under optimal conditions. For example, by adjusting the reaction temperature and solution concentration, the rate of change of supersaturation can be effectively controlled, thereby affecting the nucleation and growth process of the crystal.

[0173] Improve product quality:

[0174] Precise control of supersaturation helps produce crystals of uniform size and high purity. Higher supersaturation may result in a wider crystal size distribution, while moderate supersaturation helps form uniform crystals.

[0175] Reservoir Management:

[0176] In reservoir management, changes in supersaturation directly impact reservoir permeability and porosity. By simulating the spatiotemporal distribution of supersaturation, high-risk areas for reservoir damage, such as blockages and areas of decreased permeability, can be identified and corresponding optimization measures can be proposed.

[0177] Risk Assessment and Decision Support:

[0178] Based on the supersaturation simulation results, a risk assessment report can be generated to provide decision support for optimizing the injection-production plan. For example, by adjusting the solution pH and adding inhibitors, reservoir damage can be effectively mitigated and injection-production efficiency can be improved.

[0179] Example 13: A high-precision dynamic acquisition method for barite precipitation in reservoir modeling includes:

[0180] The key reaction parameters of the reservoir fluid-rock system (such as Table 1 to Table 5, the temperature is 25 ° C) are loaded into the three-dimensional Y-type pore throat structure model (see Figure 1 ) or three-dimensional irregular pore structure model (see Figure 2 ), using the three-dimensional pore model loaded with the key reaction parameters to carry out numerical simulation of mineral dissolution and barite precipitation process, based on the numerical simulation of mineral dissolution and barite precipitation process, calculate the amount of barite precipitation, the barite precipitation amount is calculated as follows:

[0181] Supersaturation is calculated as follows:

[0182]

[0183] Where S represents supersaturation; γ represents activity coefficient; c Ba2+ Indicates the mass concentration of barium ions, kg / m 3 ; Indicates the mass concentration of sulfate ions, kg / m 3 ;K SP represents the solubility product, mol 2 / m 6 ;

[0184] The minimum crystal size of barium sulfate is calculated as follows:

[0185]

[0186] Where L0 represents the minimum crystal size of barium sulfate, m; ν d represents the dissociation number, which describes the degree of dissociation of crystal molecules in solution; γ CL represents the interface energy, mJ / m 2 , describing the interfacial energy between the crystal and the solution; k B represents the Boltzmann constant, J / K, which is used to describe the thermal motion of microscopic particles; V m Represents the molecular volume, m 3 , reflecting the size of the space occupied by the molecules; T represents temperature, K, which affects the thermodynamics and kinetics of the crystallization process; S represents supersaturation;

[0187] The barium sulfate crystal growth rate is calculated as follows:

[0188]

[0189] G represents the growth rate of barium sulfate crystals, m / s; k a and k V are the shape factors for area and volume, respectively, used to relate the geometric properties of the particle to its size; Sh is the Sherwood number, which describes the mass transfer efficiency of solutes on the particle surface; Mw,c represents the molar mass of barium sulfate crystals, kg / mol, reflecting the chemical composition of the crystals; ρ c Indicates the density of barium sulfate crystals, kg / m 3 , describing the physical properties of crystals; D AB represents the apparent diffusion coefficient, m 2 / s, describes the diffusion ability of solute in solution; K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; L represents the size of barium sulfate crystal particles, m;

[0190] The particle number density is calculated as follows:

[0191]

[0192] Where L represents the size of the barium sulfate crystal particles, m; N represents the total number of barium sulfate crystal particles per unit volume, 1 / m3; the particle number density n can more accurately describe the distribution characteristics of the particle group;

[0193] Nucleation rate B of the barium sulfate particle source N Calculate as follows:

[0194]

[0195] Where B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s;D AB represents the apparent diffusion coefficient, m 2 / s, describes the diffusion ability of solute in solution; K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; N A represents Avogadro's constant, 1 / mol, which is used to relate the amount of a substance to the number of particles; γ CL represents the interfacial energy, mJ / m 2 , describing the interfacial energy between the crystal and the solution; k B represents the Boltzmann constant, J / K, which is used to describe the thermal motion of microscopic particles; V m Represents the molecular volume, m 3 , reflecting the size of the space occupied by the molecules; T represents the temperature, K, which affects the thermodynamics and kinetics of the crystallization process;

[0196]

[0197] m=r×t

[0198] Where m represents the amount of barite precipitation, mol / m 3 ; t represents time, s; r represents the chemical reaction rate of barium sulfate crystal formation, mol / m 3 / s;B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s; G represents the growth rate of barium sulfate crystals, m / s; n represents the particle number density, 1 / m 4 ρ C Indicates the density of barium sulfate crystals, kg / m 3 ; L0 represents the minimum crystal size of barium sulfate, m; k V The shape factor, which represents the volume, is used to relate the geometric properties of a particle to its size; ΔL i Indicates the length of the barium sulfate crystal size interval, m; L i represents the barium sulfate crystal size, m.

[0199] In this embodiment 13, the variation of supersaturation with time is shown in FIG. Figure 3 .Depend on Figure 3 As can be seen, the initial supersaturation is consumed as barium sulfate precipitates until an equilibrium state is reached where the ion concentration and solubility are equal. This reflects the dynamic behavior of the barium sulfate crystallization process. Initially, the supersaturation ratio is high. As barium sulfate precipitates, the supersaturation ratio gradually decreases and eventually stabilizes, indicating that the ion concentration and solubility in the solution have reached equilibrium. This process reveals the changing driving forces of crystal nucleation and growth.

[0200] Figure 4 The time-varying mass concentrations of the reactants (except the crystal products) and the crystal products are shown, reflecting the transformation process of the solute reactants into crystal particles.

[0201] Figure 4 The trend of sulfate (blue), barium (red), water-soluble barium sulfate (green), crystalline barium sulfate (cyan) and total concentration (purple) is shown in the figure. Figure 4 As can be seen, the concentrations of sulfate and barium gradually decrease over time, indicating that they are consumed in the formation of crystals. The concentration of water-soluble barium sulfate first increases and then decreases, indicating that it participates in the reaction as an intermediate product. The concentration of crystalline barium sulfate continues to rise, indicating that crystals gradually form and accumulate. The total concentration remains constant, verifying the conservation of mass in the system.

[0202] Figure 5A 、 Figure 5B The effect of saturation change on crystal size distribution in three-dimensional Y-type pore structure model and three-dimensional irregular pore structure model is demonstrated.

[0203] Figure 5A 、 Figure 5BThe temporal evolution of crystal size is demonstrated using a 3D Y-shaped pore structure model and a 3D irregular pore structure model. The particle density distribution is normalized by the total number of crystal particles per unit volume at each time point, allowing for comparison of crystal size distributions at different time points. Initially, the crystals are small and have a narrow distribution. As the supersaturation decreases, the crystals grow, and the size distribution broadens. Ultimately, the crystal size stabilizes, with a well-defined distribution range.

[0204] Figure 6 The velocity field, particle field, and barium sulfate concentration field in the three-dimensional Y-shaped pore throat structure model grid are displayed.

[0205] Depend on Figure 6 As can be seen, the mixed flow becomes increasingly uniform as the fluid flows along the mixing tube. The velocity field shows the flow direction and velocity distribution of the fluid, with colors ranging from light blue to dark red indicating a transition from low to high velocity. The particle field shows the distribution of particle number density, with colors ranging from light blue to dark red indicating a transition from low to high particle number density. The barium sulfate concentration field shows the concentration distribution of barium sulfate, with colors ranging from light blue to dark red indicating a transition from low to high concentration. Supersaturation initially occurs at the top of the mixing tube and gradually decreases as the fluid flows along the mixing tube, indicating that nucleation and crystal growth are ongoing.

[0206] Figure 7 The particle number density distribution at different positions in the three-dimensional Y-shaped pore throat structure model channel is shown.

[0207] Figure 7 The horizontal axis represents crystal size, and the vertical axis represents particle number density. Curves of different colors represent distributions at different locations. The crystal size distribution varies significantly with location, indicating that the meshing significantly influences crystal growth and distribution. Crystals are smaller and more narrowly distributed near the inlet. As the fluid flows, the crystals grow larger and their distribution becomes wider.

[0208] Figure 8 Chemical substance concentration and streamlines in the grid of a three-dimensional irregular pore structure model.

[0209] Figure 8 The concentration distribution and streamlines of four ions (barium, barium sulfate, hydrogen, and hydrosulfate) in the reservoir are displayed. Colors indicate concentrations, and streamlines depict the fluid flow path. The results show that the ion concentration distribution is uneven, with areas of high and low concentrations. Areas with dense streamlines indicate faster flow rates. The distribution of barium and sulfate ions may affect the formation of barium sulfate precipitation, while the concentration distribution of hydrogen ions may affect the pH of the solution.

[0210] Figure 9A 、 Figure 9BThey are the pressure field and particle field corresponding to the three-dimensional irregular pore structure model, which show the pressure field and particle field distribution in the reservoir.

[0211] Figure 9A 、 Figure 9B In the figure, the pressure field diagram shows the distribution of pressure within the reservoir, with colors ranging from red to purple indicating a decrease in pressure. The pressure is higher at the left entrance of the 3D irregular pore structure and gradually decreases as the fluid flows into the reservoir. The particle field diagram shows the distribution of particle number density within the reservoir, with colors ranging from red to light blue indicating a decrease in particle number density. The particle number density is higher at the left entrance of the 3D irregular pore structure and gradually decreases as the fluid flows into the reservoir. The distributions of the pressure field and particle field are somewhat correlated. The high pressure and high particle number density at the left entrance of the 3D irregular pore structure may be related to fluid injection and the nucleation and growth processes of crystals.

[0212] From the above, we can see that Figures 3 to 9B The processes of crystal nucleation, growth, and equilibrium are demonstrated from the perspectives of ion concentration and particle number density.

[0213] The above technical features respectively constitute the embodiments of the present invention, which have strong adaptability and implementation effect. Non-essential technical features can be added or removed according to actual needs to meet the requirements of different situations.

[0214] Table 1 Concentration parameters

[0215] illustrate Numerical unit Barium ion concentration at the inlet 0.5 <![CDATA[kmol / m 3 ]]> Inlet hydrosulfate ion concentration 0.31879 <![CDATA[kmol / m 3 ]]> Inlet chloride ion concentration 2× inlet barium ion concentration <![CDATA[kmol / m 3 ]]> Inlet hydrogen ion concentration 0.34121 <![CDATA[kmol / m 3 ]]> Inlet sulfate ion concentration 0.01121 <![CDATA[kmol / m 3 ]]>

[0216] Table 2 Physical property parameters

[0217]

[0218]

[0219] Table 3 Equilibrium constant parameters

[0220] illustrate Numerical unit Equilibrium constant k1 1.20E-02 <![CDATA[kmol / m 3 ]]> Equilibrium constant k2 185.19 <![CDATA[m 3 / kmol]]> Dissolution product Ks 9.82E-11 <![CDATA[mol 2 / l 2 ]]>

[0221] 4 Equilibrium constant parameters

[0222] illustrate Numerical unit Mixing tube diameter D 1 mm Channel length 10 mm Reynolds number 250 - Solvent viscosity 8.87E-04 Pa·s Solvent density 996.9 <![CDATA[kg / m 3 ]]> Solvent molar mass 0.018 kg / mol Average exit speed Re × solvent viscosity / (solvent density × D) m / s Average inlet velocity 2× average exit velocity m / s Fluid kinematic viscosity Solvent viscosity / solvent density <![CDATA[m 2 / s]]> Schmidt number 1 - Diffusion coefficient of substances Fluid kinematic viscosity / Schmidt number <![CDATA[m 2 / s]]> Particle diffusion coefficient Fluid kinematic viscosity / Schmidt number <![CDATA[m 2 / s]]>

[0223] Table 5 Interface and diffusion parameters

[0224] illustrate Numerical unit Interface energy 118.1 <![CDATA[mJ / m 2 ]]> Sherwood number 2 - Apparent diffusion coefficient 9.46E-10 <![CDATA[m 2 / s]]>

Claims

1. A high-precision dynamic acquisition method for barite precipitation in reservoir modeling, characterized in that: include: The key reaction parameters of the fluid-rock system of the reservoir are loaded into the three-dimensional pore model. The three-dimensional pore model loaded with the key reaction parameters is used to carry out numerical simulation of the mineral dissolution and barite precipitation process. Based on the numerical simulation of the mineral dissolution and barite precipitation process, the amount of barite precipitation is calculated. The barite precipitation amount is calculated as follows: Where m represents the amount of barite precipitation, mol / m 3 ; t represents time, s; r represents the chemical reaction rate of barium sulfate crystal formation, mol / m 3 / s;B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s; G represents the growth rate of barium sulfate crystals, m / s; n represents the particle number density, 1 / m 4 ρ c Indicates the density of barium sulfate crystals, kg / m 3 ; L0 represents the minimum crystal size of barium sulfate, m; k V Shape factor representing volume; ΔL i Indicates the length of the barium sulfate crystal size interval, m; L i represents the barium sulfate crystal size, m.

2. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 1, characterized in that Nucleation rate B of the barium sulfate particle source N Calculate as follows: Where B N represents the nucleation rate of the barium sulfate particle source, 1 / m 3 / s;D AB represents the apparent diffusion coefficient, m 2 / s;K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; N A represents Avogadro's constant, 1 / mol; γ CL represents the interfacial energy, mJ / m 2 ;k B represents the Boltzmann constant, J / K; V m Represents the molecular volume, m 3 ; T represents temperature, K.

3. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 2, characterized in that, Supersaturation is calculated as follows: Where S represents supersaturation; γ represents activity coefficient; Indicates the mass concentration of barium ions, kg / m 3 ; Indicates the mass concentration of sulfate ions, kg / m 3 ;K SP represents the solubility product, mol 2 / m 6 .

4. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 1, 2 or 3, characterized in that, The barium sulfate crystal growth rate is calculated as follows: G represents the growth rate of barium sulfate crystals, m / s; k a and k V Represent the shape factors of area and volume respectively; Sh represents Sherwood number; M w,c represents the molar mass of barium sulfate crystals, kg / mol; ρ c Indicates the density of barium sulfate crystals, kg / m 3 ;D AB represents the apparent diffusion coefficient, m 2 / s;K SP represents the solubility product, mol 2 / m 6 ; S represents supersaturation; L represents the size of barium sulfate crystal particles, m.

5. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 1, 2 or 3, characterized in that, The particle number density is calculated as follows: Where L represents the size of the barium sulfate crystal particles, m; N represents the total number of barium sulfate crystal particles per unit volume, 1 / m 3 ; n represents the particle number density, 1 / m 4 .

6. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 4 is characterized in that, The minimum crystal size of barium sulfate is calculated as follows: Where L0 represents the minimum crystal size of barium sulfate, m; ν d represents the dissociation number; γ CL represents the interfacial energy, mJ / m 2 ;k B represents the Boltzmann constant, J / K; V m Represents the molecular volume, m 3 ; T represents temperature, K; S represents supersaturation.

7. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 1, 2 or 3, characterized in that, The key reaction parameters of the fluid-rock system of the reservoir include the physical parameters of barium sulfate, equilibrium constant parameters, fluid dynamics parameters, interface and diffusion parameters, temperature parameters, and concentration parameters at the inlet of the three-dimensional pore model.

8. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 7 is characterized in that, The concentration parameters at the entrance of the three-dimensional pore model include the barium ion concentration, the hydrosulfate ion concentration, the sulfate ion concentration, the chloride ion concentration, and the hydrogen ion concentration at the entrance of the three-dimensional pore model; or / and, The physical parameters of the barium sulfate include the molar mass of barium sulfate, the density of barium sulfate, the molar mass of sulfate ions, and the molar mass of barium ions; or / and, The equilibrium constant parameters include equilibrium constant k1, equilibrium constant k2, dissolved products; or / and, The fluid dynamics parameters include the mixing tube diameter of the three-dimensional pore model, the channel length of the three-dimensional pore model, the Reynolds number, the solvent viscosity, the solvent density, the solvent molar mass, the average outlet velocity of the three-dimensional pore model, the average inlet velocity of the three-dimensional pore model, the fluid kinematic viscosity, the Schmidt number, the material diffusion coefficient, and the particle diffusion coefficient; or / and, The interface and diffusion parameters include interface energy, Sherwood number, and apparent diffusion coefficient.

9. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 1, 2 or 3, characterized in that, The three-dimensional pore model is a three-dimensional Y-shaped pore throat structure model or a three-dimensional irregular pore structure model.

10. The high-precision dynamic acquisition method of barite precipitation in reservoir modeling according to claim 8, characterized in that: The three-dimensional pore model is a three-dimensional Y-shaped pore throat structure model or a three-dimensional irregular pore structure model.

Images