CN108170937B - 自然风沙冲蚀环境建模方法 - Google Patents
自然风沙冲蚀环境建模方法 Download PDFInfo
- Publication number
- CN108170937B CN108170937B CN201711428336.8A CN201711428336A CN108170937B CN 108170937 B CN108170937 B CN 108170937B CN 201711428336 A CN201711428336 A CN 201711428336A CN 108170937 B CN108170937 B CN 108170937B
- Authority
- CN
- China
- Prior art keywords
- sand
- wind
- dust
- erosion
- data
- 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.)
- Active
Links
- 239000004576 sand Substances 0.000 title claims abstract description 60
- 230000003628 erosive effect Effects 0.000 title claims abstract description 56
- 238000000034 method Methods 0.000 title claims abstract description 8
- 239000000428 dust Substances 0.000 claims abstract description 44
- 238000004088 simulation Methods 0.000 claims abstract description 15
- 230000004907 flux Effects 0.000 claims abstract description 11
- 230000008878 coupling Effects 0.000 claims abstract description 6
- 238000010168 coupling process Methods 0.000 claims abstract description 6
- 238000005859 coupling reaction Methods 0.000 claims abstract description 6
- 238000007689 inspection Methods 0.000 claims description 9
- 238000009434 installation Methods 0.000 claims description 8
- 230000003068 static effect Effects 0.000 claims description 8
- 239000002689 soil Substances 0.000 claims description 7
- 239000000443 aerosol Substances 0.000 claims description 6
- 239000002245 particle Substances 0.000 claims description 4
- 238000012805 post-processing Methods 0.000 claims description 4
- 238000003908 quality control method Methods 0.000 claims description 3
- 239000000126 substance Substances 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 claims description 2
- 238000004062 sedimentation Methods 0.000 claims description 2
- 238000007619 statistical method Methods 0.000 claims description 2
- 239000000463 material Substances 0.000 abstract description 13
- 238000005516 engineering process Methods 0.000 abstract description 10
- 238000011160 research Methods 0.000 abstract description 6
- 238000012360 testing method Methods 0.000 abstract description 5
- 230000007613 environmental effect Effects 0.000 abstract description 4
- 238000013461 design Methods 0.000 abstract description 2
- 238000012423 maintenance Methods 0.000 abstract description 2
- 238000007667 floating Methods 0.000 description 3
- 230000004580 weight loss Effects 0.000 description 3
- 230000006399 behavior Effects 0.000 description 2
- 238000000576 coating method Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 229910021532 Calcite Inorganic materials 0.000 description 1
- 238000005299 abrasion Methods 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910001919 chlorite Inorganic materials 0.000 description 1
- 229910052619 chlorite group Inorganic materials 0.000 description 1
- QBWCMBCROVPCKQ-UHFFFAOYSA-N chlorous acid Chemical compound OCl=O QBWCMBCROVPCKQ-UHFFFAOYSA-N 0.000 description 1
- 239000011248 coating agent Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- -1 machinery Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 229910052751 metal Inorganic materials 0.000 description 1
- 239000002184 metal Substances 0.000 description 1
- 150000002739 metals Chemical class 0.000 description 1
- 239000010445 mica Substances 0.000 description 1
- 229910052618 mica group Inorganic materials 0.000 description 1
- 229910052655 plagioclase feldspar Inorganic materials 0.000 description 1
- 239000010453 quartz Substances 0.000 description 1
- VYPSYNLAJGMNEJ-UHFFFAOYSA-N silicon dioxide Inorganic materials O=[Si]=O VYPSYNLAJGMNEJ-UHFFFAOYSA-N 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Devices Affording Protection Of Roads Or Walls For Sound Insulation (AREA)
- Testing Resistance To Weather, Investigating Materials By Mechanical Methods (AREA)
Abstract
本发明适用于环境试验领域,提供了一种自然风沙冲蚀环境建模技术,包括:获取气候状态风向和风速资料;基于区域气候‑沙尘耦合模式RegCM4‑Dust,获取沙尘混合比模拟资料;计算沙尘浓度;计算风沙冲蚀通量;计算风向与研究对象的夹角,即风沙冲蚀角度。沙尘冲蚀通量和风沙冲蚀角度构成所建自然风沙冲蚀环境。本发明提供的自然风沙冲蚀环境建模技术使用数值模拟方法获取沙尘浓度数据,避免布设观测仪器,约了设备、维护、人力等大量成本,节省了观测时间,通过与实验室沙尘试验对比,可用于考核和评价材料对自然风沙冲蚀环境的适应性和耐久性,并可为抗风沙冲蚀材料设计、应用和防护提供仿真环境支撑,填补了自然风沙冲蚀环境建模仿真技术的空白。
Description
技术领域
本发明属于环境试验领域,涉及一种自然风沙冲蚀环境的建模技术。
背景技术
我国西北地区是沙尘活动高发区域,其中高风沙环境不仅对空气质量具有重要威胁,对材料的腐蚀磨损失效行为也非常严重,导致耐久性和安全性下降,并造成材料损失和设备破坏,直接影响到各种金属、机械、车辆等的使用寿命,给国民经济造成了巨大损失。沙漠地区的尘土包括砂土、浮土等,其中砂土颗粒较大,在高风速带动下,对暴露在外的各种材料产生冲击破坏作用。浮土比砂土更细更黏,主要成分是高硬度的石英和斜长石,以及方解石、绿泥石、云母等。浮土地带各种材料周围空气的含尘浓度极大,距地面越高空气含尘浓度越小而尘土中细尘所占比例越大,对车辆和工程机械的磨损更为严重。
不同风沙流速和角度对材料的破坏作用不同。例如,对于涂层,研究表明,涂层冲蚀失重量随着风沙流冲蚀速度的增大而显著增加,随着冲蚀下沙量的增加而增加,其在低角度冲蚀下的冲蚀失重量要大于高角度下的冲蚀失重量。目前,有关自然风沙环境的研究主要集中在单独的风场分布研究或单独的沙尘环境研究,将二者结合在一起对自然环境中风沙冲蚀环境的建模技术尚不多见。在环境试验领域,实验室中常用的气流挟沙喷射等方法建立的风沙环境,也仅用于研究对材料的冲蚀磨损机理等,无法真实地反映自然环境下风沙环境对材料的冲蚀行为。
发明内容
为了解决上述存在的问题,本发明的目的是,采用观测资料和数值模拟相结合技术,提供一种自然风沙冲蚀环境建模技术,旨在解决自然风沙冲蚀环境建模问题。
为达到上述目的,本发明采用的技术方案如下:自然风沙冲蚀环境建模技术,包括如下步骤:
步骤一,获取气候状态风向和风速资料:首先获取目标地区地面逐时观测资料,并对资料进行质量控制,包括极值检查、要素场一致性检查、连续性检查等,剔除缺测和不合理数据,选定最终可用资料样本;利用均值对风场资料进行统计分析,获取该地区气候状态风向和风速分布;
气候状态风向的计算为
气候状态风速的计算为
步骤二,获取沙尘混合比资料:在Linux操作系统下安装意大利国际理论物理中心研发的区域气候-沙尘耦合模式RegCM4-Dust,包括安装运行模式所需pgi、openmpi、netcdf等库函数并设置环境变量,采用./configure命令确认模式安装环境,以及采用make all命令建立模式运行环境;获取地面静态资料和模式驱动场资料,为模式提供初始和侧边界条件,地面静态资料包括地形、地表植被类型和土壤类型,模式驱动场资料包括海表温度资料、大气和陆表温度资料;设置regcm.in文件参数,包括模拟时段、三维格点数、模拟区域、地面静态数据参数、缓冲区、物理参数化方案,以及模拟沙尘环境需要的化学模块参数,包括是否考虑气溶胶反馈、干/湿移除、气溶胶类型、水可溶性、干沉降速率、粒径尺寸等;运行RegCM4-Dust模式,包括依次运行地形模块terrain、初始和侧边界模块icbc和模式主模块regcm;运行后处理模块postproc,获取月平均沙尘混合比资料;
步骤三,计算沙尘浓度,具体为:
其中,Cdust为沙尘浓度,Trac为模式输出沙尘混合比,Ts为模式输出地面热力学温度,ps为模式输出地面气压,R为干空气常数,R=287J/(K·kg);
步骤四,计算风沙冲蚀通量,具体为:
其中:mj为第j个风向风沙冲蚀通量;
步骤五,计算风向与冲蚀对象的夹角,即风沙冲蚀角度,具体为:
当冲蚀对象呈东南-西北方向放置时,则
当冲蚀对象呈东北-西南方向放置时,则
其中:θ为冲蚀角度(单位:°),α为冲蚀对象放置方向与东西方向夹角(单位:°,0≤α≤90),β为风向(单位:°,0≤β<360)。
沙尘冲蚀通量和风沙冲蚀角度构成所建自然风沙冲蚀环境。
本发明提供的自然风沙冲蚀环境建模技术,具有以下有益效果:
(1)填补了自然风沙冲蚀环境建模技术的空白;
(2)使用数值模拟方法获取沙尘浓度数据,避免布设观测仪器,节约了设备、维护、人力等大量成本,节省了观测时间;
(3)通过与实验室沙尘试验对比,可用于考核和评价材料对自然风沙冲蚀环境的适应性和耐久性;
(4)可为抗风沙冲蚀材料设计、应用和防护提供仿真环境支撑。
附图说明
以下结合附图和具体实施方式对本发明作进一步的说明。
图1为兰州除静风外的风向频率图(单位:%);
图2为兰州不同风向的平均风速图(单位:m/s);
图3为兰州月平均沙尘浓度(单位:μg/m3);
图4为兰州风沙冲蚀通量(单位:μg/(m2·s))。
具体实施方式
为了使本发明的技术方案更加清楚明白,以下结合附图及实施例,对本发明作进一步详细说明。应当理解,此处所描述的具体实施例仅仅用于解释本发明,并不用于限定本发明。
以下结合附图对本发明实施例作进一步详细说明。
步骤一,获取气候状态风向和风速资料。以兰州地区为例,从美国国家海洋和大气管理局(NOAA)获取1990-2009年地面逐3小时风向和风速观测资料,采用极值检查、要素场一致性检查、连续性检查等方法对资料进行质量控制,剔除缺测和不合理数据;按照风向间隔10°,利用公式(1)进行统计,其中静风(风速为0)出现频率为59.3%,其他风向出现频率见图1,利用公式(2)计算的平均风速见图2。
步骤二,获取沙尘混合比资料
(1)在Linux操作系统(RHEL(Red Hat Enterprise Linux)x86_64 V6.3版本)下安装意大利国际理论物理中心研发的区域气候-沙尘耦合模式RegCM4.1-Dust:
安装库函数及版本号分别为pgi V10.2、openmpi V1.4.3、netcdf V4.0.1,并设置环境变量;
采用./configure命令确认模式安装环境,包括安装目录、Bin目录、是否采用并行模式、是否安装BAND选项、是否安装CLM选项、编译器类型选择等;
采用make all命令建立模式运行环境,包括地形模块terrain、初始和侧边界模块icbc、模式主模块regcm和后处理模块postproc。
(2)获取地面静态资料和模式驱动场资料,其中地面静态资料采用全球陆地覆盖特征(GLCC)数据库,海温资料为NCEPReynolds周平均海温,大气和陆表温度资料为2.5°×2.5°NCEP再分析资料。
(3)设置regcm.in文件参数,模拟时段为1989年3月1日至2010年2月28日,其中第一年(1989年3月-1990年2月)作为模式的spin-up时间不作分析,模拟区域中心位于(112°E,37.5°N),格点水平分辨率为60km,东西向148个格点,南北向96个格点,模式垂直方向分18层,模式顶气压为50hPa,时间积分步长200s,缓冲区为12圈,物理参数化方案分别采用BATSle陆面方案、指数松弛侧边界方案、Holtslag行星边界层方案、混合对流方案、显式水汽方案和Zeng海洋通量方案,化学模块参数考虑气溶胶及其气候反馈(idirect=2)、粒径尺寸包括0.01~1.0μm、1.0~2.5μm、2.5~5.0μm和5.0~20.0μm。
(4)依次运行地形模块terrain、初始和侧边界模块icbc和模式主模块regcm,得到初始模拟结果。
(5)运行后处理模块postproc,得到月平均沙尘混合比资料。
步骤三,利用公式(3)计算沙尘浓度,得到兰州1990年3月至2009年2月1~12月月平均沙尘浓度,见图3,并进一步得到1990-2009年平均沙尘浓度为69.35μg/m3。
步骤四,利用公式(4),结合步骤一结果,计算得到兰州地区风沙冲蚀通量见图4。
步骤五,计算风向与材料的夹角,即风沙冲蚀角度,假设材料呈东南-西北方向放置,与东西方向夹角为25°,结合步骤一风向间隔(10°),利用公式(5)计算得到风沙冲蚀角度为
Claims (1)
1.自然风沙冲蚀环境建模方法,其特征在于,包括如下步骤:
步骤一,获取气候状态风向和风速资料:首先获取目标地区地面逐时观测资料,并对资料进行质量控制,包括极值检查、要素场一致性检查、连续性检查等,剔除缺测和不合理数据,选定最终可用资料样本;利用均值对风场资料进行统计分析,获取该地区气候状态风向和风速分布;
气候状态风向的计算为
气候状态风速的计算为
步骤二,获取沙尘混合比资料:在Linux操作系统下安装区域气候-沙尘耦合模式RegCM4-Dust,包括安装运行模式所需pgi、openmpi、netcdf等库函数并设置环境变量,采用./configure命令确认模式安装环境,以及采用make all命令建立模式运行环境;获取地面静态资料和模式驱动场资料,为模式提供初始和侧边界条件,地面静态资料包括地形、地表植被类型和土壤类型,模式驱动场资料包括海表温度资料、大气和陆表温度资料;设置regcm.in文件参数,包括模拟时段、三维格点数、模拟区域、地面静态数据参数、缓冲区、物理参数化方案,以及模拟沙尘环境需要的化学模块参数,包括是否考虑气溶胶反馈、干/湿移除、气溶胶类型、水可溶性、干沉降速率、粒径尺寸等;运行RegCM4-Dust模式,包括依次运行地形模块terrain、初始和侧边界模块icbc和模式主模块regcm;运行后处理模块postproc,获取月平均沙尘混合比资料;
步骤三,计算沙尘浓度,具体为:
其中,Cdust为沙尘浓度,Trac为模式输出沙尘混合比,Ts为模式输出地面热力学温度,ps为模式输出地面气压,R为干空气常数,R=287J/(K·kg);
步骤四,计算风沙冲蚀通量,具体为:
其中:mj为第j个风向风沙冲蚀通量;
步骤五,计算风向与冲蚀对象的夹角,即风沙冲蚀角度,具体为:
当冲蚀对象呈东南-西北方向放置时,则
当冲蚀对象呈东北-西南方向放置时,则
其中:θ为冲蚀角度(单位:°),α为冲蚀对象与东西方向夹角(单位:°,0≤α≤90),β为风向(单位:°,0≤β<360);
沙尘冲蚀通量和风沙冲蚀角度构成所建自然风沙冲蚀环境。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711428336.8A CN108170937B (zh) | 2017-12-26 | 2017-12-26 | 自然风沙冲蚀环境建模方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711428336.8A CN108170937B (zh) | 2017-12-26 | 2017-12-26 | 自然风沙冲蚀环境建模方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108170937A CN108170937A (zh) | 2018-06-15 |
CN108170937B true CN108170937B (zh) | 2021-03-05 |
Family
ID=62520719
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711428336.8A Active CN108170937B (zh) | 2017-12-26 | 2017-12-26 | 自然风沙冲蚀环境建模方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108170937B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113902603B (zh) * | 2021-10-13 | 2022-08-09 | 中科三清科技有限公司 | 沙尘排放通量的计算方法和装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015160927A1 (en) * | 2014-04-15 | 2015-10-22 | The University Of Akron | Methods for evaluation and estimation of external corrosion damage on buried pipelines |
CN105868501A (zh) * | 2016-04-21 | 2016-08-17 | 湘潭大学 | 热障涂层冲蚀率模型及含涂层涡轮叶片冲蚀工况模拟方法 |
CN106153290A (zh) * | 2016-08-30 | 2016-11-23 | 敦煌研究院 | 一种防水式四方位沙尘水平收集与垂直降尘采集仪 |
CN106768821A (zh) * | 2017-01-22 | 2017-05-31 | 北京师范大学 | 一种地块内土壤风蚀量的测量方法 |
CN106982100A (zh) * | 2017-03-31 | 2017-07-25 | 北京润科通用技术有限公司 | 一种戈壁滩环境空地信道建模方法及装置 |
-
2017
- 2017-12-26 CN CN201711428336.8A patent/CN108170937B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2015160927A1 (en) * | 2014-04-15 | 2015-10-22 | The University Of Akron | Methods for evaluation and estimation of external corrosion damage on buried pipelines |
CN105868501A (zh) * | 2016-04-21 | 2016-08-17 | 湘潭大学 | 热障涂层冲蚀率模型及含涂层涡轮叶片冲蚀工况模拟方法 |
CN106153290A (zh) * | 2016-08-30 | 2016-11-23 | 敦煌研究院 | 一种防水式四方位沙尘水平收集与垂直降尘采集仪 |
CN106768821A (zh) * | 2017-01-22 | 2017-05-31 | 北京师范大学 | 一种地块内土壤风蚀量的测量方法 |
CN106982100A (zh) * | 2017-03-31 | 2017-07-25 | 北京润科通用技术有限公司 | 一种戈壁滩环境空地信道建模方法及装置 |
Non-Patent Citations (2)
Title |
---|
"Impact of desert dust on the east Asian summer monsoon";宿兴涛等;《热带气象学报(英文版)》;20160229;page252-264 * |
"RegCm4模式对云南及周边地区干旱化趋势的预估";王美丽等;《高原气象》;20150630;第34卷(第3期);第707-714页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108170937A (zh) | 2018-06-15 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Bowman | Large‐scale isentropic mixing properties of the Antarctic polar vortex from analyzed winds | |
Woodward | Modeling the atmospheric life cycle and radiative impact of mineral dust in the Hadley Centre climate model | |
Xue et al. | The dispersion of silver iodide particles from ground-based generators over complex terrain. Part II: WRF large-eddy simulations versus observations | |
CN105067975A (zh) | 一种分析雾霾对输变电设备外绝缘影响的试验方法 | |
Cramer | 14 Using plant functional types in a global vegetation model | |
Jin et al. | Diagnostic analysis of wintertime PM2. 5 pollution in the North China Plain: The impacts of regional transport and atmospheric boundary layer variation | |
CN108170937B (zh) | 自然风沙冲蚀环境建模方法 | |
Woodage et al. | UK HiGEM: simulations of desert dust and biomass burning aerosols with a high-resolution atmospheric GCM | |
CN112326474A (zh) | 腐蚀-疲劳协同加载的寿命加速试验方法 | |
Park et al. | Analyses of high aerosol concentration events (dense haze/mist) occurred in East Asia during 10-16 January 2013 using the data simulated by the Aerosol Modeling System | |
Coe et al. | The water balance of northern Africa during the mid-Holocene: an evaluation of the 6 ka BP PMIP simulations | |
Godowitch et al. | Investigating the impact on modeled ozone concentrations using meteorological fields from WRF with an updated four–dimensional data assimilation approach | |
Wang et al. | Influence of meteorological reanalysis field on air quality modeling in the Yangtze River Delta, China | |
Diop et al. | The use of ALADIN model and MERRA-2 reanalysis to represent dust seasonal dry deposition from 2006 to 2010 in Senegal, West Africa | |
Rungratanaubon et al. | Vertical variation of nitrogen oxide (NOx) concentration using a backward air mass trajectories model in an urban area of Bangkok, Thailand | |
Navascues | Analysis of 2 meter temperature and relative humidity | |
Gantner et al. | The impact of land-surface parameter properties and resolution on the simulated cloud-topped atmospheric boundary layer | |
Chen et al. | Application of LIDAR technique and MM5-CMAQ modeling approach for the assessment of winter PM 10 air pollution: a case study in Beijing, China | |
Yuwono et al. | Generation of total suspended particulate (TSP) in ambient air from four soil types in Indonesia | |
Fratini et al. | Evaluation of a wind erosion model in a desert area of northern Asia by eddy covariance | |
Weathers et al. | Deposition | |
Vongmahadlek et al. | Applicability of RAMS for a simulation to provide inputs to an air quality model: Modeling evaluation and sensitivity test | |
Prakash | A new physics-based model for guiding bat carcass surveys at wind turbines | |
Omar et al. | Spatial Analysis of Predicting Soil Physicochemical Properties Using Satellite Images: A Case Study in Erbil Governorate | |
Varble et al. | Cloud, Aerosol, and Complex Terrain Interactions (CACTI) Preliminary Science Plan |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |