直流电阻率2.5维有限元数值模拟

从稳定电流场满足的基本微分方程出发,推导了地面点源电场电位边值问题对应的等价变分问题.用网格离散积分区域时,以四边形为最小研究单元.针对电场在电源附近衰减快,变化梯度大,数值模拟结果在电源附近位置有较大误差问题,采用边界/奇异校正算法改善模拟结果在电源附近的奇异性问题,同时减少边界条件的影响,提高模拟结果精度.从波数域电位求取空间域电位,利用最优化离散波数进行反付氏变换.模型算例结果验证了该算法的可行性与有效性.     

基于有限单元法的二维／三维大地电磁正演模拟策略

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

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

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

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

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

通道控制的岩浆型矿床:以加拿大沃伊西湾铜镍矿床为例

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

荒漠土壤微生物碳垂直分布规律对有机碳库的表征作用

以古尔班通古特沙漠南缘原始盐漠为研究对象,测定不同深度的土壤有机碳和土壤微生物碳含量,以分析它们之间的响应关系。结果表明:(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种有机碳库的储量。     

地形畸变对偶极-偶极电阻率测深拟断面图的影响

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

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

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。结合雨洪法和形态调查法结果分析,此次泥石流重现周期为百年一遇。结合泥石流发育趋势,文章认为流域仍有可能暴发大型泥石流,提出了上游固坡、中游调控、下游排导的防治建议。     

改进的小波分频重构算法在石油地震勘探中的应用

在地震资料处理过程中,为了提高信号的信噪比以及保留数据的可靠频段,通常采用频率域滤波的方法.但是,传统的频率域滤波方法是基于傅里叶变换,在窗口与窗口之间的接合部位会出现假异常.且在滤除噪声的同时,很可能同时也将有用信息一起滤除掉. 本文针对上述问题,将小波域分频处理与重构的方法应用到石油地震勘探数据处理中.在提高信号信噪比的同时又可兼顾其它频段信息,使处理出的地震剖面信噪比提高,频带有所拓宽,从而达到最佳重构.文中先后给出了小波变换及其逆变换的数值计算方法以及最佳小波和最佳重构系数的求取方法.最后通过实例分析,其结果表明,该方法在理论上是正确的,实际应用效果也是成功的.