大气污染物SOx输送方程的尺度分析

文中用量纲分析方法分析了控制大气污染物输送方程的动力学特性。提出了 6种反映大气污染物输送过程中各种动力、物理和化学过程相对重要性的动力学参数。并以硫氧化物为例进行了各参数量级大小分析 ,尤其对 3种起汇作用 (去除机理 )的机理 (化学转化、干沉降和湿清除 )在不同尺度大气污染过程中的作用进行了详细分析。结果表明 :在一般情况下 ,SO2 的气相化学作用小于干沉降和湿清除的作用 ;干沉降作用很依赖于模式最底层厚度的选取 ;在有降雨时湿清除作用一般较大。文中还对大气污染物SOx 输送方程各项特征量的量级作了分析对比 ,得出了控制不同尺度大气污染物输送过程的零级近似方程和一级近似方程 ,并指出了这些方程的基本特征     

平面相控阵脉冲压缩的雷达气象方程

在低发射功率条件下,相控阵多普勒天气雷达为了增加探测距离和提高分辨率,需要采用脉冲压缩技术.由于采用相控阵技术,波束扫描过程中,在平面阵非法线方向上会产生波束展宽,引起波束内功率下降,需要在计算雷达气象方程时弥补功率下降的误差.平面相控阵天气雷达在探测远距离目标时,使用宽脉冲发射,经脉冲压缩处理以提高分辨率,因此需要建立适用的雷达气象方程.首先讨论了脉冲压缩低峰值功率的平面相控阵多普勒天气雷达的结构、线性调频脉冲压缩的距离(多普勒测量)、调频波形的耦合问题以及目标距离的测量;然后给出了适用于脉冲压缩的平面相控阵雷达气象方程,此方程同时也适用于一维线性阵的脉冲压缩的天气雷达.     

拉萨光合有效辐射的变化趋势及其估算方程的建立

利用2005年1月-2010年12月中国生态研究网络拉萨农田站的观测数据,分析了拉萨地区光合有效辐射Qp的变化特征.结果表明,拉萨地区的Qp与太阳总辐射Rs具有相同的月变化和日变化特征,Rs与Qp的平均日累计值分别为20.2 MJ·m-2和37.2 mol·m-2,Qp呈逐年递减趋势.Rs与Qp的月变化表现为6月最大、1月最小,其最大、最小值分别为25.4 MJ·m-2、14.0 M J·m-2和48.7 mol·m-2、25.2 mol·m-2.日变化表现为正午大、早晚小.光合有效辐射系数(Qp/Rs)的变化范围在1.3～3.2 mol·MJ-1之间,年平均值为1.89 mol·MJ-1,呈逐年递减趋势,2010年比2005年降低了6.7％.Qp/Rs的日变化表现为早晚大、正午小,季节变化表现为夏季最大、春秋季次之、冬季最小,而且夏、冬两季的季节变化较为明显,冬季比夏季低10.6％.通过分析Qp与晴空指数和太阳天顶角之间的关系建立了拉萨地区Qp的估算方程,相对误差在4.49％左右.     

安阳市区11～3月雾生成短期预报

统计了安阳市区1970-1989年气候资料，分析了雾生成前一天的天气形势和影响系统，并以此确定了雾生成的起始场条件。在满足起始场条件的基础上，选择与雾生成关系密切的相对湿度、气层稳定度和风向风速等因子，建立预报方程。     

热液成矿作用的一般动力学方程

分析了热液成矿作用的基本特点和研究现状，推导出热液成矿作用动力学体系的一般方程，并着重讨论了方程中表征热液与围岩之间化学反应的速率项，除目前一般采用动力学实验实测反应速率常数外，作者根据局域平衡原理，认为成矿作用可视为体系内各点随着温度和压力的变化不断重建平衡的过程，可用温度和压力变化时成矿反应的平衡常数的变化速度表征成矿反应的动力学行为，成矿化学反应的速率还可用主要成矿物质的溶解的变化率来表征。     

改进时间分数阶模型模拟非Fick溶质运移

