Pollution monitoring method and device based on stable isotope, computer equipment and storage medium

By constructing a pollution monitoring method using stable isotopes, obtaining volatilization, dissolution, and adsorption parameters, and training and evaluating models, the problems of low monitoring accuracy and spatial resolution of high-density non-aqueous liquid pollutants were solved, and efficient monitoring of pollutant distribution and flux was achieved.

CN121978271APending Publication Date: 2026-05-05ZJU HANGZHOU GLOBAL SCI & TECH INNOVATION CENT
View PDF 0 Cites 0 Cited by

Patent Information

Authority / Receiving Office
CN · China
Patent Type
Applications(China)
Current Assignee / Owner
ZJU HANGZHOU GLOBAL SCI & TECH INNOVATION CENT
Filing Date
2025-12-26
Publication Date
2026-05-05

AI Technical Summary

Technical Problem

Existing technologies for investigating high-density non-aqueous liquid DNAPLs contamination suffer from high costs, poor spatial coverage, and low accuracy and spatial resolution in in-situ characterization and fine delineation of DNAPLs multiphase states.

Method used

By obtaining the volatile fractionation, dissolution fractionation, and adsorption fractionation parameters of the stable isotopes of the target pollutant, an initial evaluation model is constructed. The target evaluation model is then trained using sample fractionation data to obtain local fractionation data and determine the fractionation data and mass flux of the pollutant in different phases.

Benefits of technology

It enables efficient monitoring of high-density non-aqueous liquid pollutants, accurately characterizes their distribution and mass flux in different phases, and solves the problem of accurate monitoring of multiphase transformation and migration processes.

✦ Generated by Eureka AI based on patent content.

Smart Images

  • Figure CN121978271A_ABST
    Figure CN121978271A_ABST
Patent Text Reader

Abstract

The invention relates to a pollution monitoring method and device based on stable isotopes, computer equipment and a storage medium. The method comprises the following steps: acquiring a volatilization fractionation parameter, a dissolution fractionation parameter and an adsorption fractionation parameter of a target pollutant and sample fractionation data corresponding to the target pollutant; constructing an initial evaluation model according to the volatilization fractionation parameters, the dissolution fractionation parameters and the adsorption fractionation parameters; training the initial evaluation model according to the sample fractionation data, and determining a target evaluation model; acquiring local fractionation data of the target pollutant; evaluating local fractionation data according to the target evaluation model, and determining fractionation data of a second target phase state of a target pollutant in a target environment; and determining the mass flux of the target pollutant in the target environment according to the local fractionation data and the fractionation data of the second target phase state, thereby efficiently monitoring the high-density non-aqueous phase liquid pollutant. The method can be adopted.
Need to check novelty before this filing date? Find Prior Art

Description

Technical Field

[0001] This application relates to the field of pollution monitoring technology, and in particular to a pollution monitoring method, apparatus, computer equipment and storage medium based on stable isotopes. Background Technology

[0002] Currently, the investigation of contamination of high-density non-aqueous liquid DNAPLs in a site mainly relies on traditional invasive drilling sampling and experimental analysis, which has disadvantages such as high cost and poor spatial coverage. At the same time, due to the heterogeneity of the site medium and the complexity of the multiphase and multi-interface migration and transformation processes of DNAPLs, the in-situ characterization and fine characterization techniques of DNAPLs multiphase have problems with low accuracy and spatial resolution.

[0003] There is an urgent need for a solution that can efficiently monitor high-density non-aqueous liquid pollutants. Summary of the Invention

[0004] Therefore, it is necessary to provide a stable isotope-based pollution monitoring method, device, computer equipment, and storage medium that can efficiently monitor high-density non-aqueous liquid pollutants, addressing the aforementioned technical problems.

[0005] In a first aspect, this application provides a pollution monitoring method based on stable isotopes. The method includes: acquiring volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters, and sample fractionation data corresponding to the stable isotopes of a target pollutant; the target pollutant is a high-density non-aqueous liquid pollutant; the sample fractionation data includes the molar quantity and stable isotope abundance of the target pollutant in the NAPL phase, aqueous phase, gas phase, and adsorption phase under stable conditions in any environment; constructing an initial evaluation model based on the volatile fractionation parameters, dissolution fractionation parameters, and adsorption fractionation parameters; and training the initial evaluation model based on the sample fractionation data. The process involves: determining a target assessment model; acquiring partial fractionation data of the target pollutant; the partial fractionation data being the fractionation data of the target pollutant in a first target phase in the target environment; the first target phase being any one of the NAPL phase, aqueous phase, gas phase, or adsorbed phase; evaluating the partial fractionation data according to the target assessment model to determine the fractionation data of the target pollutant in a second target phase in the target environment; the second target phase being any phase other than the first target phase; and determining the mass flux of the target pollutant in the target environment based on the partial fractionation data and the fractionation data of the second target phase.

[0006] In one embodiment, obtaining the volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters, and sample fractionation data corresponding to the stable isotopes of the target pollutant includes: obtaining sample fractionation data corresponding to the target pollutant; conducting a volatilization experiment to determine the volatile fractionation parameters of the stable isotopes of the target pollutant; conducting a dissolution experiment to determine the dissolution fractionation parameters of the stable isotopes of the target pollutant; and conducting an adsorption experiment to determine the adsorption fractionation parameters of the stable isotopes of the target pollutant.

[0007] In one embodiment, the step of conducting a volatilization experiment on the target pollutant to determine the volatilization fractionation parameters of the stable isotopes of the target pollutant includes: conducting a volatilization experiment based on a preset volatilization environment to determine the volatile isotope ratio and volatilization ratio of the target pollutant in the NAPL phase after volatilization; determining a first volatile isotope abundance based on the volatile isotope ratio and the corresponding natural isotope ratio of the target pollutant; determining a second volatile isotope abundance based on the preset isotope ratio and the corresponding natural isotope ratio of the target pollutant; the preset isotope ratio includes the isotope ratio of the target pollutant in the NAPL phase before volatilization; determining a volatilization fractionation factor of the target pollutant based on the first volatile isotope abundance, the second volatile isotope abundance, and the volatilization ratio, and using the volatilization fractionation factor as a volatilization fractionation parameter.

[0008] In one embodiment, the step of conducting a dissolution experiment on the target pollutant to determine the dissolution fractionation parameters of the stable isotopes of the target pollutant includes: conducting a dissolution experiment based on a preset dissolution environment to determine the first dissolved isotope ratio of the target pollutant in the NAPL phase and the second dissolved isotope ratio of the target pollutant in the aqueous phase after dissolution; determining the first dissolved isotope abundance based on the first dissolved isotope ratio and the corresponding natural isotope ratio of the target pollutant; determining the second dissolved isotope abundance based on the second dissolved isotope ratio and the corresponding natural isotope ratio of the target pollutant; determining the dissolution fractionation factor of the target pollutant based on the first dissolved isotope abundance and the second dissolved isotope abundance, and using the dissolution fractionation factor as the dissolution fractionation parameter.

