基于有限单元法的二维／三维大地电磁正演模拟策略   总被引：1，自引：0，他引：1  

徐凌华  童孝忠  柳建新  郭荣文  王浩 《物探化探计算技术》2009,31(5):421-425

对于二维和三维大地电磁正演问题,有限单元法最后形成了一个线性方程组KX=p。方程组中的K是大型稀疏的带状对称复系数矩阵,其条件数远大于1,为严重病态矩阵,求解其对应方程组会遇到很多困难。不完全LU分解处理的BICGSTAB算法,可用于该线性方程组的求解,并且具有速度快,精度高,稳定性好等优点。为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下,设计了非均匀网格剖分。在程序编制中,因只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间。通过对二维模型和三维模型电磁响应的计算,验证了该算法的正确性。  相似文献   

一种新型CSAMT观测参数的计算方法     

佟铁钢  柳建新 《物探化探计算技术》2009,31(3):210-212

常规的可控源电磁法理论在计算视电阻率公式上,由于主要采用了其电磁场,所以难以直接反映全频域视电阻率的值,也不能直观地显现出地下介质的地质构造。这里采用水平偶极子激发的电磁场,提出了电场的全域精确表达式,并可直接计算出大地电阻率。通过与卡尼亚电阻率计算的对比结果表明：该方法的结果在远区等价于卡尼亚电阻率,在近区和过渡带则明显地改善了卡尼亚电阻率的非波场区场畸变,从而能更好地接近基底的真电阻率,更形象地反映了地下介质的垂向电性变化。  相似文献   

拖曳锚安装缆绳在海床土中的反悬链形态研究     

柳成林  刘海笑  赵燕兵  王宸 《岩土力学》2013,34(11):3221-3228

由于在抗拔承载力和深水安装中的优势,拖曳锚在深水系泊系统中具有良好的发展前景。在拖曳锚安装过程中,安装缆绳和锚之间存在复杂的相互作用,安装缆绳在海床土中的反悬链形态将直接影响锚的嵌入和运动,对反悬链形态的研究有助于提高锚的嵌入性能、准确预测锚的运动轨迹以及解决锚在安装中的精确定位问题。基于嵌入缆单元在海床土中的力学模型,推导出适用于黏性土和无黏性土的反悬链方程,利用该方程可求出安装缆绳在海床土中的反悬链形态;通过在缆绳上布置倾角传感器并运用圆弧递推方法,设计出可实时测量嵌入缆形态的试验技术,为验证试验技术的合理性,在空气中分别选取3种典型缆绳形态进行模拟,结果表明：模拟形态与真实形态吻合良好,并具有较高精度;利用实时测量技术开展模型试验,对反悬链方程进行了考察,验证了其模拟嵌入缆反悬链形态的精度,并获取了对拖曳锚安装过程中反悬链形态变化规律的直观认识。  相似文献   

基于DDS技术的航空重力仪电容测微正弦激励源设计     

李文萍  柳林涛  吴鹏飞 《大地测量与地球动力学》2013,(Z2):84-87

介绍基于直接数字合成技术(DDS)研制的正弦激励源系统工作原理,运用matlab-simulink仿真模块仿真正弦激励源数据,结果表明,幅值稳定性可达到10-6~10-7,频率稳定性达到10-4,该系统能解决以往正弦电路中幅值稳定度和高频时频率稳定度差的问题,提高航空重力仪的观测精度。  相似文献   

基于全区视电阻率的线源FCSEM二维正演模拟     

李杰  柳建新  童孝忠  刘颖 《地球物理学进展》2012,27(1):343-348

频率域可控源电磁是在大地电磁测深的基础上发展起来的一种人工源电磁测深法,其二维电磁响应的计算须采用数值模拟方法.本文以Matlab为程序编译工具,采用双二次插值的有限单元法,推导出相应的计算公式.为了模拟无穷远边界及满足计算机的内存需求,在保证计算精度的情况下设计了非均匀网格剖分.在程序编制中,只存储有限元系数矩阵的非零元素,大大减少了正演计算的时间.针对频率域可控源电磁法中卡尼亚电阻率在过渡区和近区畸变的问题,给出了全区视电阻率的迭代公式,并对典型的一维层状模型以及简单二维模型进行了计算.过渡区和近区数据经过校正后,可以正确反映出模型的地电特征,证明了线源下近区勘探的可能性.  相似文献   

通道控制的岩浆型矿床:以加拿大沃伊西湾铜镍矿床为例          下载免费PDF全文

定立  辛江  杨永强  黄磊  柳建平  侯玉树  王晓亮  周聪  高知睿 《地质科技通报》2013,32(5):58

加拿大沃伊西湾铜镍矿床是20世纪末加拿大发现的最大矿床,而且近年还有新增储量的发现。该矿床储量不仅在世
界同类矿床中排在前几名,而且是岩浆通道模式最经典的矿床。岩体控制矿体,该基性侵入体位于加拿大拉布拉多半岛(Labrador)
北部中元古代奈恩深层岩套(Nainplutonicsuite)中,岩体主要由西帝普斯岩体和东帝普斯岩体组成,它们由一条隐伏的岩浆
通道所连接,矿体就位于岩浆通道中,形成典型岩浆通道成矿系统。岩浆通道的膨大部分形成了奥韦尔德矿床。对该矿成矿特征
及其成矿模式的总结,对于寻找与岩浆通道成矿有关的矿床具有指导意义。   相似文献   

荒漠土壤微生物碳垂直分布规律对有机碳库的表征作用   总被引：1，自引：1，他引：0  

房飞  胡玉昆  公延明  杨秀娟  柳妍妍 《中国沙漠》2013,33(3):777-781

以古尔班通古特沙漠南缘原始盐漠为研究对象,测定不同深度的土壤有机碳和土壤微生物碳含量,以分析它们之间的响应关系。结果表明:(1)在土壤垂直剖面上,土壤微生物碳（SMC）含量与有机碳（SOC）含量呈现极显著正线性相关（R2=0.63,p=0.0003）。（2）SMC出现了2个明显的改变界面（20 cm,80 cm）,0~20、20~80、80~500 cm值分别为:2.24~3.06、0.19~0.72、0.0017~0.0097 mg·kg-1;0~20 cm和20~80 cm的SMC差异极显著（p<0.0001）,20~80 cm和80~500 cm的SMC差异显著（p<0.05）。（3）对应于SMC的土壤层划分,SOC在0~20 cm、20~80 cm和80~500 cm同样具有一定的分层性。（4）我们把具有不同微生物活性的有机碳层分别定义为活性、惰性、稳定性有机碳库,土壤垂直剖面上微生物碳的分布很好地表征了土壤中活性、惰性、稳定性有机碳库的分布;通过对这3种碳库所在土层进行合理划分,可以定量分析土壤中3种有机碳库的储量。  相似文献   

海拉尔盆地贝尔凹陷伸展断层转折褶皱作用及其对沉积作用的制约     

刘志宏  梅梅  柳行军  吴相梅  万传彪  林东成  高军义 《吉林大学学报(地球科学版)》2012,42(5):1330-1337

海拉尔盆地是叠置于内蒙－大兴安岭古生代碰撞造山带之上的中新生代盆地,总体上呈北东向展布,由3坳、2隆5个一级构造单元组成,自西向东依次为扎贲诺尔坳陷、嵯岗隆起、贝尔湖坳陷、巴彦山隆起和呼和湖坳陷, 其中贝尔凹陷是贝尔湖坳陷南部的1个二级构造单元。贝尔凹陷的控陷断层为1个上凸、2个上凹断层转折组合构成的犁式正断层,在深部发育伸展双重构造。在犁式正断层上盘由于上凸断层转折产生的活动轴面与断层面为同旋向剪切,而上凹断层转折产生的活动轴面与断层面为反旋向剪切,在断层上盘形成2个背斜中间夹1个向斜的构造组合。由活动轴面与固定轴面之间滚动面的宽度确定控陷断层F₁在早白垩世南屯期晚期的水平伸展量为6 850 m。贝尔坳陷伸展断层转折褶皱的几何学和运动学特征,反映了断层滑移速率、构造沉降速率和沉积速率对凹陷内沉积地层变形特征的影响和构造活动对沉积作用和油气地质条件的制约。贝尔凹陷南屯组上段内部角度不整合的形成与变形作用中的剥蚀作用和沉积间断无关,而是由于半地堑从略欠补偿充填向过补偿充填状态转变时,沉积速率相对于构造沉降速率显著增加所致。  相似文献   

地形畸变对偶极-偶极电阻率测深拟断面图的影响   总被引：4，自引：3，他引：1       下载免费PDF全文

吴小平  柳建新  段无悔  蓝燕 《地球物理学进展》2010,25(2):635-641

偶极－偶极剖面方法是一种非常流行的物探方法,但就地形对该方法的实质性影响及修正研究得比较少.我们就几种典型的地形畸变情况进行了正演计算,并绘制了拟断面图.通过对这些图形分析,可以看到在地形畸变处的拟断面图不能反映地下地质体的情况,且有效拟断面图是不连续的,需要分段处理,分级解释.  相似文献   

四川木里县项脚沟“7·5”特大型泥石流特征及发展趋势分析     

杨华铨  柳金峰  孙昊  赵万玉  张文涛 《中国地质灾害与防治学报》2024,(1):100-107

2021年7月5日,凉山州木里县项脚沟暴发特大型泥石流灾害。文章通过野外调查和特征参数计算,研究了本次泥石流的形成条件、暴发过程和暴发特征。根据实地考察判断,泥石流主要在森林火灾、短时强降雨和沟道地形条件等因素共同作用下暴发,为暴雨径流冲刷引发的火后泥石流。暴雨形成的洪水冲刷坡面、侵蚀沟道,导致沟道两岸坍塌滑坡,堵溃效应明显,泥石流规模扩大。项脚沟过火面积达74.61%,高烈度火烧区面积达57.98%,泥石流临界降雨强度为77.84 mm/h,累计降雨达141.60 mm。根据泥石流特征参数计算结果,此次泥石流密度范围为1.83～1.93 g/cm³,属黏性泥石流,主沟下游出口流速为7.22 m/s,峰值流量759.08 m³/s。结合雨洪法和形态调查法结果分析,此次泥石流重现周期为百年一遇。结合泥石流发育趋势,文章认为流域仍有可能暴发大型泥石流,提出了上游固坡、中游调控、下游排导的防治建议。  相似文献