通过将经典时间分数阶对流-弥散方程的等待时间分布函数的尾部修改为指数型,推导出了改进时间分数阶对流-弥散方程,并提出有效的时空算子分裂数值求解方法。对两个理想算例和一个实际算例进行计算,结果表明,改进的时间分数阶对流-弥散方程继承了时间分数阶对流-弥散方程能模拟穿透曲线幂率型拖尾分布的优点,还可模拟穿透曲线尾部由幂率型转换到指数型的过程;特征时间λ、分数阶指数γ和两相容量比例系数β共同决定了运移行为。改进的新模型可以区分非均质介质中流动相和非流动相中的溶质浓度, 更细微地模拟非Fick溶质运移行为。     

平顶山地下盐穴储气库泥岩夹层稳定性评价

泥岩夹层在盐溶建腔阶段难以被溶蚀,盐层溶采后仍赋存于腔内。针对平顶山地下储气库在储气运行阶段,未被溶 蚀的夹层因岩体流变性出现垮塌而影响储气库安全运行的难题,提出一种基于夹层垮塌时间的稳定性评价方法。首先,根 据弹性板壳理论建立夹层力学模型,解得夹层局部破坏的极限应变;其次,根据夹层岩样蠕变实验,建立蠕变本构模型。 最后,结合夹层力学模型与蠕变本构模型,建立垮塌时间计算模型,并利用基于PSO 算法的Matlab 程序对平顶山储气库夹 层不同预设厚度的垮塌时间值计算,对临界垮塌厚度值进行反演计算。结果表明：随夹层厚度增加,其计算垮塌时间增 大,夹层稳定性提高;厚度相同的夹层,其跨径越大,计算垮塌时间值越小,计算临界垮塌厚度值越大,即稳定性越差; 对小于临界垮塌厚度的夹层,其在设计使用年限内垮塌失稳风险较高,反之则垮塌失稳风险较小。     

三角形网格下求解二维浅水方程的和谐Godunov格式

为保证计算格式的和谐性,通过特殊的底坡源项处理技术,在三角形网格上建立了求解二维浅水流动方程的具有空间二阶精度的Godunov格式。应用准确Riemann解求解法向数值通量,用改正的干底Riemann解处理动边界问题。经典型算例和钱塘江河口涌潮计算验证,表明模型健全,分辨率高,具有较大的推广应用价值。     

三维交叉裂隙渗流传质特性数值模拟

对交叉裂隙渗流传质特性的定量描述是研究整个裂隙网络渗透传质特性的基础。为真实模拟水流及溶质在三维交叉裂隙中的运移过程,首先通过三维轮廓仪获取天然岩石裂隙表面的形貌数据,再应用三维重构技术生成相应的三维交叉裂隙模型,随后求解Navier-Stokes方程,假定溶质运移满足Fick定律,模拟水流和溶质在三维交叉裂隙中的运移过程。通过对比粗糙裂隙模型与平行平板模型的模拟结果发现：粗糙度对流体的分布及流动状态存在显著的影响;不同进、出口工况下的流体流动及溶质运移状态亦表明：裂隙交叉的几何形貌会显著地影响溶质混合行为。这些结果表明,目前被广泛采用的平行平板模型在评估岩体内特别是交叉口的物质运移特性时将导致较大的偏差,在将来的研究中有必要针对裂隙交叉口的几何特征建立修正的模型以提高评估的准确性。     

应用最优控制理论自动率定二维浅水方程的糙率参数

应用偏微分方程最优化控制理论,导出了二维浅水方程糙率自动率定的最优化系统,实现了二维水力学模型计算中分布式糙率参数的自动率定算法。结合实际应用的可操作性,按区块设定糙率,对恒定流和非恒定流两种流态分别进行了常数糙率和随水深变化的糙率自动率定实验。数值实验表明,在观测数据包含足够的信息时,该方法不仅能够准确地自动率定出相关的糙率参数,还能成功地将糙率随水深的变化关系式识别出来。但是,糙率水深变化关系识别迭代收敛速度显著变慢,因此,建议采用相应的实现对策来考虑变糙率的洪水预测问题。