[0009] In one embodiment, the step of conducting an adsorption experiment on the target pollutant to determine the adsorption fractionation parameters of the stable isotopes of the target pollutant includes: conducting an adsorption experiment based on a preset adsorption environment to determine the first adsorbed isotope ratio and adsorption ratio of the target pollutant in the target pollutant solution after adsorption; determining the first adsorbed isotope abundance based on the first adsorbed isotope ratio and the corresponding natural isotope ratio of the target pollutant; determining the second adsorbed isotope abundance based on the first adsorbed isotope abundance, the adsorption ratio, and the preset isotope abundance; the preset isotope abundance includes the isotope abundance of the target pollutant in the target pollutant solution before adsorption; and determining the adsorption fractionation factor of the target pollutant based on the first adsorbed isotope abundance and the second adsorbed isotope abundance, and using the adsorption fractionation factor as the adsorption fractionation parameter.

[0010] In one embodiment, training the initial evaluation model based on the sample fractionation data to determine the target evaluation model includes: fitting the first constant, the second constant, and the third constant in the initial evaluation model based on the sample fractionation data to determine the target evaluation model; the first constant is the dissolution rate constant corresponding to the target pollutant, the second constant is the volatilization rate constant corresponding to the target pollutant, and the third constant is the adsorption rate constant corresponding to the target pollutant.

[0011] In one embodiment, the fractionation data of the target pollutant in the second target phase of the target environment includes the molar quantity and isotopic abundance of the target pollutant in the second target phase of the target environment. Determining the mass flux of the target pollutant in the target environment based on the partial fractionation data and the fractionation data of the second target phase includes: determining the mass flux of the target pollutant in the target environment during the dissolution process based on a first preset flux calculation formula, according to the first constant and the molar quantity of the target pollutant in the NAPL phase of the target environment; and determining the mass flux of the target pollutant in the target environment during the volatilization process based on a second preset flux calculation formula, according to the second constant and the molar quantity of the target pollutant in the NAPL phase of the target environment.

[0012] Secondly, this application also provides a pollution monitoring device based on stable isotopes. The device includes:

[0013] The first acquisition module is used to acquire the volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters, and sample fractionation data of the stable isotopes of the target pollutant; the target pollutant is a high-density non-aqueous liquid pollutant; the sample fractionation data is the molar quantity and stable isotope abundance of the NAPL phase, aqueous phase, gas phase, and adsorption phase of the target pollutant under any environment and stable conditions.

[0014] A construction module is used to construct an initial evaluation model based on the volatile fractionation parameters, dissolution fractionation parameters, and adsorption fractionation parameters;

[0015] The first determining module is used to train the initial evaluation model based on the sample fractionation data to determine the target evaluation model;

[0016] The second acquisition module is used to acquire partial fractionation data of the target pollutant; the partial fractionation data is the fractionation data of the target pollutant in the first target phase in the target environment; the first target phase is any one of the NAPL phase, aqueous phase, gas phase or adsorbed phase;

[0017] The second determining module is used to evaluate the partial fractionation data according to the target evaluation model and determine the fractionation data of the second target phase of the target pollutant in the target environment; the second target phase is the other phases besides the first target phase.

[0018] The third determining module is used to determine the mass flux of the target pollutant in the target environment based on the partial fractionation data and the fractionation data of the second target phase.

[0019] Thirdly, this application also provides a computer device, the computer device including a memory and a processor, the memory storing a computer program, and the processor executing the computer program to implement any of the methods in the first aspect above.

[0020] Fourthly, this application also provides a computer-readable storage medium having a computer program stored thereon, which, when executed by a processor, implements any of the methods in the first aspect described above.

[0021] Fifthly, this application also provides a computer program product, which includes a computer program that, when executed by a processor, implements any of the methods described in the first aspect above.

[0022] The aforementioned pollution monitoring method, apparatus, computer equipment, and storage medium based on stable isotopes acquire the volatile fractionation parameters, dissolved fractionation parameters, adsorption fractionation parameters, and corresponding sample fractionation data of the target pollutant; construct an initial evaluation model based on the volatile fractionation parameters, dissolved fractionation parameters, and adsorption fractionation parameters; train the initial evaluation model based on the sample fractionation data to determine the target evaluation model; acquire partial fractionation data of the target pollutant; the partial fractionation data is the fractionation data of the target pollutant in the first target phase in the target environment; evaluate the partial fractionation data based on the target evaluation model to determine the fractionation data of the target pollutant in the second target phase in the target environment; and determine the mass flux of the target pollutant in the target environment based on the partial fractionation data and the fractionation data of the second target phase, thereby achieving efficient monitoring of high-density non-aqueous liquid pollutants. Attached Figure Description

[0023] Figure 1 This is an application environment diagram of a pollution monitoring method based on stable isotopes in one embodiment;

[0024] Figure 2 This is a flowchart illustrating a pollution monitoring method based on stable isotopes in one embodiment;

[0025] Figure 3 This is a schematic diagram of the process for obtaining the volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters, and sample fractionation data corresponding to the target pollutant in one embodiment.

[0026] Figure 4 This is a structural block diagram of a pollution monitoring device based on stable isotopes in one embodiment;

[0027] Figure 5 This is an internal structural diagram of a computer device in one embodiment. Detailed Implementation

[0028] To make the objectives, technical solutions, and advantages of this application clearer, the following detailed description is provided in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the scope of this application.

[0029] The pollution monitoring method based on stable isotopes provided in this application can be applied to, for example... Figure 1In the application environment shown, terminal 102 communicates with server 104 via a network. A data storage system can store the data that server 104 needs to process. The data storage system can be integrated onto server 104 or located in the cloud or on other network servers. Server 104 is used to execute the pollution monitoring method based on stable isotopes. Terminal 102 can be, but is not limited to, various personal computers, laptops, smartphones, tablets, IoT devices, and portable wearable devices. IoT devices can include smart speakers, smart TVs, smart air conditioners, smart vehicle devices, etc. Portable wearable devices can include smartwatches, smart bracelets, head-mounted devices, etc. Server 104 can be implemented using a standalone server or a server cluster consisting of multiple servers.

[0030] To address the aforementioned problems, in one embodiment of this application, such as Figure 2 As shown, a pollution monitoring method based on stable isotopes is provided, including the following steps:

[0031] Step 201: Obtain the volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters of the stable isotopes of the target pollutant, as well as the sample fractionation data corresponding to the target pollutant.

[0032] The target pollutant is a high-density non-aqueous liquid pollutant, such as o-dichlorobenzene. Sample fractionation data includes the molar amounts and stable isotopic abundances of the target pollutant in the NAPL phase, aqueous phase, gas phase, and adsorbed phase under stable conditions in any environment. Molar amount refers to the amount of substance of the target pollutant. Isotopic abundance is the ratio of the isotopic ratio of the target pollutant in a specific sample or system to the isotopic ratio of the target pollutant in nature. The NAPL phase is the target pollutant in its independent, water-insoluble liquid form. The aqueous phase is the state in which the target pollutant is dissolved in water. The gas phase is the state in which the target pollutant volatilizes into the air. The adsorbed phase is the state in which the target pollutant is adsorbed into the soil. The volatile fractionation parameter is the volatile fractionation factor of the target pollutant in any environment, used to characterize the volatilization of the target pollutant in any environment. The dissolution fractionation parameter is the dissolution fractionation factor of the target pollutant in any environment, used to characterize the dissolution of the target pollutant in any environment. The adsorption fractionation parameter is the adsorption fractionation factor of the target pollutant in any environment. The adsorption fractionation factor is used to characterize the adsorption of the target pollutant in any environment.

[0033] Obtain sample fractionation data corresponding to the target pollutant, and conduct volatilization, dissolution and adsorption experiments on the target pollutant to determine the volatilization fractionation parameters, dissolution fractionation parameters and adsorption fractionation parameters of the target pollutant.

[0034] It should be noted that the sample fractionation data may be in one or more sets.

[0035] Step 202: Construct an initial evaluation model based on the volatile fractionation parameters, dissolution fractionation parameters, and adsorption fractionation parameters.

[0036] The initial assessment model is an initial model used to assess the target pollutant.

[0037] Based on the principle of mass conservation, an initial evaluation model is constructed according to the volatile fractionation parameters, dissolution fractionation parameters, and adsorption fractionation parameters.

[0038] For example, the initial evaluation model is:

[0039]

[0040] Where Mn, Mw, Mg, and Ms represent the molar quantities of the target pollutant in the NAPL, aqueous, gaseous, and adsorbed phases, respectively; Rn, Rw, Rg, and Rs represent the isotopic ratios of the target pollutant in the NAPL, aqueous, gaseous, and adsorbed phases, respectively; δn, δw, δg, and δs represent the isotopic abundances of the target pollutant in the NAPL, aqueous, gaseous, and adsorbed phases, respectively; kdiss is the first constant; kvol is the second constant; and ksorb is the third constant. The isotopic ratio R is obtained by dividing the isotopic abundance of the target pollutant by 1000, adding 1, and then multiplying the result by the natural isotopic ratio.

[0041] Step 203: Train the initial evaluation model based on the sample fractionation data to determine the target evaluation model.

[0042] The target assessment model is the final model used to assess the target pollutants.

[0043] Based on the sample fractionation data, the first constant, the second constant, and the third constant in the initial evaluation model are fitted to determine the target evaluation model.

[0044] Wherein, the first constant is the dissolution rate constant corresponding to the target pollutant, the second constant is the volatilization rate constant corresponding to the target pollutant, and the third constant is the adsorption rate constant corresponding to the target pollutant.

[0045] The molar quantities and isotopic abundances of the NAPL phase, aqueous phase, gas phase, and adsorbed phase corresponding to the sample fractionation data are input into the initial evaluation model. Simultaneous solutions are then obtained based on sample fractionation data under various stability conditions to fit the first, second, and third constants. The initial evaluation model with the first, second, and third constants obtained is the target evaluation model.

[0046] It should be noted that when there are multiple sets of sample fractionation data, the average values ​​of the first, second, and third constants obtained from each fitting will be calculated to obtain the final first, second, and third constants.

[0047] It should be noted that, in this embodiment, the initial evaluation model and the target evaluation model are in the form of multiple formulas.

[0048] Step 204: Obtain partial fractionation data of the target pollutant.

[0049] The partial fractionation data refers to the fractionation data of the target pollutant in the first target phase in the target environment, including the molar quantity and isotopic abundance of the target pollutant in the first target phase in the target environment. The target environment is the actual environment that needs to be monitored. The first target phase can be any one of the following phases: NAPL phase, aqueous phase, gas phase, or adsorbed phase.

[0050] Step 205: Evaluate the local fractionation data according to the target assessment model to determine the fractionation data of the second target phase of the target pollutant in the target environment.

[0051] The second target phase is any phase other than the first target phase.

[0052] The fractionation data is input into the target assessment model to determine the fractionation data of the second target phase of the target pollutant in the target environment.

[0053] The fractionation data of the target pollutant in the second target phase in the target environment includes the molar quantity and isotopic abundance of the target pollutant in the second target phase in the target environment.

[0054] For example, when the first target phase is the NAPL phase, the second target phase is the aqueous phase, the gas phase, and the adsorbed phase.

[0055] Step 206: Determine the mass flux of the target pollutant in the target environment based on the partial fractionation data and the fractionation data of the second target phase.

[0056] Based on the preset flux calculation formula, and according to the first constant, the second constant, and the molar quantity of the target pollutant in the NAPL phase in the target environment, the mass flux of the target pollutant in the dissolution process and the mass flux of the target pollutant in the volatilization process in the target environment are determined.

[0057] The mass flux of a target pollutant during its dissolution process in the target environment is the number of moles of the target pollutant passing through a unit area per unit time during the dissolution process in the target environment. The mass flux of a target pollutant during its volatilization process in the target environment is the number of moles of the target pollutant passing through a unit area per unit time during the volatilization process in the target environment.

[0058] The aforementioned pollution monitoring method based on stable isotopes involves acquiring the volatile fractionation parameters, dissolved fractionation parameters, adsorption fractionation parameters, and corresponding sample fractionation data of the target pollutant; constructing an initial assessment model based on the volatile fractionation parameters, dissolved fractionation parameters, and adsorption fractionation parameters; training the initial assessment model using the sample fractionation data to determine the target assessment model; acquiring partial fractionation data of the target pollutant; the partial fractionation data representing the fractionation data of the target pollutant in the first target phase in the target environment; evaluating the partial fractionation data using the target assessment model to determine the fractionation data of the target pollutant in the second target phase in the target environment; and determining the mass flux of the target pollutant in the target environment based on the partial fractionation data and the fractionation data of the second target phase. This achieves efficient monitoring of high-density non-aqueous liquid pollutants.

[0059] It should be noted that during the migration and transformation of high-density non-aqueous liquid pollutants in groundwater, the dissolved phase, adsorbed phase, and gas phase are in a dynamic equilibrium process. A change in any one phase will affect the other two. Utilizing stable isotope tracing technology helps to accurately characterize the phase transformation and migration processes of high-density non-aqueous liquid pollutants. By sampling and analyzing the changes in the stable isotope ratios of each phase of high-density non-aqueous liquid pollutants, the concentration and mass flux of pollutants in each phase at water-soil and air-soil interfaces can be monitored.

[0060] It should be noted that in physicochemical processes, such as volatilization, dissolution, adsorption, and diffusion, the isotopic composition of pollutants changes before and after the reaction due to differences in kinetic or thermodynamic effects, a phenomenon known as isotopic fractionation. For example, during volatilization, lighter isotopes (¹²C, ¹H) are more readily volatilized than heavier isotopes (¹³C, ²D), leading to an enrichment of lighter isotopes in the gas phase and heavier isotopes in the residual DNAPL. Therefore, by monitoring changes in the deuterium-to-hydrogen ratio (δ²H) in the gas, aqueous, and DNAPL phases, the intensity and pathway of these processes can theoretically be deduced.

