首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 250 毫秒
1.
本文构造矩形有限元的概率算法,应用这种方法可以直接求出有限元解在一个或少数几个网格点的近似值。  相似文献   

2.
倾斜断层Vzz、Vzzz的联合反问题解   总被引:1,自引:0,他引:1  
顾心如 《物探与化探》1994,18(5):353-362
本文利用两个上延高度的Vzz极值距和其中之一高度的Vzzz的零值点获得了倾斜断层的位置、埋深、断距、倾向、倾角以及断层正逆性质等的反问题解。  相似文献   

3.
采用三角网格剖分的有限元法,研究了2.5维航空瞬变电磁法正演模拟问题。利用时频变换数值方法将时间域电磁场转换到拉氏域,再利用傅里叶变换将三维问题降维变为2.5维问题,然后由有限元法求解得到拉氏域二维电磁场,逆拉氏变换后得到时间域航空瞬变响应。为了回避正演模拟中总感应磁场在场源处的奇异性问题,采用异常场算法,场源响应通过在微分方程中施加背景电磁场实现。由于瞬变电磁信号具有较大的动态范围,而且需要经过两次正、逆拉氏变换和傅里叶变换,每个环节的计算精度和速度要严格控制在较高的水平上,否则积累误差会非常大。模型计算表明均匀大地和层状大地模型解析解与数值解吻合很好。这证明该算法是正确可行的,可作为研究二维复杂地质体的方法手段。  相似文献   

4.
利用G-S逆拉氏变换法计算瞬变测深正演问题   总被引:3,自引:0,他引:3  
求解时间域电磁法问题,通常有两种途径。第一,直接求解时间域问题;第二,先解相应的频率域问题,然后变换为时间域问题。本文采用第二种方法,利用Gaver-stohfest逆拉氏变换法将频率测深正演问题转化为对应的瞬变测深正演问题。实际计算表明,该方法是将频率域问题转化为时间域问题的有效方法。  相似文献   

5.
叠前弹性波逆时深度偏移及波场分离技术探讨   总被引:3,自引:0,他引:3  
这里探讨了Sun和McMechan提出的叠前弹性波标量逆时深度偏移方法,即在地表附近,对地表接收到的弹性波波场分量进行波场逆时延拓,然后分别对在波场延拓过程中通过波场分离得到的纯纵、横波反射波场进行偏移。在模型中,每个网格点的纵、横波成像时间,为震源到该网格点的纵波初至时间。利用纵波速度模型和声波传播方程的有限差分解,对反射纵波进行逆时延拓和成像,利用横波速度模型和声波传播方程的有限差分解,对转换横波进行逆时延拓和成像,并在偏移前进行极性校正。然后,在Sun和McMechan偏移成像方法的基础上,提出了直接偏移成像方法,即对地表接收到的弹性波波场分量在整个模型范围内进行波场逆时延拓。在波场延拓过程中,对符合成像条件的网格点进行波场分离、偏移成像和转换波极性校正。最后,通过数值模拟对二种偏移成像方法的结果进行了比较。  相似文献   

6.
三角网格模型剖分方法的研究与应用   总被引:1,自引:0,他引:1  
数值波场正演在地震学和勘探地球物理学领域正得到越来越广泛的应用,用一个好的方法来建立一个复杂的地质模型则显得尤为关键。三角网格模型相对于层状模型和矩形网格模型在反映地质界面形态、三角网格剖分个数调节及射线追踪速度等方面具有明显优势。利用三角网格剖分的建模方法,在VC2008编程环境下实现了复杂地质模型的描述。通过对复杂模型的建模测试,表明采用三角网格剖分方法,可解决以下模型难题:①对复杂模型界面的描述;②对复杂逆断层的描述;③对封闭块体的描述;④对模型弹性参数进行描述。  相似文献   

7.
大地电磁(MT)数值模拟中通常使用有限单元法,通过伽辽金(Galerkin)法将微分方程转化为与其等价的泛函形式,对泛函求取极值并在单元上定义插值基函数,得到节点上电磁场值的线性方程组,最终形成大型复对称稀疏矩阵。要达到较高的有限元计算精度,一般采用密集的网格或高次插值的方法,这样做大大的减慢了正演的速度。结合两者的优点利用三次插值和h-型自适应相结合的有限元法来实现MT的正演算法。首先从一个粗网格出发并利用三次插值,通过后验误差估计方法局部加密网格,在计算量较小的情况获得较高的计算精度。这种方法可以针对目标区域和介质分界面发生突变处进行网格加密,不需要全局加密网格。最后通过对国际标准模型COMMEMI-2D1的模拟,分别比较二次插值与三次插值的自适应网格数量和数值模拟结果,证明了三次插值自适应有限元算法的可行性。  相似文献   

