首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 140 毫秒
1.
激发极化数据的最小二乘二维反演方法   总被引:12,自引:1,他引:11  
给出了一种电阻雍激发极化法数据的最小二乘二维反演方法,该方法与以往这类方法的主要差别在于网格单元中电导率和极化率参数呈线性变化而不是均匀值,从而使电经和极化率的反演结果更准确,更容易用等值线成图,方法在二维有限元正演计算中,采用三角单元,使实测数据反演前不需要进行地形改正;在目标函数中加入最简单模型以及背景场等先验信息,既民反演问题的多解性,又使反结果更接近实际情况,还采用了前人用电位函数与模型参  相似文献   

2.
Zohdy方法在三维电阻率测深数据反演中的应用   总被引:2,自引:0,他引:2  
三维电阻率测深数据采用Zohdy方法近似反演,正演计算采用有限单元法,反演初始模型由测量视电阻率数据给出,通过比较实测视电阻率值和预测模型计算的视电阻率值对数差来修改模型网格电阻率.为了解决任意电极距测深数据的反演,采用大、小双网格剖分,大网格反映地下电性分布情况,小网格用于实际有限元正演计算.文章对加5%随机噪声的模型理论视电阻率测深断面数据和实测视电阻率测深断面数据进行反演.结果表明该反演方法需要的内存小、计算量少,并能够取得很好的反演效果.  相似文献   

3.
有限元法在大地电磁场正演计算中的应用及改进   总被引:1,自引:0,他引:1  
有限单元法是一种灵活的、行之有效的解偏微分方程的数值方法,它被广泛地应用于各种科学技术问题中.近年来,也被引用到大地电磁场的正演计算中、大地电磁测深资料的二维自动反演解释,除要求正演计算出二维地质模型条件下大地电磁场沿走向分量的分布外,还需要计算这些分量在地面处沿垂直方向的偏导数,以求出辅助的场分量在地表处的值,从而算出地面阻抗值.另外,还要求能给出计算地面阻抗相对模型参数的偏导数的公式,以便能用计算机自动完成二维反演解释过程.这后面两点确定了有限元法在大地电磁测深应用中的特点.Rodi,W.L提出采用矩形有限单元进行计算的方案,这除了便于计算和编制程序外,还由于用不完全的二次式作为插值公式,提高了计算场分量垂直偏导数和辅助场分量的精确程度和效率,适应大地电磁测深正演计算中的上述特殊要求.这种方法的缺点是矩形  相似文献   

4.
本文对药园山铜矿的控岩、控矿构造进行了讨论,认为东西向隐伏断裂与南北向构造复合,为成矿母岩构成了通道,新华夏系构造是主导控矿构造.利用有限单元法反演了控矿构造应力场和应变能量场,计算结果与地质事实基本相符.  相似文献   

5.
有限单元法三维电阻率最小二乘反演中存在问题的研究   总被引:6,自引:1,他引:5  
针对有限单元法在三维电阻率最小二乘反演中关于深部单元的反演能力、反演稳定性和速度问题进行了讨论与改进.首先构造了全新的"体积因子"作为先验信息加入到电阻率反演中,使地质体赋存区域的电阻率变化更符合实际情况,提高了深度较大的地质体的分辨能力.反演计算中对实测电阻率数据和模型正演数据分别取对数,并改进迭代修正量的计算方法,将以前的加减修正量改为与修正量指数之积,这样,不但可以改进电阻率模型出现负值的缺陷,提高反演计算的稳定性,而且还可减少迭代次数;针对反演时间过长的缺陷,将三维全域反演改为三维局域反演,使每次迭代的时间由4 h减少到目前的15 min;对模型的反演计算表明,以上改进是非常有效的.  相似文献   

6.
电阻率层析成像反演中若干问题的讨论   总被引:1,自引:0,他引:1  
孙洪星  李大海 《江苏地质》1996,20(3):181-184
以双异性点电流源电流场二维地电断面为例,着重讨论了电阻率层析成像反演中有限单元法数值模拟边界条件的选择,网格部分以及模拟精度等问题,同时也讨论了成像反演中最优化反演方法及拟合参数的评价等。文中提出的网格部分方法和选择的边界条件可有效地提高计算速度和数值模拟的精度。  相似文献   

7.
一、引言有限单元法是根据边界条件进行应力应变正演计算最常用而又有效的数值方法。在地球动力学研究中,我们常常希望利用测得有限点的应力或应变值等来确定影响应力应变分布的动力机制,这就是地球动力学问题的反演。在线弹性有限元反演中,可以利用叠加原理获得力学方程组,解该方程组即可得到反演结果。但利用叠加原理反演的物理量种类有限,尤其该方法不能推广到非线性有限单元模型;另外,既使在边界力的反演中,由于观测资料的欠缺,使得叠加原理得到的常  相似文献   

8.
赵永贵 《中国岩溶》1988,7(3):241-246
这里提出一个新的水文地质参数的数值反演方法。这个反演方法是以有限单元和边界单元的交连耦合法(CFEBEM)以及最小二乘法(LSM)为基础构成的。CFEBEM方法同有限单元法(FEM)一样灵活,同边界单元法(BEM)一样精确。应用该反演方法不但可以确定水文地质参数的分布,同时也可获得地下水流场的计算结果。文中还研究了三种反演迭代技术。在太原地下水系统研究中的应用说明这种方法是十分有效的。   相似文献   

9.
对于野外实测数据的反演解释,一维反演仍然占据着重要地位。提出一种使用探测深度和等对数域剖分的方式来剖分层厚、只进行电阻率单参数反演的方案。反演时,先使用MT反演来建立初始模型,之后采用MATLAB自带的parfor循环和最优化工具箱来进行最优化计算,使用解析法来计算雅克比矩阵,从而大大提高了计算速度。对于实际地电模型,充分考虑到已有的先验信息,通过使用井震约束反演有效降低了反演多解性,使反演结果更加接近真实的地层情况。  相似文献   

10.
构造应力场边界载荷反演的有限元逆逼近法   总被引:2,自引:0,他引:2  
付玉华  王兴明  袁海平 《岩土力学》2009,30(6):1850-1855
构造应力场问题的边界条件、加载方式是地应力场分析的难题,由于构造运动未知和地质构造复杂,无法直接求解。边界力反演是进行构造应力场研究的重要手段。从弹性力学基本方程出发,针对油气储层地应力场反问题,结合对区域构造应力场认识,利用阻尼最小二乘法建立了应力场反演的约束模型,提出了一种用于反演构造应力场边界力的反分析方法,该方法根据部分实测点的应力值进行边界力反演,并与有限元方法和有限差分法相结合,以边界力为参量逼近目标函数。根据研究区少量地应力实测资料采用该方法进行最优化计算,使得计算应力场与实测应力达到最优拟合,以弥补目前的回归反演方法和边界载荷调整法的不足来提高计算精度。算例表明了该方法的有效性,可用于解决构造应力场边界力的反演问题。  相似文献   

11.
提出了计算渗流场流量的虚拟单元法,推导了虚拟单元法有限元方程,并探讨了虚拟单元法的有限元实现方法,建立了虚拟单元法求解渗流流量的理论和方法体系。二维径向渗流和平底坝基算例比较分析表明,虚拟单元法具有计算精度高及有限元实现简单并容易编写通用程序等优点。虚拟单元法与中线法计算流量具有相同的精度,但在应用上虚拟单元法有一定优势。  相似文献   

12.
基于FLAC3D的饱和-非饱和渗流分析   总被引:1,自引:0,他引:1  
李毅  伍嘉  李坤 《岩土力学》2012,33(2):617-622
当使用基于有限差分法的FLAC3D计算渗流场时,自由面位置过高,非饱和区渗透系数被定义成常数,这与实际情况不符。运用变单元渗透系数法,对其进行二次开发,并提出了修正的变单元渗透系数法,对非饱和区的渗透系数进行较为严格的理论修正,得出准确的连续光滑的自由面,求解出饱和-非饱和渗流场。利用有限元法与该方法进行对比,论证了该方法的可行性和有效性。提出了附着接触面的方法,解决了三维复杂边界中寻找溢出点的困难;并以实际的三维复杂工程为例,计算了某边坡的渗流场分布情况,经验证,结果准确。研究成果对研究饱和-非饱和渗流以及该软件中与渗流有关的计算具有有指导意义。  相似文献   

13.
土质边坡非饱和渗流场与应力场耦合数值分析   总被引:2,自引:0,他引:2  
阐述了坡面径流-非饱和渗流分析与应力场耦合分析的重要意义和作用;介绍了在非饱和渗流分析中具有重要作用的土-水特征曲线,包括考虑固结压力条件下的土-水特征曲线,并介绍了部分试验成果以及坡面径流-非饱和渗流耦合分析的计算方法;提出坡面径流-非饱和渗流分析与应力场的耦合计算方法,并编制相应的有限元程序;通过算例对比了耦合与非耦合情况下应力场、渗流场的差异,说明耦合分析更符合实际,且所提出的方法是切实可行的。  相似文献   

14.
当滑体或潜在滑体内有倾斜的相对弱面或端部条块外侧面倾斜且有水压力作用时,采用一般条分法进行稳定性分析应比采用垂直条分法更为恰当。本文研究了对一般条分法也适用的渗流作用简化计算方法,首先给出土条渗透力的合理计算式,再根据水的面力与体力的关系给出土条周边孔隙水压力计算方法。同已有的计算方法相比,本文方法的合理性显著提高。  相似文献   