[0061] The essence of moving from multi-process fractionation models to flux calculation models is to use isotopes as "process tracers" to constrain and correct traditional concentration-based flux models. The most core innovation in this embodiment is the establishment of a quantitative conversion relationship between "isotope fractionation signals" and "interfacial mass transfer rates," thereby solving the industry pain point of inaccurate flux assessment caused by pool-like distribution in DNAPL-contaminated sites.

[0062] It should be noted that the core of this embodiment is that, previously, no stable isotope method had been used to assess the multiphase distribution of target pollutants in soil and groundwater. Previously, methods relied solely on drilling and sampling, assessing only the phases encountered, leaving the phase distribution of the target pollutant underground as a black box. The stable isotope-based pollution monitoring method provided in this embodiment can distinguish the distribution and mass flux of different phases of the target pollutant.

[0063] In other embodiments of this application, such as Figure 3 As shown, the acquisition of stable isotope volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters of the target pollutant, and sample fractionation data corresponding to the target pollutant includes:

[0064] Step 301: Obtain sample fractionation data corresponding to the target pollutant.

[0065] The sample fractionation data are the molar quantities and isotopic abundances of the target pollutant in the NAPL phase, aqueous phase, gas phase, and adsorbed phase under any stable environment and conditions.

[0066] Step 302: Conduct a volatilization experiment on the target pollutant to determine the volatilization fractionation parameters of the stable isotopes of the target pollutant.

[0067] The volatile fractionation parameter is the volatile fractionation factor of the target pollutant in any environment. The volatile fractionation factor is used to characterize the volatilization of the target pollutant in any environment.

[0068] Based on a pre-defined volatilization environment, volatilization experiments are conducted to determine the volatile isotope ratios and volatilization proportions of the target pollutant in the NAPL phase after volatilization. Then, the first volatile isotope abundance is determined based on the volatile isotope ratios and the corresponding natural isotope ratios of the target pollutant. The second volatile isotope abundance is determined based on the pre-defined isotope ratios and the corresponding natural isotope ratios of the target pollutant. Finally, the volatilization fractionation factor of the target pollutant is determined based on the first volatile isotope abundance, the second volatile isotope abundance, and the volatilization proportion. The volatilization fractionation factor is used as the volatilization fractionation parameter.

[0069] It should be noted that the evaporation rate was obtained using traditional testing methods.

[0070] The volatile environment is a pre-defined volatile environment. In this embodiment, it is a sealed glass bottle with a sampling port, containing only the NAPL phase and the gas phase. The target pollutants include isotopic pollutants and common pollutants. The isotopic pollutants and common pollutants have the same chemical formula, but the number of neutrons in the hydrogen atoms composing them differs. For example, the hydrogen atoms in the isotopic pollutants are... The hydrogen atoms of the pollutants are The volatile isotope ratio is the ratio between the molar number of the target pollutant in the NAPL phase after volatilization and the molar number of the common pollutant. The volatilization ratio is the ratio between the molar number of the target pollutant in the NAPL phase after volatilization and the molar number of the target pollutant in the APL phase before volatilization.

[0071] The natural isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant and the molar number of common pollutants in nature. The first volatile isotope abundance is the part-of-thousands difference between the volatile isotope ratio and the natural isotope ratio. The preset isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant and the molar number of common pollutants in the NAPL phase before volatilization. The second volatile isotope abundance is the part-of-thousands difference between the preset isotope ratio and the natural isotope ratio.

[0072] Step 303: Conduct a dissolution experiment on the target pollutant to determine the dissolution and fractionation parameters of the stable isotopes of the target pollutant.

[0073] The dissolution fractionation parameter is the dissolution fractionation factor of the target pollutant in any environment. The dissolution fractionation factor is used to characterize the dissolution of the target pollutant in any environment.

[0074] Based on the preset dissolution environment, dissolution experiments were conducted to determine the first dissolved isotope ratio of the target pollutant in the NAPL phase and the second dissolved isotope ratio of the target pollutant in the aqueous phase after dissolution. Then, based on the first dissolved isotope ratio and the corresponding natural isotope ratio of the target pollutant, the abundance of the first dissolved isotope was determined. Then, based on the second dissolved isotope ratio and the corresponding natural isotope ratio of the target pollutant, the abundance of the second dissolved isotope was determined. Finally, based on the abundance of the first and second dissolved isotopes, the dissolution fractionation factor of the target pollutant was determined, and the dissolution fractionation factor was used as the dissolution fractionation parameter.

[0075] In this embodiment, the dissolution environment is a sealed reaction bottle, such as a water-sealed bottle, that is almost filled with water and has a very small headspace. The sealed reaction bottle contains a magnetic stir bar and contains a NAPL phase and an aqueous phase.

[0076] The first dissolved isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant phase and the molar number of common pollutants in the dissolved NAPL phase. The second dissolved isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant phase and the molar number of common pollutants in the dissolved aqueous phase. The natural isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant phase and the molar number of common pollutants in nature. The first dissolved isotope abundance is the isotope abundance corresponding to the first dissolved isotope ratio and the natural isotope ratio. The second dissolved isotope abundance is the isotope abundance corresponding to the second dissolved isotope ratio and the natural isotope ratio.

[0077] Step 304: Conduct adsorption experiments on the target pollutant to determine the adsorption and fractionation parameters of the stable isotopes of the target pollutant.

[0078] The adsorption fractionation parameter is the adsorption fractionation factor of the target pollutant in any environment. The adsorption fractionation factor is used to characterize the adsorption of the target pollutant in any environment.

[0079] Based on the preset adsorption environment, adsorption experiments are conducted to determine the first adsorbed isotope ratio and adsorption ratio of the target pollutant in the target pollutant solution after adsorption. Then, the first adsorbed isotope abundance is determined according to the first adsorbed isotope ratio and the corresponding natural isotope ratio of the target pollutant. Next, the second adsorbed isotope abundance is determined according to the first adsorbed isotope abundance, the adsorption ratio, and the preset isotope abundance. Finally, the adsorption fractionation factor of the target pollutant is determined according to the first and second adsorbed isotope abundances, and the adsorption fractionation factor is used as the adsorption fractionation parameter.

[0080] The preset isotopic abundance is the isotopic abundance of the target pollutant solution before adsorption.

[0081] In this embodiment, the adsorption environment is a set of centrifuge tubes containing soil medium and a dissolved o-dichlorobenzene solution with a preset isotopic abundance.

[0082] The first adsorption isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant solution after adsorption and the molar number of common pollutants. The adsorption ratio is the ratio between the molar number of target pollutants in the adsorbed phase after adsorption and the molar number of target pollutants in the target pollutant solution before adsorption. The natural isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant in nature and the molar number of common pollutants. The first adsorption isotope abundance is the isotope abundance between the first adsorption isotope ratio and the natural isotope ratio. The second adsorption isotope abundance is the isotope abundance between the isotope ratio of the target pollutant in the adsorbed phase after adsorption and the natural isotope ratio.

[0083] In other embodiments of this application, a volatilization experiment is conducted on the target pollutant to determine the volatilization fractionation parameters of the stable isotopes of the target pollutant, including:

[0084] Step 1: Based on the preset volatilization environment, conduct a volatilization experiment to determine the volatilization isotope ratio and volatilization proportion of the target pollutants in the NAPL phase after volatilization.

[0085] It should be noted that the volatile isotope ratio and volatile proportion were obtained through traditional detection methods.

[0086] The evaporation environment is a pre-defined evaporation environment, which in this embodiment is a sealed glass bottle with a sampling port, containing only the NAPL phase and the gas phase. The dissolution experiment involves placing the glass bottle at a constant temperature until the environment inside the glass bottle stabilizes.

[0087] The target pollutants include isotopic pollutants and ordinary pollutants. The isotopic pollutants and ordinary pollutants have the same chemical formula, but differ in the number of neutrons in the hydrogen atoms that make up the isotopic pollutants and ordinary pollutants. For example, the hydrogen atoms in the isotopic pollutants are... The hydrogen atoms of the pollutants are The volatile isotope ratio is the ratio between the molar number of the target pollutant in the NAPL phase after volatilization and the molar number of the common pollutant. The volatilization ratio is the ratio between the molar number of the target pollutant in the NAPL phase after volatilization and the molar number of the target pollutant in the APL phase before volatilization.

[0088] Step 2: Determine the abundance of the first volatile isotope based on the ratio of volatile isotopes and the ratio of natural isotopes corresponding to the target pollutant.

[0089] The natural isotope ratio is the ratio between the molar number of isotopic pollutants in a target pollutant and the molar number of common pollutants in nature. The first volatile isotope abundance is the difference in parts per thousand between the volatile isotope ratio and the natural isotope ratio.

[0090] Step 3: Determine the abundance of the second volatile isotope based on the preset isotope ratio and the natural isotope ratio corresponding to the target pollutant.

[0091] The preset isotope ratio includes the isotope ratio of the target pollutant in the NAPL phase before volatilization. The preset isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant in the NAPL phase before volatilization and the molar number of ordinary pollutants. The second volatile isotope abundance is the isotope abundance between the preset isotope ratio and the natural isotope ratio.

[0092] Step 4: Determine the volatile fractionation factor of the target pollutant based on the abundance of the first volatile isotope, the abundance of the second volatile isotope, and the volatile ratio, and use the volatile fractionation factor as the volatile fractionation parameter.

[0093] The volatile fractionation parameter is the volatile fractionation factor of the target pollutant in any environment. The volatile fractionation factor is used to characterize the volatilization of the target pollutant in any environment.

[0094] For example, the formula for calculating the volatile fractionation factor of a target pollutant is as follows:

[0095]

[0096] The abundance of the first volatile isotope, The abundance of the second volatile isotope is given by f, where f is the volatile proportion. It is the volatile fractionation factor.

[0097] In other embodiments of this application, a dissolution experiment is conducted on the target pollutant to determine the dissolution and fractionation parameters of the stable isotopes of the target pollutant, including:

[0098] Step 1: Based on the preset dissolution environment, conduct a dissolution experiment to determine the first dissolved isotope ratio of the target pollutant in the NAPL phase and the second dissolved isotope ratio of the target pollutant in the aqueous phase after dissolution.

[0099] It should be noted that the first dissolved isotope ratio and the second dissolved isotope ratio of the target pollutant in the aqueous phase were obtained by traditional detection methods.

[0100] In this embodiment, the dissolution environment is a nearly water-filled, sealed reaction flask with minimal headspace, such as a water-sealed flask, containing a magnetic stir bar. The sealed reaction flask contains both the NAPL phase and the aqueous phase. The dissolution experiment involves high-speed stirring to ensure sufficient contact between the NAPL phase and the aqueous phase until the environment within the sealed reaction flask stabilizes.

[0101] The first dissolved isotope ratio is the ratio between the molar number of the target pollutant and the molar number of common pollutants in the dissolved NAPL phase. The second dissolved isotope ratio is the ratio between the molar number of the target pollutant and the molar number of common pollutants in the dissolved aqueous phase.

[0102] It should be noted that the target pollutant has low solubility in water, so solid phase microextraction and other techniques are needed to enrich the dissolved o-dichlorobenzene in the water in order to calculate the second dissolved isotope ratio.

[0103] Step 2: Determine the abundance of the first dissolved isotope based on the first dissolved isotope ratio and the natural isotope ratio corresponding to the target pollutant.

[0104] The natural isotope ratio is the ratio between the molar number of isotopic pollutants in a target pollutant and the molar number of common pollutants in nature. The first dissolved isotope abundance is the difference in per mille between the first dissolved isotope ratio and the natural isotope ratio.

[0105] Step 3: Determine the abundance of the second dissolved isotope based on the second dissolved isotope ratio and the natural isotope ratio corresponding to the target pollutant.

[0106] The second dissolved isotope abundance is the isotope abundance corresponding to the second dissolved isotope ratio and the natural isotope ratio.

[0107] Step 4: Determine the solubility fractionation factor of the target pollutant based on the first and second solubility isotope abundances, and use the solubility fractionation factor as the solubility fractionation parameter.

[0108] The dissolution fractionation parameter is the dissolution fractionation factor of the target pollutant in any environment. The dissolution fractionation factor is used to characterize the dissolution of the target pollutant in any environment.

[0109] For example, the formula for calculating the solubility fractionation factor of a target pollutant is as follows:

[0110]

[0111] The first dissolved isotope abundance, The second dissolved isotope abundance, This is the solubility fractionation factor.

[0112] In other embodiments of this application, adsorption experiments are conducted on the target pollutant to determine the adsorption fractionation parameters of the stable isotopes of the target pollutant, including:

[0113] Step 1: Based on the preset adsorption environment, conduct an adsorption experiment to determine the first adsorbed isotope ratio and adsorption ratio of the target pollutant in the target pollutant solution after adsorption.

[0114] It should be noted that the first adsorbed isotope ratio and adsorption ratio were obtained using traditional detection methods.

[0115] In this embodiment, the adsorption environment is a set of centrifuge tubes containing soil medium and a dissolved o-dichlorobenzene solution with a preset isotopic abundance. The adsorption experiment is conducted by oscillation at a constant temperature until the environment inside the centrifuge tubes stabilizes and adsorption equilibrium is reached.

[0116] The first adsorption isotope ratio is the ratio between the molar number of isotopic pollutants in the target pollutant solution after adsorption and the molar number of ordinary pollutants. The adsorption ratio is the ratio between the molar number of target pollutants in the adsorbed phase after adsorption and the molar number of target pollutants in the target pollutant solution before adsorption.

[0117] Step 2: Determine the abundance of the first adsorbed isotope based on the first adsorbed isotope ratio and the natural isotope ratio corresponding to the target pollutant.

