钱塘江河口盐度数值模拟

强潮河口盐水入侵对饮用水源地危害极大。基于平面二维水动力盐度模型, 对典型强潮河口—钱塘江的水动力及盐水入侵过程进行了数值模拟研究。结果表明枯水径流时盐度变化与潮位过程曲线类似, 潮差对盐度大小影响显著, 径流量的增加将逐渐减小其相似程度。当流量增加到一定程度后, 继续增加的一定径流量所产生的抑咸效果减弱, 水资源有效利用率降低, 此时允许水源地盐度超标并改从蓄淡避咸水库取水可有效节约水资源。盐度平面分布显示, 盐水入侵在强潮河口弯道处受涨潮流主流线影响明显, 靠近主流线一岸的盐度大于对岸, 单从盐水入侵角度考虑, 强潮河口弯道段的取水口应设置在远离涨潮流主流线一岸。钱塘江河口盐度数值模拟对于研究减轻盐水入侵对水源地危害的措施具有指导意义。     

非开挖HDD地表变形数值模拟及DCRP识别技术研究

采用非开挖水平定向钻进技术近地表施工时,由于地层及钻井液压力等原因易产生地表变形。结合美国梅萨市一处通信管道施工实例,对0-0.9 MPa钻井液压力作用下的地表变形程度进行了数值分析,得到地表最大隆起量为26 mm。采用DCRP技术及水准仪方法对该变形区域钻井液压力作用前后的地表进行位移检测,得到地表最大隆起量分别为28 mm和30 mm,表明数值模拟分析结果同现场检测结果比较吻合。同时,在现场所布置的26个检测点中,两种方法检测到的隆起量最大相差2.7 mm,表明其精度比较接近。此外,由于DCRP技术属于非接触式三维检测,受施工现场条件及人员操作限制少,故具有较为理想的应用前景。     

汶川地震引发高速远程滑坡运动机理数值模拟研究——以谢家店子滑坡为例

针对汶川地震引发的谢家店子滑坡,在现场调查分析的基础上,建立了二维离散元数值模拟模型,采用2D-Block软件对其进行了全过程的数值模拟研究,并通过对跟踪块体的深入分析,研究了相应地质体在不同阶段下的运动特征。模拟结果表明,谢家店子滑坡经历了剧动启程抛掷阶段、快速撞击飞行阶段、铲刮减速碎屑流阶段及堆积掩埋阶段。为了揭示地震引发高速滑坡的发生规律,分别研究了地震震级、斜坡地形和斜坡上岩块的尺寸对高速滑坡启动和运动过程的影响规律。地震震级对边坡的启动、变形、破坏和运动有很大的影响,地震震级越大,滑体启动的加速度和速度也就越大,从而易形成高速远程滑坡。斜坡体本身的地形地貌对滑体运动也有较大影响。在震级和岩石力学参数一定的条件下,斜坡上岩块的大小对其启动、变形和运动过程有一定影响,随着岩块的增大,滑体运动的每个阶段历时都在减小,但当岩体十分破碎时,滑体虽然能够运动,但是很难发生抛掷。将地震滑坡的启动机理概括为积累变形效应、振荡启程效应和振荡加速效应。     

煤层瓦斯卸压抽放动态过程的气-固耦合模型研究

根据瓦斯渗流与煤体变形的基本理论,引入煤体孔隙变形与透气性演化的耦合作用方程,建立了考虑煤层吸附、解吸作用的含瓦斯煤岩固-气耦合作用模型。应用该模型模拟研究了不同压力影响下瓦斯抽放过程中煤层透气性的演化和抽放孔周围瓦斯压力的变化规律,其结果对于深入理解瓦斯抽放作用机制并采取相应的瓦斯预防和控制措施等具有重要的理论和实践意义。     

多步分解法初步研究石英砂中甲烷水合物p-T稳定条件

