Delving into three-dimensional structure of the West Luzon Eddy in a regional ocean model

The three-dimensional structure and associated dynamics of the prominent cold (cyclonic) West Luzon Eddy (WLE) were investigated by a high-resolution regional ocean model. The WLE was horizontally and vertically heterogeneous, exhibiting asymmetric structures in the circulation, vorticity, vertical motion and energy distributions within the eddy. The asymmetry was mainly attributed to the existence of an eddy dipole formed by a coexisting warm (anti-cyclonic) eddy to the south of the WLE. Analysis of the momentum balance revealed that the coexistence of two eddies intensified barotropic pressure gradients in the southern WLE to locally enhance the eastward jet. The positive (negative) vorticity of the jet strengthened (weakened) the eddy in the southern sector (periphery), which, together with the formation of a subsurface density front, intensified (suppressed) the corresponding upward motion and cooling. The baroclinic pressure gradients opposed the dominant barotropic components and spun down the eddy at greater depths with stronger weakening in the southern sector near the front. Asymmetric energy distributions showed that larger mean kinetic energy (MKE) and eddy available potential energy (EAPE) were stored in the southern sector of the WLE. While the larger MKE was directly linked with the stronger barotropic currents, the larger EAPE in the southern WLE was formed by baroclinic energy conversions due to a strong density gradient at the front.     

Numerical modeling of water yield of mine in Yangzhuang Iron Mine,Anhui Province of China

This study develops a three-dimensional heterogeneous numerical model to simulate the water inrush process and predict the water yield for mineral exploration in Yangzhuang Iron Mine in Anhui Province. To identify the hydrogeological parameters of the aquifer in the study area, the model was calibrated and validated using the observed heads through the integrated trial-and-error and automated techniques. Also, the sensitivity analysis of the model was performed to evaluate the uncertainty associated with the calibrated model. According to the mine construction plan at different mining levels of -500 m, -600 m, and -700 m, the calibrated model was then applied to predict the water yields dependent on the different mining levels. As indicated by the prediction results, the numerical simulation model can systematically describe the groundwater system in the mining area and determine the source of water inrush in this iron mine. In conclusion, numerical analyses carried out in this study can provide guidance to decision-makers in balancing the iron ore mining and mine dewatering in the future.     

衰减夹层GPR模拟的时频域全局反射误差

全局反射误差分析是深入研究探地雷达（GPR）吸收边界条件吸收效率的有力工具.基于常规完全匹配层（PML）的标准交错网格有限差分算法必须满足严格的CFL条件限制,即在单位时间步长内,不容许电磁波传播的距离超过单元网格尺寸.为了提高主区域所有网格节点的计算效率,并有效地吸收传播后期出现的低频隐失波,提出基于非分裂递归卷积完全匹配层（UCPML）的旋转交错网格（RSG）GPR正演算法.该算法采用不同的网格交错策略,并在边界条件中引入了吸收低频隐失波的自由可变因子,使得该算法允许选取较大的时间步长,提高了计算效率,并且实现了对低频隐失波的高效吸收.本文首先给出了RSG差分格式,推导了基于UCPML的RSG差分更新方程,探讨了数值色散的稳定性条件,然后以绕射现象严重的衰减夹层数值模拟为例,分别从波场快照、单道波记录、时间域反射误差（TDRE）、频率域反射误差（FDRE）四个方面分析了UCPML与常规PML的全局反射误差,说明了基于UCPML和RSG的GPR正演算法能更有效地吸收低频隐失波.

     

库车坳陷克拉2气田异常流体压力演化史

根据克拉2气田现今实测地层压力分析了异常高压发育的特征,并通过流体包裹体测试求取包裹体形成时期的古压力。计算结果表明:下白垩统和古近系库姆格列木群砂岩储层在距今8～6Ma(康村组沉积期末)时,古压力系数大致为1.09～1.18,发育过剩压力为4.6～9.6MPa的异常高压。在此基础上,以古、今压力数据为约束条件,利用数值模拟技术定量恢复出下白垩统与侏罗系异常流体压力的演化历史。研究发现:下白垩统储层流体压力经历了沉积型异常流体压力形成发展(白垩纪-库车组沉积期末)和沉积型异常流体压力萎缩—构造挤压型异常流体压力孕育(库车组沉积期末至今)两个阶段。康村组沉积早期沉积型超压开始出现,库车组沉积期末达到顶峰,过剩压力达40～50MPa;此后,在沉积型超压萎缩—构造挤压型超压孕育过程的耦合下,过剩压力曾降至25～30MPa,更新世起超压略有升高。侏罗系流体压力演化较为简单,吉迪克组沉积期开始出现超压,至库车组沉积期末达到高峰,此后异常高压逐渐萎缩。在克拉2气田的主要成藏时期(距今3～1Ma),烃源岩较储层具有更高的异常压力,源—储压力差为10～30MPa。     

不同边界层参数化方案在暴雨数值模拟中的对比分析

本文利用非静力平衡的中尺度模式(MM5),通过对一次华南锋前暖区暴雨的数值模拟,比较了3种边界层参数化方案的作用。结果表明：从24小时累积降水的分布来看,MRF方案的预报与实况对比效果最好,其次是Blaekadar高分辨率PBL模式,最差是no PBI,方案。分析MRF方案的输出结果,发现对应每个时刻的累积雨团上空,925hPa都有较强的辐合辐散偶极子与之相配合,而在300hPa则有较强的辐散。不同边界层方案模拟的辐散场有一定的差别,导致了降水的不同预报效果。     

三维变分同化中不同物理量对台风预报的影响

采用NMC方法统计了2006年9月的背景误差协方差,利用Bogus资料对台风进行了初始化,并对2006年“桑美”台风进行了同化试验。结果发现,同化不同的Bogus资料,所得到的台风初始场各不相同,对台风预报的影响也各不相同。对海平面气压同化,可以使台风强度明显加强,形成成熟的暖心结构;基于风速同化,对台风强度的改变较弱,对暖心结构的改进不是很明显。在同化海平面气压和风速的基础上,针对相对湿度的同化,在一定程度上可以改善台风强度预报,有利于提高台风路径的预报精度。     

中国长三角地区PM2.5时空分布数值模拟研究

利用第三代空气质量预报模式LOTOS-EUROS(Long Term Ozone Simulation-European Operational Smog)对2018年中国长三角地区细颗粒物(PM2.5)浓度的时空分布进行数值模拟,通过对比模拟结果与地面观测值,验证模式对PM2.5长期特征模拟的合理性并探讨长三角地区P...     

磁层亚暴扩展相的数值模拟

本文用多步隐格式方法,求解包含电阻的磁流体力学方程组,对磁层亚暴扩展相作数值模拟.计算结果展现了亚暴扩展相期间磁尾变化的主要特征.它表明,球向和尾向等离子体流与准稳态重联有关.等离子体团的尾向喷发,致使中性线尾侧电流片内的密度降低,等离子体片变薄.中性线近地侧的等离子体团朝着地球运动,并合并于地球附近的重联区内.     

Reflected and transmitted waves in a channel with side porous mattresses

The problem of wave propagation and wave damping in a channel with side porous mattresses of arbitrary shape protruding from the walls is studied. The solution was achieved by applying 3-D boundary element method and was employed to study wave field in the channel and to analyze the effect of the geometry of the mattresses and physical and hydraulic properties of porous material on wave damping. The results show that wave damping in the channel strongly depends on wave parameters, especially, on the wave number. Wave reflection and transmission decrease with increasing the wave number. The results also show that the wave field in the channel strongly depends on the geometry of the mattresses as well as on physical and hydraulic properties of porous material used to build these wave dampers. The geometry of the mattresses and physical and hydraulic properties of porous material have a moderate effect on wave reflection and a significant effect on wave transmission. The results show that wave transmission down the channel decreases with increasing the length and thickness of the mattresses. Moreover, wave transmission decreases with increasing the porosity and damping properties of porous media used to build the mattresses. The analysis shows that porous mattresses protruding from the channel walls are very efficient in damping water waves propagating down the channel and may be built in channels to reduce high waves and achieve desired wave conditions. Theoretical results are in reasonable agreement with experimental data.     

1966年以来华北地区一系列七级大震破裂过程的数值模拟

用非连续变形分析方法(DDA+FEM)数值模拟，在华北地区各地块相互制约的块体系统环境中，地块边界断层上发生1966年邢台地震、 1969年渤海地震、 1975年海城地震、 1976年唐山地震等不同类型破裂模式大震的破裂过程. 数值模拟结果给出每个大震释放的主应力场，最大剪应力变化等值线图，地震前后位移变化矢量图，发震断层滑移随时间变化以及走滑错距和应力降等震源参数. 这些结果与地震的震源机制，用地震波资料研究得到的震源参数,宏观等震线,地表观测的水平位移矢量图基本一致. 其中1969年渤海地震正交破裂模式的结果与宏观等震线及小震分布图像更接近. 1976年唐山地震复杂震源模式与该震早期余震分布图像更相符. 表明用DDA+FEM方法进行数值模拟研究，能较好地模拟地震破裂过程.