8.
地下工程建设过程中,地下水分布以及岩体渗透参数确定至关重要。基于对解决欠定反分析问题的目标函数分析,通过利用适应度函数、地质统计原理和变差函数等优化,得到了解决地层渗透系数一类水文地质问题的罚函数,为解决欠定的反分析问题提供了一个优化评判准则,并结合粒子群算法建立优化反分析数学模型;提出利用压水试验钻孔水头变化反分析得到地层渗透系数的新方法。对南京地铁上元门车站,进行现场压水试验得到钻孔水头高度,利用经过罚函数优化的反分析数值计算获得地层渗透系数,对压水试验区域的电阻率跨孔CT探测验证了该方法区域上的正确性,该区域的钻孔原位渗透系数测试验证了该方法在数值上的准确性。结果表明,经过罚函数优化的反分析计算的准确率达到90%,说明该方法有助于得到全面的地层水文地质信息,对后续的治理具有指导作用,希望能对类似工程有一定的借鉴意义。  相似文献   

9.
时间域有限元算法已被广泛应用于航空电磁三维正演模拟当中,然而由于航空电磁测区面积大,且采样密集,造成正演计算量巨大,传统的串行算法已经无法满足计算效率要求,为此,开展了并行加速算法研究以解决计算效率不足的问题。基于航空电磁系统的影响范围有限,采用局部网格技术将计算任务划分成多个子网格,即每个发射源一套网格,各网格的正演计算相互独立,不存在数据依赖性,具有很好的可并行性;利用MPI技术对多个子网格正演任务进行分配,在各个进程上进行并行计算;针对每个正演子网格,在进行时间域有限元算法正演模拟过程中,采用OpenMP技术对单元矩阵进行并行计算。典型地电模型的数值模拟结果表明本文开发的MPI+OpenMP并行正演算法可以有效提高正演速度,最高加速比可达10倍。  相似文献   

10.
流域虚拟仿真中水沙模拟时空过程三维可视化   总被引:1,自引:1,他引:0       下载免费PDF全文
将计算机图形学领域的粒子系统引入数字流域研究,研发了基于通用网格的水沙模拟时空过程三维可视化系统.该系统建立数据结构描述粒子信息和计算或者虚拟网格单元信息传递关系;采用网格控制粒子随机生成;提出局部搜索定位方法来快速确定粒子位置;通过反距离加权平方插值提取水流信息;利用自适应龙格-库塔方法实现粒子的动态变化;设计了附有更多水流属性的水体粒子;研究了流场变化过程中粒子分布不均的问题.将上述方法应用于三峡与葛洲坝间河道水流的动态可视化模拟,验证了该系统的有效性和实用性.  相似文献   

11.
李兰 《水科学进展》1999,10(1):7-13
根据逆边界逆动态控制理论,将河流水污染动态控制问题提为逆边界逆动态混合控制问题。针对多个或单个污染源排放浓度和排放总量计算,提出了一维对流-扩散方程逆控制的精确算法。该方法与现行最优控制方法相比,其优点是充分考虑了河流沿程的稀释混合容量,并能充分考虑水质动态标准和社会经济变化等因素,可获得动态控制精确解的近似解。  相似文献   

12.
讨论了一类二阶抛物型方程反问题的数值解法。应用拟解法的思想,把原问题分解为一系列适定的正问题和一个不适定的线性代数方程组。对于相应的正问题,证明了解连续依赖于初始分布,由此得到了在时刻的稳定性估计。使用古典欧拉差分格式求解正问题和用截断奇异值分解法求解病态方程组。数值结果显示数值解与理论解吻合很好。  相似文献   

13.
The paper presents the results of a joint inversion of magnetotelluric sounding (MTS) curves taken at several stations without preliminary selection and normalization and distorted by the presence of lateral electrical-conductivity inhomogeneities in the medium. In the calculations, we used synthetic MTS data for a three-dimensional model. Preparation and interpretation of data are carried out by the Trefftz method using a numerical model of the field and MTS curve distortions. To solve the inverse problem and optimize the subsurface model, we used a nonlinear least-squares method and an iterative process with calculation of the sensitivity matrix and its singular decomposition. The target functional is determined by the discrepancies between the model and synthetic experimental apparent-resistivity curves corresponding to the elements of the impedance tensors on the lateral diagonal. The reliability of the reconstructed subsurface model is characterized by the dispersion of the deviations of its parameters from the parameters of the known model used for the preparation of synthetic experimental data. The joint consideration of distorted apparent-resistivity curves at several stations increases the reliability of interpretation results. The obtained solution to the inverse problem is approximate and can be used as a starting model for more complex algorithms and programs.  相似文献   

14.
线性方程组的正则化解法及其应用   总被引:3,自引:0,他引:3  
在重力资料反演解释中,正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子,求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中,要求对所求的解答预先有一个初步的推断解释。为此,要求对所解释的问题先进行一番研究,并对综合利用各种方法和手段所获得的先验资料进行对比,提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面,与正则化法和其它线性方程解法的结果作对比,说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。  相似文献   

15.
焦荣昌 《物探与化探》1994,18(3):200-208
在重力资料反演解释中,正则化法解线性方程组是一个行之有效的方法。这种方法的特点是引入正则化因子,求泛函极小和构组一族正则化近似并从中寻求最佳解。在求解过程中,要求对所求的解答预先有一个初步的推断解释。为此,要求对所解释的问题先进行一番研究,并对综合利用各种方法和手段所获得的先验资料进行对比,提出比较合适的基望值。本文给出了一个参数为已知的理论模型体的横剖面,与正则化法和其它线性方程解法的结果作对比,说明正则化法的有效性。文章也给出了利用正则化法对横穿吕宋岛、马尼拉海沟的一段地质地球物理剖面以及中国南黄海海域内─地质地球物理调查剖面进行反演解释的应用效果。  相似文献   

16.
岩土工程优化反分析是一个典型的复杂非线性函数优化问题,采用全局优化算法是解决这个问题的理想途径。结合ABAQUS有限元软件,提出遗传算法与有限元联合反演法,将有限元程序作为一个单独模块嵌入到遗传算法程序中,以测点的实测值与计算值建立误差函数,编制了遗传算法反演分析程序。并给出应用实例验证了该法的有效性,表明该方法可应用于岩土工程中的反演分析工作。  相似文献   

17.
稀疏脉冲反演实际上就是利用反褶积原理, 从带有噪声的地震道中计算出具有稀疏分布特征的反射系数的振幅和时间.稀疏脉冲反演是非线性优化问题, 通常都是把非线性优化问题转化为线性优化问题, 然后用线性优化算法求解.以范数约束为基础, 提出L1-L2范数联合约束求解的方法.该算法采用了目前国际流行的内点算法, 与传统的优化算法相比, 这种算法具有精度高, 速度快的优点.通过研究人工模型和南海某油田实际数据, 表明该算法的计算结果和测井记录吻合好, 提高了地震分辨率, 目的层段分辨率优于8m.利用反射系数剖面预测的储层厚度和开发井吻合很好, 大大地降低了资源量计算的风险和油田开发的不确定性.   相似文献   

18.
基于广义逆和奇异值分解理论,研究核磁共振T2谱反演的截断法和阻尼法。首先给出反演问题的M-P广义逆解;然后对系数矩阵进行奇异值分解求其M-P广义逆,对奇异值进行截断或者加阻尼,保证解的稳定性兼顾其分辨率;最后利用迭代技术实现解的非负约束。数值仿真实验和岩心T2谱反演表明:对于双峰谱,在信噪比SNR≥10时,截断法和阻尼法反演T2谱的分辨率和稳定性都很高;截断法截断因子选择范围小,在SNR<10时解的分辨率和稳定性差;阻尼法阻尼因子选择范围大,在SNR<10时解的稳定性好,但短组分分辨率低,长组分分辨率高。对于三峰谱,当SNR≥20时,截断法和阻尼法反演的T2谱光滑连续,具有三峰结构;当SNR<20时两者均不能反演出与构造谱接近的三峰谱。因此,2种方法可应用于SNR≥10的双峰谱和SNR≥20的三峰谱反演;从分辨率和稳定性、平滑噪声及截断和阻尼因子的选取上,阻尼法略优于截断法。  相似文献   

19.
预应力锚索抗滑桩加固边坡(滑坡)的稳定性分析   总被引:1,自引:0,他引:1       下载免费PDF全文
预应力锚索抗滑桩作为一种实用有效的支挡工程措施己在边坡整治工程中得到了广泛地应用,然而其理论研究远远滞后于工程应用。本文以极限平衡理论为基础,假定滑动面为圆弧面或近似于圆弧面,采用Bishop(毕肖普法)法,并考虑竖向条间剪力对滑移面抗剪强度的影响,建立了滑坡最危险滑移面的搜索模型,并推导了稳定系数的求解公式,运用遗传算法对支护结构最危险滑移面进行动态搜索和确定,计算其稳定性系数,避免了在优化设计中陷入局部最小值的缺点,更具有有效性和准确性。  相似文献   

20.
在油气田开发过程中,微震监测是获得水力压裂引起裂缝分布的一种较为有效的方法。微震的定位成像与裂缝解释需要利用有效微震信号位置,而微震信号具有低信噪比的特点,传统信号拾取方法无法有效实现较低信噪比条件下初至时刻的准确拾取。本文提出一种基于时频谱熵的初至拾取新方法,该方法首先通过S变换获取含噪信号的时频谱;然后对谱内各个采样点沿频率方向进行分帧操作,并计算每帧频段内的近似负熵值,以最小近似负熵值作为该谱点的负熵值;最后沿时间方向比较各谱点的负熵值,最小值对应的时刻即为初至时刻。本文利用不同信噪比的合成地震数据对该方法进行效果验证,并与长短时窗能量比(STA/LTA)法进行拾取结果对比,结果表明:信噪比在-5 dB时,两种方法拾取效果都很好;信噪比在-10 dB时,时频谱熵法拾取效果更好。时频谱熵法更适合低信噪比情况下的信号初至拾取。  相似文献   

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

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