CN114021493A - Complex terrain CFD calculation domain modeling method and geological sequestration carbon dioxide leakage simulation and risk assessment system - Google Patents

Complex terrain CFD calculation domain modeling method and geological sequestration carbon dioxide leakage simulation and risk assessment system Download PDF

Info

Publication number
CN114021493A
CN114021493A CN202111268181.2A CN202111268181A CN114021493A CN 114021493 A CN114021493 A CN 114021493A CN 202111268181 A CN202111268181 A CN 202111268181A CN 114021493 A CN114021493 A CN 114021493A
Authority
CN
China
Prior art keywords
model
dimensional
command
terrain
module
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202111268181.2A
Other languages
Chinese (zh)
Inventor
张一梅
栗帅
李鱼
郭文瑾
郭晓倩
吴百苗
林千果
梁希
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
North China Electric Power University
Original Assignee
North China Electric Power University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by North China Electric Power University filed Critical North China Electric Power University
Priority to CN202111268181.2A priority Critical patent/CN114021493A/en
Publication of CN114021493A publication Critical patent/CN114021493A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/28—Design optimisation, verification or simulation using fluid dynamics, e.g. using Navier-Stokes equations or computational fluid dynamics [CFD]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0635—Risk analysis of enterprise or organisation activities
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00—Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/14—Force analysis or force optimisation, e.g. static or dynamic forces

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Business, Economics & Management (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Human Resources & Organizations (AREA)
  • Evolutionary Computation (AREA)
  • General Engineering & Computer Science (AREA)
  • Strategic Management (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Economics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Computing Systems (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Development Economics (AREA)
  • Fluid Mechanics (AREA)
  • Educational Administration (AREA)
  • Mathematical Physics (AREA)
  • Game Theory and Decision Science (AREA)
  • Algebra (AREA)
  • Marketing (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Tourism & Hospitality (AREA)
  • General Business, Economics & Management (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention belongs to the technical field of carbon peak and carbon neutralization, and particularly relates to a CFD (computational fluid dynamics) computational domain modeling method for complex terrain and geological sequestration CO2A leakage simulation and risk assessment system. According to the modeling method, three-dimensional modeling of complex terrains such as loess hilly ridge regions can be realized simply through data conversion. The simulation and risk assessment system comprises a CO2Exposure concentration acquisition module, CO2An exposure reference concentration module, a health risk evaluation module and a health risk space-time distribution visualization module, a model based on the modelingThe planned risk assessment system has good visualization characteristics, and CCUS potential leakage risk assessment with rapidness, high operability, high precision and high visualization height is realized.

Description

Complex terrain CFD calculation domain modeling method and geological sequestration carbon dioxide leakage simulation and risk assessment system
Technical Field
The invention belongs to the technical field of carbon peak and carbon neutralization, and particularly relates to a CFD (computational fluid dynamics) computational domain modeling method for complex terrain and geological sequestration CO2A leakage simulation and risk assessment system.
Background
Carbon capture, utilization and storage (CCUS) is an economical and effective measure against global climate change, but is vulnerable to public panic and rejection of CCUS due to lack of understanding of the hazards of potential leakage, and is one of the major obstacles to large-scale popularization of CCUS. An oil field in the northwest region, one of the leading demonstration bases of the CCUS in China, belongs to the typical landforms of loess hills and hills, comprehensively describes routes, the range and the hazard level of potential carbon dioxide leakage in a demonstration engineering region, and has important significance for the popularization of the CCUS technology and the realization of the strategic targets of carbon neutralization countries.
69.1% of landforms in China are mountainous areas, oil reservoir storage areas where CCUS pilot demonstration engineering is located mostly belong to the landforms of loess hills and hills, and the complex landforms can obviously influence the diffusion of leaked gas, so that the data uncertainty of a risk assessment concentration source is increased. The risk assessment of the CCUS leakage of the loess hills and the ridge has the following problems: CO 22No threshold data and evaluation model are available; the gas diffusion analog computation domain under complex terrain is difficult to model; the risk characterization approach is poorly understood by the public. Thus developing CO2Health risk evaluation method, complex terrain calculation domain method, risk characterization method closer to actual situation, safe and efficient implementation of CCUS pilot engineering, large-scale popularization and application and carbon-neutral stateThe realization of the strategy goal is of great significance.
The current CFD simulation and risk assessment of gas diffusion for complex terrain/topography has several approaches: (1) importing SURFER software based on the scanning point cloud data and outputting the SURFER software to ANSYS GAMBIT software, and editing point-by-point/line/surface/body in the ANSYS GAMBIT software to generate a mountain terrain simplified to a certain extent; (2) establishing a three-dimensional building model by using SKETCH UP in an auxiliary manner, butting the three-dimensional building model with ICEM to generate a CFD (computational fluid dynamics) calculation domain, and visualizing a diffusion simulation result in an ARCGIS (autoregressive integrated circuit) by adopting a coordinate data conversion method; (3) method for carrying out CCUS sequestration CO by adopting risk matrix method2Semi-quantitative risk assessment of potential leakage and hazard risks.
Although the above way realizes CCUS sequestration of CO to a certain extent2May be CO2The leakage diffusion simulation provides reference, but the preprocessing software GAMBIT software is stopped in a new ANSYS system; the single SKETCH UP software can only draw a small number of buildings or software with low complexity; the risk matrix method is based on the factor weight matrix, can provide risk levels under simple conditions, but is too rich in subjective judgment, limited in risk classification range, incapable of making proper response to actual environmental parameters and extremely low in accuracy. In addition, the methods have not yet realized the CFD three-dimensional calculation domain modeling and the geological sequestration of CO in complex terrains such as loess hills and hills2And (4) simulating leakage and evaluating risk.
Therefore, development of a convenient series of COs2A health risk evaluation method, a complex terrain calculation domain method and a risk characterization method with high visualization and universality are adopted, multiple leakage sources and leakage scenes are identified from the CCUS engineering overall angle system, CFD three-dimensional calculation domain modeling and geological sequestration CO are constructed2The leakage simulation and risk assessment are important measures for solving the current CCUS risk management problem in China.
Disclosure of Invention
The invention aims to provide a method for evaluating and managing the risk of loess hills and hills by utilizing the CCUS in the conventional complex landforms, particularly the problems of complex diffusion simulation modeling, less objective risk evaluation standard, low risk characterization visualization and low popularity in CCUS sealing leakage risk evaluation and managementComplex terrain CFD calculation domain modeling method and geological sequestration CO2According to the leakage simulation and risk assessment system, three-dimensional modeling of complex terrains such as loess hills and hills regions can be achieved simply through data conversion, the modeling and risk assessment system based on the modeling has good visualization characteristics, and CCUS potential leakage risk assessment with high speed, high operability, high precision and high visualization height is achieved.
The technical scheme of the invention is as follows: a complex terrain CFD computational domain modeling method comprises the following steps:
(1) importing digital terrain elevation model image source data of a region to be researched into Global Mapper software, setting a model coordinate system and a projection standard according to the longitude and latitude of the region and a coordinate system of the image source data, and utilizing a drawing command to define the peripheral boundary of the research region; generating a three-dimensional contour model of the area to be researched by utilizing an analysis command and a TIN-based contour generation command; outputting a primary digital graph of the three-dimensional contour line by using a derivation mode and a vector radar format mode;
(2) opening the primary digital graph of the three-dimensional contour line derived in the step (1) by using AUTOCAD software, and viewing and deleting the single points of the unclosed line and the non-highest/lowest elevation by using a three-dimensional view command and a dynamic constraint rotation command; generating a smooth and closed three-dimensional contour line digital graph by utilizing an output command;
(3) leading in the three-dimensional contour line digital graph led out in the step (2) by using SketchUp software, and converting the three-dimensional contour line digital graph into a terrain three-dimensional curved surface model by using a 'sandbox' and 'curved surface fluctuation' command; carrying out simple three-dimensional modeling according to a preset building by adopting Auto CAD software, and deriving a three-dimensional model of the building;
(4) importing the terrain three-dimensional curved surface model and the building three-dimensional model in the step (3) by using Rhinoseros software, and cutting a non-calculation domain part of the three-dimensional model by using a Boolean division command, a Boolean difference set command and a Boolean intersection command to realize seamless butt joint of the building three-dimensional model and the terrain three-dimensional curved surface model so as to generate a real three-dimensional curved surface model containing the terrain and the building; and generating four-directional boundaries of the research area by using a 'suspension line' command, sealing the top of the model by using a 'capping' command to form a hollow 6-surface combined model, and materializing the hollow 6-surface combined model by using an 'entity' command to generate a three-dimensional entity calculation domain model.
The data of the digital terrain elevation model image source of the area to be researched in the step (1) adopts ASTER GDEM remote sensing satellite data V3 version, the spatial resolution of the data is 30m, and the format is tiff; the model coordinate System and the projection standard are ink card tray projection, and a World Geodetic System1984 projection coordinate System and a BeiJing80 coordinate System are adopted; the three-dimensional contour line model is generated based on the same elevation base point.
And (4) drawing a planar wire frame model according to the actual coordinate position of the building in the step (3) and stretching to form the three-dimensional model of the building, wherein the size of the three-dimensional model of the building is consistent with the actual measurement size of the building.
The primary digital graph of the three-dimensional contour line is in a dxf format; the three-dimensional contour line digital graph is in a dwg format; the terrain three-dimensional curved surface model is in a skp format; building a three-dimensional model of the building in a dxf format; the three-dimensional entity calculation domain model is in prasolid format.
Geological sequestration CO based on obtained three-dimensional entity computational domain model2The leakage simulation and risk assessment system comprises a CO2Exposure concentration acquisition module, CO2The system comprises an exposure reference concentration module, a health risk evaluation module and a health risk space-time distribution visualization module; wherein CO is2The exposure concentration acquisition module is used for acquiring CO of each receptor point in different leakage scenes on the basis of the three-dimensional entity calculation domain model2Exposure concentration; CO 22The exposure reference concentration module is based on toxicological analysis for obtaining CO2An exposure concentration reference value; the health risk assessment module is used for acquiring CO of each space point in different leakage scenes2Exposing point cloud data of a health risk coefficient; the health risk space-time distribution visualization module is used for acquiring a highly visualized image integrating terrain, landform, land use and receptor health risk information.
The CO is2The exposure concentration acquisition module is used for firstly calculating the three-dimensional entity into a domainThe model is led into ANSYS ICEM to divide grids, and then the divided grids are led into ANSYS FLUENT to be subjected to gas diffusion simulation, so that the CO of each receptor point in different leakage scenes is obtained2The exposure concentration.
The divided grids adopt hexahedral grids, and space cells in the range of 200m around the building are encrypted by using ANSYS ICEM software local encryption commands; the gas diffusion simulation of ANSYS FLUENT uses a k-epsilon turbulence model.
The CO is2Exposure reference concentration module CO obtained based on toxicology analysis2Exposure concentration baseline values are shown in table 1:
TABLE 1CO2Reference value of exposure concentration
Figure BDA0003327664360000041
The health risk assessment module adopts a quotient method and calculates CO in batches through EXCEL software2Exposing the receptor site CO acquired by the concentration acquisition module2Exposure concentration and CO2Exposing the reference concentration Module resulting CO2The ratio of the exposure concentration reference values is obtained, namely the CO of each space point under different leakage scenes2Exposing point cloud data of the health risk coefficients.
The health risk space-time distribution visualization module uses SURFER software to perform XY coordinates and CO on a receptor point space plane2Performing Critical interpolation on the exposure health risk coefficient to obtain different CO2The receptor health risk spatial distribution map caused by a leakage scene is exported into a kml format file by adopting an export command, a high-degree visual image integrating the information such as terrain, landform, land use, receptor health risk and the like is obtained by opening the kml format file by using GOOGLE EARTH software, the risk distribution details and the spatial distribution trend can be browsed by utilizing a flight command, and further the geological storage CO of the loess hills and hills landform can be researched2Leakage potential hazard risk; wherein the kml file coordinate system is WGS 1984.
The invention has the beneficial effects that: aiming at the problem of difficult modeling of the conventional complex terrain three-dimensional calculation domain, the complex terrain three-dimensional modeling method is established based on four types of general commercial software, i.e. GLOBAL MAPPER, AUTOCAD, SKETCH UP and RHINOCEROS, only 1-2 steps of operation are needed for single software, the complex terrain can be rapidly and accurately modeled through simple graph editing, and the method has high applicability and popularization.
Geological sequestration CO based on three-dimensional entity computational domain model2Leakage simulation and risk assessment system for current CO2The problem that a publicly available exposure concentration threshold value and a risk calculation model are not available for human health risk assessment is solved, and based on toxicological analysis, the CO with different hazard grades is established by adopting an international literature retrieval and data analysis method2Threshold data of exposure concentration and exposure duration, and establishing CO by adopting a health risk assessment accepted quotient method based on the threshold data2Exposing a human health risk assessment mathematical model, judging the risk level of a receptor through simple concentration comparison, having reliable data and simple calculation, and being used for CCUS risk assessment and other CO risk assessment2Assessment of exposure potential.
The assessment system also aims at the problem that the large-area risk visualization of the engineering field adopts a mode of taking an engineering drawing as a base drawing and simply superposing the spatial distribution of risks, so that the risk assessment conclusion of residents is vague, and the panic psychology of the residents to the potential leakage irresponsible cannot be effectively relieved.
In conclusion, the modeling method can realize rapid three-dimensional modeling on complex terrains such as loess hilly ridge regions and the like; quantitative evaluation of the simulation and risk evaluation system setEstimating geological sequestration of CO2The leakage risk and the visual representation of the risk level are integrated, the rapid, high-operability, high-precision and high-visualization CCUS potential leakage risk assessment is realized, the engineering design and the risk prevention and control are facilitated, the identification and cognition degree of residents on the potential risk can be improved, the implementation of carbon neutralization measures in China and the realization of double-carbon targets are greatly promoted, and the method has high applicability and popularization.
Drawings
FIG. 1 is an overall framework for module and model development involved in an embodiment of the present invention.
Fig. 2 is a diagram of a topographic image data map box in the sealing area and a predetermined range around the sealing area according to an embodiment of the present invention.
Fig. 3 is a smooth closed dwg format three-dimensional contour digital image according to an embodiment of the invention.
FIG. 4 is a skp format three-dimensional surface model according to an embodiment of the present invention.
FIG. 5 is a three-dimensional surface model of a terrain or structure according to an embodiment of the present invention.
FIG. 6 is an ANSYS FLUENT receptor site CO of an embodiment of the present invention2The density image is exposed.
FIG. 7 is a spatiotemporal visualization image of risk index in GOOGLE EARTH according to an embodiment of the present invention.
Detailed Description
The present invention will be described in detail below with reference to the accompanying drawings.
The complex terrain CFD calculation domain modeling method comprises the following steps:
(1) and acquiring terrain image data in the sealed area and a nearby preset range, wherein the terrain image source data is ASTER GDEM remote sensing satellite data V3 version, the spatial resolution of the data is 30m, and the format is tiff, and referring to fig. 2. Importing the terrain image source data into Global Mapper (Blue Marble Geographic Inc.v.19) software, setting a model coordinate system and a projection standard according to the longitude and latitude of a region and a coordinate system of the image source data, and using a drawing command to define the peripheral boundary of a research area; generating a three-dimensional contour model of the area to be researched by utilizing an analysis command and a TIN-based contour generation command; outputting a primary digital graph of the three-dimensional contour line by using a derivation mode and a vector radar format mode; the model coordinate System and the projection standard are ink card tray projection, and a World Geodetic System1984 projection coordinate System and a BeiJing80 coordinate System are adopted; the three-dimensional contour line model is generated based on the same elevation base point.
(2) Opening the primary digital graph of the three-dimensional contour line in the dxf format by using AUTOCAD software, and viewing and deleting single points of unclosed lines and non-highest/lowest elevations by using a three-dimensional view command and a dynamic constraint rotation command; a smooth, closed three-dimensional contour digital graph is generated using the "output" command, see fig. 3.
(3) Importing the three-dimensional contour digital graph in the dwg format by using SketchUp software (TrimbleInc.v.2018) in the step (2), and converting the three-dimensional contour digital graph into a terrain three-dimensional curved surface model in the skp format by using a 'sandbox' and 'curved surface fluctuation' command, referring to FIG. 4; and carrying out simple three-dimensional modeling according to a preset building by adopting Auto CAD software, and deriving a building three-dimensional model in a dxf format, wherein the building three-dimensional model is formed by drawing a planar wire frame model according to the actual coordinate position of the building and stretching, and the size of the building three-dimensional model is consistent with the actual measurement size of the building.
(4) Introducing a terrain three-dimensional curved surface model in the skp format and a building three-dimensional solid model in the dxf format in the step (3) by using Rhinoseros (v6.0) software, and cutting a non-calculation domain part of the three-dimensional solid model by using a Boolean division command, a Boolean difference set command and a Boolean intersection command to realize seamless butt joint of the building three-dimensional solid model and the terrain three-dimensional curved surface model so as to generate a real three-dimensional curved surface model containing the terrain and the building, and referring to FIG. 5; and generating four-directional boundaries of the research area by using a 'suspension line' command, closing the top of the model by using a 'capping' command to form a hollow 6-surface combined model, and materializing the hollow 6-surface combined model by using an 'entity' command to generate a three-dimensional entity computational domain model in a parasolid format.
Geological sequestration CO based on obtained three-dimensional entity computational domain model2A leak simulation and risk assessment system comprising a CO2ExposingConcentration acquisition module, CO2The system comprises an exposure reference concentration module, a health risk evaluation module and a health risk space-time distribution visualization module.
The CO is2An exposure concentration acquisition module: firstly, introducing the three-dimensional entity computational domain model in the parasolid format into ANSYS ICEM to divide grids, and then introducing the divided grids into ANSYS FLUENT to perform gas diffusion simulation, thereby obtaining CO of each receptor point in different leakage scenes2Exposure concentration, see fig. 6. The division grid adopts hexahedral grid, and ANSYS ICEM software 'local encryption' command is utilized to encrypt the 200m range space cells around the building; the gas diffusion simulation of ANSYS FLUENT uses a k-epsilon turbulence model.
The CO is2An exposure reference concentration module: obtaining CO based on literature experimental data arrangement and toxicology analysis2Exposure concentration baseline values.
The experimental data of the literature is data published in authoritative periodicals at home and abroad during 2001-2020, and the lists of the references directly cited are shown in the following tables 2 and 3:
TABLE 2 CO2Exposure concentration baseline and reference
Figure BDA0003327664360000071
Table 3 reference details
Figure BDA0003327664360000072
Figure BDA0003327664360000081
The health risk assessment module: calculating CO in batches by EXCEL software by using a quotient method2Exposing the receptor site CO acquired by the concentration acquisition module2Exposure concentration and CO2Exposing the reference concentration Module resulting CO2Ratio of reference values of exposure concentration, i.e. obtainingObtaining CO of each space point under different leakage scenes2Exposing point cloud data of the health risk coefficients.
The health risk spatial-temporal distribution visualization module: using SURFER software to carry out X-Y coordinate and CO on the space plane of the receptor point2Performing Critical interpolation on the exposure health risk coefficient to obtain different CO2The receptor health risk spatial distribution map caused by a leakage scene is exported into a kml format file by adopting an export command, a high-degree visual image integrating the information such as terrain, landform, land use, receptor health risk and the like is obtained by opening the kml format file by using GOOGLE EARTH software, the risk distribution details and the spatial distribution trend can be browsed by utilizing a flight command, and further the geological storage CO of the loess hills and hills landform can be researched2Leakage potential hazard risk, see fig. 7; wherein the kml file coordinate system is WGS 1984.

Claims (10)

1. A CFD computational domain modeling method for complex terrain is characterized by comprising the following steps:
(1) importing digital terrain elevation model image source data of a region to be researched into Global Mapper software, setting a model coordinate system and a projection standard according to the longitude and latitude of the region and a coordinate system of the image source data, and utilizing a drawing command to define the peripheral boundary of the research region; generating a three-dimensional contour model of the area to be researched by utilizing an analysis command and a TIN-based contour generation command; outputting a primary digital graph of the three-dimensional contour line by using a derivation mode and a vector radar format mode;
(2) opening the primary digital graph of the three-dimensional contour line derived in the step (1) by using AUTOCAD software, and viewing and deleting the single points of the unclosed line and the non-highest/lowest elevation by using a three-dimensional view command and a dynamic constraint rotation command; generating a smooth and closed three-dimensional contour line digital graph by utilizing an output command;
(3) leading in the three-dimensional contour line digital graph led out in the step (2) by using SketchUp software, and converting the three-dimensional contour line digital graph into a terrain three-dimensional curved surface model by using a 'sandbox' and 'curved surface fluctuation' command; carrying out simple three-dimensional modeling according to a preset building by adopting Auto CAD software, and deriving a three-dimensional model of the building;
(4) importing the terrain three-dimensional curved surface model and the building three-dimensional model in the step (3) by using Rhinoseros software, and cutting a non-calculation domain part of the three-dimensional model by using a Boolean division command, a Boolean difference set command and a Boolean intersection command to realize seamless butt joint of the building three-dimensional model and the terrain three-dimensional curved surface model so as to generate a real three-dimensional curved surface model containing the terrain and the building; and generating four-directional boundaries of the research area by using a 'suspension line' command, sealing the top of the model by using a 'capping' command to form a hollow 6-surface combined model, and materializing the hollow 6-surface combined model by using an 'entity' command to generate a three-dimensional entity calculation domain model.
2. The complex terrain CFD computational domain modeling method of claim 1, wherein in the step (1), the digital terrain elevation model image source data of the area to be researched adopts ASTER GDEM remote sensing satellite data V3 version, the data has a spatial resolution of 30m and a format of tiff; the model coordinate System and the projection standard are ink card tray projection, and a World Geodetic System1984 projection coordinate System and a BeiJing80 coordinate System are adopted; the three-dimensional contour line model is generated based on the same elevation base point.
3. The complex terrain CFD computational domain modeling method according to claim 1, wherein the building three-dimensional stereo model in the step (3) is formed by drawing a planar wire frame model according to the actual coordinate position of the building and stretching, and the dimension of the model is consistent with the measured dimension of the building.
4. The complex terrain CFD computational domain modeling method of claim 1, wherein the primary digital graph of the three-dimensional contour is in a dxf format; the three-dimensional contour line digital graph is in a dwg format; the terrain three-dimensional curved surface model is in a skp format; building a three-dimensional model of the building in a dxf format; the three-dimensional entity calculation domain model is in prasolid format.
5. Geological sequestration CO based on the three-dimensional entity computational domain model obtained in claim 12A leak simulation and risk assessment system, characterized in that the system comprises a CO2Exposure concentration acquisition module, CO2The system comprises an exposure reference concentration module, a health risk evaluation module and a health risk space-time distribution visualization module; wherein CO is2The exposure concentration acquisition module is used for acquiring CO of each receptor point in different leakage scenes on the basis of the three-dimensional entity calculation domain model2Exposure concentration; CO 22The exposure reference concentration module is based on toxicological analysis for obtaining CO2An exposure concentration reference value; the health risk assessment module is used for acquiring CO of each space point in different leakage scenes2Exposing point cloud data of a health risk coefficient; the health risk space-time distribution visualization module is used for acquiring a highly visualized image integrating terrain, landform, land use and receptor health risk information.
6. Geological sequestration of CO according to claim 52A leak simulation and risk assessment system, wherein the CO2The exposure concentration acquisition module is used for guiding the three-dimensional entity computational domain model into ANSYS ICEM to divide grids, and then guiding the divided grids into ANSYS FLUENT to perform gas diffusion simulation, so that CO of each receptor point in different leakage scenes is acquired2The exposure concentration.
7. Geological sequestration of CO according to claim 62The leakage simulation and risk assessment system is characterized in that the divided grids adopt hexahedral grids, and space cells in a range of 200m around a building are encrypted by using ANSYS ICEM software local encryption commands; the gas diffusion simulation of ANSYS FLUENT uses a k-epsilon turbulence model.
8. Geological sequestration of CO according to claim 52A leak simulation and risk assessment system, wherein the CO2Exposure reference concentration module CO obtained based on toxicology analysis2The exposure concentration reference values are shown in Table 1：
TABLE 1CO2Reference value of exposure concentration
Figure FDA0003327664350000021
9. Geological sequestration of CO according to claim 52The leakage simulation and risk assessment system is characterized in that the health risk assessment module calculates CO in batches by EXCEL software by adopting a quotient method2Exposing the receptor site CO acquired by the concentration acquisition module2Exposure concentration and CO2Exposing the reference concentration Module resulting CO2The ratio of the exposure concentration reference values is obtained, namely the CO of each space point under different leakage scenes2Exposing point cloud data of the health risk coefficients.
10. Geological sequestration of CO according to claim 52The leakage simulation and risk assessment system is characterized in that the health risk space-time distribution visualization module uses SURFER software to carry out XY coordinates and CO on a spatial plane of a receptor point2Performing Critical interpolation on the exposure health risk coefficient to obtain different CO2The receptor health risk spatial distribution map caused by a leakage scene is exported into a kml format file by adopting an export command, a high-degree visual image integrating the information such as terrain, landform, land use, receptor health risk and the like is obtained by opening the kml format file by using GOOGLE EARTH software, the risk distribution details and the spatial distribution trend can be browsed by utilizing a flight command, and further the geological storage CO of the loess hills and hills landform can be researched2Leakage potential hazard risk; wherein the kml file coordinate system is WGS 1984.
CN202111268181.2A 2021-10-29 2021-10-29 Complex terrain CFD calculation domain modeling method and geological sequestration carbon dioxide leakage simulation and risk assessment system Pending CN114021493A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111268181.2A CN114021493A (en) 2021-10-29 2021-10-29 Complex terrain CFD calculation domain modeling method and geological sequestration carbon dioxide leakage simulation and risk assessment system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111268181.2A CN114021493A (en) 2021-10-29 2021-10-29 Complex terrain CFD calculation domain modeling method and geological sequestration carbon dioxide leakage simulation and risk assessment system

Publications (1)

Publication Number Publication Date
CN114021493A true CN114021493A (en) 2022-02-08

Family

ID=80058519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111268181.2A Pending CN114021493A (en) 2021-10-29 2021-10-29 Complex terrain CFD calculation domain modeling method and geological sequestration carbon dioxide leakage simulation and risk assessment system

Country Status (1)

Country Link
CN (1) CN114021493A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114595648A (en) * 2022-05-11 2022-06-07 浙江浙能航天氢能技术有限公司 Design method and system for inert gas explosion suppression scheme of hydrogenation station
CN115994496A (en) * 2022-12-02 2023-04-21 上海师范大学 Urban park high resolution atmospheric CO 2 Numerical simulation method of concentration three-dimensional field

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114595648A (en) * 2022-05-11 2022-06-07 浙江浙能航天氢能技术有限公司 Design method and system for inert gas explosion suppression scheme of hydrogenation station
CN115994496A (en) * 2022-12-02 2023-04-21 上海师范大学 Urban park high resolution atmospheric CO 2 Numerical simulation method of concentration three-dimensional field
CN115994496B (en) * 2022-12-02 2024-03-12 上海师范大学 Urban park high resolution atmospheric CO 2 Numerical simulation method of concentration three-dimensional field

Similar Documents

Publication Publication Date Title
CN114021493A (en) Complex terrain CFD calculation domain modeling method and geological sequestration carbon dioxide leakage simulation and risk assessment system
CN104361616B (en) A kind of topography and geomorphology file acquisition method assessed for wind energy turbine set wind-resources
CN108763825B (en) Numerical simulation method for simulating wind field of complex terrain
CN108009314B (en) Archaeological excavation data-based three-dimensional reconstruction method for early topography of site area
Yan et al. Integration of 3D objects and terrain for 3D modelling supporting the digital twin
CN104766366A (en) Method for establishing three-dimensional virtual reality demonstration
CN102693328A (en) Application of ARCVIEW GIS (geographic information system) 3.2 technology in mountainous land building design
CN105261067A (en) Overground and underground integrated modeling method based on true 3D volumetric display technique and system
Khayyal et al. Creation and spatial analysis of 3D city modeling based on GIS data
Malinverni et al. GIS-based smart cartography using 3D modeling
Urech Point-cloud modeling: Exploring a site-specific approach for landscape design
Xu et al. Methods for the construction of DEMs of artificial slopes considering morphological features and semantic information
Araújo et al. An approach to simulate wind fields around an urban environment for wind energy application
Stendardo et al. Urban design strategies for the upcycling of urban infrastructure residual pockets: 3D city modelling from Open Data and low-cost rapid mapping tools
CN103714494B (en) Power circuit method for building up and device
Koukofikis et al. An integration of urban spatial data with energy simulation to produce X3D city models: the case of Landkreis Ludwigsburg
Lu et al. Computer aided design system based on 3D GIS for park design
Prasertsuk et al. 3D Data Visualization and Analysis Tools for AI Ready City: Space Syntax and Social Media Data
Trubina et al. Environmental assessment of urban areas: approaches based on geospatial data integration
Liu et al. An improved LOS method for implementing visibility analysis of 3D complex landscapes
Ting et al. Application of ground-based 3D laser scanning technology in engineering surveying
CN107833278A (en) Terrain simulation method, apparatus and electronic equipment
Hao et al. Accuracy analysis of earthwork calculation based on triangulated irregular network (TIN)
CN117874893B (en) Chemical pipeline visualization method based on BIM and UE5
Dimen et al. 3D Modeling of the Historical Buildings and Archeological Sites

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination