Blowout diffusion simulation method and device for sulfur-containing gas well under complex terrain and meteorological conditions and storage medium

Abstract

The invention discloses a method and equipment for simulating blowout diffusion of a sulfur-containing gas well under complex terrain and meteorological conditions and a storage medium. The method comprises the following steps: collecting meteorological data, terrain data and diffusion source data; judging a boundary layer state according to the meteorological data; calculating scale parameters according to the boundary layer state, the meteorological data and the topographic data; calculating a mixing height, a profile equation and diffusion parameters according to the scale parameters; according to the mixing height, the profile equation and the diffusion parameters, simulation calculation of the sulfur-containing gas diffusion concentration distribution after blowout is conducted, and the total diffusion concentration is obtained. Aiming at the blowout diffusion problem of the sulfur-containing gas well under complex terrain and meteorological conditions, the sulfur-containing gas diffusion rule after out-of-control blowout can be accurately simulated in time, quickly and effectively, a basis can be provided for a safety and environmental protection department to timely formulate effective emergency measures, and a basis is provided for well site arrangement of the sulfur-containing gas well and formulation of safety prevention and control measures.

Description

Technical Field

[0001] The invention relates to a method, equipment and storage medium for simulating the blowout diffusion of a sulfur-containing gas well under complex terrain and meteorological conditions, and belongs to the technical field. Background Art

[0002] The Sichuan Basin has a natural gas reserve of 20.25 trillion cubic meters and has become the main battlefield and ballast stone for natural gas development in my country. However, the production of sour natural gas in the Sichuan Basin accounts for about 80% of the total production, which not only contains hydrocarbon gases such as methane and ethane, but also highly toxic and corrosive gases such as hydrogen sulfide and sulfur dioxide. Blowout accidents may occur during drilling, completion, well testing, well repair and even production. Once a blowout occurs, the sour gas will leak from the wellhead to the atmosphere, forming a high-concentration toxic gas cloud, which will spread with the wind direction and wind speed, posing a serious threat to the health of humans and the ecological environment around it. Accurately simulating the diffusion law of sour gas after a blowout plays an important guiding role in timely, rapid, accurate and effective grasp of the distribution of sour gas, and in the formulation and implementation of emergency and rescue measures. It has become a problem that the safety and environmental protection departments need to solve urgently. However, the diffusion of sour gas is greatly affected by terrain and meteorological conditions. Fluid dynamics simulation software and atmospheric diffusion simulation software are mainly used for simulation at home and abroad. It is difficult to consider the influence of real terrain, meteorology and other factors, and the calculation efficiency is low and the calculation time is long. Summary of the invention

[0003] The purpose of the present invention is to provide a method, device and storage medium for simulating the blowout diffusion of a sulfur-containing gas well under complex terrain and meteorological conditions in order to solve the problems existing in the prior art.

[0004] The present invention provides a technical solution to solve the above technical problems: a method for simulating the diffusion of a blowout of a sulfur-containing gas well under complex terrain and meteorological conditions, comprising the following steps:

[0005] Step 1: Collect meteorological data, topographic data, and diffusion source data;

[0006] Step 2: Determine the boundary layer state based on meteorological data;

[0007] Step 3: Calculate the scale parameters based on the boundary layer state, meteorological data, and terrain data;

[0008] Step 4: Calculate the mixing height, profile equation and diffusion parameter based on the scale parameter;

[0009] Step 5: simulate and calculate the diffusion concentration distribution of the sulfur-containing gas after the blowout based on the mixing height, profile equation and diffusion parameters to obtain the total diffusion concentration.

[0010] A further technical solution is that the meteorological data includes wind speed, wind direction, and temperature; the terrain data includes outlet height and outlet diameter; and the diffusion source data includes hydrogen sulfide content and hydrogen sulfide flow rate.

[0011] A further technical solution is that in step 2, according to the temperature and the cloud coefficient, combined with the following formula, it is determined whether the boundary layer state is in the convective boundary layer or the stable boundary layer;

[0012]

[0013] Where: c 1 =5.31×10 -13 Wm -2 K -6 ;c 2 =60Wm -2 ; σ sb is the Stefin Boltzman constant; T is the temperature; r{φ} is the reflectivity; φ crit is the solar elevation angle; n is the cloud factor.

[0014] A further technical solution is that the scale parameters include Monin-Obuhov length, convection velocity scale, and temperature scale;

[0015] The calculation formula of the Monin-Obuhov length is:

[0016]

[0017] Where: L is the Monin-Obuhov length; ρ is the air density; g is the gravitational acceleration; c p is the specific heat of air at constant pressure; T ref is the reference temperature; H is the height of the gas; k is the von-Karman constant, which is 0.4; u * is the ground friction speed;

[0018] The calculation formula of the convection velocity scale is:

[0019]

[0020] Where: z ic is the tropospheric mixing height; w * is the convection velocity scale;

[0021] The temperature scale is calculated as:

[0022]

[0023] Where: θ * is the temperature scale.

[0024] A further technical solution is that in step 3, when the boundary layer state is under the convective boundary layer condition, the calculation formula of the ground friction speed is:

[0025]

[0026] Where: k is the von-Karman constant, which is 0.4; u ref is the wind speed at reference height; ref is the standard measurement height for measuring wind speed; 0 is the rough length; is the ground friction velocity under flow boundary layer conditions; Ψ m is the stability correction function;

[0027] When the boundary layer state is under stable boundary layer conditions, the calculation formula of the ground friction velocity is:

[0028]

[0029] Where: C D is the traction coefficient; u is the wind speed; is the ground friction velocity under stable boundary layer conditions.

[0030] A further technical solution is that the calculation formula in step 4 includes:

[0031] Mixing height z i The expression is:

[0032] z i =max[z ic ,z im ]

[0033] Where: z ic is the convective mixing height; z im is the mechanical mixing height;

[0034] The wind speed profile equation is expressed as:

[0035]

[0036] The temperature gradient profile equation is expressed as:

[0037]

[0038] Where: z is the measurement height; mx is the starting depth where the temperature gradient begins to decay; z iθ The characteristic depth to describe the temperature gradient decay characteristics;

[0039] Vertical diffusion parameter σ z The calculation formula is:

[0040]

[0041] Where: zc is the vertical convection diffusion parameter; σ zm is the mechanical vertical diffusion parameter;

[0042] Horizontal diffusion parameter σ y The calculation formula is:

[0043]

[0044] Where: yc is the horizontal convection diffusion parameter; σ ym is the mechanical level diffusion parameter.

[0045] A further technical solution is that the calculation formula in step 5 is:

[0046] C c {x r ,y r ,z r =C d {x r ,y r ,z r}+C r {x r ,y r ,z r}+C p {x r ,y r ,z r}

[0047] in,

[0048]

[0049] Where: C c is the total concentration of sulfur-containing gas; C d is the concentration of airflow directly diffused to the ground; C r is the mass concentration of the indirect source diffused by the updraft to the top of the mixed layer; Q is the source intensity; C p is the mass concentration of the penetrating source penetrating into the upper stable layer of the mixed layer; u is the wind speed at the top of the leakage source; F y is the horizontal distribution function; σ y is the horizontal diffusion parameter; y is the crosswind distance; z i is the height of the mixing layer; σ zj is the vertical diffusion parameter of the indirect source; σ zp is the vertical diffusion parameter of the penetrating source; dj is the total height of the direct source plume; rj is the total height of the indirect source plume;j is the weight coefficient of the rising and falling parts of the plume; f p is the plume density under convective conditions; h ep is the penetration source height; z ieff is the height of the reflecting surface in the stable layer.

[0050] A simulation system for blowout diffusion of a sulfur-containing gas well, comprising:

[0051] Data preprocessing module, which preprocesses the meteorological, topographical, diffusion source and other data in the database or imported from outside into the data required for diffusion simulation calculation;

[0052] The diffusion simulation calculation module calculates the temporal and spatial distribution of the leaked gas based on the input current terrain, weather and leakage conditions and by using a blowout diffusion simulation method for sulfur-containing gas wells under complex terrain and weather conditions;

[0053] The calculation result processing module processes the calculated temporal and spatial distribution law of the leaked gas, including data extraction, scalar number processing and vector data generation;

[0054] The visualization module uses the GIS system to perform 2D and 3D visualization of the temporal and spatial distribution of leaked gas, allowing for a clearer and more direct observation of the diffusion results of leaked gas.

[0055] An electronic device, comprising:

[0056] Memory and processor;

[0057] The memory can be used to store data and executable instructions received by the computer, and the processor is used to execute the executable instructions received by the computer, process data and control other components, thereby realizing a method for simulating the blowout diffusion of sulfur-containing gas wells under complex terrain and meteorological conditions.

[0058] A computer-readable storage medium is a physical material or device for storing data or instructions. When the computer-executable instructions are executed by a processor, a method for simulating the blowout diffusion of a sulfur-containing gas well under complex terrain and meteorological conditions is implemented.

[0059] The beneficial effects of the present invention are as follows: the method can timely, rapidly and effectively simulate the diffusion law of sulfur-containing gas after the blowout is out of control, aiming at the problem of blowout diffusion of sulfur-containing gas wells under complex terrain and meteorological conditions, and can provide a basis for the safety and environmental protection departments to formulate effective emergency measures in a timely manner, and provide a basis for the layout of sulfur-containing gas well sites and the formulation of safety prevention and control measures. It can also be used for safety analysis and evaluation of sulfur-containing natural gas exploration and development projects, and has an important role and significance for ensuring the safety of people's lives and protecting the environment. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] Figure 1 It is a schematic diagram of the diffusion simulation process of the present invention;

[0061] Figure 2 To display the simulation results on the three-dimensional terrain;

[0062] Figure 3 It is a structural schematic diagram of a simulation system for blowout diffusion of a sulfur-containing gas well;

[0063] Figure 4 The figure is a schematic diagram of the structure of an electronic device. DETAILED DESCRIPTION

[0064] The technical solution of the present invention will be described clearly and completely below in conjunction with the accompanying drawings. Obviously, the described embodiments are only part of the embodiments of the invention, not all of them. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0065] The present invention provides a method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and meteorological conditions, comprising the following steps:

[0066] Step 1: Collect meteorological data, topographic data, and diffusion source data, including wind speed, wind direction, outlet height, natural gas flow rate, outlet diameter, temperature, hydrogen sulfide content, hydrogen sulfide flow rate and other basic data;

[0067] Step 2: According to parameters such as temperature and cloud coefficient, combined with formula (1), determine whether the boundary layer state is CBL (convective boundary layer) or SBL (stable boundary layer);

[0068]

[0069] Where: c 1 =5.31×10 -13 Wm -2 K -6 ;c 2 =60Wm -2 ; σ sb is the Stefin Boltzman constant; T is the temperature; r{φ} is the reflectivity; φ crit is the sun elevation angle; n is the cloud factor;

[0070] Generally, under cloudless conditions, the average solar elevation angle at the change point from stable to convective is 13°, and on cloudy days it is 23°;

[0071] That is, in cloudless conditions, the solar elevation angle φ crit When ≤13°, the boundary layer is in SBL (stable boundary layer); the solar elevation angle φ critWhen >13°, the boundary layer state is CBL (convective boundary layer);

[0072] In cloudy conditions, the sun's elevation angle φ crit When ≤23°, the boundary layer is in SBL (stable boundary layer); the solar elevation angle φ crit When >23°, the boundary layer state is CBL (convective boundary layer);

[0073] Step 3: Calculate the scale parameters based on the boundary layer state, meteorological data, and terrain data;

[0074] The scale parameters include Monin-Obuhov length, convection velocity scale, and temperature scale;

[0075] The calculation formula of the Monin-Obuhov length is:

[0076]

[0077] Where: L is the Monin-Obuhov length; ρ is the air density; g is the gravitational acceleration; c p is the specific heat of air at constant pressure; T ref is the reference temperature; H is the height of the gas; k is the von-Karman constant, which is 0.4; u * is the ground friction speed;

[0078] The calculation formula of the convection velocity scale is:

[0079]

[0080] Where: z ic is the tropospheric mixing height; w * is the convection velocity scale;

[0081] The temperature scale is calculated as:

[0082]

[0083] Where: θ * is the temperature scale.

[0084] The ground friction speed u * Different calculation formulas are selected according to different boundary layer states; specifically, when the boundary layer state is under the convective boundary layer condition, the calculation formula for the ground friction velocity is:

[0085]

[0086] in,

[0087]

[0088] μ=(1-16z ref / L) 1 / 4 (8)

[0089] μ 0 =(1-16z 0 / L) 1 / 4 (9)