设计加工了一套研究海底沉积物水合物稳定条件实验装置,采用多步加热分解的方法对不同孔隙度的石英砂沉积物甲烷水合物相态平衡点(p-T稳定条件)进行了初步研究。为校验实验装置和实验方法,首先对甲烷水合物在十二烷基硫酸钠(SDS)溶液中的生成和分解过程进行了研究,获得的实验数据与文献数据吻合得很好。在此基础上,以4种不同粒径(96~180μm、180~380μm、212~380μm和380~830μm)的石英砂为介质,分别研究了甲烷水合物的p-T稳定条件。结果表明,实验测得的数据与纯水溶液中的结果基本一致,说明在实验所用的不同粒径石英砂构成的孔隙尺度内,毛细管作用对水合物相态点的影响可以忽略。     

水平长导线源瞬变电磁三维积分方程模拟

以体积分方程法模拟三维电磁正演问题的原理为基础,利用数值滤波方法实现了时间域水平长导线源的瞬变电磁三维问题的模拟,对电阻率、极化率、时间常数和频率相关系数的变化分别进行了三维数值模拟。为了提高计算的精度,变换过程中采用双精度数值滤波方法。利用偶极源的解析解和数值解的对比来检验数值滤波方法的可靠性。数值模拟结果表明,此种数值滤波方法精度高,结果稳定。高精度数值滤波结果为瞬变电磁法在石油、矿产、地热以及其他方面的应用提供了理论指导。     

鄂西恩施地区非饱和土降雨入渗规律

在分析鄂西恩施地区降雨特征的基础上,选取6组典型降雨条件,设计并进行了室内人工降雨土柱非饱和入渗试验。通过实测土体体积含水率,揭示了鄂西恩施地区非饱和土降雨入渗规律,获得了降雨入渗影响区体积含水率变化规律、入渗前锋运移规律以及降雨强度和历时对入渗的影响规律。通过GeoStudio软件对室内降雨试验进行了数值模拟,数值模拟结果与试验结果基本相同,进一步验证了降雨入渗规律。同时依据实测结果,分析了数值模拟的不足之处,为通过数值模拟方法研究非饱和土降雨入渗奠定了基础。     

松辽盆地徐家围子三维构造应力场数值模拟研究

本文采用有限元数值模拟的方法,以徐家围子断陷深层的重点勘探区域为研究对象,对深部构造应力场进行仿真。结果表明:研究区在断陷期的应力场以张性为主,局部表现为挤压性质;挤压比较强烈的区域多位于现今的构造高位置,徐西斜坡、丰乐低隆起和宋站凸起等地方相对明显;营城组沉积期间不一定发生过区域性的挤压运动,在左旋的张性环境下也可以形成挤压构造。     

速度误差对时间偏移和深度偏移的影响分析

地震数值计算表明:地震速度模型出现误差,对时间偏移和深度偏移有不同的影响,且深度偏移结果的误差往往大于时间偏移的误差,影响了深度偏移方法的广泛应用。根本原因在于时间偏移使用的是均方根速度,而深度偏移使用的是层速度,因此在速度模型相同的情况下,层速度的变化远大于均方根速度的变化。地震数值模拟分析同时说明,深度偏移结果的误差还取决于地下介质的结构和岩性。当界面二边速度差异比较大时,速度摄动值δ引起深度偏移结果的误差,远大于界面二边速度平缓变化的情况。深度成像的误差等于速度摄动值δ和界面二边速度对比度vi+1/vi的乘积。     

长江中下游江湖关系演变趋势数值模拟

以长江中下游防洪系统为对象,概述了在大型复杂防洪系统洪水行为数值模拟基础上,成功地将长江中下游洪水演进数学模型转化为面向长江防洪系统防洪规划方案评估需求的长江中下游江湖水沙演变的数学模型.为适应防洪规划方案论证涉及江湖水沙相互制衡相互关联客观情况,建立了面向江湖水沙关系及其演变的数学模型.针对长江中下游江湖水沙运动特点,在水沙数值模拟的范围内侧重对下荆江河道冲刷、荆江三口分流分沙模式、洞庭湖泥沙淤积、江湖耦合等环节进行了讨论,提出了合理可行的数值处理方法.模拟结果较好反映了江湖水沙演变规律,主要成果已应用于长江中下游防洪规划和防汛调度方案中.