基于地层信息插值的FLAC3D简易建模方法

为建立一种简易的FLAC3D建模方法,降低FLAC3D前处理时的建模难度,本文基于Python编程语言开发了建模文件生成程序。该程序以岩体各层位的高程信息为基础,采用Kriging插值生成离散网格点数据;并采用离散网格点自动生成建模所需要的网格数据文件;将该文件导入FLAC3D软件中,即可生成与真实钻孔数据和地层信息吻合的三维模型。程序提供了3种FLAC3D内置单元模型可供选择。该方法采用Kriging插值算法,可生成精确的三维地层模型数据。与FLAC3D软件自身的脚本建模语言相比,该方法简单快速。本文使用3类应用数据进行建模,验证提出的插值建模方法可以针对复杂的地形和地质信息生成精确的三维模型。研究成果为采用FLAC3D进行复杂地质区域数值模拟计算提供了一种简单、高效、准确的建模工具。     

Critical conditions for natural gas charging and delineation of effective gas source rocks for tight sandstone reservoirs

Although it has been shown that the potential of tight‐sand gas resources is large, the research into the mechanisms of hydrocarbon charging of tight sandstone reservoirs has been relatively sparse. Researchers have found that there is a force balance during hydrocarbon charging, but discriminant models still have not been established. Based on the force balance conditions observed during gas migration from source rocks to tight sandstone reservoirs, a calculation formula was established. A formula for identifying effective source rocks was developed with the gas expulsion intensity as the discrimination parameter. The critical gas expulsion intensity under conditions of various burial depths, temperatures, and pressures can be obtained using the calculation formula. This method was applied in the Jurassic tight sandstone reservoirs of the eastern Kuqa Depression, Tarim Basin, and it was calculated that the critical expulsion intensity range from 6.05 × 10⁸ m³/km² to 10.07 × 10⁸ m³/km². The critical gas charging force first increases with depth and later decreases with greater depths. The distribution range of effective gas source rocks and total expelled gas volume can be determined based on this threshold. This method provides new insight into and method for predicting favourable tight‐sand gas‐bearing regions and estimating their resource potentials. Copyright © 2014 John Wiley & Sons, Ltd.     

气溶胶对雷暴云起电以及闪电发生率影响的数值模拟

本文利用二维耦合气溶胶模块的雷暴云起电模式,结合一次南京雷暴个例,进行250 m分辨率雷暴云起电模拟实验,探讨了气溶胶浓度对雷暴云空间电荷分布以及闪电发生率的影响。在这个气溶胶模块中,假定一个三模态的气溶胶对数分布,考虑了气溶胶活化过程。结果显示:(1)随着气溶胶浓度增大,雷暴云电荷结构保持为三极型。(2)当气溶胶浓度从50 cm-3增加至1000 cm-3时,水成物粒子浓度上升,雷暴云电荷量和闪电发生率增加明显。(3)气溶胶浓度在1000~3000 cm-3范围时,云水竞争限制了冰晶的增长,导致雷暴云上部主正电荷堆电荷量降低。云滴和霰粒子浓度缓慢上升促进中部主负电荷堆和底部次正电荷堆电荷量继续增大。闪电发生率保持稳定。(4)当气溶胶浓度大于3000 cm-3时,水成物粒子浓度稳定,云内的电荷量以及闪电发生率保持为一定量级。     

闪电和固、液态降水关系差异的数值模拟

为了进一步认识闪电和固、液态降水的关系,本文利用三维雷暴云动力-电耦合数值模式,通过设置敏感性试验组,模拟了一次雷暴过程,分析雷暴中闪电和降水的特征,以及闪电和固、液态降水对垂直风速的依赖关系,探讨闪电与固、液态降水的时空分布关系和单次闪电表征的降水量(RPF:rainyields per flash)。结果表明:对流云降水中,液态降水占主要部分,但固态降水比液态降水对于垂直风速的依赖性更强。随着对流的增强,固态降水在总降水中占的比重越来越大。首次放电时间不断提前,闪电峰值落后垂直风速峰值,总闪数一开始随对流的增强而增加,对流一旦增强到一定程度,总闪数则逐渐减小。固态降水和液态降水的开始时间和峰值时间均随着对流的增强而不断提前,而液态降水出现时间和峰值时间均提前于固态降水。雷暴云首次放电的时间滞后于液态降水,而闪电峰值提前固态降水峰值或与固态降水峰值同时产生。雷暴云中的放电活动集中在强降水区域前缘的较弱降水区,强降水区对应的闪电较少,对流的增强会使降水区域面积、降水量和降水强度增加。由于液态降水总量远大于固态降水总量,固、液态RPF的数值相差达到一个量级,但单位时间内固态降水和液态降水增加的速率相近。在单位时间内闪电次数越多,RPF则越小,而固态RPF和闪电次数的线性相关性明显好于液态RPF,所以利用固态降水可以更好地预报闪电。这些结果有助于进一步认识闪电和降水的关系,并可为闪电预报提供新的思路。     

青海两次多单体降雹过程的雹谱分布特征

利用三维冰雹分档数值模式研究了1979年7月20日和8月5日青海西宁两次多单体降雹发生发展过程中云内微物理结构及雹谱演变特征。结果表明:1)云水、雨水的比含水量随雹云的发展呈先增长后减少的特点,两次个例的冰雹主要通过过冷雨滴冻结而形成,其增长方式是凝华及碰并过冷云水。2)雹谱谱宽与上升气流有关,个例一中倾斜上升气流有利于冰雹循环增长,雹谱较宽;个例二中没有倾斜上升气流,不利于冰雹循环增长,雹谱较窄。     

一次江淮强暴雨过程的湿有效能量及其收支特征

利用高分辨率的WRF模式模拟结果,采用基于格点的湿有效能量计算方案,对一次江淮梅雨期强暴雨发生发展过程中湿有效能量的时空演变特征进行分析,并从定性和定量的角度探讨了能量方程转换项、平流项和垂直输送项对强暴雨过程中湿有效能量的输送和积聚作用。结果表明,强暴雨过程中湿有效能量的时空特征与强暴雨发生发展具有良好的对应关系,对流层低层800 h Pa湿有效能量40×104J·h Pa-1·m-2的等值线范围和该等值线伸展至500 h Pa附近可作为判断强暴雨发生的必要条件。暴雨发生前2~3 h的能量快速积聚及其对暴雨区移向的引导,对强暴雨预报具有良好的指示作用。湿有效能量的水平和垂直输送及转换确保了能量的积聚和对流层中层能量的增加,为强暴雨的发生和维持提供了充足的能量。     

2008年1月贵州冻雨的数值模拟和层结结构分析

针对2008年初发生在贵州地区的严重冻雨过程,分别从环流背景、低空急流和水汽输送条件等方面分析了准静止锋维持的原因,并选取本次灾害最严重的第3次过程为典型个例,利用WRF模式针对准静止锋影响下的贵州冻雨进行数值模拟来研究冻雨的发生机制。模拟结果较好地反映出高低空环流形势场特征,强雨雪降水带的走向、落区,以及地面温度的分布,均与观测基本吻合。通过分析高分辨率模式的模拟结果,揭示了准静止锋上贵州地区冻雨的层结结构特征及云物质在冻雨区的分布特征。研究结果表明,贵州中部的冻雨区除一般的三层结构(包含冰晶层、暖层和冷层)外,还具有典型的两层结构特征,即:高空的固体降水粒子稀少,900~600 h Pa深厚的逆温层和0℃以上的暖层使中低空存在大量液态粒子,下落的液滴经过近地面的浅薄冷层,形成大量过冷却雨滴,而后降落至地面迅速冻结。     

建筑物高度对上行闪电触发以及传播影响的数值模拟

为了探讨建筑物高度对单个上行闪电触发以及传播的影响,设定了一个固定的背景电场,并结合自行触发的上行闪电随机放电参数化方案,进行了二维高分辨率上行闪电放电的模拟试验。结果表明:(1)上行闪电在初始阶段分支比较少;发展到离地面2 km左右后,闪电开始出现大量的分支,闪电通道开始出现明显的分叉:一部分通道继续向高电荷密度中心垂直传播,另一部分通道绕过高电荷密度中心,向外水平传播;模拟的上行闪电只能垂直传播到4 km处的负电荷中心,不能穿过0电势线向上方的正电荷区传播。(2)建筑物高度对上行闪电的触发起了关键作用,建筑物越高,越容易触发上行闪电。(3)建筑物高度对上行闪电传播具有一定的反作用,随着建筑物高度增高,模拟出的上行闪电的水平和垂直传播距离都有所减小,通道的分形维数变小,通道传播的总长度也逐渐减小。     

一种新型高度地形追随坐标在GRAPES区域模式中的实现:理想试验与比较研究

将一种新的高度地形追随坐标(Klemp坐标)引入了GRAPES区域模式,并与传统追随坐标(Gal-Chen坐标)和平缓地形追随坐标(SLEVE,Smooth Level Vertical coordinate)进行了比较。对不同坐标下气压梯度力的计算误差通过理想静止大气试验进行了评估,结果表明:与Gal-Chen坐标和SLEVE坐标相比,Klemp坐标有效地减小了气压梯度力的计算误差。理想重力波模拟试验表明,Klemp坐标下对重力波的模拟相比其他两种坐标也更接近于解析解。模式进一步采用了Mahrer气压梯度计算方案减少了计算误差,并提高了模式的精度和稳定性。实际个例试验与理想试验的结论相似。     

2015年4月大气环流和天气分析

2015年4月环流特征如下：极涡为偶极型环流,极涡主体位于格陵兰西侧与加拿大之间,另一中心位于亚洲的东北部;两个极涡中心较常年平均偏低4~8 dagpm,乌拉尔山高压脊偏高8~12 dagpm;副热带高压面积偏大、西脊点偏西,南支槽位置与强度接近常年平均;4月全国平均气温11.6℃,较常年同期（11℃）偏高0.6℃;全国平均降水量43 mm,接近常年同期（44.7 mm）,我国中东部长江以北地区降水偏多,江南、华南降水偏少。4月全国强对流天气过程频发,1—4、19—20和28—29日出现了范围较大的强对流天气。北方地区出现2次沙尘天气。