正演反演现状
毕设粘声波正演模拟
02
有限元法：另一种常用的数值方法，通过将声波传播问题转化
为求解偏微分方程的问题，利用有限元方法进行求解。
边界元法：一种用于求解边界问题的数值方法，通过将问题转
03
化为求解积分方程的问题，利用边界元方法进行求解。
03
毕设粘声波正演模拟的实现过 程
毕设粘声波正演模拟的流程设计
确定模拟目标
明确模拟的目的和需求，例如模拟某地区的 地震波场或声波场。
提高了对粘声波正演模拟的认识
通过实践操作，我们深入了解了粘声波正演模拟的原理、方法和应用，提高了对该领域 的认识和理解。
积累了宝贵的实践经验
在毕设过程中，我们遇到了许多问题和挑战，通过不断尝试和探索，积累了丰富的实践 经验，为今后的学习和工作打下了坚实的基础。
展望
进一步优化模拟算
法
虽然我们已经实现了粘声波正演 模拟，但还可以进一步优化算法 ，提高模拟的精度和效率。
1 2
模型改进
对模型进行改进，使其更加符合实际物理现象， 提高模拟精度。
参数优化
对模型参数进行优化，减少误差，提高模拟结果 的可信度。
3
算法改进
对算法进行改进，提高计算效率和精度，减少误 差。
毕设粘声波正演模拟的未来发展趋势
01
智能化
利用人工智能技术，实现自动化 、智能化的模拟过程，提高模拟 效率和精度。
发展趋势：随着计算 机技术和数值模拟方 法的不断发展，粘声 波正演模拟将更加精 细化、高效化和智能 化。未来研究方向包 括
改进粘声波正演模拟 算法，提高模拟精度 和效率；
结合实际勘探数据， 开展多参数、多尺度 模拟研究；
利用人工智能和机器 学习等技术，实现自 动化和智能化模拟。
地球科学中的正反问题
一、什么是正演问题和反演问题在地球科学中，有两大问题是离不开的，正演问题和反演问题。
由物理定律根据给定物理模型的参数计算出数据的问题是正演问题。
而由观测数据通过适当的方法计算物理模型参数来重建物理模型的问题是反演问题。
由卫星云图预报天气、由遥感影像估计粮食产量都是正演问题。
从思路上而言，正演问题比较简单。
如果给定物理模型的系数，由物理定律能够计算出与观测数据相比对的理论数据。
在模型比较精确的情况下，正演一般能够获得比较好的效果。
当然，反演问题也在多个领域有应用，这里可以给出很多实例，比如太阳的内部结构探测、储油层厚度的估计、莫霍面深度的推断、核幔边界形态的分析等等。
由于我们不单对模型系数不清楚，甚至有时对物理模型本身都不甚清楚，所以我们可以断言反演比正演问题将面临更多更大的困难。
根据百度百科，正演问题（direct problem）定义：在地球物理磁法勘探的理论研究中，根据磁性体的形状、产状和磁性数据，通过理论计算、模拟计算或模型实验等方法，得到磁异常的理论数值或理论曲线，统称为正演问题。
反演问题（inversed problem）在磁法勘探理论研究和解释磁测成果时，根据磁异常特征，确定磁性体的形状、产状及其磁性等，称为“反演问题”。
这个概念给的范围太狭隘，就简单的地磁勘探而已，所以仅作为参考。
二、哪个先提出来现在有一个逻辑问题，是先有正演问题还是先有反演问题？似乎直观上先有前者，然而我认为，对大多数问题，尤其是系统复杂的问题，应当是先有后者。
科学研究的先驱们没有今天的人有这么好的条件，不可能通过课堂学习系统地掌握成体系的知识，也没有条件去图书馆查阅资料，更不用说利用检索工具搜集信息了，他们掌握的资料和信息是极其的贫乏的。
当先驱们涉足新的研究领域时，是没有经验可循的，也没有什么物理模型可以利用。
他们看到的是规律或者说模型所呈现出来的现象，他们的任务是找出规律、建立模型，这个任务本身就是反演问题。
全波形反演研究现状及发展趋势
文章编号 &333%&((&#$3&($3&%33<<%3<收稿日期 $3&#%&$%3+%改回日期 $3&(%3&%3#&作者简介 杨勤勇#&+,('$"男"博士"教授级高级工程师"长期从事物探技术研究和物探科技管理工作&基金项目 国家科技重大专项#$3&&[\353&(%33&%33$$资助&"全波形反演研究现状及发展趋势杨勤勇 胡光辉 王立歆中国石油化工股份有限公司石油物探技术研究院 江苏南京$&&&3#摘要 全波形反演技术是当前勘探地球物理领域的研究热点之一&新一轮的全波形反演研究触及了波场模拟+梯度估计+数据预处理+目标泛函的选择等诸多深层次的内容"逐渐将全波形反演对实际观测系统+地震数据等要求的局限性以及其具备的潜力揭示出来&通过对全波形反演理论和技术研究进展及应用现状的全面调研"介绍了全波形反演算法研究从时间域到频率域+再到混合域和拉普拉斯域的发展进程%阐明了全波形反演技术已实现海上地震资料应用但对陆上地震资料还没有真正成功实例的应用现状%着重分析了陆上资料全波形反演应用的瓶颈主要在于观测系统的限制+低频数据的缺失+数据预处理面临的挑战以及近地表条件和激发接收的影响等%指出了发展分步骤+分尺度的反演方法和反演策略以及多种手段的有效联合是实现陆上资料全波形反演的有效途径&关键词 全波形反演%研究进展%应用现状%实用化瓶颈%前景展望A %#(&34#+,+!V47I I C 4&333%&((&4$3&(43&43&&中图分类号 1,#&4(文献标识码 .3.4.56704/5/>45;22.?.918<.;//6.;21C C >99F 5?.C 16<:;?.64:1;_B C O T 7C Q ?C O "R G^G B C O 0G 7"Y B C O U7L 7C #"'$#N *7E *#N ,/@'7(3B *@*(C 7,W $@!'!0!*"K ($<'$%$&&&3#"+,'$($)B 4/657/(Z G @@M B N :K ?;F7C N :;I 7?C #Z Y D $7I ?C :?K E 0:0?E 7I I G :7C E 0:O :?H 0Q I 78I :L H@?;B E 7?C 4)0:C :M;?G C J I E G J Q 7C N ?@N ::I 7CF B C Q J ::H :;7I I G :I "I G 80B IM B N :K 7:@J I 7F G @B E 7?C "O ;B J 7:C E :I E 7F B E 7?C "J B E B H ;:%H ;?8:I I 7C O "?-%V :8E 7N :K G C 8E 7?C B C J I ??C 4)0?I :K G ;E 0:;I E G J 7:I ;:N :B @E 0:@7F E B E 7?C I ?K O :?F :E ;QB C J I :7I F 78J B E B "B I M :@@B I E 0:H ?E 7:C E 7B @?K Z Y D "E 0;?G O 0B C ?N :;B @@7C N :I E 7O B E 7?C ?C E 0:Z Y D E 0:?;Q "E :80C ?@?O Q B C J B H H @78B E 7?C I E B E G I "M :7C E ;?%J G 8:Z Y D J :N :@?HF :C E "K ;?FE 7F :J ?F B 7C E ?K ;:XG :C 8Q J ?F B 7C B C JU BH @B 8:J ?F B 7C ":I H :87B @@Q E ?E 0:0Q -;7JJ ?%F B 7C 4D E 7I I 0?M C E 0B E Z Y D 0B I -::C B @;:B J Q I G 88:I IK G @@Q G I :J 7C E 0:F B ;7C :J B E B "0?M :N :;7E 7I I E 7@@B 80B @@:C O :7C @B C J J B E B B H H @78B E 7?C %-?E E @:C :89I ?K E 0:@B C J J B E B B H H @78B E 7?C B ;:B C B @Q /:J "I G 80B I @7F 7E B E 7?C ?K O :?F :E ;Q"@B 89?K @?M %K ;:X G :C 8Q 8?F H ?C :C E ?K J B E B "E 0:80B @@:C O :?K J B E B H ;:%H ;?8:I I 7C O "E 0:C :B ;I G ;K B 8:8?C J 7E 7?C I B C J I ?G ;8:%;:%8:7N :;@7F 7E B E 7?C ?K @B C J B 8X G 7I 7E 7?C %Z 7C B @@Q 7E 7I 8?C 8@G J :J E 0B E I E :H %-Q %I E :H 07:;B ;8078B @F G @E 7%I 8B @:I E ;B E :O 7:I B C J 8?F -7C B E 7?C ?K J 7K K :;:C EF :E 0?J IM ?G @J -:B N B 7@B -@:E ?;:B @7/:E 0:@B C J J B E B Z Y D B H H@78B E 7?C 4D .EF 1624(Z Y D ";:I :B ;80H ;?O ;:I I "B H H @78B E 7?C I E B EG I "Z Y DH ;B 8E 78B @-?E E @:C :89I "H ;?I H :8E ""纵观石油工业的发展历程"其实是一部勘探地球物理技术的发展史&每一次勘探地球物理技术的进步都带来了油气产量的飞跃&随着勘探开发程度的不断深入"对高精度成像+储层描述和精细地震地质解释提出了更高的要求"亟需地球物理技术"尤其是精确的反演成像技术的进步&全波形反演#K G @@M B N :K ?;F7C N :;I 7?C "Z Y D$理论和技术以其高精度+多参数建模的能力"吸引了越来越多的勘探地球物理学家的目光"已成为当前勘探地球物理领域的研究热点&它利用叠前地震波场中的运动学和动力学信息"具有揭示复杂地质背景下构造细节及岩性的潜在能力&但是"在数学上它是一个高度病态的非线性问题)&*"需要解决其多解性及收敛性问题&从地球物理角度看"它涉及模型的参数<<第5#卷第&期$3&(年&月石"油"物"探^*>1R _W D !.U16>W 1*!)D "^Z >61*)6>U *]P=?@45#""?4&A B C 4"$3&(化+误差泛函的建立+数据预处理+波场的数值模拟+子波的估计等研究内容)$*&在某种程度上可以说"全波形反演理论和技术是地球物理领域的一项终极技术&因此"该项技术的发展必然是长期的+逐步完善的过程&全波形反演是c B Q:I估计理论在勘探地球物理领域的一个应用范例&它可以描述为一个基于地震全波场模拟的数据拟合过程"其使用了地震记录中的全波形信息)&*"而不像其它传统的方法仅使用地震波形中的部分信息#如旅行时层析成像等技术$&全波形反演的实现是在正则化约束下通过更新迭代初始模型进而减小计算数据和观测数据之间的误差"逐步逼近真实模型的过程)#*&理论上"全波形反演已被证明为一种建立高精度速度模型的有效手段)(%5*&全波形反演理论是一套完美的体系"但它对数据+初始模型+地震波正演+激发子波是有理论假设的&某些数值模型验证全波形反演过程可以得到几乎和真实模型完全匹配的反演结果"但是在复杂介质情况下"即使理论模型数据的Z Y D反演结果也不收敛"更不用说实际地震数据了&Z Y D反演结果很难收敛到正确的结果上"核心问题在于地震波波场#尤其是反射地震波场$与反演参数之间的关系是严重非线性的"误差泛函存在非常多的局部极值点"因此需要很好的初始模型来降低误差泛函的非线性性&然而"在复杂介质和低信噪比情况下"正确的初始模型的获取本身就十分困难%其次"简单的地震波正演模拟算子不能模拟实测地震波场中复杂的波现象%第三"地震子波的空变特征更加重了地震波场与反演参数之间的非线性关系"这是陆上地震数据进行Z Y D测试更为困难的基本原因&更深入地看"地震波在地下介质中传播"不同的波现象与不同的介质成分紧密关联&譬如潜波#S7N7C O Y B N:$可以认为是一种地表观测的透射波"利用S7N7C O Y B N:可以反演背景速度#初始速度模型$&但是"S7N7C O Y B N:要在深层介质中传播"一定是低频的&另一方面"还需要长偏移距才能观测到中深层的S7N7C O Y B N:&这是Z Y D需要低频长偏移距数据的基本原因&缺乏了低频长偏移距数据"背景速度就需要用其它的方式获取&与此同时"Z Y D的实用化必然与计算效率相关&在一般的迭代算法中"Z Y D的一个迭代步骤至少需要全部炮集的#次地震波正演计算&据此可以看出"Z Y D的计算量多么巨大9总之"尽管Z Y D理论和技术看起来很完美"但是它的理论限制和应用瓶颈很多&在海上地震资料处理中"全波形速度反演的精度比目前的层析速度反演要高&陆上资料的全波形反演应用主要受限于子波空变+信噪比低"经典意义下的全波形反演还没有很成功的应用实例&但是"无论如何"当前的Z Y D反演技术主要功能还是估计比较精确的背景速度"与逆时偏移#;:N:;I:%E7F:F7O;B E7?C"6)P$一起实现精确的反射地震波成像"距离精确的储层参数反演还有较长的路要走&深入分析Z Y D理论及其受限制原因"提出满足其基本假设的数据预处理方法和实现Z Y D的更合理的流程"是将Z Y D导向实用化所必须的举措&我们在全面调研的基础上"介绍了全波形反演技术的研究进展+应用现状+存在的瓶颈及发展趋势"着重分析了陆上地震资料应用的困难及可能的解决方案"最后展望了全波形反演技术的实用化进程&&"全波形反演技术研究进展$3世纪'3年代")B;B C E?@B)&",*提出了基于广义最小二乘的时间域全波形反演"这一方法的产生推动了全波形反演的发展"对其后全波形反演理论体系的发展产生了深远影响&)B;B C E?@B借鉴了共轭状态法求梯度的思想"通过炮点正传波场与检波点残差逆传波场的互相关估计出梯度方向"避开了Z;:80:E导数的直接计算"使得二维时间域全波形反演的实现成为可能&据此也可以看出"高精度+高效率地模拟地震波传播的波场是全波形反演计算的核心问题)<*&=7;7:G L)'%+*首先将二阶有限差分交错网格方法运用于地震波场模拟%之后"U:N B C J:;)&3*将该方法推广至四阶"由于其在精度及效率上与全波形反演的要求吻合)&&*"四阶有限差分交错网格法一直被广泛采用&此外"`?I@?K K等)&$%&#*提出的伪谱法和P B;K G;E)&(*提出的有限元法也得到了充分的应用和发展&在反演策略上"由于时间域全波形反演可以灵活地对数据进行必要的预处理"以及灵活选取所需的特征波等特点而受到关注&但时间域反演同时反演所有频率成分"增大了问题的非线性性"因此"c G C9I等)&5*提出了时间域的多尺度反演"这种算法通过对资料的处理将问题分解为不同的空间尺度"增加了问题求解的稳定性"降低其非线性程度&$3世纪+3年代"1;B E E)&,%&<*将)B;B C E?@B的理'<石"油"物"探第5#卷论发展到了频率域"由此奠定了频率域全波形反演发展的基础"并极大地推动了全波形反演的实用化进程&在此之前"P B;K G;E)&(*就明确指出"对于处理多震源问题"频率域的有限元或有限差分法是最有效的数值离散化手段&但后来被>H:;E?等)&'*证明频率域波场模拟在处理三维较大模型时有很大的局限性"尤其是对内存的需求&频率域全波形反演仅需几个离散的频率即可完成模型的高精度重建"因此"对每道地震记录而言"几个相应的傅里叶级数取代时间域整个时间序列"大大节省了存储空间%其次"频率域全波形反演直接在频率域求解"反演过程中直接处理频率域的解"因此"容易实现从低频到高频的多尺度反演%此外"频率域求解过程中容易加入吸收因子等参数)&&*&所以"近年来频率域全波形反演备受关注&频率域正演的发展为频率域反演的实现奠定了基础&其中具代表性的有>H:;E?等)&'*提出的基于三维$<点加权平均算子的粘声波正演法%c;?I I7:;等)&+*提出的省资源有限体积法%*E7:C C:等)$3*提出的不连续的^B@:;% 97C方法&由于三维频率域正演对内存的超大要求"一般的机群很难满足"使三维频率域全波形反演的发展受到了很大的限制"也因此促进了频率域反演联合时间域正演的混合域反演方法的发展)$&*&混合域反演联合了频率域反演和时间域正演的优点"利用傅里叶变换在波场模拟过程中直接求解频率域的解"不增加额外的计算量"同时可以获得多个频率域的波场信息&这种算法既节约了存储空间+方便实现多尺度算法"又可以灵活地对数据进行预处理"且不需要超大内存空间"更适应现有机群设施"在计算资源满足的条件下适用于大规模并行计算"最大程度地提高计算效率)&&*&所以"混合域反演方法的诞生极大程度地促进了三维全波形反演在实际资料中的应用进程&低频数据与初始模型的耦合是全波形反演在实际资料应用中遇到的最大瓶颈&低频信息的缺失"使得常规建模手段难以满足全波形反演对初始模型精度的要求&因此"无论是时间域反演还是频率域反演"往往造成模型更新迭代过程中的错误收敛&为了解决实际应用中的这一问题"W07C 等)$$%$(*利用拉普拉斯域对频率不敏感的特性发展了拉普拉斯域的全波形反演&遗憾的是拉普拉斯域的全波形反演只能恢复模型中的长波长信息"并不能获得像频率域全波形反演那么高的精度&但有了准确的长波长速度分量"作为频率域全波形反演的初始模型"是帮助频率域全波形反演绕过低频信息的有效手段&因此"`7F等)$5*又发展了拉普拉斯域联合频率域的全波形反演&近年来"全波形反演得到了越来越多的关注"一些优化算法的提出更促进了全波形反演的实用化进程&全波形反演从最初的时间域已经逐步发展到了频率域+拉普拉斯域"值得一提的是混合域算法的提出对全波形反演的实用化进程有着巨大的推进作用&在当前计算能力快速发展的支撑下"全波形反演的应用已经从二维走向了三维"从模型验证阶段逐渐走向了实用化阶段&$"全波形反演的应用现状全波形反演的应用最早出现于$3世纪'3年代"^B G E07:;等)$,*和P?;B)$<*实现了二维地震资料的全波形反演&这些应用实例证明了全波形反演是一种高精度的建模手段"它具有精细刻画地下构造及岩性的能力"但同时也指出了全波形反演对初始模型的严重依赖"这一0病态1性使得全波形反演在缺少低频信息的情况下很难取得成功&受计算资源的限制"当时Z Y D技术的应用始终局限在二维情况&之后"地球物理学家对全波形反演的应用进行了更深入的研究"一些新方法被应用到全波形反演中"并出现了很多成功应用的案例)$'%#&*&三维全波形反演在$3世纪+3年代中后期陆续出现)#$%##*"当然计算水平的提高在三维全波形反演应用中扮演了重要角色&$3&3年"W7;O G:等)#(*率先对挪威北海油田#=B@0B@@地区$>c!数据实现了三维全波形反演"这一成果极大地鼓舞了全波形反演的研究热潮&=B@0B@@地区由于浅层气云的覆盖"一直是建模+成像的难点"甚至被称为勘探的0盲区1&全波形反演不仅对该气云形态进行了准确的描述"对其周边充气的断裂构造也进行了精细的刻画&由此也发展了直接对高精度速度体进行地质解释的地震资料解释新思路&浅层建模一直是建模的难点"从=B@0B@@地区反演结果来看"全波形反演完全有能力实现浅层的高精度建模"其对该地区浅表层古河道的精细刻画"达到了常规建模手段无法达到的精度&之后"全波形反演对海上三维实际资料的应用陆续出现)#5%#<*&遗憾的是"这一技术仍然停留在只能应用于海上地震资料阶段"陆上资料的应用还存在很大挑战"主要是无法提供满足全波形反演要求精度的初始速度模型及足够低频的观测数据&+<第&期杨勤勇等4全波形反演研究现状及发展趋势$3&$年"壳牌公司与东方地球物理公司合作实现了二维陆上资料的全波形反演)#'*&虽然这一研究是基于低频大偏移距的特殊观测系统"在实际地震资料采集中很难大规模实施"但这一结果也证明了陆上资料全波形反演应用的可行性&从&h5R/的初始频率开始反演"减少了全波形反演对初始模型的依赖&大偏移距的勘探"保证了记录中的全波场信息&因此"从一个均匀递增的背景速度场出发"也得出了很好的反演效果&低频数据对恢复长波长速度分量起到了重要作用&也再一次证明了低频数据的重要性&在$3&#年美国勘探地球物理学家学会#W*^$年会上"出现了很多不同区块的全波形反演实际资料应用实例&此外"拉普拉斯域的全波形反演也得到了发展&1Q G C等)#+*发展了拉普拉斯域的三维弹性波全波形反演&与直接解法不同"他们采取和频率域近似的迭代解法"使该方法更适用于大尺度三维全波形反演&针对海洋深水环境"U::等)(3*提出了一种剥离直达波的拉普拉斯域全波形反演"与常规的对数正态的拉普拉斯域全波形反演相比"剥离直达波的全波形反演方法更适合海洋深水环境&此外"拉普拉斯混合频率域的全波形反演也在实际资料应用中得到了验证)$(*&在缺少低频信息的情况下"利用拉普拉斯域全波形反演对频率不敏感的特性首先恢复长波长信息"继而以此为初始模型利用频率域全波形反演恢复模型的短波长分量"实现模型的高精度重建&国内对全波形反演的研究起步较晚"技术水平较国外还有一定的差距"目前国内还没有看到全波形反演成功应用的实例&针对国内陆上探区特点"中国石油化工股份有限公司石油物探技术研究院开展了针对目标层的全波形反演研究"并取得了阶段性成果)(&*"该结果展示了全波形反演对中浅层建模的积极作用"对成像剖面有较大改善"构造与测井信息更加吻合&目前"无论是海上资料还是陆上资料的全波形反演"都使用层析反演的手段来获取初始速度模型"但结果表明"这一速度模型精度并不能满足全波形反演的要求&因此"W07C提出的拉普拉斯域全波形反演为经典意义下的全波形反演提供初始速度模型备受关注&此外"特征波的波形反演也为初始模型建模带来了新的曙光&虽然我们已看到了陆上资料全波形反演成功应用的实例"但这些大多需要特殊的观测系统和低频勘探作为保障"又或者与经典意义下的全波形反演不同"仅实现了特征波的全波形反演&所以"能否实现陆上地震资料的全波形反演"更进一步地说"能不能实现常规观测系统下陆上资料的全波形反演"将是极具挑战性的研究课题&#"陆上地震资料应用的瓶颈及发展动向""理论上"全波形反演处理应用于海上地震资料或陆上地震资料时并无差异&无论是波场模拟+梯度求取还是优化算法都是同一套理论体系&那么"为什么海上地震资料可以正确收敛到一个高精度的速度模型"而陆上地震资料的应用却一直受到限制呢6王华忠等)($*从概率论的观点分析了波形反演的本质"并指出在假设观测噪声为高斯白噪的情况下"c B Q:I估计可以在最小二乘意义下实现%分析了全波形反演难以实现的根本原因在于数据空间向参数空间映射的强非线性关系以及介质模型的复杂性和描述地震波场物理传播过程的正演算子的复杂性&胡光辉等)#*通过模型验证分析了全波形反演的应用现状并指出了其在陆上资料应用中的困难&对比海上地震资料与陆上地震资料的特点"从数据域的角度出发"我们认为陆上地震资料全波形反演的应用瓶颈主要在于(&$全波形反演要求全方位角+大偏移距的观测系统&海洋资料>c!数据很好地适应了这一要求"=B@0B@@油田的成功应用也证实了这一点&而陆上地震资料的常规三维勘探多为滚动采集"这对全方位角的要求有一定的限制&但最为严重的是"传统的观测系统设计多基于反射波勘探"偏移距较小"波路径正交程度差"这在很大程度上限制了全波场信息的获取&$$低频信息的缺失&海上数据去除海洋噪声干扰"一般最低可用频带为#43%#45R/&而传统的陆上检波器所能响应的频带有效范围往往在,R/以上"即使数字检波器可以响应低频信息"而面波的干扰也使得该频段不可用&此外"常规的建模手段又很难较准确地恢复模型的长波长信息&因此"没有低频数据与初始模型很好地耦合"使这一0病态1的反演问题极易陷入局部极小"不能获得全局的最优解&#$数据预处理面临挑战&全波形反演除了运用地震波的走时信息"最重要的是它考虑了波传播的动力学特征"因此在数据拟合过程中"其对噪声非常敏感&这就对数据预处理工作提出了更高的3'石"油"物"探第5#卷要求"既要消除噪声干扰"又不能破坏波的动力学特征&尤其在声波近似的全波形反演过程中"面波作为干扰波需要清除"而在去除面波的过程中"往往损害了地震记录的低频信息&甚至由于其线性特征"同时损害了初至及浅层折射等有明显线性特征的有效波信息&($与海上气枪震源相比"陆上炸药震源稳定性较差"受药量+激发岩性等因素影响"随机干扰严重&地面检波器受地表条件影响也制约了陆上资料的品质"且检波器与地面的耦合性远远不能做到水上检波器与海水这种一致性的耦合&而且与海水这一均质的速度体相比"近地表介质极其复杂"许多不确定因素导致波动方程并不能完全模拟全部波形"因此在构建误差函数时累积了更多的错误与误差"进而影响其收敛方向&全波形反演是一套完美的数学理论体系"已逐渐从模型试算走向实际资料的应用"并已在海上资料的实际应用中获得了成功&陆上资料的成功应用虽然还存在种种困难"但地球物理学家们已经注意到了其中的种种挑战"意识到陆上资料全波形反演应用的限制已经不属其理论范畴&所以"当前全波形反演的发展已经从理论研究转向了推动其实用化的研究"以及为全波形反演实用化的各种准备工作的研究&针对陆上资料特点"\G等)(#*发展了基于反射波的波形反演研究"并在实际资料应用中取得了较好的效果"也使特征波的全波形反演受到了人们的关注&c;?I I7:;等)((*改进了反射波全波形反演的目标函数"使用伴随状态法的全波形反演策略实现了反射波这一特征波的全波形反演&特征波全波形反演"使用某一或某一些特征波的走时和波形信息逐步逼近全波形反演的结果及精度"是解决陆上资料全波形反演问题的有效途径&这一方法也在实际资料应用中得到了验证)(5%(,*&此外"针对陆上资料低频信息缺失问题的研究也取得了新的进展)(<%('*&通过低频补偿+频率外推以及改变目标函数)(+%53*对频率的敏感性等方式"减小低频信息缺失带来的不利影响"完成缺少低频信息下的全波形反演&("全波形反演应用前景展望全波形反演技术是地球物理领域一个新的研究热点"近年来得到了快速的发展&在某些海上实际资料的应用中得到了优于现有其它速度建模方法的结果"陆上地震资料全波形反演的应用还受到种种限制&随着地震采集技术水平的不断提高"全波场+宽频带+炮检对等观测手段的实现"以及预处理保真技术的发展"陆上全波形反演的应用将逐步走向实用化&我们认为"经典的全波形反演方法的实用化是非常困难的"需要发展一种反演流程"通过分步骤+分尺度+多种手段联合来实现全波形反演"大幅度提高目前的速度估计和成像的精度&王华忠等)($*分析指出"速度场的反演利用特征波的波形反演#!Y D$要优于全波形反演%需要合理地增加相位信息在泛函中的比例以及考虑模型的先验信息& Z Y D反演的根本问题在于地震波场与反演参数之间的非线性性"在于实测地震波场的概率特征并不完全符合高斯分布"在于地震波正演不能很好地模拟实测波场"在于反射波振幅不只受速度改变的影响&针对这些根本问题"首先要弄清楚地震波场中的不同成分与速度参数的不同成分之间的关系"尽量让这些关系紧密化和线性化&因此"将地震波场分解成特征波场"不同特征波场成分反演速度场的不同成分"综合应用波场中的旅行时+相位+振幅信息"构建实用化的全波形反演技术的反演流程"是将全波形反演技术应用于陆上地震数据当前需要解决的问题&\G等)(#*发展的反射波全波形反演就是这种特征波反演思想的体现&此外"通过设置目标函数对减小非线性问题"提取波场特征量信息以及对不同摄动尺度的响应扮演着重要角色&董良国等)53*分析了不同的目标函数随不同物性参数的不同摄动尺度的变化性态"并给出了模型验证结果"对选择合理的反演方法以及反演策略具有重要意义&不同目标函数的非线性程度不同"对目标函数形态的有效分析是研究陆上资料全波形反演应用的恰当举措&拉普拉斯域全波形反演以及地震道包络构造的目标函数正是利用了它们在数据空间和模型参数之间的相对较低的非线性程度&所以"分步骤+分尺度反演方法和反演策略以及多种手段的有效联合是实现陆上地震资料全波形反演的有效途径&我们相信"在不久的将来"全波形反演基本思想和方法一定会对弹性参数反演+多参数反演!建模+高精度成像以及储层参数估计与储层描述产生积极而又深远的影响&参"考"文"献)&*")B;B C E?@B.4D C N:;I7?C?K I:7I F78;:K@:8E7?CJ B E B7CE0:B8?G I E78B H H;?L7F B E7?C)A*4^:?H0Q I78I"&+'("(+#'$(&$5+%&$,,&'第&期杨勤勇等4全波形反演研究现状及发展趋势。
地球物理学中的电磁场正演与反演
第22卷 第4期地 球 物 理 学 进 展Vol.22 No.42007年8月(页码:1181～1194)PRO GRESS IN GEOP H YSICSAug. 2007地球物理学中的电磁场正演与反演汤井田, 任政勇, 化希瑞(中南大学信息物理工程学院,长沙410083)摘 要 本文在近年来众多的地球物理研究者的研究基础上,总结了当前地电磁模型正反演已有成果,定量分析了各种主要正反演的性能测试,指出不同正反演法的优点、缺点以及应用范围局限,提出了各种方法的发展趋势以及当前计算地球物理领域的核心内容,指出了计算地球物理领域的数值模拟发展方向.关键词 有限差分,有限元,积分方程,线性迭代,蒙特卡罗,电磁模型,正演,反演中图分类号 P318,P319 文献标识码 A 文章编号 100422903(2007)0421181214The forw ard modeling and inversion ingeophysical electrom agnetic f ieldTAN G Jing 2tian , REN Zheng 2yong , HUA Xi 2rui(S chool of I nf o -p hysics and Geomatics Engineering ,Changsha ,410083)Abstract Based on the excellent achievements by many geophysical researchers at present ,this paper has analyzed performances of the most used forward and inversion methods in electromagnetic quantitatively ,and then ,has pointed the merits and faults of this algorithms.So ,with the quantitative analysis and testing ,the development trend of a 2bove forward and inversion in geophysical electromagnetic filed is pointed and also the core of computational geophys 2ics is listed.At the end ,we have clearly listed the development trend of the numerical simulation in computational ge 2ophysics.K eyw ords finite difference method ,finite element method ,integral equation method ,linear iterate method ,monte 2carlo method ,electromagnetic simulation ,forward model ,inversion收稿日期 2007204210; 修回日期 2007206220.基金项目 国家863计划(2006AA06Z105,2007AA06Z134)项目资助.作者简介 汤井田,博士,博士生导师,中南大学教授,中国地球物理学会会员,美国勘探地球物理学家协会(SEG )会员.主要从事电磁场理论和应用、地球物理信号处理及反演成像等研究.(Email :jttang @ ).0 引 言在地球物理学,电磁场的复杂性决定了地球物理模型的复杂性,一般而言,地球物理模型无法以解析法得到解析解[1],因此,数值模拟方法在地球物理学中得到了广泛地应用,并以此,地球物理学家得到了许多经典的地电模型的电磁场分布数据.借助于这些电磁场分布数据,结合地电模型结构,我们可以初步建立电磁场数据与模型之间的对应关系.但对于复杂的模型来说,其电磁场分布也非常之复杂,在这种情况下,模型与电磁场数据之间的关系变得十分复杂,因此,需要一种高效的、准确的方法来建立模型与电磁模型之间的关系,这种方法即称为电磁模型的反演[2].目前而言,反演主要集中在完全2D 、3D 非线性模型上[3～5],在其中,3D 电磁场数值模拟是3D 电磁反演的核心引擎,因此反演与正演是相得益彰,互相促进的.限于篇幅,本文只讨论广泛应用于地球物理电磁场正演的有限差分[6～8],有限元[9,10],积分方程法[11,12]等,基于此的方法变种,如微分2积分法[13]等不具体讨论;对于反演来说,只讨论线性迭代法[14]、蒙特卡洛反演方法[3].而其它一些变种如微分2积分方法[13]等不具体论述.本文第一部分,讨论有限差分、有限元和积分方程法,分析其现有应用效果,其优点与缺点,基于此分析其发展趋势;第二部分详细论述反演算法的应用以及发展趋势,集中讨论线性化迭代法,蒙特卡洛地 球 物 理 学 进 展22卷非线性全局最优化方法等,分析其优点与缺点,并讨论解决当前阻碍其发展的解决方法,指出非线性反演的在电磁模型中发展趋势.表1 符号的意义T able1 The meaning of symbolsV Laplace算子ε介电常数μ磁导率σ介质电导率ω～角频率j ext外加电流σ～=σ-iωt复电导率e-iωt时间依赖常数E电场强度H磁感应强度E0初始电场强度H0初始磁感应强度r空间坐标V s体积A系统矩阵D空间维数X节点值向量B右边向量φ目标函数m目标模型δm模型增量λ罚函数因子J n×m灵敏矩阵H n×m海森矩阵i,n,k不清索引计数β调整因子M=N M模型集x模型参数w,v模型个体1 电磁场正演分析电磁正演模型的宏观控制方程为Maxwell方程,就其在频率领域的形式为[2]:Δ×H=σ～E+j ext,Δ×E=iωμH.(1)求解(1)式,便可获得H和E.对于绝大多数模型,(1)式只能够通常数值方法来求解,下面列举主要数值方法最新进展.1.1 有限差分法[4,7]有限差分(Finite-deference met hod,FDM)是最为古老的数值计算方法之一,其被用于应用地球物理邻域始于20世纪60代(Yee1966[7];Jones and Pascoe,1972[15];Dey and Morrison,1979[16]; Madden and Mackie,1989[17]),特别进入90年代,交错式样网格被广泛用于地电磁场的分析中来,使有限差分法步入全盛时期(Smit h and Booker, 1991[8];Mackie et al,1993,1994[18,19];Wang and Hohmann,1993[20];Weaver,1994[21];Newman and Alumbaugh,1995,1997[22,23];Smit h,1996a, b[24,25];Varent sov,1999[26];Champagne etal, 1999[27];Xiong et al.,2000[28];Fomenko and Mogi, 2002[29];Newman and Alumbaugh,2002[30]).有限差分的基本原理为:方程(Ⅰ)控制的模型被分为规则的网格,其规模为M=N x×N y×N z,N i为直角笛卡尔坐标系的坐标轴方向的节点距,电磁与磁场被离散到节点,并导致一些关于电磁场节点值的线性方程组,A FD X=B,A FD为3M×3M的复数、对称、大型、稀疏矩阵,X为3M长的各节点电场或磁场的三方向值的向量,B 为由j ext等激励和边界条件生成的长度为3M的向量.同上可知,有限差分的最大不足之处为,它要求模型能够被剖分成规则的单元如四边形,六面体等,严重制约了其在复杂地球物理模型中的应用;最大优点在于能够非常好的处理内部介质中电磁性差异引出的磁场与电磁不连续现象,这是由交错网格的基本性质决定.目前来说,作为电磁数值模拟方法的主导者,有限差分法(FD)正处于各向同性介质模型转向各向异性介质模型的升级(Weidelt,1999[31]; Weiss and Newman2002,2003[32,33]);正处于频率域电磁模型的模拟向时间域电磁模型模拟的空间转换,并借助于并行技术求解(Wang and Hohmann, 1993[20],Wang and Tripp,1996[34],Haber et al, 2002[35];Commer and Newman,2004[36]).1.2 有限单元法[4,8,9]有限单元法(Finite element met hod,FEM)并未广泛地被应用到地电磁场数值模拟计算当中来, FEM利用节点值与节点基函数来形成整个电磁场的分布.不同于FDM,FEM是基于电磁场的积分形28114期汤井田,等:地球物理学中的电磁场正演与反演式,它是由电磁场的微分形式通过Green等定理变换而得,通常也称有限法的解为微分形式的弱解.同于FDM,FEM最也形成大型,对称,复数,稀疏矩阵,A FE X=B.不同于FDM,FEM并不一定要求模型能够被剖分成规则单元,如三角形与六面体单元(其被理论与实践证明可以无限度精确地模拟地球物理模型),因此,FEM能够求解FDM不能够求解的复杂地球物理模型,并被应用于实际中(Reddy,1977[37] Coggen,1976[9];Pridmore,1981[38]Pasulsen, 1988[39]Boyce,1992[40]Livelybrooks,1993[41]Lager and Mur,1998[42]Sugeng,1999[43]unorbi,1999[44] Ratz,1999[45]Ellis,1999[46]Haber,1999[47]Zyserman and Santos,2000[48]Badea,2001[49]Mit subhata and Uchida,2004[50].由上可知,FEM不仅能够处理FDM能处理的简单模型,更能够处理复杂的模型,因此,FEM能够作为地电磁场数值模拟的通用者. FEM显然肯定一些不足之处:对于复杂的模型,其结果不能给人以绝对的信服,其解没有相应的误差分析,并且这种分析是非常之必要.FEM的发展趋势:(1)对复杂的模型给予相应的精确的误差分布,难以肯定结果的真实可靠性[24～30];(2)基于势理论的成长,电磁场借助于矢量势与(或)标量势的方程系统能够完美的代表电磁场分布,有限元求解这些系统是一种大势所趋[44,49,50].(3)虽然FD能够处理内部边界电磁场不连续现象,但是基于节点的有限元法不能处理此理解,从而给结果带来误差,基于边的矢量有限元能够很好的处理节点有限元的不足[43,50],因此,随着对误差的要求越来越小,矢量有限元将会越来越多的应用到地电磁场的分析中来.1.3 积分方程法[4,11,12]积分方程法实现了均匀导电半空间三维大地电磁响应的数值模拟.即求取张量格林函数积分时,采用二次剖分算法解决计算中奇异值问题,对于含有贝塞尔函数的积分项,利用结合连分式展开的高斯求积代替常规的快速汉克尔变换方法,确保了张量格林函数的正确计算并提高了计算精度.最后通过数值模拟结果的对比及模型试算验证了算法的正确性.积分方程法(Integral equation met hod,IE)把Maxwell方程变成Fredholm积分方程(Raiche, 1974[11,12])E(r)=E0(r)+∫V s G(r,r′)(σ～-σ～0)E(r′)dr′,(2)(2)式为电场表达式,此方程即为著名的散射方程(Scattering Equation,SE).(2)式中,E0(r)通常为已经项,G为3×3的Green函数在1D参考介质中矩阵,V s为(σ～-σ～0)不为0处的体积.通过离散化方程(2),产生线性方程组,AIEX=B为复数、密实矩阵.由此可见,IEM的主要优点为线性方程的维数相对FDM、FEM要小的多,可以快速求解模型;不足之处为,解的精度严重依赖于AIE的精确度,但一般来讲,AIE的精确无法得出有限保证,并且其本身也是一项十分耗时的工作.但是由于其速度快的优点,特别是在3D电磁模型计算中,被广泛地应用(Ting and Hohmann,1981[51];Wannamaker, 1984[52];Newman and HOhmann,1988[53];Hohm2 ann,1988[54];Wannamaker,1991[55];Dmit riev and Newmeyanova,1992[56];Xiong,1992[57];Xiong and Tripp,1995[58];Kauf man and Eaton,2001[59]).由于其速度快的原因,IE的发展趋势为求解三维大型、超大型基本电磁模型上面,由此可见,IE是所有电磁场数值模型中的效率快速者.积分方程法主要优点为,1.积分方程法只须对异常体进行剖分和求积,不涉及微分方法中的吸收边界等复杂问题,在三维电磁数值模拟研究中具有快速、方便等特点,与有限元和有限差分法相比,这种方法在模拟有限大小三维体电磁响应时更为有效,计算速度快,占用内存少因而积分方程法近年来受到人们的关注和重视,并取得较快的发展. 2.由于计算机的迅速发展,对异常体进行三维网格剖分和数值求积已变得越来越方便.同样的问题,用计算机计算的时间比以前大大降低.三维电磁响应数值模拟不再是“昂贵”和“费时”,从而可以成为一种廉价、快速、能推广的解释技术.1.4 频谱Lancsoz分解法[4]频谱Lancsoz分解法(Spectral Lancsoz Decompo2 sition Method,SLDM)是一种频率中非常有效的数值模拟方法(Druskin and Knizhernam,1994[60]; Druskin,1999[61]).特别是有模型多频率情况下的首先者,因为SLDM在求解多频模型所需时间与其它方法如FDM、FEM、IDM求解单频模型所需时间相当.SLDM由于其在多频模型模拟上的优点,算得上电磁场模型模拟中的高效者.目前而看,SDLM正转向各向异性模型的模拟(Wang and Fang,2001[62]),3811地 球 物 理 学 进 展22卷Davydycheva(2003)[63]提出了特别的电导率平均法与最优化网格法来减小网格大小与数目,从而加速了SDLM的速度,使其效率更上一层.综观上述各种数值模型方法,正演各种数值方法不外乎把地球物理模拟转化为复数,大型的线性方程组.因而如何快速、准确地求解此线性方程成为重中之重,在数据表明,此线性方程的求解时间约为总求解时间的80%[2].通常来说,由FEM、FDM、ID、SDL M等法生成的线性方程的条件数(Condi2tion Number,CN)非常之大(109-1012,Tamarch2enko,1999[64]),而求解速度与CN成正比,因此十分之有必要减小线性方程式的CN,从而加速成了方程组的预条件处理器(p reconditioners)的发展.在IEM方面,通常利用M ID E(modified iterative2dissipative met hod)来加快方程的收敛速度(Sing2er,1995[65];Pankratov,1995,1997[66][67];Singerand Fainberg,1995,1997[68][69];Avdeev and Zha2nov,2002[70]),通常与FDM法(Newman andAlumbaugh,2002[30])相对比,足见M IDE在ID中的作用,表2列出了IE与FD方法中各种预处理器的性能.表2 各种预处理器的性能,模型为三维感应测井(引用Avdeev(2002)[30])IE测试平台为PC P2350MH z,FD测试平台为IBM R S-6000590工作站T able2 The performance testing of differentpreconditioners,testing on3D induction loggingmodel(cited from Avdeev(2002)[30]).testing platform is PC P2350MH z for IE andIBM R S-6000590for FD正演方法网格大小N x×N y×N z=M频率(k Hz)预处理器迭代次数运行时间A(s)IE 31×31×32=30752101600MIDM72950500056332810L IN172121FD435334160J acobi60005686 4353345000J acobi12001101对于FDM、FEM、SLDM来说,最通常用预处理器则为J acobi,SSOR与不完全L U分解器(例如,M=25×22×21=11550,N bicgstab=396;T CPU= 18min在P31-Ghz PC上,Mit suhata and Uchida, 2004).另外,还有低感应数法(Low induction num2 ber,IN,Newman and Alumbaugh,2002[18])与多重网格预处理器等,表3、4列出L IN与J acobi处理器的测试性能.表3 IE法中的L IN与Jacobi处理器的测试性能,模型为3D感应测井模型的结果统计(采用Avdeev,2002[30]),本次Jacobi测试平台为P350MH z,LIN平台为IBM RS-6000590工作站T able3 The perform ance testing of L IN and Jacobi on IE method,testing models is3D induction logging models(cited from Avdeev,2002[30])Jacobi is tested on PC P2350MH z, L IN is tested on IBM RS26000590w orkstation预处理器迭代次数相对残差J acobi1 1.00E-035 2.00E-11L IN1 1.10E-011009.40E-051000 1.30E-10由上表各表定量分析可知,经预处理过的线性方程组不仅在收敛速度上加快,而且在精确度上也有所提高.因此,寻找最优的预处理器是今后地电模型电磁正演的发展趋势之一.2 电磁模型反演反演领域十分活跃,目前反演存在三个主要问题:(1)理论表明反演的收敛速度严重依赖于正演模型的精确,但目前正演的准确度仍然无法得以保证(Zhdanov,2000[70];Torres2Verdin and Ha2 bashy,2002[71];Zhang,2003[72]).(2)反演问题通常规模较大,通常需要在成千上万的节点上反演成千上万的参数.就目前而言,计算机速度较难以提供如此之动力.(3)地球物理模型的反演通常是非线性的、病态的,这有增加了数值模拟上的困难,结果很难以收敛到精确解,只可以把误差控制在一定的范围之内.非线性成倍增加了反演的计算负担,使反演很难在完全现实的状态中完成.(4)反演存在非唯一性、非稳定性,要解决此困难,通常要包括稳定罚顶(Stabilizing Penalty Func2 tion,SPF,Tikhonov and Arenin,1977[73]);通常SPF依赖于先念信息,可影响解的平稳性、精确性等等(Part niaguine and Zhdanov,1999[74];Sasaki, 2004[75];Heber,2005[76]).因此,选取合理的SPF 在反演过程是十分重要(Farquharson and Olden2 burg,1998[77]).因此,完全反演将会是十分活跃的领域,以下为48114期汤井田,等:地球物理学中的电磁场正演与反演当前反演的主要方法和最新进展.2.1 线性化迭代法线性化迭代法(linear iterator met hod,L IM)地电磁模型反演算法中最为古老的方法(Eato n,1989[78]),在其产生的10之中,发展较为缓慢.非约束非线性最优化(Unconst rained nonlinear optimi2zation,Nocedal and Wright,1999[79])思想的引入使得L IM得到快速发展,数学理论的完善更是推动了L IM的进步.L IM的标准迭代公式可表示为:φ(m,λ)=φd (m)+λR(m)→m,λmin,(3)一般来讲,要求解(3)式的最小值值问题,可应用非线性牛顿迭代性(nonlinear Newto n2type itera2 tions,NN I;如,Newton Iterations,N I、Gauss2 Newton Iterations,CN T、quasi2Newton Iterations, QN I)求解模型空间参数.一旦(3)式得到了满足,得是反演具休来说,在每的最优模型.L IM算法描述如下:Step1:初始化模型参数。
测井约束地震反演的现状和发展
5科技资讯科技资讯S I N &T N OLOGY I N FORM TI ON 2008N O.01SCI ENC E &TEC HNO LO GY I N FO RM A TI ON 工业技术地震反演技术是伴随着地震技术在油田勘探开发中的不断深入应用而发展起来的,是20世纪80年代兴起的一门新的方法技术,该方法能利用观测数据恢复地下地质结构和岩石性质。
地震反演技术一直是地震勘探中的一项核心技术,其目的是用地震反射资料,反推地下的波阻抗或速度的分布,估算储层参数,并进行储层预测和油藏描述,为油气勘探提供可靠的基础资料。
现行的地震技术难于满足地质任务对地层纵横向分辨率的要求,测井技术虽然具有很高的纵向分辨率,但不能满足地质任务对井间地层参数变化的要求。
为此我们利用地震和测井资料的互补特征,研究测井约束地震反演技术,以取得井间地层物性参数的分布信息。
1测井约束地震反演发展过程测井约束地震反演技术是一种基于模型的反演技术,一般要求求解一个最优化问题。
测井约束地震反演充分利用测井的低频－高频成分和丰富的地震中频信息,以地震剖面所过井位的声波测井资料和地震层位解释结果作为约束条件[1~2],通过迭代反演对地质模型进行不断反复修改,使合成地震纪录资料与实际地震资料尽可能逼近,最终模型就是反演结果。
由此可见,这种反演实际上是通过正演来完成的。
反演的精度和分辨率与初始模型给定有很大的关系,也与正演合成方法、钻井数量、井位分布以及模型修改量确定的方法有关,同时也取决于地震测井资料处理和解释工作的精细程度。
借助于这种反演方法,人们可以获得精确描述所观测到的数据集的地下模型。
现行测井资料约束的波阻抗反演技术起源于合成声波测井技术。
80年代后期,Se ym our 等[2]提出了利用地震剖面所过井位的声波测井资料作为约束条件,正、反演结合进行迭代,求取地下波阻抗的方法。
由于这种方法利用了测井资料的高频信息,大幅度拓宽了地震资料的频带,地震剖而的视分辨率得到了很大的提高,因而这种方法表现出了强劲的发展势头,成为国内外众多公司竞相发展的对象。
地震波阻抗反演方法综述
地震波阻抗反演方法综述、地震反演技术研究现状地震反演方法是一门综合运用数学、物理、计算机科学等学科发展起来的新技术新方法，每当数学方法、物理理论有了新的认识和发展时，就会有新的地震反演技术、方法的提出。
随着计算机技术的不断发展、硬件设施的不断升级，这些方法技术得到了实践验证和提升，反过来地震反演技术运用中出现的新问题、新思路又不断促使数学方法、地球物理学理论的再次发展。
时至今日，地震反演技术仍然是一个不断发展、不断成熟、不断丰富着的领域。
反演是正演的逆过程，在地震勘探中正演是已知地下的地质构造情况、岩性物性分布情况，根据地震波传播规律和适当的数学计算方法模拟地震波在地下传播以及接收地震波传输到地表信息的过程。
地球物理反演就是使用已知的地震波传播规律和计算方法，将地表接收到的地震数据通过逆向运算，预测地下构造情况、岩性物性分布情况的过程。
地震波阻抗正演是对反演的理论基础和实现手段。
1959 年美国人Edwin Laurentine Drake 在宾夕法尼亚州开凿的第一口钻井揭开了世界石油工业的序幕。
从刚开始的查看地质露头、寻找构造高点寻找石油，到通过地震剖面的亮点技术寻找石油，再到现在运用多种科学技术手段进行油气资源的预测，石油勘探经历了一个飞速的发展历程。
声波阻抗（AI ）是介质密度和波在介质中传播速度的乘积，它能够反映地下地质的岩性信息。
声波阻抗反演技术是20 世纪70 年代加拿大Roy Lindseth 博士提出的，通过反演能够将反映地层界面信息的地震数据变为反映岩性变化的波阻抗（或速度）信息。
由于波阻抗与地下岩石的密度、速度等信息紧密联系，又可以直接与已知地质、钻井测井信息对比，因此广泛应用于储层的预测和油藏描述中，深受石油工作者的喜爱。
70 年代后期，从地震道提取声波资料的合成声波技术得到了快速发展，以此为基础发展的基于模型的一维有井波阻抗反演技术，提高了反演结果的可靠性。
进入80 年代，Cooke 等人将数学中的广义线性方法运用于地震资料反演，提出了广义线性地震反演。
地震反演技术回顾与展望
地震反演技术回顾与展望一、概述地震反演技术，作为地球物理学领域的重要分支，一直以来在油气资源勘探、地质构造解析以及地震灾害预测等方面发挥着关键作用。
该技术利用地震波在地下介质中传播的信息，通过反演算法处理地震数据，进而推导出地下岩层的物理属性，如速度、密度等。
这些属性信息对于深入了解地下构造、识别油气藏以及评估地震风险具有不可估量的价值。
随着科技的不断进步，地震反演技术也经历了从简单到复杂、从粗放到精细的发展历程。
早期的地震反演方法主要基于射线理论或波动方程的一阶近似，这些方法虽然计算效率高，但精度相对较低，难以满足复杂地质条件下的勘探需求。
随着计算机技术的发展，基于全波形反演、多属性联合反演等高精度反演方法逐渐得到应用，这些方法能够更准确地刻画地下介质的物理属性，为油气勘探等领域提供了更为可靠的依据。
地震反演技术仍面临诸多挑战。
一方面，地震数据的采集和处理过程中不可避免地存在噪声干扰和信号衰减等问题，这些问题会严重影响反演的准确性和稳定性。
另一方面，地下介质的复杂性以及地震波传播的多路径效应也给反演工作带来了极大的困难。
如何在保证计算效率的同时提高反演的精度和稳定性，是当前地震反演技术研究的热点和难点。
展望未来，随着计算机技术的持续进步和人工智能等新技术的应用，地震反演技术有望实现更大的突破。
一方面，高性能计算技术的发展将为地震反演提供更为强大的计算支持，使得更复杂的反演算法得以实施。
另一方面，人工智能技术的应用将有助于提高地震数据的处理效率和反演的准确性，例如通过深度学习等方法对地震数据进行智能降噪和增强，以及通过机器学习等方法优化反演算法等。
随着多源多尺度地球物理数据的融合利用以及大数据、云计算等技术的引入，地震反演技术有望进一步拓展其应用领域和深化其研究内涵。
地震反演技术作为地球物理学领域的重要技术手段，在油气勘探、地质构造解析以及地震灾害预测等方面具有广泛的应用前景。
面对当前的挑战和未来的机遇，地震反演技术的研究和发展需要不断创新和突破，以更好地服务于人类社会的可持续发展。
正演和反演的区别
正演和反演的区别
在不同学科正反演定义有微小差别,但大体上一致的.举个几个的例子：
在地球物理中,已经知道地球介质的性质,如（地震波传播速度等）,求地震波的走时（即地震波在地球中的传播时间等）这是正演；又如,已知地下介质结构,物质特性等,求波速,重力值,电磁学,地热学上的等一些指标因素,这叫正演.
反过来,已经知道如地震波的传播速度,让求地球介质的性质等,这叫反演；又如,已经知道某地的重力异常值,反求该地区地下的物质特性等,是反演.
在遥感影像解译中,反演就是不知道影像上的地面物体是什么,而根据光谱信息等反求地面物体.
总之,一般来说,在实际应用中,反演用的比较多,因为往往人们都是想根据手里的资料反求未知的东西.但是正演是反演的基础,只有有了大量正演知识,你才能正确的反演.。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
正演方法
一、解析方法研究
大斜度井和水平井的随钻电磁波正演可归结为电磁波场量的求解。
在忽略井眼泥浆和侵入带影响的情况下,倾斜井眼周围的感应电磁场可以分解成水平磁偶极子(HMD) 和垂直磁偶极子(VMD)源产生的场的叠加。
Hardman R H和Shen L C通过建立双界面三层水平层状地层模型，研究任意方向的磁偶极子在此模型下的电磁场分布，再通过数值方法求得响应。
Kcnnedy W D等利用赫兹矢量求解Maxwell方程, 得到了任意地层中的电磁场分布,而且从理论上解释了井斜对响应信号的影响,。
肖加奇和张庚骥通过研究有耗层状介质中任意方向磁偶极子电磁波的辐射与传播, 得到一种求解任意层中、任意方向磁偶极子电磁波场量的解析方法。
由于井眼的倾斜, 电磁场呈非轴对称分布。
将发射线圈近似为磁偶极子, 根据电磁场叠加原理, VMD 只产生TE 波, 而HMD 同时产生TE 波和TM 波, 可以将这种电磁波分解为TE V、波、TE H波和TM H波，且彼此间相互独立互不藕合。
利用Green函数分别求解它们在有源和无源区的解析解。
可利用FHT(快速汉克尔变化)或FFHT（快速傅里叶一汉克尔变换方法）对其解析解进行了数值计算。
其木苏荣等在利用与肖加奇类似的地层模型, 应用层状介质中TE 波和TM 波分解技术以及广义反射、透射系数和振幅的递推算法推导出了倾斜地层中电磁场分布的解析解, 并
给出了电磁波振幅的递推公式。
二、数值模拟方法
对于普通直井，随钻正演方法有限元素法（FEM）、有限差分法（FDM）、数值模式匹配法（NMM）、积分方程法（IEM）、逐次逼近法（SAM）、积分变换法（ITM）等。
对于大斜度井和水平井，正演方法主要有以下几种：
快速傅里叶-汉克尔变换法
针对最简单的纯粹的倾斜地层模型, 最有效的途径是使用快速傅里叶-汉克尔变换算法（FFHT）。
随钻电磁波测发射线圈简化为磁偶极子，考虑层状各向同性介质，当忽略井眼的存在时，正演计算问题就简化成了偶极天线电磁波在层状介质中的传播问题，这种问题可以通过汉克尔积分变换来求解。
有限差分法
Yee K S提出了基于交错网格的有限差分方法。
在交错网格中, 电场在单元棱边上采样, 磁场在单元侧面中心采样。
采用交错网格的优点在于, 对于微分Maxwell方程组成立的梯度、旋度和散度关系式对其差分模拟仍然成立。
Wang Tsili和Hohmann利用有限差分方法求取了三维电磁场在时间域中的解。
针对Maxwell方程, 采用交错网格技术求解, 时间步长采用Du Fort-Frankel方法。
边界条件主要包括电磁场的法向分量连续和电流磁流法向连续。
利用该条件定解出电磁场具体表达式,。
避免了导数计算, 并且减小了数值误差, 也不必计算大型的转置矩阵, 因此, 很容易在计算机上实现。
Newman G A等人研究了地层呈横向各向异性时三维感应测井响应问题。
他们利用交错网格有限差分方法近似求解电场的矢量方程, 并采用Krylov子空间迭代法和一种新的预处理方法来提高运算速度和精度。
这种新的预处理方法将电场分解为旋度和散度两部分, 以构造近似逆算子。
这种预处理方法比简单的雅可比预处理计算速度提高了一个数量级。
数值模式匹配法
W.C.Chew提出的数值模式匹配法（NMM）特别适用于求解二维轴对称非均匀介质问题。
该方法的核心是将二维的数值计算转化为一维的解析递推和另一维的数值计算，在轴向上用解析解，在径向上用有限元素法（FEM）求解，同时给出了广特征值问题中矩阵元素的解析表达式，与完全采用数值方法的有限元素法相比，计算速度大大提高。
谱分解法
Druskin V和Knizhnerman提出的谱分解法（SLDM）, 是目前比较实用的一种方法。
该方法建立在时间域和频域中的全局Krylov子空间逼近基础上, 利用交错网格技术分离模型空间。
选用共扼梯度法迭代出上述子空间, 但计算量太大, 通过Fourier变化, 可将空间域转换成频域, 降低了计算量。
SLDM方法的优点是不必每次都计算Krylov子空间的正交基和迭代系数矩阵, 不过它的缺点是不能利用预条件技术进行加速, 收敛速度受到限制。
逆谱分解法
Druskin V L等, 改进了谱分解法提出了逆谱分解法。
因为在中,SLDM不能进行预条件加速, 收敛速度受到限制, 于是Druskin等人提出可以利用A’的伪逆A’-1构造子Krylov子空间K m（此方法称为SLDMINV）。
计算A’-1时, 转变成求解个标量方程（包括个3个Laplace方程和个邪方程, 和1个Possion方程），这在某种程度上提高了运算速度。
在求解Possion方程时, 利用预条件共扼梯度迭代法, 既加快速度, 又提高了结果的收敛性。
汪功礼等人在Druskin的基础上, 利用Krylov子空间不变性求解离散得到的大型稀疏复对称线性方程组。
在构造Krylov子空间时, 使用其系数矩阵的伪逆以改善迭代的收敛性。
迭代中,使用不完全Cholesky分解共扼梯度法求解4个三维方程以得到新的Lanczos向量。
经过迭代不超过20次可得到理想结果。
三、软件
Ansys有限元软件
Ansys程序是一个功能强大的设计分析及优化软件包，其特点有：数据统一，使用统一的数据库来存储模型数据及求解结果，实现前后处理、分析求解及多场分析的数据统一；强大的建模能力具备三维建模能力，仅靠GUI的（图形界面）就可建立各种复杂的几何模型；强大的求解功能，提供了数种求解器，用户可以根据分析要求选择合适的求解器；强大的非线性分析功能，具有强大的非线性分析功能，可进行几何非线性、材料非线性及状态非线性分析；智能网格划分，具有智能网格划分功能，根据模型的特点自动生成有限元网格。
良好的优化功能；良好的用户开发环境。
FEPG系统
FEPG系统采用组件化的程序设计和人工智能技术，根据有限元方法的数学原理及其内在规律，以类似于数学公式推理的方式，由计算机自动产生有限元程序，使用者只需编写微分方程表达式和算法表达式，就可自动生成fortran源代码。
该系统突破了目前通用有限元程序只用于特定领域和特定问题的限制，适合于各种领域的各种工程与科学计算问题，免去了大量繁琐的有限元编程劳动，并保证了程序的正确性和统一性。
它可根据具体问题自动产生合适需要的有限元程序，可用来分析电磁场的多方面的问题，使得复杂的电磁场问题得到了大大的简化。
COMSOL软件
COMSOL Multiphysies是全球第一款真正的多物理场耦合的数值仿真软件，它是以有限元法为基础，通过求解偏微分方程(单场)或偏微分方程组(多场)来实现真实物理现象的仿真。
用数学方法求解真实世界的物理现象，COMSOL Multiphysies以高效的计算性能和杰出的多
场双向直接耦合分析能力实现了高度精确的数值仿真。
COMSOL Mutiphysies软件系统中的电磁场软件包提供三种应用模式：导电介质中电流场(emdc)模式、似稳场(emqa)模式和电磁波(rfw)模式。
emdc应用模式适合各种普通电测井、侧向聚焦电测井在导电媒质中电流场的仿真计算。
emqa应用模式满足感应测井、聚焦感应类测井和阵列感应测井响应的建模和仿真计算要求。
rfw应用模式满足电磁波传播电阻率测井、电磁波传播介电测井仪器测量响应的建模和仿真计算要求。
本文采用该软件分析随钻感应仪器在三维地层中的响应特性。
反演方法
反演方法分为两大类线性反演和非线性反演。
线性反演方法包括广义线性反演、Born近似反演、Rytov近似反演、最速下降法;共扼梯度法；牛顿法、高斯牛顿法；脉冲谱方法、广义脉冲谱方法等。
这类方法都是将反问题局部近似为一个线性问题,其优点是收敛速度快，不足之处在于它们都是局部收敛的方法,迭代结果明显依赖于初始猜测的选取,容易陷入大量局部极值的陷阱,很难捕获到令人满意的全局最优解。
线性化反演方法
根据是否需要线性化,是否具有全局搜索能力,非线性反演又可以分为:线性化(或拟线性化)反演和完全非线性反。
完全非线性反演不进行问题的局部线性近似,而是将反演问题转化为一个非线性优化问题(或非线性算子方程问题)，通过各种途径直接求解，实现从数据空间到模型空间的映。
许多学者认为完全非线性反演是解决非线性反问题的根本方法。
完全非线性反演方法包括穷举法、MonteCarlo法、混沌优化法、模拟退火法及遗传算法等等。