三维地电场数值计算中的LANCZOS迭代过程及其改进算法

针对三维地电场正演数值计算过程中形成的超大规模稀疏线性方程组,在分析此类线性方程组的一般解法基础上,着重阐述一种适宜求解此类方程组的Lanczos迭代过程与算法原理。同时,当地下介质的电性差异较大时,形成系数矩阵A的条件数就很大,可对算法进行适当改进。讨论采用不完全Cholesky分解方法进行预条件处理,经过条件数改善后,形成新的线性方程组系数矩阵,就会变为一个近似的单位矩阵。经改进后的Lanczos算法,将提高数值计算稳定性,从而加快迭代收敛速度,为提高反演质量提供基础。     

深部软弱岩体置换体对岩质边坡支护效果分析

大岗山水电站右岸边坡位于大渡河断裂和磨西断裂所切割的断块上,地震发生频繁且烈度高,岩体卸荷强烈,构造稳定性差,裂隙发育集中。已有研究成果显示出边坡开挖过程中锚索支护的重要性,也有较多学者对抗剪洞支护的设计进行了较全面的研究。在已有的锚索支护基础上,基于卸荷理论对设计抗剪洞的再次支护效果和和采用拟静力法模拟的抗震效果进行评估,通过不同高程的关键点的位移响应,较系统地论证抗剪洞设置的必要性。通过研究应力的变化发现锚索支护段和抗剪洞支护段之间形成较为明显的拉裂区域。对工程后期监测提出建议。     

解决土地碎片化的制度框架设计

为解决土地碎片化问题,采用文献统计和比较分析等研究方法,在对土地碎片化现状、原因与问题分析的基础上,通过总结国内外实践经验,提出了解决土地碎片化的制度框架。研究结果表明,土地碎片化的产生既有自然原因也有制度原因,所带来的问题主要包括:1田埂面积多而导致的耕地资源浪费;2地界纠纷多而导致新技术新产品难以推广;3谈判成本高导致农地资源无法合理进行配置;4机械水平低导致劳动生产率、土地生产率和成本产值率也随之下降等。必须要通过明晰产权关系、构建市场平台、创建新型农业经营体系、推广"良田良种良法"等政策措施来逐步破解中国土地碎片化难题。     

岩石蠕变模型研究进展及若干问题探讨

岩石蠕变是岩土工程变形失稳的主要原因之一。近年来蠕变研究正处于一个探索阶段,本文从四个方面综述了蠕变模型的研究进展。研究发现,在岩石蠕变的三个阶段中利用经典本构模型均很难描述加速蠕变阶段,研究者们通过新的元件或者改进的非线性黏弹塑性本构模型可以很好的模拟岩石蠕变实际曲线;基于损伤理论的岩石蠕变模型是近年来发展的主要方向,可以很好的解决岩石微观裂纹所带来的蠕变;随着岩石深部工程的发展,岩体受到周围实际环境下的影响是不可忽略的,从而研究含水量的变化与水力和其它应力耦合下的岩石蠕变也是今后的重点。最后指出,由于试验仪器的原因,高温高压和各向异性下的岩石蠕变模型研究进行的还不是很多,是今后岩石蠕变研究的难点。     

基于卸荷和抗震分析理论的大岗山坝肩边坡锚固优化研究

基于卸荷岩体的理论和方法,对大岗山水电站坝肩边坡进行研究。以大型通用有限元程序ADINA为计算平台,对比分析坝肩边坡卸荷后,不同地震加速度及不同锚固长度和锚固吨位在卸荷条件下关键点的位移和应力以及塑性区的变化规律,认为边坡开挖变形由两部分组成,一部分为卸荷回弹,另一部分是由于大吨位锚索施加后产生的朝向坡内的压缩变形。研究表明,由于锚索穿过V岩体底部界线,在可能的滑坡剪出口形成强大的阻滑键,增加锚索长度或吨位均可有效地减小坝肩塑性区,地震加速度为0.1g和0.15g时,在各加固吨位下均可能形成塑性区的贯通;在不同的地震加速度下,受不同锚固长度和锚固力影响,不同关键点处水平向位移呈现线性变化和非线性变化,非线性关系中的拐点应重点关注;若要减少岩体塑性区,相同比例下增大锚固长度效果更好。为了有效减少在高地震烈度时边坡塑性区的发育,建议采用增加120%的锚固长度。     

用改进的光滑约束最小二乘正交分解法实现电阻率三维反演

对三维电阻率反演问题进行了深入研究,提供了一种利用地表观测数据实现三维反演的实用算法.该方法应用有限差分求正演解,并通过对粗糙度矩阵元素进行适当改进,使之适用于各种情况下粗糙度矩阵的求取,进而建立在模型的总粗糙度极小条件下的反演方程.对反演方程采用收敛速度快且稳定的最小二乘正交分解（LSQR）法进行迭代求解,在迭代求解过程中只需利用偏导数矩阵和其转置矩阵乘以一个向量的结果,回避了直接求偏导数矩阵的繁琐计算,节省了内存,加快了反演的计算速度.不同的计算实例表明上述方法是求解大规模三维电阻率反演问题的有效方法.     

饱和砂岩对周期性循环载荷的动态响应

在单轴循环加载实验中,改变平均应力和正弦波应力振幅的条件下,研究了饱和砂岩的衰减、杨氏模量和泊松比的动态响应。获得随平均应力增大杨氏模量增大,随应力振幅增大应力～应变滞回曲线向应变增大方向移动,随应力振幅减小滞回圈包含在峰值应力的滞回圈内,反映出饱和岩石记忆特性。饱和岩石随应变振幅的增大衰减增大,模量呈线性下降,泊松比呈非线性增大。在低频段,这些物理力学参数的弹性非线性特征,对研究天然地震和储层勘探,以及它们的资料解释具有现实意义。     

三维电阻率正演计算中的Lanczos迭代算法

在三维电阻率的正反演计算中,快速、准确的正演计算是反演的关键。而正演计算往往涉及到求解大型线性方程组Ax=b的问题,通过Lanczos迭代构造出对称三对角阵方程组,并采用正交分解法进行求解,与传统算法相比,此算法占用内存少、收敛速度快、且稳定;针对大型稀疏矩阵的特点,采用简单地记录矩阵的非零元素值及其所在行、列值的方法,来存储大型稀疏矩阵,可大大节省机器内存,提高运算速度。通过理论分析和点电源三维地电场计算实例,阐述该法是地电三维正演计算的有效方法。     

LANCZOS迭代算法及其在三维地电场数值模拟计算中的应用研究

在三维电阻率的正演计算中往往涉及到快速、准确求解大型线性方程纽Ax=b的问题。通过采用有限差分法来构造出求解点电源三维地电场的大型稀疏对称线性方程组。并引入Lanczos迭代技术，构造出三对角阵方程组，然后采用正交分解法进行求解，它是Krylov子空间方法中的一种。与传统迭代算法相比，它占用内存少，收敛速度快且稳定。针对大型稀疏矩阵及MATLAB语言的特点，采用简单记录矩阵的非零元素值及其所在行、列值的方法存储大型稀疏矩阵，可大大节省机器内存，提高运算速度。理论分析和计算实例显示，此算法是地电三维正演计算的有效方法，为下一步的反演计算打好基础。     

循环载荷下饱和砂岩的疲劳损伤

通过对大庆石油饱和长石砂岩的单轴循环试验和桔红色砂岩的三轴循环破坏试验，获得了应变随时间变化曲线(类似一条常载荷的蠕变曲线)，并研究分析了岩石的模量随循环数增加而下降，衰减随循环数的增加而增大，损伤随循环数的增加而增大直至破坏的变化规律。岩石的这种损伤演化研究对增加石油的采收率和岩土工程的稳定性研究具有现实意义。