[0090] Where: k is the von-Karman constant, which is 0.4; u ref is the wind speed at reference height; ref is the standard measurement height for measuring wind speed; 0 is the rough length; is the ground friction velocity under flow boundary layer conditions; Ψ m is the stability correction function;

[0091] When the boundary layer state is under stable boundary layer conditions, the calculation formula of the ground friction velocity is:

[0092]

[0093] Where: C D is the traction coefficient; u is the wind speed; is the ground friction velocity under stable boundary layer conditions;

[0094] Step 4: Calculate the mixing height, profile equation and diffusion parameter based on the scale parameter;

[0095] The mixing height z i The expression is:

[0096] z i =max[z ic ,z im ] (11)

[0097]

[0098] Where: z ic is the convective mixing height; z im is the mechanical mixing height;

[0099] The wind speed profile equation is expressed as:

[0100]

[0101] The temperature gradient profile equation is expressed as:

[0102]

[0103] Where: z is the measurement height; mxis the starting depth where the temperature gradient begins to decay, and the present invention takes 100m; iθ is the characteristic depth describing the temperature gradient attenuation characteristics, z iθ =max{z im , 100m}, z im is the mechanical mixing height;

[0104] Vertical diffusion parameter σ z The calculation formula is:

[0105]

[0106] Where: zc is the vertical convection diffusion parameter; σ zm is the mechanical vertical diffusion parameter;

[0107] Horizontal diffusion parameter σ y The calculation formula is:

[0108]

[0109] Where: yc is the horizontal convection diffusion parameter; σ ym is the mechanical level diffusion parameter;

[0110] Step 5: simulate and calculate the diffusion concentration distribution of the sulfur-containing gas after the blowout according to the mixing height, profile equation and diffusion parameters to obtain the total diffusion concentration;