15.
李树忱  王兆清  袁超 《岩土力学》2013,34(7):1867-1873
岩土体的渗透破坏、地下工程的防渗设计等无不与渗流计算有关。针对渗流自由面问题,提出一种重心拉格朗日插值的配点型无网格方法。由于渗流自由面问题的求解区域是不规则区域,该方法通过将不规则求解区域嵌入一个正则矩形区域,在正则区域上采用重心拉格朗日插值近似未知函数,利用配点法离散渗流问题的控制方程,将重心拉格朗日插值的微分矩阵离散成代数方程表达的矩阵形式。将自由面上的边界条件通过重心拉格朗日插值离散,通过置换方程法和附加方程法施加边界条件,利用正则区域上的重心插值配点法,通过迭代确定最终自由面的位置。数值算例表明所提出的无网格方法对于求解渗流自由面问题的正确性和高精度。  相似文献   

16.
拼缝渗漏水是盾构隧道最常见的病害表现形式,对衬砌的服役性能有重要影响。在前人研究成果的基础上,提出基于防水密封垫失效的拼缝有效开度假定,进而推导出基于拼缝有效开度的盾构隧道渗流量及衬砌等效渗透系数的实用计算公式。以国内较为通用的常规断面盾构隧道为例,利用上述实用公式估算其渗流量,并与基于块体离散元的数值模拟结果相互验证。再以3种不同几何形态的盾构隧道为原型,完成了135种不同工况下的渗流计算,探讨了有效开度、衬砌外水压及衬砌几何形态,对盾构隧道渗流量及衬砌等效渗透系数的影响。结果表明:拼缝有效开度对渗流量及等效渗透系数的影响十分显著,未考虑拼缝间密封垫失效的传统计算方法可能高估了盾构隧道衬砌的防水性能。本文所述实用计算公式,具有表达清晰、计算简便的特点,可为盾构隧道衬砌的防水设计及后续服役性能评价提供定性和定量的参考。  相似文献   

17.
地下工程渗流排水孔数值模拟的隐式复合单元法   总被引:1,自引:0,他引:1  
倪绍虎  肖明 《岩土力学》2008,29(6):1659-1664
根据渗流基本原理,提出了排水孔在地下工程中渗流数值模拟的隐式复合单元法。将排水孔隐含于模型网格中,得到排水孔在模型中的数据信息。判断所有排水孔与单元的相对位置,通过修正排水孔所穿过单元的渗透传导矩阵来模拟排水孔的强排水效果,并通过实例对排水子结构法和隐式复合单元法进行了比较。运用编制的三维有限元计算程序对某地下水电站厂房渗流场及排水孔进行了模拟,结果表明,采用隐式复合单元法对排水孔进行模拟是可行的。  相似文献   

18.
基于非饱和土固结理论的有限元强度折减法   总被引:1,自引:0,他引:1  
周桂云  李同春 《岩土力学》2008,29(4):1133-1137
库水位下降使库岸坡体内产生饱和-非饱和非稳定渗流,非稳定渗流作用是边坡失稳的重要原因之一。为了分析非稳定渗流对边坡稳定性的影响,首先推导了固、液两相孔隙介质的固结方程,考虑了土体渗流与变形的耦合作用。在此基础上结合有限元强度折减法求解边坡稳定安全系数,将渗流、变形及稳定分析采用一套统一的有限元方法。并通过算例分析了库水位骤降情况下,坡体的渗透系数、水位降落比对稳定安全系数的影响,计算结果表明,所提的理论和方法是有效可行的,为饱和-非饱和非稳定渗流作用下边坡稳定问题的分析提供了实用工具。  相似文献   

19.
李杰  李强 《地质力学学报》2001,7(3):271-277
在水电、石油、交通等部门的地下工程中,经常会遇到地下水的渗流问题,如果处理不当,可能诱发应力恶化而影响结构稳定。本文以非饱和介质的渗流模型为基础,实现了地下洞室开挖过程中的非稳定渗流场解析。
按分期施工方案,在保持渗流场变化连续性的前提下,对各施工阶段进行了准稳定渗流分析,并在此基础上建立了地下工程非稳定渗流数值模型。其中,对于洞室系统,基于水量平衡原理,采用了子结构模型;对于防渗帷幕,采用了隔水层模型;对于排水井列,采用了杆单元模型。
根据上述模型,采用有限差分法作为数值模拟方法,编写了地下工程非稳定渗流场分析的有限元程序,实现了复杂渗流区域的渗流场的动态模拟。工程实例计算分析表明,用该模型计算的结果是合理的。   相似文献   

20.
降雨条件下岩土饱和-非饱和渗流分析   总被引:4,自引:1,他引:4  
刘建军  熊俊  何翔 《岩土力学》2004,25(Z2):559-563
基于Buckley-Leverett两相渗流方程,提出了新的岩土饱和-非饱和渗流数学模型,利用有限差分方法给出了隐式压力显示饱和度数值求解方法,编制了饱和-非饱和渗流计算程序.结合工程实例,模拟了降雨入渗条件下边坡岩体渗流场孔隙压力变化和含水饱和度变化,模拟结果验证了所提出的模型对饱和-非饱和渗流的有效性.  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号