[0118] The natural isotope ratio is the ratio between the molar number of isotopic pollutants in a target pollutant and the molar number of common pollutants in nature. The first adsorbed isotope abundance is the difference in parts per thousand between the first adsorbed isotope ratio and the natural isotope ratio.

[0119] Step 3: Determine the second adsorbed isotope abundance based on the first adsorbed isotope abundance, adsorption ratio, and preset isotope abundance.

[0120] The second adsorption isotope abundance is the isotope abundance between the isotope ratio of the target pollutant in the adsorbed phase after adsorption and the natural isotope ratio. The preset isotope abundance includes the parts-per-thousands difference of the target pollutant in the target pollutant solution before adsorption.

[0121] For example, the formula for determining the abundance of the second adsorbed isotope is as follows:

[0122]

[0123] To preset isotope abundance, The abundance of the first adsorbed isotope. The abundance of the second adsorbed isotope. This represents the adsorption ratio.

[0124] Step 4: Determine the adsorption fractionation factor of the target pollutant based on the abundance of the first and second adsorption isotopes, and use the adsorption fractionation factor as the adsorption fractionation parameter.

[0125] The adsorption fractionation parameter is the adsorption fractionation factor of the target pollutant in any environment. The adsorption fractionation factor is used to characterize the adsorption of the target pollutant in any environment.

[0126] For example, the formula for calculating the adsorption fractionation factor of the target pollutant is as follows:

[0127]

[0128] The abundance of the second adsorbed isotope. The abundance of the first adsorbed isotope. It is the adsorption fractionation factor.

[0129] In other embodiments of this application, training an initial evaluation model based on sample fractionation data to determine a target evaluation model includes:

[0130] Based on the sample fractionation data, the first, second, and third constants in the initial evaluation model are fitted to determine the target evaluation model.

[0131] The first constant is the rate constant for the dissolution of the target pollutant, the second constant is the rate constant for the volatilization of the target pollutant, and the third constant is the rate constant for the adsorption of the target pollutant.

[0132] The sample fractionation data is substituted into the initial evaluation model to calculate the first constant, the second constant, and the third constant, thereby updating the initial evaluation model and finally obtaining the target evaluation model.

[0133] It should be noted that by adjusting the first, second, and third constants, the curves of the phase concentrations and δ values ​​over time evaluated by the target assessment model correspond to the experimental data, thereby improving the capability of the target assessment model. This lays the foundation for subsequent evaluation of local fractionation data based on the target assessment model to determine the fractionation data of the second target phase of the target pollutant in the target environment.

[0134] In other embodiments of this application, the fractionation data of the target pollutant in the second target phase in the target environment includes the molar quantity and isotopic abundance of the target pollutant in the second target phase in the target environment.

[0135] Based on partial fractionation data and fractionation data of the second target phase, the mass flux of the target pollutant in the target environment is determined to include:

[0136] Step 1: Based on the first preset flux calculation formula, determine the mass flux of the target pollutant in the dissolution process in the target environment according to the first constant and the molar quantity of the target pollutant in the NAPL phase in the target environment.

[0137] The first preset flux calculation formula is a pre-defined formula used to calculate the mass flux of the target pollutant during the dissolution process in the target environment. The mass flux of the target pollutant during the dissolution process in the target environment is the number of moles of the target pollutant passing through a unit area per unit time during the dissolution process in the target environment.

[0138] For example, the first preset flux calculation formula is:

[0139]

[0140] Let be the mass flux of the target pollutant during the dissolution process in the target environment, k_vol be the second constant, M_n be the molar quantity of the target pollutant in the NAPL phase, A be the interfacial area of ​​the target pollutant during the dissolution process in the target environment, and V be the control volume of the target pollutant during the dissolution process in the target environment. Here, the interfacial area A and the control volume V are known quantities.

[0141] Step 2: Based on the second preset flux calculation formula, determine the mass flux of the target pollutant in the volatilization process in the target environment according to the second constant and the molar quantity of the target pollutant in the NAPL phase in the target environment.

[0142] The second preset flux calculation formula is a pre-defined formula used to calculate the mass flux of the target pollutant during the volatilization process in the target environment. The mass flux of the target pollutant during the volatilization process in the target environment is the number of moles of the target pollutant passing through a unit area per unit time during the volatilization process in the target environment.

[0143] For example, the second preset flux calculation formula is:

[0144]

[0145] in, Let M_n be the mass flux of the target pollutant during its volatilization process in the target environment, k_diss be the first constant, M_n be the molar quantity of the target pollutant in the NAPL phase, A be the interfacial area of ​​the target pollutant during its volatilization process in the target environment, and V be the control volume of the target pollutant during its volatilization process in the target environment. Here, the interfacial area A and the control volume V are known quantities.

[0146] It should be understood that although the steps in the flowcharts of the embodiments described above are shown sequentially according to the arrows, these steps are not necessarily executed in the order indicated by the arrows. Unless explicitly stated herein, there is no strict order restriction on the execution of these steps, and they can be executed in other orders. Moreover, at least some steps in the flowcharts of the embodiments described above may include multiple steps or multiple stages. These steps or stages are not necessarily completed at the same time, but can be executed at different times. The execution order of these steps or stages is not necessarily sequential, but can be performed alternately or in turn with other steps or at least some of the steps or stages of other steps.

[0147] Based on the same inventive concept, this application also provides a stable isotope-based pollution monitoring device for implementing the above-described stable isotope-based pollution monitoring method. The solution provided by this device is similar to the solution described in the above-described method; therefore, the specific limitations of one or more stable isotope-based pollution monitoring device embodiments provided below can be found in the limitations of the stable isotope-based pollution monitoring method described above, and will not be repeated here.

[0148] In one embodiment of this application, such as Figure 4 As shown, a pollution monitoring device based on stable isotopes is provided, comprising:

[0149] The first acquisition module 100 is used to acquire the volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters, and sample fractionation data of the stable isotopes of the target pollutant; the target pollutant is a high-density non-aqueous liquid pollutant; the sample fractionation data is the molar quantity and stable isotope abundance of the NAPL phase, aqueous phase, gas phase, and adsorption phase of the target pollutant under stable conditions in any environment.

[0150] The construction module 200 is used to construct an initial evaluation model based on the volatile fractionation parameters, dissolution fractionation parameters, and adsorption fractionation parameters;

[0151] The first determining module 300 is used to train the initial evaluation model based on the sample fractionation data and determine the target evaluation model;

[0152] The second acquisition module 400 is used to acquire partial fractionation data of the target pollutant; the partial fractionation data is the fractionation data of the target pollutant in the first target phase in the target environment; the first target phase is any one of the NAPL phase, aqueous phase, gas phase or adsorbed phase;

[0153] The second determining module 500 is used to evaluate the partial fractionation data according to the target evaluation model and determine the fractionation data of the second target phase of the target pollutant in the target environment; the second target phase is the other phases besides the first target phase.

[0154] The third determining module 600 is used to determine the mass flux of the target pollutant in the target environment based on the partial fractionation data and the fractionation data of the second target phase.

[0155] In other embodiments of this application, the first acquisition module 100 is further configured to acquire sample fractionation data corresponding to the target pollutant; conduct a volatilization experiment for the target pollutant to determine the volatilization fractionation parameters of the stable isotopes of the target pollutant; conduct a dissolution experiment for the target pollutant to determine the dissolution fractionation parameters of the stable isotopes of the target pollutant; and conduct an adsorption experiment for the target pollutant to determine the adsorption fractionation parameters of the stable isotopes of the target pollutant.

[0156] In other embodiments of this application, the first acquisition module 100 is further configured to conduct a volatilization experiment based on a preset volatilization environment to determine the volatile isotope ratio and volatilization ratio of the target pollutant in the NAPL phase after volatilization; determine a first volatile isotope abundance based on the volatile isotope ratio and the natural isotope ratio corresponding to the target pollutant; determine a second volatile isotope abundance based on the preset isotope ratio and the natural isotope ratio corresponding to the target pollutant; the preset isotope ratio includes the isotope ratio of the target pollutant in the NAPL phase before volatilization; determine the volatilization fractionation factor of the target pollutant based on the first volatile isotope abundance, the second volatile isotope abundance, and the volatilization ratio, and use the volatilization fractionation factor as a volatilization fractionation parameter.

[0157] In other embodiments of this application, the first acquisition module 100 is further configured to conduct a dissolution experiment based on a preset dissolution environment, determine the first dissolved isotope ratio of the target pollutant in the NAPL phase and the second dissolved isotope ratio of the target pollutant in the aqueous phase after dissolution; determine the first dissolved isotope abundance based on the first dissolved isotope ratio and the natural isotope ratio corresponding to the target pollutant; determine the second dissolved isotope abundance based on the second dissolved isotope ratio and the natural isotope ratio corresponding to the target pollutant; determine the dissolution fractionation factor of the target pollutant based on the first dissolved isotope abundance and the second dissolved isotope abundance, and use the dissolution fractionation factor as a dissolution fractionation parameter.

[0158] In other embodiments of this application, the first acquisition module 100 is further configured to perform an adsorption experiment based on a preset adsorption environment to determine the first adsorbed isotope ratio and adsorption ratio of the target pollutant in the target pollutant solution after adsorption; determine the first adsorbed isotope abundance based on the first adsorbed isotope ratio and the natural isotope ratio corresponding to the target pollutant; determine the second adsorbed isotope abundance based on the first adsorbed isotope abundance, the adsorption ratio, and the preset isotope abundance; the preset isotope abundance includes the isotope abundance of the target pollutant in the target pollutant solution before adsorption; determine the adsorption fractionation factor of the target pollutant based on the first adsorbed isotope abundance and the second adsorbed isotope abundance, and use the adsorption fractionation factor as an adsorption fractionation parameter.

[0159] In other embodiments of this application, the first determining module 300 is further configured to fit the first constant, the second constant, and the third constant in the initial evaluation model based on the sample fractionation data to determine the target evaluation model; the first constant is the dissolution rate constant corresponding to the target pollutant, the second constant is the volatilization rate constant corresponding to the target pollutant, and the third constant is the adsorption rate constant corresponding to the target pollutant.

[0160] In other embodiments of this application, the third determining module 600 is further configured to determine the mass flux of the target pollutant in the dissolution process in the target environment based on the first preset flux calculation formula, according to the first constant and the molar quantity of the target pollutant in the NAPL phase in the target environment; and to determine the mass flux of the target pollutant in the volatilization process in the target environment based on the second preset flux calculation formula, according to the second constant and the molar quantity of the target pollutant in the NAPL phase in the target environment.

[0161] The modules in the aforementioned pollution monitoring device based on stable isotopes can be implemented entirely or partially through software, hardware, or a combination thereof. These modules can be embedded in the processor of a computer device in hardware form or independent of it, or stored in the memory of the computer device in software form, so that the processor can call and execute the corresponding operations of each module.

[0162] In one embodiment of this application, a computer device is provided, which may be a server, and its internal structure diagram may be as follows. Figure 5 As shown, the computer device includes a processor, memory, and a network interface connected via a system bus. The processor provides computational and control capabilities. The memory includes non-volatile storage media and internal memory. The non-volatile storage media stores the operating system, computer programs, and a database. The internal memory provides an environment for the operation of the operating system and computer programs in the non-volatile storage media. The database stores all relevant data for implementing a pollution monitoring method based on stable isotopes. The network interface communicates with external terminals via a network connection. When executed by the processor, the computer program implements a pollution monitoring method based on stable isotopes.

[0163] Those skilled in the art will understand that Figure 5 The structure shown is merely a block diagram of a portion of the structure related to the present application and does not constitute a limitation on the computer device to which the present application is applied. Specific computer devices may include more or fewer components than those shown in the figure, or combine certain components, or have different component arrangements.

[0164] In one embodiment of this application, a computer device is provided, including a memory and a processor. The memory stores a computer program, and the processor executes the computer program to implement the steps of the pollution monitoring method based on stable isotopes in the above embodiments.

[0165] In one embodiment of this application, a computer-readable storage medium is provided, on which a computer program is stored, the computer program being executed by a processor to implement the steps of the pollution monitoring method based on stable isotopes in the above-described method embodiments.

[0166] In one embodiment of this application, a computer program product is provided, including a computer program that, when executed by a processor, implements the steps of the pollution monitoring method based on stable isotopes in the above-described method embodiments.

[0167] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data used for analysis, data stored, data displayed, etc.) involved in this application are all information and data authorized by the user or fully authorized by all parties.

[0168] Those skilled in the art will understand that all or part of the processes in the methods of the above embodiments can be implemented by a computer program instructing related hardware. The computer program can be stored in a non-volatile computer-readable storage medium. When executed, the computer program can include the processes of the embodiments of the above methods. Any references to memory, databases, or other media used in the embodiments provided in this application can include at least one of non-volatile and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. By way of illustration and not limitation, RAM can take many forms, such as Static Random Access Memory (SRAM) or Dynamic Random Access Memory (DRAM). The databases involved in the embodiments provided in this application may include at least one type of relational database and non-relational database. Non-relational databases may include, but are not limited to, blockchain-based distributed databases. The processors involved in the embodiments provided in this application may be general-purpose processors, central processing units, graphics processing units, digital signal processors, programmable logic devices, quantum computing-based data processing logic devices, etc., and are not limited to these.

[0169] The technical features of the above embodiments can be combined in any way. For the sake of brevity, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this specification.

[0170] The embodiments described above are merely illustrative of several implementation methods of this application, and while the descriptions are specific and detailed, they should not be construed as limiting the scope of this application. It should be noted that those skilled in the art can make various modifications and improvements without departing from the concept of this application, and these all fall within the protection scope of this application. Therefore, the protection scope of this application should be determined by the appended claims.

Claims

1. A pollution monitoring method based on stable isotopes, characterized in that, The method includes: The volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters, and sample fractionation data of the stable isotopes of the target pollutant are obtained; the target pollutant is a high-density non-aqueous liquid pollutant; the sample fractionation data are the molar quantities and stable isotope abundances of the NAPL phase, aqueous phase, gas phase, and adsorption phase of the target pollutant under any environment and stable conditions. An initial evaluation model is constructed based on the volatile fractionation parameters, dissolution fractionation parameters, and adsorption fractionation parameters. The initial evaluation model is trained based on the sample fractionation data to determine the target evaluation model; Obtain partial fractionation data of the target pollutant; the partial fractionation data are the fractionation data of the target pollutant in the first target phase in the target environment; the first target phase is any one of the NAPL phase, aqueous phase, gas phase or adsorbed phase; The partial fractionation data are evaluated based on the target assessment model to determine the fractionation data of the target pollutant in the second target phase in the target environment; the second target phase is the other phases besides the first target phase. Based on the partial fractionation data and the fractionation data of the second target phase, the mass flux of the target pollutant in the target environment is determined.

2. The pollution monitoring method based on stable isotopes according to claim 1, characterized in that, The acquisition of stable isotope volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters of the target pollutant, and sample fractionation data corresponding to the target pollutant includes: Obtain sample fractionation data corresponding to the target pollutant; For the target pollutant, a volatilization experiment was conducted to determine the volatilization fractionation parameters of the stable isotopes of the target pollutant. Dissolution experiments were conducted on the target pollutant to determine the dissolution and fractionation parameters of the stable isotopes of the target pollutant. Adsorption experiments were conducted on the target pollutant to determine the adsorption and fractionation parameters of the stable isotopes of the target pollutant.

3. The pollution monitoring method based on stable isotopes according to claim 2, characterized in that, The process of conducting volatilization experiments on the target pollutant to determine the volatilization fractionation parameters of the stable isotopes of the target pollutant includes: Based on the preset volatilization environment, volatilization experiments were conducted to determine the volatilization isotope ratio and volatilization proportion of the target pollutants in the NAPL phase after volatilization. The first volatile isotope abundance is determined based on the volatile isotope ratio and the natural isotope ratio corresponding to the target pollutant. The second volatile isotope abundance is determined based on the preset isotope ratio and the natural isotope ratio corresponding to the target pollutant; the preset isotope ratio includes the isotope ratio of the target pollutant in the NAPL phase before volatilization. Based on the first volatile isotope abundance, the second volatile isotope abundance, and the volatile ratio, the volatile fractionation factor of the target pollutant is determined, and the volatile fractionation factor is used as the volatile fractionation parameter.

4. The pollution monitoring method based on stable isotopes according to claim 2, characterized in that, The process of conducting dissolution experiments on the target pollutant to determine the dissolution and fractionation parameters of its stable isotopes includes: Based on the preset dissolution environment, a dissolution experiment was conducted to determine the first dissolution isotope ratio of the target pollutant in the NAPL phase and the second dissolution isotope ratio of the target pollutant in the aqueous phase after dissolution. The abundance of the first dissolved isotope is determined based on the first dissolved isotope ratio and the natural isotope ratio corresponding to the target pollutant. The abundance of the second dissolved isotope is determined based on the second dissolved isotope ratio and the natural isotope ratio corresponding to the target pollutant. Based on the first and second dissolved isotope abundances, the solubility fractionation factor of the target pollutant is determined, and the solubility fractionation factor is used as the solubility fractionation parameter.

5. The pollution monitoring method based on stable isotopes according to claim 2, characterized in that, The adsorption experiments conducted on the target pollutant to determine the adsorption and fractionation parameters of its stable isotopes include: Based on the preset adsorption environment, an adsorption experiment was conducted to determine the first adsorbed isotope ratio and adsorption ratio of the target pollutant in the target pollutant solution after adsorption. The abundance of the first adsorbed isotope is determined based on the first adsorbed isotope ratio and the natural isotope ratio corresponding to the target pollutant. The second adsorbed isotope abundance is determined based on the first adsorbed isotope abundance, the adsorption ratio, and the preset isotope abundance; the preset isotope abundance includes the isotope abundance of the target pollutant in the target pollutant solution before adsorption. Based on the abundance of the first adsorbed isotope and the abundance of the second adsorbed isotope, the adsorption fractionation factor of the target pollutant is determined, and the adsorption fractionation factor is used as the adsorption fractionation parameter.

6. The pollution monitoring method based on stable isotopes according to claim 1, characterized in that, The step of training the initial evaluation model based on the sample fractionation data to determine the target evaluation model includes: Based on the sample fractionation data, the first constant, the second constant, and the third constant in the initial evaluation model are fitted to determine the target evaluation model; the first constant is the dissolution rate constant corresponding to the target pollutant, the second constant is the volatilization rate constant corresponding to the target pollutant, and the third constant is the adsorption rate constant corresponding to the target pollutant.

7. The pollution monitoring method based on stable isotopes according to claim 6, characterized in that, The fractionation data of the target pollutant in the second target phase in the target environment includes the molar quantity and isotopic abundance of the target pollutant in the second target phase in the target environment. Based on the partial fractionation data and the fractionation data of the second target phase, the mass flux of the target pollutant in the target environment is determined, including: Based on the first preset flux calculation formula, the mass flux of the target pollutant in the target environment during the dissolution process is determined according to the first constant and the molar amount of the target pollutant in the NAPL phase in the target environment. Based on the second preset flux calculation formula, and according to the second constant and the molar quantity of the target pollutant in the NAPL phase in the target environment, the mass flux of the target pollutant in the volatilization process in the target environment is determined.

8. A pollution monitoring device based on stable isotopes, characterized in that, The device includes: The first acquisition module is used to acquire the volatile fractionation parameters, dissolution fractionation parameters, adsorption fractionation parameters, and sample fractionation data of the stable isotopes of the target pollutant; the target pollutant is a high-density non-aqueous liquid pollutant; the sample fractionation data is the molar quantity and stable isotope abundance of the NAPL phase, aqueous phase, gas phase, and adsorption phase of the target pollutant under any environment and stable conditions. A construction module is used to construct an initial evaluation model based on the volatile fractionation parameters, dissolution fractionation parameters, and adsorption fractionation parameters; The first determining module is used to train the initial evaluation model based on the sample fractionation data to determine the target evaluation model; The second acquisition module is used to acquire partial fractionation data of the target pollutant; the partial fractionation data is the fractionation data of the target pollutant in the first target phase in the target environment; the first target phase is any one of the NAPL phase, aqueous phase, gas phase or adsorbed phase; The second determining module is used to evaluate the partial fractionation data according to the target evaluation model and determine the fractionation data of the second target phase of the target pollutant in the target environment; the second target phase is the other phases besides the first target phase. The third determining module is used to determine the mass flux of the target pollutant in the target environment based on the partial fractionation data and the fractionation data of the second target phase.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, characterized in that, When the processor executes the computer program, it implements the steps of the method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed by a processor, it implements the steps of the method according to any one of claims 1 to 7.