大地电磁测深三维正演模拟及异常特征研究

为了实现大地电磁测深三维高效快速正演算法,本文研究了大地电磁测深积分方程三维正演模拟的理论与方法.首先详细地推导了Fredholm型积分方程;然后给出积分方程的离散方法、并矢格林函数积分的求解方法以及大型矩阵方程的存储与求解方法;最后,通过理论地电模型的计算检验了方法和程序编制的正确性,并展示了几种理论地电模型的数值计算结果,分析其异常特征.算例结果表明,该正演算法在CPU主频为1.81GHz、内存为1GB的计算机上能正确地模拟多个异常体的异常特征;地表电性不均匀体会使地表视电阻率曲线产生十分严重的畸变;TE模式相对TM模式纵向分辨率比较好,而TM模式的横向分辨率较好;随着异常体埋深的增加,反映异常的分辨能力进一步降低.     

山地冰川流动模型探讨

近30 a来冰川动力学模型有了快速发展, 在南极、 格林兰冰盖预测中取得一系列重要成果, 对山地冰川的研究也初见端倪. 从冰川流动的力学过程出发, 利用本构方程、 理想冰川假设、 浅冰层近似(Shallow ice approximation)假设完整地推导了理想冰川流动的物理过程, 揭示了冰川流动的机理, 建立了气候变化和冰川自身重力引起的理想冰川物质和能量再分配的温度耦合三维流动模型. 结合山地冰川的冰床形态, 将理想冰川与实际冰川相结合, 使理想冰川流动模型更好地近似山地冰川的流动.     

具有精确色散性的非线性波浪水平二维数学模型

建立具有色散性的水平二维非线性波浪方程,方程的非线性近似到了三阶。方程以波面升高和自由表面速度势表达的微分-积分型数学方程,给出方程的数值求解方法和算例,对方程积分项的处理给出了计算方法。计算结果与Boussinesq方程模型和缓坡方程模型的对应计算结果进行了对比。     

一类具有梯度项的超定边值问题中解的对称性

利用经典的平行平面移动法研究一类具有梯度项的拉普拉斯方程超定问题中解的对称性,得到此类超定边值问题解和区域对称的充分条件.结果发现,解和区域的对称性依赖于非齐次项关于空间变量的连续性、对称性及关于梯度变量偏导数的连续性.     

印度尼西亚方程的改进及其在大庆长垣低孔渗储层的应用

针对大庆长垣F油层储层具有孔隙度小、渗透率低、泥质含量高、微孔隙发育、束缚水饱和度高等特点,通过理论和实验分析长垣F油层低孔、低渗储层的主要导电因素,得出泥质附加导电和束缚水导电是该区F油层的低孔、低渗泥质砂岩储层的重要导电因素.考虑束缚水和可动水导电路径不同,基于三水导电模型的推导原理,建立混合水等效电导率方程;依据该区矿化度较低和含泥高的储层特征,对适用于淡水高泥储层的印度尼西亚方程进行改进,用混合水等效电导率代替原方程中的地层水电导率,建立长垣地区F油层低孔低渗泥质砂岩含水饱和度模型.利用岩心分析数据和常规测井资料,采用多元回归法确定长垣F油层饱和度模型中的泥质体积分数、孔隙度、束缚水饱和度.利用岩电实验数据,采用最优化方法确定长垣F油层饱和度模型中的饱和度指数、胶结指数等参数.经密闭取心井含水饱和度对比和试油结果验证,建立的含水饱和度模型可用于F油层的低孔低渗泥质砂岩储层定量评价.     

基于CTRW的聚驱后示踪剂产出特征分析

应用传统的对流扩散方程模拟聚驱后示踪剂产出曲线存在不适应性,将产生质量不守恒或模拟不出实际产出曲线问题.根据连续时间随机游走( CTRW)理论,利用CTRW模型对含聚合物的示踪剂测试实验的产出曲线进行模拟.结果表明:CTRW模型模拟的示踪剂产出曲线克服了基于传统的对流扩散方程模拟示踪剂产出曲线时所遇到的问题;CTRW模...     

深水悬臂桥墩动力特性及地震响应分析

利用ADINA中的势流体单元对一已建成的深水高墩建模,桥梁上部结构对桥墩的影响以墩顶集中质量的形式体现,分别对不同水深条件下桥墩-水体系的动力特性和地震动响应进行了系统的求解;为了检验Morison方程在深水桥墩动力分析中的有效性,将基于忽略速度力项的Morison方程的计算结果与基于势流体的结果进行了比较。研究结果表明：①水体对结构的自振频率折减影响随着水深的增大而增大,基于Morison方程得到的结构频率的影响大于势流体模型;②水深对墩顶位移的影响比加速度明显,墩底剪力受水深的影响大于其对墩底弯矩的影响;③考虑上部集中质量的影响,结构的响应会不同程度地减小;④相对于势流模型,忽略速度项的Morison方程所得到的结构动力反应略显保守。     

渗透作用下多孔介质中循环浓度污染物的迁移过程研究

根据饱和多孔介质中考虑释放效应的一维渗透作用的可溶性污染物的迁移控制方程,通过Laplace变换和Fourier变换及其逆变换求得相应的通解形式。根据半无限体表面可溶性污染物点源注入情形下的基本解,通过积分方法得到半无限体表面圆形区域上作用循环浓度污染物后,在多孔介质内部污染物迁移过程的求解方法。作为一个典型算例,对渗透作用下循环浓度污染源引起的多孔介质中的迁移过程进行分析。算例表明,当污染源浓度为周期变化时,多孔介质内部污染物的浓度随时间增长由不稳定的周期变化过程逐渐过渡到稳定的循环过程,而某一深度处污染物浓度的相位则相应滞后。此时,稳定后的污染物浓度周期与污染源浓度周期相同。实际上,随多孔介质表面污染源浓度的周期变化,在靠近多孔介质表面一定深度范围内的污染物浓度在深度方向也呈增大或减小的交替变化过程。另一方面,在污染源循环变化过程中污染物不断向深处推进,而最终其影响范围限定在某一深度内。     

基于有限差分原理的预应力锚索抗滑桩改进计算方法

针对现有预应力锚索抗滑桩计算模型及其相应计算理论存在的问题,基于有限差分法的原理,提出了改进的锚索桩计算模型并进行了详细的理论推导。根据实际的施工和受力过程,该模型包括预应力施加和滑坡推力作用两个阶段,通过滑坡推力分布形式的改进和桩锚协调变形方程的修正来真实反映预应力锚索抗滑桩的主动式受力特点。同时利用Matlab编制成相应的计算程序进行算例验证,并同目前普遍采用的计算理论进行对比分析。结果表明,该改进计算方法及其相应的计算程序充分考虑了预应力施加这一锚索桩特有的过程所带来的影响,且在计算中不必将锚索拉力等效为在桩顶处施加的剪力和弯矩。另外还可对桩前滑面以上存在任意厚度滑体和任意复杂地基条件的情况进行合理地设计计算,相对于目前已有的计算理论,该改进方法更加合理、可靠。     

修正Van Genuchten模型拟合干旱、半干旱地区遗址土颗分曲线的可靠性研究

不均匀系数、曲率系数在土力学中的应用很广,但对精确求解d_(10)、d_(30)、d_(60) 3个参数的方法研究较少。提出了一种修正Van Genuchten(MVG)方程描述颗分曲线。以Matlab编程语言中的cftool为拟合工具,以青海黄土、青海遗址土、甘肃遗址土和新疆遗址土等33组土为研究对象,对比拟合方程描述遗址土颗分曲线的效果,分析MVG模型中3个参数a、m、n对曲线的影响。结果表明:MVG模型在拟合33组土的颗分曲线时均具有较好的效果,相关系数2R(29)0.98;参数a导致曲线偏移并可以推测土粒中的最大粒径,参数m影响土中的各颗粒组分的比例,参数n影响曲线的斜率。提出的MVG方程能够精确地求解d_(10)、d_(30)、d_(60) 3个参数。