[0111] C c {x r ,y r ,z r =C d {x r ,y r ,z r}+C r {x r ,y r ,z r}+C p {x r ,y r ,z r} (twenty two)

[0112] in,

[0113]

[0114] Where: C c is the total concentration of sulfur-containing gas; C d is the concentration of airflow directly diffused to the ground; C r is the mass concentration of the indirect source diffused by the updraft to the top of the mixed layer; Q is the source intensity; Cp is the mass concentration of the penetrating source penetrating into the upper stable layer of the mixed layer; u is the wind speed at the top of the leakage source; F y is the horizontal distribution function; σ y is the horizontal diffusion parameter; y is the crosswind distance; z i is the height of the mixing layer; σ zj is the vertical diffusion parameter of the indirect source; σ zp is the vertical diffusion parameter of the penetrating source; dj is the total height of the direct source plume; rj is the total height of the indirect source plume; j is the weight coefficient of the rising and falling parts of the plume; f p is the plume density under convective conditions; h ep is the penetration source height; z ieff is the height of the reflecting surface in the stable layer.

[0115] Step 6. Establish a digital elevation model that represents the three-dimensional characteristics of the terrain to realize the import and processing of complex terrain; study the visualization extraction, generation and display methods of diffusion simulation results, establish a three-dimensional terrain and diffusion calculation data visualization model, and realize the three-dimensional GIS spatial display of diffusion simulation results.

[0116] Example

[0117] According to the actual situation of the target area, basic data such as wind speed, wind direction, outlet height, natural gas flow rate, outlet diameter, temperature, hydrogen sulfide content, and hydrogen sulfide flow were obtained, and an example calculation of the diffusion concentration distribution after a short-term sulfur-containing gas leakage was performed. The specific steps are as follows:

[0118] 1. According to the given parameters such as temperature and cloud coefficient, combined with formula (1), determine whether the boundary layer state is CBL or SBL;

[0119] 2. After determining the state of the boundary layer, if it is a convective boundary layer, the meteorological and topographic data such as wind speed, roughness length, leakage source height, etc. are substituted into equations (2) to (9) to calculate the scale parameters; if it is a stable boundary layer, it is substituted into equations (2) to (4) and (10) to calculate;

[0120] 3. After calculating the scale parameters, substitute them into equations (11) to (13) to calculate the mixing height; substitute them into equation (14) to obtain the wind speed profile equation; substitute them into equation (15) to obtain the temperature gradient profile equation; substitute them into equations (16) to (18) to obtain the vertical diffusion parameters; substitute them into equations (19) to (21) to obtain the horizontal diffusion parameters.

[0121] 4. After obtaining the profile equations and diffusion parameters, the total diffusion concentration is calculated by equations (22) to (26);

[0122] The simulation results are shown in Table 1;

[0123] Table 1

[0124] 11 hours measurement point concentration 18-hour measurement point concentration Measured <![CDATA[34.56mg / m 3 ]]> <![CDATA[69.12mg / m 3 ]]> Calculation results <![CDATA[40.832mg / m 3 ]]> <![CDATA[76.368mg / m 3 ]]> error 18％ 10％

[0125] The embodiment of the present application also provides a simulation system for blowout diffusion of a sulfur-containing gas well, including:

[0126] Data preprocessing module 1, preprocessing the meteorological, topographical, diffusion source and other data in the database or imported from outside into the data required for diffusion simulation calculation;

[0127] Diffusion simulation calculation module 2 calculates the temporal and spatial distribution of leaked gas based on the input current terrain, weather and leakage conditions and by using a blowout diffusion simulation method for sulfur-containing gas wells under complex terrain and weather conditions;

[0128] The calculation result processing module 3 processes the calculated temporal and spatial distribution law of the leaked gas, including data extraction, scalar number processing and vector data generation;

[0129] Visualization module 4 uses the GIS system to perform 2D and 3D visualization of the temporal and spatial distribution of leaked gas, so as to observe the diffusion results of leaked gas more clearly and directly.

[0130] The embodiment of the present application also provides an electronic device, including:

[0131] Memory 5 and processor 6;

[0132] The memory 5 can be used to store data and executable instructions received by the computer, and the processor 6 is used to execute the executable instructions received by the computer, process data and control other components, thereby realizing a method for simulating the blowout diffusion of sulfur-containing gas wells under complex terrain and meteorological conditions.

[0133] The embodiment of the present application also provides a computer-readable storage medium, which is a physical material or device for storing data or instructions. When the computer-executable instructions are executed by a processor, a method for simulating the blowout diffusion of a sulfur-containing gas well under complex terrain and meteorological conditions is implemented.

[0134] The above description is not intended to limit the present invention in any form. Although the present invention has been disclosed through the above embodiments, it is not intended to limit the present invention. Any technician familiar with the profession can make some changes or modifications to equivalent embodiments of equivalent changes using the technical contents disclosed above without departing from the scope of the technical solution of the present invention. However, any simple modification, equivalent change and modification made to the above embodiments based on the technical essence of the present invention without departing from the content of the technical solution of the present invention still falls within the scope of the technical solution of the present invention.

Claims

1. A method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and meteorological conditions, characterized in that: The following steps are involved: Step 1: Collect meteorological data, topographic data, and diffusion source data; Step 2: Determine the boundary layer state based on meteorological data; Step 3: Calculate the scale parameters based on the boundary layer state, meteorological data, and terrain data; Step 4: Calculate the mixing height, profile equation and diffusion parameter based on the scale parameter; Step 5: simulate and calculate the diffusion concentration distribution of the sulfur-containing gas after the blowout based on the mixing height, profile equation and diffusion parameters to obtain the total diffusion concentration.

2. The method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and meteorological conditions according to claim 1, characterized in that: The meteorological data include wind speed, wind direction, and temperature; the terrain data include outlet height and outlet diameter; and the diffusion source data include hydrogen sulfide content and hydrogen sulfide flow rate.

3. The method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and meteorological conditions according to claim 2, characterized in that: In the step 2, based on the temperature and the cloud coefficient, the following formula is combined to determine whether the boundary layer state is in the convective boundary layer or the stable boundary layer; Where: c1 = 5.31 × 10 -13 Wm -2 K -6 ; c2=60Wm -2 ; σ sb is the Stefin Boltzman constant; T is the temperature; r{φ} is the reflectivity; φ crit is the solar elevation angle; n is the cloud factor.

4. The method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and meteorological conditions according to claim 1, characterized in that: The scale parameters include Monin-Obuhov length, convection velocity scale, and temperature scale; The calculation formula of the Monin-Obuhov length is: Where: L is the Monin-Obuhov length; ρ is the air density; g is the gravitational acceleration; c p is the specific heat of air at constant pressure; T ref is the reference temperature; H is the height of the gas; k is the von-Karman constant, which is 0.4; u * is the ground friction speed; The calculation formula of the convection velocity scale is: Where: z ic is the tropospheric mixing height; w * is the convection velocity scale; The temperature scale is calculated as: Where: θ * is the temperature scale.

5. The method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and meteorological conditions according to claim 4, characterized in that: In step 3, when the boundary layer state is under the convective boundary layer condition, the calculation formula of the ground friction speed is: Where: k is the von-Karman constant, which is 0.4; u ref is the wind speed at reference height; ref is the standard measurement height when measuring wind speed; z0 is the roughness length; is the ground friction velocity under flow boundary layer conditions; Ψ m is the stability correction function; When the boundary layer state is under stable boundary layer conditions, the calculation formula of the ground friction velocity is: Where: C D is the traction coefficient; u is the wind speed; is the ground friction velocity under stable boundary layer conditions.

6. The method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and meteorological conditions according to claim 1, characterized in that: The calculation formula in step 4 includes: Mixing height z i The expression is: With i =max[z ic ,With im ] Where: z ic is the convective mixing height; z im is the mechanical mixing height; The wind speed profile equation is expressed as: The temperature gradient profile equation is expressed as: Where: z is the measurement height; mx is the starting depth where the temperature gradient begins to decay; z iθ The characteristic depth to describe the temperature gradient decay characteristics; Vertical diffusion parameter σ z The calculation formula is: Where: zc is the vertical convection diffusion parameter; σ zm is the mechanical vertical diffusion parameter; Horizontal diffusion parameter σ y The calculation formula is: Where: yc is the horizontal convection diffusion parameter; σ ym is the mechanical level diffusion parameter.

7. The method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and meteorological conditions according to claim 1, characterized in that: The calculation formula in step 5 is: C c {x r ,y r ,z r }=C d {x r ,y r ,z r }+C r {x r ,y r ,z r }+C p {x r ,y r ,z r } in, Where: C c is the total concentration of sulfur-containing gas; C d is the concentration of airflow directly diffused to the ground; C r is the mass concentration of the indirect source diffused by the updraft to the top of the mixed layer; Q is the source intensity; C p is the mass concentration of the penetrating source penetrating into the upper stable layer of the mixed layer; u is the wind speed at the top of the leakage source; F y is the horizontal distribution function; σ y is the horizontal diffusion parameter; y is the crosswind distance; z i is the height of the mixing layer; σ zj is the vertical diffusion parameter of the indirect source; σ zp is the vertical diffusion parameter of the penetrating source; dj is the total height of the direct source plume; rj is the total height of the indirect source plume; j is the weight coefficient of the rising and falling parts of the plume; f p is the plume density under convective conditions; h ep is the penetration source height; z ieff is the height of the reflecting surface in the stable layer.

8. A simulation system for the diffusion of blowout in a sulfur-containing gas well, characterized in that: include: Data preprocessing module, which preprocesses the meteorological, topographical, diffusion source and other data in the database or imported from outside into the data required for diffusion simulation calculation; A diffusion simulation calculation module, which calculates the temporal and spatial distribution of the leaked gas based on the input current terrain, weather and leakage conditions and using a method for simulating the diffusion of a sulfur-containing gas well blowout under complex terrain and weather conditions as described in any one of claims 1 to 7; The calculation result processing module processes the calculated temporal and spatial distribution law of the leaked gas, including data extraction, scalar number processing and vector data generation; The visualization module uses the GIS system to perform 2D and 3D visualization of the temporal and spatial distribution of leaked gas, allowing for a clearer and more direct observation of the diffusion results of leaked gas.

9. An electronic device, characterized in that: include: Memory and processor; The memory can be used to store data and executable instructions received by the computer, and the processor is used to execute the executable instructions received by the computer, process data and control other components, which can realize the process of a method for simulating the blowout diffusion of sulfur-containing gas wells under complex terrain and meteorological conditions as described in any one of claims 1-7.

10. A computer-readable storage medium is a physical material or device for storing data or instructions, characterized in that: When the computer executable instructions are executed by the processor, the process of simulating the blowout diffusion of a sulfur-containing gas well under complex terrain and meteorological conditions as described in any one of claims 1 to 7 is realized.

Images