共查询到18条相似文献，搜索用时 218 毫秒
1.
频谱法是瞬变电磁法正演计算中常用的计算方法之一,其关键在于内层含贝塞尔函数的积分计算。对于重叠回线装置的瞬变电磁法正演的内层积分,以往采用在积分区间内寻找贝塞尔函数的零点分布,依次在相邻零点之间采用一般的数值积分法在求得各自的积分值后再叠加的方法进行计算。这种方法精度较高,但效率低。利用贝塞尔函数的大宗量渐近特性,提出一种新的计算方法。计算结果表明,该方法计算效率高、方法简单、精度可靠。 相似文献
2.
多层介质点电流源磁场响应计算 总被引：1，自引：0，他引：1
介绍了点电流源层状介质模型磁电阻率法响应的模拟方法。由电流磁场的旋度公式出发,推导了多层介质中点电流源磁场的表达式,其为一个关于空间径向坐标和贝塞尔函数空间波数的二重积分。利用贝塞尔函数的性质可将该二重积分化为标准的贝塞尔积分,然后采用直接数值积分的方法计算该积分。在此过程中通过连分式展开加快积分求和的收敛速度。利用上述方法对海底天然气水合物模型点电流源产生的磁场进行了计算。模型计算结果表明,海底天然气水合物的存在能引起磁场较为明显的相对变化。 相似文献
3.
利用传输线理论、Sommerfeld积分快速计算以及最小平方拟合技术研究建立多分量感应测井数据的一种新的快速参数化迭代反演算法,同时重构水平层状横向同性地层的纵、横向电阻率以及水平层界面深度。首先,通过Fourier变换与传输线理论给出频率波数域中电磁场并矢Green函数在各个地层中的解析解,并利用三次样条插值和贝塞尔函数递推公式建立Sommerfeld积分的半解析算法,快速计算多分量感应的测井响应。然后在此基础上,利用摄动理论建立磁场并矢Green函数与模型向量间变化关系的摄动方程,并将摄动方程中各个积分转化为Sommerfeld积分,实现正演模拟的同时用半解析算法快速确定多分量感应测井响应的Fréchet导数。最后,利用归一化处理和奇异值分解技术,同时反演所有地层的纵、横向电阻率和层界面深度,实现输入数据和反演模型的模拟数据优化拟合。理论模型的数值结果验证了该反演算法的有效性及抗噪性。 相似文献
4.
导线环垂直激发场高精度数值模拟 总被引：1，自引：0，他引：1
导线环垂直激发场是导线环在地下产生的激发场中垂直于地磁场的分量,在地面核磁共振找水方法中,垂直激发场在地下的分布直接决定地面核磁共振响应,探讨影响垂直激发场地下分布的诸因素进而研究这些因素对地面核磁共振响应的作用,对核磁共振找水正确野外施工具有一定的指导意义,有助于采集更高质量的野外数据.利用经典Chave算法对含有贝塞尔函数的强振荡积分核进行了数值积分试验,将其计算结果和对应的解析解进行了对比,试验结果显示该算法用于强震荡函数无穷积分的计算具有良好的收敛性、较高的精度和较快的计算速度,采用Chave算法对含有贝塞尔函数的强振荡函数进行积分是可行的.在此基础上,进一步计算并数值模拟了导线环垂直激发场及其相位,分别从地磁倾角、磁偏角、地下介质导电性、深度等方面探讨了影响垂直激发场及其相位分布的主要因素. 相似文献
5.
最近,国内一些单位在物探中应用PC—1500袖珍计算机做了不少工作。在电测深方面文献[1]、[2]提供了电阻率测深及激发极化测深的方法和程序。本文给出了在PC—1500机上实现的电测深曲线T函数自动拟合反演解释方法和程序,供大家参考。一、理论基础众所周知,在水平层状大地表面上对称四极排列的视电阻率表达式为 p_s=r~2 intergral from n=0 to ∞ T_1(λ)J_1(λr)λdλ (1)式中r=AB/2,为对称四极排列的电极距;J_1(λr)是一阶贝塞尔函数;T_1(λ)为视电阻率转换函数(简称T函数),可由下列递推公式求得: Tn(λ)=pn, 相似文献
6.
光学、电磁学和地震学都需要计算Hankel变换。除少数核函数的Hankel变换有解析式,绝大多数的通常采用数值积分的方法对其近似计算。虽然数字滤波法是常规的计算方法,但由于滤波系数的不统一,导致计算结果并不完全一致,甚至在某些情况结果是错误的。为了解决这个问题,采用直接数值积分的方法计算Hankel变换。具体过程是:首先将Hankel变换分解为积分限从零到Bessel函数的第一个零点值的积分段和后续各相邻零点值作为积分限的积分段的和,然后用高斯求积计算每段积分。对高斯求积结果组成的序列采用连分式求和,无论是计算结果的精确性还是收敛速度都要优于直接求和。这里首先详细阐述了联用高斯求积与连分式求和计算Hankel变换的基本原理和相关算法,然后用Matlab编程对比了直接求和与连分式求和的计算效果,最后将此方法应用于电偶源在均匀半空间激发的地下场计算的两个实例中。应用表明:将高斯求积与连分式求和进行联用是计算Hankel变换的一种有效方法,可以广泛应用于地球物理应用中。 相似文献
7.
水平井开发特殊类型油气藏有其独到优势,针对双重孔隙介质水平井复杂渗流,以Warren-Root双重介质模型为基础,建立了水平井三维不稳定渗流模型,采用Laplace变换、有限Fourier余弦变换及叠加积分等方法对模型进行求解,获得定产条件下无量纲井底压力解,再利用Duhamel褶积原理将定产条件压力解变换为井底定压条件下的产量解;通过重新定义无量纲量建立McCray型产量积分平均递减方程,结合Stehfest数值反演,获得无量纲产量积分平均递减曲线及其导数曲线图版。计算结果表明,产量积分平均导数曲线上出现双重介质特有的“凹子”, 应用该类图版拟合生产数据可方便进行产量预测及储层动态参数求取。 相似文献
8.
数值积分是基于节点计算的数值分析方法的关键技术之一。针对完全无网格法的计算要求，提出了一种能自动适应任意计算域上各种节点分布方式的数值积分算法。该算法能随计算点的位置不同，自动确定积分域及积分域内的求积点，计算出数值积分结果，给出了临空边坡滑移体和圆孔方板2个实例。结果表明，按文中算法编制的程序具有良好的边界再生能力和较高的计算精度，可作为基于节点数值分析方法的一个程序模块，用于实际计算。底层算法思想，使得程序可方便地转换为其他计算语言。 相似文献
9.
概率积分法是煤矿开采沉陷预测的重要方法,可适用于缓倾斜、倾斜煤层开采地表移动分布曲线符合正态分布规律的地表移动和变形计算。针对原方法对任意形状工作面存在的不足,在原有走向和倾向概率积分函数的基础上,以走向和倾向为基准划分积分函数区域,在该区域内剖分非结构化的三角形单元作为二重积分的基本单位;通过坐标旋转变换的二重积分换元法,转化为以直线为上、下限的二重积分,计算旋转后的影响半径,对新的二重积分采用复合辛普森数值积分公式,积分得出地表任意点在基本三角形开采单元下的积分下沉值;最后,通过叠加计算完成任意点、任意形状工作面概率积分法的沉陷预测计算。算法在GIS系统下实现,实例应用效果良好,可为“三下”采煤提供地表移动预测计算的支持。 相似文献
10.
应用迭代法的困难之处在于：一是迭代公式是否存在导数运算；二是初始值选定是否影响迭代公式的收敛性；三是收敛快速和达到需要精度等问题。我们获得了求方程根不用计算导数的平方收敛迭代公式，并设计了求根的大范围收敛算法，编写了C^ 语言程序，进行了算法和数值分析。与其它算法比较，该算法具有无导数计算、初值任意选定、平方快速收敛、大范围收敛和双精度控制(根的精度和函数值的精度控制)等优点。 相似文献
11.
体积分方程法模拟复电阻率三维体电磁响应 总被引：1，自引：0，他引：1
利用体积分方程法计算了均匀半空间中复电阻率(激电)效应和电磁效应同时存在时的三维体响应。在计算中,对于需要计算三重积分的张量格林函数电荷项一次部分,应用一种差分近似的方法求解,这种方法在保证计算精度的同时更加便于计算机实现;采用二次剖分的算法解决了计算张量格林函数时的奇异值问题;计算含有贝塞尔函数的积分项时,利用一种结合连分式展开的高斯求积代替常规的快速汉克尔变换方法。验证了计算结果并分析了三维体复电阻率模型(Cole-Cole)参数对正演结果的影响,为三维体复电阻率及其参数反演提供了正演依据。 相似文献
12.
NK8310螯合树脂分离富集地质样品中痕量金银铂钯及其测定 总被引：14，自引：3，他引：14
研究了硫脲螯合树脂（NK8310）分类富集地质样品中痕量Au、Ag、Pt和Pd的实验条件。在φ=10％的王水介质中，[AuCI4]^-、[AgCI2]^-、[PtCI6]^2-和[PdCI4]^2-定量吸附于树脂上并与大量贱金属分离；用5g/L硫脲－0.12mol/L HCI溶液洗脱Au、Ag、Pt和Pd，回收率为97％－104％。用硫镍矿管理样以及国家一级标准物质进行分析验证，分析结果与推荐值及标准值吻合，表明NK8310螯合树脂适用于地质样品中Au、Ag、Pt和Pd的分离富集。 相似文献
13.
一般的商业软件只提供比较简单的算法,能够对矢量数据几何纠正实现简单的、全局意义上的处理。研究以Delaunay算法为基础,提出了基于二维平面剖分的、具有局部意义的矢量数据几何纠正算法。在Visual Studio2005开发环境下,使用C#2.0与C 分别实现相关功能,完成了算法的实现。文中提到的解决方法在东北地区2005年土地利用数据库更新项目的空间数据集成问题中得到了成功应用。 相似文献
14.
动力响应问题的求解对于地基在外荷载作用下引起的弹性波动问题研究有重要的意义。本文提出了一种求解横观各向同性层状地基在施加时间简谐荷载作用下任意点的应力响应的算法。此算法利用傅里叶变换将广义平面应变问题频率-空间域的动力方程转化到频率-波数域内,结合对偶变量的引入,利用高精度的精细积分算法对状态方程进行求解,在得到频率-波数域内的位移响应的基础上,利用傅里叶逆变换得到任意点的动应力响应。简谐荷载不仅可以施加在地基表面,而且可以施加在地基内部。对比算例验证了本文算法的准确性,同时对地基各向异性特性、激励频率和阻尼比对动应力响应的影响进行了参数分析,为工程实际提供可靠的数值依据。 相似文献
15.
地球自转速度在长期的减慢,而内核自转 又较其外部的地球部分转动快,这是由日 、月对地球的潮汐摩擦引起的同一动力学背景下导致的相关现象。通过能量守衡的定量分析 ,得出地球潮汐耗散能的11% 即可导致内核自转较其外部的地球部分自转快的结论。通过推 理比较和定量分析,提出了岩石圈层、主地幔圈层与内核三个圈层的差异旋转模型：岩石圈 自转速度<主地慢自转速度<内核自转速度。它可能是导致大陆向西飘移-板块构造及地磁产 生的重要机制之一。 相似文献
16.
Snow avalanches,which are widely and frequently developed at high elevations,seriously threatens the built traffic corridors in the Tibetan Plateau. Susceptibility evaluation of snow avalanche via machine learning model with a high forecast accuracy can be appled to quickly and effectively assess the regional avalanche risk. This paper took the central Shaluli Mountain region as the study area,in which the snow avalanche inventory was established through remote sensing interpretation and field investigation verification. We quantitatively extracted 17 evaluation factors via GIS-based analysis,and these factors were selected through the variance expansion factor（VIF）. Four machine learning models containing SVM,DT,MLP and KNN were used to compile the susceptibility index map of snow avalanches,and kappa coefficient and ROC curve were used to verify the accuracy. The results suggested that the susceptibility indexes obtained from SVM,DT,MLP and KNN were in the range of［0,0. 964］,［0,815］,［0,0. 995］and［0,1］,respectively. The accuracy test results show that these four models all have good prediction accuracy. Among them,the SVM model is the best. The results also indicated that the areas with the high snow avalanche susceptibility mainly distributed in Genie Mountain and Rigong Mountain,most of which were above the planation surface of the Tibetan Plateau. The average altitude of the extremely high snow-avalanche-prone areas is 4 939 m,while the average altitude of the high snow avalanche-prone areas is 4 859 m. The snow avalanche has low perniciousness on the Sichuan-Tibet Highway and the Sichuan-Tibet Railway in the study area. This study can provide theoretical basis and method reference for disaster prevention and mitigation of snow avalanche along Sichuan-Tibet Railway and other major projects across Shaluli Mountains region. © 2022 Science Press (China). 相似文献
17.
用广角度X射线散射 (WAXS)和由它得到的径向分布函数 (RDF)研究了偏岭石的结构。偏岭石是一种以非晶组份为主的部分有序的中间序性物质。构成偏岭石的结构单元是四面体 [MO4 ](M =Si,Al)和八面体 [AlO6](O =O,OH)。偏岭石结构中近邻四面体与四面体共顶连接构成M-O层,八面体与八面体共棱或共面连接构成Al-O层;而四面体与八面体通过不规则的共顶连接构成近程层状结构。它的近程结构将不随序性或结晶性变化。它的微结构及其随热处理温度的变化主要决定于结构基元之间的关系和序性。偏岭石是结构不同于变高岭石的一种粘土. 相似文献
18.
广域视电阻率与卡尼亚视电阻率相比,能更好地表征地下介质的电性变化规律,一般常使用迭代法和逆插值法求解,这里使用这两种算法对E—Ex方式广域视电阻率进行了计算.结果表明,使用这两种算法所得广域视电阻率,都正确反映了地下介质的电性变化规律,逆插值的影响因素较少,计算过程稳定.其不足之处在于要求插值目标函数单调,迭代法步骤简单,并存在迭代发散的情况.在实际应用时,可根据广域视电阻率公式的特点,选择适当的方法进行计算. 相似文献