企业级的网络设计方案的规划与测试
中海工业（江苏）有限公司作为一个大型的造船业基地，负责船舶工业的各项任务，而各项任务作为一个整体又需要各种数据信息的交流整合，因此公司的信息化建设是对于各项业务正常稳定进行的通讯保证，信息化网络的设计必须充分考虑到数字化使用和智能化管理的要求，遵循功能性、实用性、兼容性、前瞻性、安全性、扩展性和经济性的原则。素,并对大型企计算机企业网的建设提出了相应的方法和策略。在总结以往网络设计方案的基础上,提出了一种基于交换机组网模式和路由器组网模式相结合的网络设计方案,通过对核心层网络设计、汇聚层网络设计和接入层网络设计的方案设计的总结,使得网络设计方案能够将两个地域上分开的部分,采用通过各自的核心交换机连接的方式,保证线路和设备接口的稳定,并且保证数据能安全和可靠地传输。在计算机和通信內络系统的结构设计攴网络规划设计时,常常要进行大量的仿真分析,本文介绍的就是一种先进的网络仿真软件OPNET。文中介绍了 OPNET的组成、特点及其网络仿真的方法与步骤,并且给出了 OPNET在局域网规划和设计中的应用企业网,网络规划和设计,络仿真, OPNET第页第页第页附件四上海交通大学学位论文原创性声明本人郑重声明:所呈交的学位论文,是本人在导师的指导下,独立进行研究工作所取得的成果。除文中已经注明引用的内容外,本论文不包含任何其他个人或集体已经发表或撰写过的作品成果。对本文的研究做出重要贡献的个人和集体,均已在文中以明确方式标明。本人完全意识到本声明的法律结果由本人承担。学位论文作者签名;支训口期:208年2月必日附件五上海交通大学学位论文版权使用授权书本学位论文作者完全了解学校有关保留、使用学位论文的规定,同意学校保留并向国家有关部门或机构送交论文的复印件和电子版,允许论文被查阅和借阅。本人授权上海交通大学可以将本学位论文的全部或部分内容编入有关数据库进行检索,可以采用影印、缩印或扫描等复制手段保存和汇编本学位论文保密口,在年解密后适用本授权书。本学位论文属于不保密M(请在以上方框内打学位论文作者签名,士X指导教师签名:龙多日期:208年10月28日日期>1年0月28为适应企业信息化的发展,满足日益增长的通讯需求和网络的稳定运行,如今的大型企业网络建设比传统企业网络建设提出更高的要求,主要表现在如下几个方间:)现代大型企业网络需要高带宽和高性能,以满足企业员工工日益增长的通讯需求。随着计算机技术的多样化发展,基于网络的各种应用日益增多,如今的企业网络已经发展成为一个多业务载平台,它不仅要继续承载企业的办公自动化和浏览等简单的数据业务,还要承载涉及个业生产运营的各种业务应用系统数据,以及带宽和时延都要求很高的电话、视频会议等多媒体业务,因此数据流量将大大增加,尤其是对核心网络的数据交换能力提出前所未有的要求)现代大型企业网终需要更高的可靠性和实时性,用以保障企业生产运营的止常进行。随着企业各种业务需要在计算机网络上进行操作,如何保障企业网络无屮断运行已经成为保证企业正常的生产运营的关键。现代人型企业网络在可靠性设计方面主要应从三方面考虑:首先是设备级可靠性设计,这里不仅要考察网络设备是否实现了关键部件的冗余备份,还要从网络设备整体设计架构、处理引擎种类笭多方面去考察;其次是业务的可靠性设计,这里要注意网络设备在故障倒换过程中是香对业务的正常运行有影响;再次是链路的可靠性设计,以人网的链路安仝来自于它的多路径选择,所以在企业网络建设时要考虑网络设备是否能够提供冇效的链路自愈手段和快速重路由协议的支持现代大型企业网络需要提供完善的端到端保障,以满足企业网多业务承载的需求。大型企业网络承载业务的不断增多,单纯的提高带宽并不能够有效的保縫数据交换的畅通无阻,所以如今的大型企业网络建设必须要考虑到网络应能够智能的识別应用事件的紧急和重要程度,如视频、音频、数据流,同时能够调度网终中的资源,保证重要和紧急业务的带宽、时延、优先级和无阻塞的传送,实现对业务的合理调度才是一个大型企业网络提供高效实用服务的保障。)现代大犁企业网络应提供更完善的网络安全解决方案,以阻击病毒和黑客的攻击,减少企业的经济损失传统个业网络的安全措施主要是通过部署防火墙、、杀毒软件以及配合交换机或路由器的访问列表来实现对于病毒和黑客攻击的防御,但实践证明这些第页被动的防御措施并不能有效的解决企业网络的全问题。在企业网络已经成为公司生产运营的重要组成部分的今天,现代企业网络必须要有一整套从用户接入控制,病毒报文识别到主动抑制的一系列安全控制于段,才能有效的保证企业网终的稳定运行。网络设计人员通常有能力创建出一个符合要求的企业网络,但当问题出现的时候,他们却不能采用构建网络时的思维来解决这些问题。每一次对网络的升级、打补」以及修改都会进一步增加网络的复杂性。这些问题可能导致网络难以被人理解,也不易于故障排除。随着时间的推移,可能导致网络不如预期的那样运行良好,随着网络规模的不断增长而不能很好的扩展,并且不能达到客户的需求解决这一问题的方法就是采用规范化、系统化的于段来构建企业网络。自顶向下的网络设计采用了结构化系统分析的思想。结构化系统分析的主要目标是能够更准确的描述用户需求,但在实际工作中,用户的需求常常被设计人员忽略或理解错误。另一个目标是将项目分解成更易被维护与修改的模块,使其便于管理。结构化系统分析有下列特性网络按照自顶向下的顺序设计。在设计期间,现有的许多技术和模型可以被用来刻划现有的系统、新的用户需求和将来的系统结构等特征。●理解薮据流量、数据类型以及能够访问或修改数据的进程。理解企业屮各阶层用户的定位和需求。●逻辑模型必须在物理模型之前被开发岀来。逻辑模型描述了功能划分和系统结构的基本组成模块。物理模型描述了设备的特定的技术和实现。对亍大型的网络设计项目,模块化显得十分必要。设计应该依功能划分,从而使项目更容易管理。在本文中,网络设计主要分成四个阶段,并以周而复始的方式执行1)需求分析。在这个阶段内,网终分析人员通过与用户和技术人员当面交流来获取对于新的或升级系统所需的商业和技术目标的一种理解。然后归纳岀当前网络的状况,包括逻辑和物理拓扑结构以及网络性能。这个阶段的最后一步是分析现在和将来的网络通信量,包括流量以及负载、协议行为和服务质量(QoS)的要求。2)逻辑设计。这个阶段需要设计新的或者升级网络的逻辑拓扑结构,网络层的编址,命名,交换机路由选择协议。逻辑设计也包括安全规划,网络管理设计,以及对服务提供商能够满足的广域网和远程访问的需求作初步调查。3)物理设计。在物理设计阶段,实垗逻辑设计的具体的技术和产品被选择第页出来。此外,在逻辑设计阶段开始对服务提供商的调查,必须在这个阶段完成。4)测试、优化、和记录文档。自顶向下网络设计的最后步骤是编写和实施测试计划,建造一个原型或者试验环境,优化网络设计,把所有网终设计的提议以及工作进程绘制成文档。当用户反馈的网络监测建议对网络加以提高或者需要增加新的应用时,这些主要网络设计阶段将不断周而复始地运行。图1-1显示了网络设计和实施循环周期。分析需求监测优化逻车设计网络性能实施及测试网终物埋设计测试,优化并记录设计图1-1网络设计和实施循环周期首先,作为该工程项目的参与人员之一,在阅读了大量企业网设计与规划方面的资料与文献后,逐渐认清了目前我国大型企业网络建设的研究背景及研究意义,然后和导师及同事们分析探讨了影响大型企业计算机企业网建设的主要因素,并对大型企业计算机企业网的建设提出了相应的方法和策略,在整个设计过程中釆用结构化的网络设计方法,将这个设计过程分为需求分析、逻辑设计、物理设计、测试等四步骤。其次,这个项目的客户是一家船舶制造单位,在详细研究了客户的招标文件中技术方的需求后,根据其业务流程表现出了其特定的网络建设的需求,我方在和客户充分的交流后,将客户的技术需求转化为各项网络技术目标,并得到了客户的认可再次,多次实地考察,结合实际情况和客户的需要,提出适合于该企业的组网模式,兼顾网络性能、安仝性、可管理性、可扩展性以及成本控制。在总结以往网络设计方案的基础上,提出了一种基于交换机组网模式和路由器组网模式相结合的网络设计方案,通过对核心层网络设计、汇聚层网络设计和接入层网络设第页
- 2020-12-05下载
- 积分：1
三维孔洞储层建模及其地震波场正演模拟
三维孔洞储层建模及其地震波场正演模拟，理论讲解很透彻，分析思路清晰1290地球物理学进展26卷预测结果,即从具有确定性资料的控制点(如井点)解释.同时,利用该过程中产生的6冂井的时深关出发,推测出点间(如井间)确定的、唯一的储层参系,用三角剖分网格建立了速度模型并经过井点校数随机模拟是从一个随机函数z(v)中抽取多个可正,实现了对工区构造框架的时深转换通过以波阻能实现,即人工合成反映Z()空间分布的可供选择抗为协变量的孔隙度属性模拟,借助于三维可视化的、等概率的高分辨率实现技术,我们可以大致看到孔洞储集体的形态、分布、对于该工区来说,三维地震资料分辨率较高,对规模及连通性(图1).根据孔隙度发育情况,我们将孔洞储层已经有一定的反映(常表现为低频、不连续强储层分为孔洞欠发育(致密)、孔洞较发育(较致密)振幅反射).通过岩石物理分析又发现孔洞储层低速、和孔洞发育(较疏松)三种类型并分别设计了各自相低密,常规波阻抗反演能够刈其几何形态、空间接触应的弹性参数,同时以模拟出的孔洞形态约束弹性关系定量表征.因此本文将波阻抗数据体作为协变波正演模拟时孔洞体的空间分布量地震属性纵向等值法),采用确定性的协克里金2波动方程正演模拟原理插值算法,对孔洞储层的物性参数进行了三维建模反演所起到的作用,是通过归一化的测井曲线对碳酸盐岩岩溶风化壳孔洞型油气藏属于一种典原始地震数据进行校正,使数据在空间上得到了有型的、复杂的非均质范畴,可以视为由准均匀介质中效的平衡,从而使孔洞反映的更清楚;二是在地质建呈不规则分布的、大小和形状各异的低速体共同组模过程中通过宏观控制,充分利用空间变量的相关成的非层状储集体.在地震剖面上看到的储集体的性,克服低频模型的不足提高属性模拟的分辨率波应是这些低速体的散射(绕射)波.若利用常规波2【V动方程正演模拟方法所使用的均匀介质中的声波方N()程或弹性波方程,难以得到具有复杂非均匀性的孔cline洞型油气藏的地震波场响应2.因此,本文采用非均匀横向各向同性弹性介质中的弹性波波动方程进l()行正演模拟计算,取z轴为垂直对称轴,它可以表示为如下的一阶方程组:a0awta(λ+2naU)+AW(1)7(1(λ1+21)2+λax图1地质模型孔洞储集体俯视图Fig 1 Top view of the cavity reservoiμ(ain geological model其中:(U(x,z,t),W(x,z,t)是速度向量;B(x,z)是密度ρ=g(x,z)的倒数,或者叫疏度;r建模过程中最大的难点是求取准确的速度场,τ(x,z,t),za=x(x,z;t),n=rn(x,z,t)是应木文首先收集整理了工区内6口探井、评价井的钻力张量.A,P/和A1,p分别为水平和垂直方向上井分层数据及多种测井曲线(电阻率、声波时差、密的拉梅系数;为一新的弹性常数可见,在横向各度、自然伽玛等),对其进行了归一化和环境校正,并向同性弹性介质中,独立的弹性常数有五个,它们是制作了合成记录.通过与井旁地震道对比准确标定密度、在垂直方向上的纵、横波速度及纵、横波的各了前中生界侵蚀顶面(15)石炭系双峰灰岩顶面向异性系数,即:(T5。),中下奥陶统顶面(T7)、下奥陶统蓬莱坝组顶面(T7),它们都是区域性的波峰反射在此基础+2uL, Vsi-A pL上采用25m*25m测线密度对该区块6.25km2的3U地震数据体T5°、T5、T74、T78层位进行了精细√λ1+2,CSV闵小刚,等:三维孔洞储层建模其地震波场演模拟1291在具体的有限差分解法上,除了规则网格外,非均匀介质模型的弹性波动方程正演模拟特别是种较为先进的交错网格(图2)最早由 madariaga提当每一个波长中的网格点数多于10个时, Levander出, Virtex在模拟各向同性介质屮SH波和P(1988)2的结果显示,网格色散与网格各向异性均Sⅴ波时也使用了这种差分网格,其精度为o(△2十可忽略不计△x2),在不增加计算工作量和存储容间的前提下,假设U,W分别为介质在x,z两个方向上的速和常规差分网格相比局部精度提高了4倍,且收敛度分量的离散量,R,T,H分别为rxr=和τx的离速度也较快. Levander2又将这种差分网格的精度散量,Lo,M,L1,M1和M2分别为y,kM,A⊥?P1和提高到o(△t2+△x). Crase2则发展了精度可达任g的离散量在(1)式中,各导数项均用中心差分来意阶的高阶交错网格法,但其计算量和内存要求比代替,在如图2所示的一个交错网格中,U,B在节低阶有限差分法大幅度增加.本文使用的是 Virieux点1处计算;W,B在节点2处计算;R,T,M,L,(1986)1的交错网格差分公式,其差分精度为和M1,L1在节点3处计算;H,M2在节点4处计算(△2+△x2), Ikelle l t和 Yung$ K(199)21说这样()式离散为4明该算法可以糈确、稳定地应用于任总复杂变化的=U+B,(R年,-R…)十B,(I1+-H),wH,n-v++B++△(r}一rn)+By2(T+-+),ry=对++(n+2M4)+,△m+-RW(2)+T+,;+(L1+2M1)△tW+U2)+M2△t鲁←z以将震源函数直接赋在rx和n的节点上来模拟震源,即Soure,t)=R(t)t_(source_x, soure_x, t)=R(t.此外,在震源没有激发之前地下介质内部所质点都是静止的,包括质点振动速度为零和所受应力为零.因此,初始条件为图2一个交错树格Fig. 2 A st0,r(x,z,t)=0(t≤0)(3)对于自山表面边界条件,本文采用了模型空间其中,上标k为时间t的离散量,下标i,分别为x的上部加空气的条件,然后再采用吸收边界条件把和z的离散量.△,△x,△z分别为t,x,z的步长空气上边界的弹性波吸收掉,对于空气的下界面,则鉴于 Ricker子波对地震波的分辨率较其它子作正常的分界面来认识,从而获得和实际应用中波函数高,因此,震源选用 Ricker子波,其形式为所采用的地表放炮、地表接收达到一致的效果.R(t)=[1-2Lmf(t-to)] Jexp[-(rf(E to))2]有限差分法在求解波动方程时,会产生不期望式中f表示子波的主频,t为子波持续时间,t为f的数值频散或称网格频散,导致数值模拟结果分辨的函数,在模拟地下激发的地震波时,有限差分交错辛降低2所谓数值频散实质上是一种因离散化求网格中的正应力x和x=是在同一节点上赋值的,解波动方程而产生的伪波动,这种频散既不同于波而vr和vn在节点处的数值并没有参与计算,因此可动方程本身引起的频散,也不同于因波传播的速度1292地球物理学进展26卷频率和角度变化而引起的频散,它是有限差分方法果我们在这里仅分为三种类型:孔涧欠发育(致求解波动方程时所固有的本质特征,无法避免.为了密)、孔洞较发育(较致密)和孔洞发育〔较疏松).消除这种数值频散,前人进行了大量研究,他们的结论是基木一致的,即为了消除数值频散,在使用二阶表1地层框架内各层物性参数有限差分方法时,每个功率对应的波长至少必须使Table I The properties ot each layer用11个网格点,面四阶有限差分则可用二阶差分网in stratigraphie framewor格点数的一半.木文采用的稳定性条件,即计算稳定p(nu/s) v(m/s) (kg/m3)的离散参数区域为151:r4G界面2500三叠系)以Lmd2m≤1(2m-1)fT50界面下石炭系顶)~T46l730≤Ld2m≤T56界面(2m-1)!(巴楚组顶)~T5023102350其中,T74界囿(下奥陶系顶)~T56±8002470T78界面(蓬芠坝组颠)灬T7460002650界面以3702此外,在做波动方程的模型计算时,由于只能在对于试验工区的每条线,其长度均为1625m个有限区域进行,而弹性波在其计算边界上能量为了侏证该区域内均为满叠、孔洞的绕射波收敛以衰诚并不为零,从而产生很强的边界反射,这是模型及边界吸收较为干净,我们在模型的左边延长了计算时所不希望的,需要做人工吸收戌衰减处理,计1200m,右边延长了1575m(延长部分的地层接触算吸收边界的方法有许多种,一般情况下网格周围系并不代表真实情况),即模型总长度为1、4km,的耗散采用质点的速度和应力值乘上一个小于1深度范围为4000~6500m每条线均采用了同样的的因子来平滑的衰减;另种可能性是在网格周围观测系统,具体为:采用零相位对称雷克子波作为震使用低Q值来实现吸收作用,但是后者的吸收效果源(主频40玎z),单边放炮(共20炮,每炮128个检不如前者的吸收效果好,因此本文采用的是第一种波器接收)炮间距160m,检波器间距20m,8次叠方法,具体实珧时釆用了〔 eran等的吸收边界条加,最小偏移距0m,最大偏移距2540m,记录长件实现边界吸收1.6s,Δt=2ms,第一炮的坐标位置为(-1200,0)exp[-a2(I-i)2],1≤i≤1.基于差分稳定条件,取模型中最小介质速度2500m/s其中,I为给定的吸收边界带总节点数;i为吸收边为参考,得到的计算参数为:网格剖分间隔3m界内的节点编号;a为衰减系数,其值的选定与1的3m,时间延拓步长为0.27ms,每个波长(62.5m)大小密切相关,且对吸收效果的影响很大本文中Ⅰ内有20.8个网格.我们一共对33条线进行了正演取为40(即围绕计算区域,再向外设置宽度为40个模拟,图3展示了较为典型的 inline2585线(位于研网格的条形吸收区域)a=0.305/40,i取从0~40究工区的中心部位,地层接触关系以及孔洞体的分节点号.在条形吸收区域中的每个网格结点处,对全布相对比较复杂),从中可以大致看出二维正演模拟部的5个波场量(U,W,R,T,H),在每计算一个时的普遍情况与孔洞体波场响应特征的一般规律问步长后,都做少量的波场减表2展示∫该条线上各孔洞体的几何及物性参3模型计算数,其中④号属于欠发育(充填致密物)类型,①③⑥号属于较发育(充填较致密物)类型,②⑤号属于发在正演之前,我们统计了工区的速度、密度资育(充填较疏松物)类型.此外,建模过程中,我们还料,为了简化模型,并使得孔洞体的地震响应特征更考虑了线与线之间地层起伏渐变、孔洞大小渐变孔具有针对性,我们采用了背景为常速介质、蜜度参数洞物性参数渐变的过程,即所有建模因素都渐变由( arner公式计算的思路(表1).对于孔洞储集的而不是突变,最终保证了三维地震数据体的连体,根据钴井揭示和前面提到的孔隙度属性模拟结续性4期闵小刚,等;三维孔洞储层建模及其地震波场正演模拟1293表2各孔洞储集体的几何及其物性参数最大振幅,且绕射曲率与反射曲率相同,表明二者具Table 2 The geometry and propcrty parameters有不同的传播速度;每个绕射波可分为左右、上下f each cavity reservoir正、反向绕射分支,正向绕射分攴的相位与反射P孔润体尚度宽度vVP波相同,反向绕射分支的相位反转180°,与反射P(m)(m)(m/s)(m/s)(kg/m3)充填物屮心距界面(m)波的相反17396500029002503较致密105弹性波正演模拟生成的炮域合成记录被导人10113480027822470较疏松6FOCUS软件进行常规处理,包括速度谱拾取、动82784500029002500较致密85校、切除、增益、滤波、叠加、偏移和变面积、变密度显①575520030222530致密示等.由于在观测系统中只设计了8次覆盖,为了增⑤18115480027822470较疏松104加速度谱拾取精度,本文采用了由相邻的7个CDP2714850029002500较致密86道集混合构成一个超道集的办法,隔10个CDP拾图4是该模型在590ms时的波场快照,其波场取一个速度文件,并在拾取前先作常规NMO校正清晰,网格频散小,边界吸收较干净这表明,在求解切除,使得原始道集记录能量更强、信噪比史高二维弹性波动方程时,将差分解法和交错网格技术图5、图6分别是TK610井、TK623井所在位置处相结合,通过较好地使用吸收边界和稳定性条件可CDP道集记录及其速度谱,从图中可见各个反射界以显著削弱数值频散,有效地提扃计算精度.同时面的同相轴清晰可辨,对应的能量团集中,而在合成在保证一定的精度前提下,可以采用铰大的空间网记录上T7界面下孔洞所在位置处都有一明显的格间距,提高计算效率.从图巾还可看出,孔洞绕射同相轴,能量团也比较集中,由于TK610升比波和反射波在绕射点处相切,在切点处绕射波具有TK623井孔洞储集体更为发育(尽管二者振幅相1200-80004008001200160020002400280040004505500图3主测线2585地质模型Fig 3 Geclcgical model ol inline258-12004004008001200160045000.10.3图4主测线2585在590ms时波场快照Fig. Snapshot of wave field at 590ms in inlinc25851294地球物理学进展26卷Sg224-230CDP49 SE QNO250030003500400045000.240.60.60.80.8TE1.01. 01.21.2141.4图5TK610井所在位置处CDP道集记录及其速度谱ig. 5 The CDP gather and velocity spectrum at well TK610Sgl58-1640.2ONO250030003500400045000.0.40.60.60.8081.0:1.01212623(2565图6TK623井所在位置处CDP道集记录及共速度谱Fig 6 The CDP gather and velocity spectrum at well TK623当,在地质模型设计时均认为是充填较疏松物,但相消),使得T7界面断断续续,并在该界面下出现TK610井比TK623井在目的层段的厚度要大,横些“短反射”通过仔细分析,我们发现“短反射”中向展布范围也更宽测试产能更高),在合成记录上较强者出现的时间,与孔洞位置相对应.从该模型的孔洞对应的同相轴振幅更强、波形更连续,速度谱上偏移剖面上(图b)可以看到,所有的孔洞体均得到能量团也更强、更集中比铰好的偏移成像,并表现为负正负三个相位的图7是处理完后的叠加和偏移剖面.从叠加剖波形.但鉴于反射波地震勘探的纵向分辨率(大于面上(图a)可以比较清楚的看到孔洞体顶、底的两1/4波长),所有能检测出的孔洞或孔洞组合在叠加组强反射,但是二者之间出现具有绕射特征的弱波剖面上都叠合在T74界面下第一个波峰轴上,在偏代替了成层的背景,这些绕射波的相互下涉(相长、移剖而上都体现在T7界面下第一个黑椭圆体上,4期闵小刚,等:三维孔洞储层建模及其地震波场正演模拟1295601001401802202603003401001401802202603003400.00.00.20.2040.40.60.60.81014露9.926c+089.926e1081022e+091022e+09图?主测线2585对应的叠加剖而(a)和偏移剖面(b)Fig. 7 The stacking section (a) and migration section (b) of inline2585至」其下的“串珠”是孔洞的假象(孔洞组合与围岩(b)之间的多次波及绕射波经偏移归位后形成较强短反L2560L2580L2600射).由于T74界面反射波与沿纵横向有一定分布的孔洞(比较明显的是①、⑤号)的绕射波叠加,使得30001300040孔洞所在位置处T74界面反射波能量变弱,而孔洞底部与下覆围岩之间的正极性反射由于受T7界200600面反射波的负值性续至波叠加,也变得较弱.此外,B40080在构造高点上(④⑥号孔洞体所在位置,④号更为明显),由于孔洞引起的绕射与隆起引起的回转波的相6003600100互丁涉,T56和T7界面不连续,甚至在其间出现空白反射,而实际资料也有这种情况.这说明,对于塔3800800200和油田碳酸盐岩孔洞储集体这类特殊的油气储层来400040001400400说,在解释时遇到层位问断时,不能轻易地开断层,而应该综合考虑构造、孔洞绕射等地震波场特征.这图8联络测线2795实际剖面(a)和正演剖面(b)对比也是塔河油出勘探开发实践中发现“表层弱反射、内Fig 8 Comparison of the actual section (a) and幕强反射”地震特征对应有利储层的一个佐证forward modeling section (b)in crossline2795依据以上思路与工作流程,我们得到了33条沿主测线方向的二维偏移剖面,在并成三维体之前,为而正演模拟釆取的是8次叠加、道间距10m),正演了尽量消除线与线之问因地层起伏造成的不闭合,剖面较好的反映了实际情况.这不仅体现在层位的我们采取先把33条线的速度文件并成三维体,整体形态、分布比较相似(由于速度取了平均,各层的厚平滑后两用每条线对应的、平滑后的速度对其原始度不一致,但不影响我们的主要的,即对孔洞体地共中心点数据进行动校、叠加、偏移的办法,得到33震响应特征的分析),更重要的是,我们所设计的孔条新的二维偏移剖面,再并成一新的三维体,此外,泂体,其位置、形状规模、振幅强弱均与实际地震资由于正演模拟数据体线间距为50m,道间距为料具有相当好的对应关系,这表明我们在止演模拟10m,其空间采样率比实际资料低,本文编制了相和处理时的设计思想和参数选取原则是合理的,这应的算法在频率域对其进行插值,使线间距加密到结果也为我们进行后续工作提供了比较好的数据25m图8是联络测线2795在时间域的实际剖面源由于实质上是2.5维,不是基于面元的真三维,(a)和正演剖面(b)对比,排除二者在采集时的一些所以沿联络测线的剖面上同相轴有抖动现象,这是差异(如实际三维采集资料为24次叠加道间距25m,不可避免的)1296地球物理学进展26卷4结论与建议2]谢桂生,刘洪,赵连功,伪谱法地震波正演模拟的多线程并行计算[冂.地球物理学进展,2005,20(1);17~23.本文从三维角度,建立了与实际资料比较吻合Xie G S, Liu H, Zhao L G. Parallel Algorit hm based on the的孔洞储层模型,并进行了弹性波正演模拟,总结了multithread Technique for pseudospectal modeling of seismic地震响应规律,主要结论如下:wave[J]. Progress in Geophysics(in Chinese), 2005, 20(1)1)结合地震资料建立储层地质模型能够有效[3]刘财,张智,邯志刚,等.线性粘弹体屮地震渡场伪谱法模拟地降低储层模型的不确定性,提高建模精度.同时利技术[门].地球物理学进展。:005,20(3),640~644,用协克里金技术,用波阻抗反演的确定性信息约束Liu C, Zhang Z, Shao Z G, et aL. Pseudo-spectral forward储层的平面非均质性,可以实现孔隙度属性的确定modeling nf seismic wave in linear viscoelasic solid [J]P1性建模),2005,20(3):640~644.4」张智,刘财,邵志刚,伪谱法在常Q粘弹介质地震彼场模拟(2)在求解二维弹性波动方程时,将差分解法和中的应用效果[].地球物理学进展,2005,20(4):945交错网格技术相结合,通过较好地使用吸收边界和949,稳定性条件可以显著削弱数值频散,有效地提高计Zhang Z, Liu C, Shao G. The application of pseudo-spectral算精度.同时,在保证一定的精度前提下,可以采用forward modeling of seismic wave field in constant Q较大的空间网格间距,提高计算效率该方法具有广viscoelastic medium [J]. Progress in Geophysics, 2005,20(4)945~949泛的适用性5]盖良国,马在出,曹景忠,等.一阶弹性波方程交错网格高阶(3)孔洞储集体在偏移剖面上表现为负-正-负差分解法[冂].地球物理学报,200,43(3):411-~419三个相位的波形,但只能确定奥陶系风化面下第Dong LG, Maz T, Cao j Z, et al. A staggered-grid high个负相位是孔洞的发育位置,其下的“串珠”是孔洞order difference method of one-order elastic wave equation]的很象.风化面反射波与沿纵横向有一定分布的孔Chinese J. Geophys. (in Chinese),2000,43(3):411-419洞体的绕射波叠加,使得孔洞所在位置处风化面反[6]董艮国,马在田,曹景忠,一阶弹性波方程交错网格高阶差分解法稳定性研究[门].地球物理学报,200,43(6):856~射能量变弱,而孔洞底部与下覆闱岩之间的正极性反射由于受风化面透射波的负值性续至波叠加,也Dong L G, Ma Z T, Cao J Z. a study on stability of the变得较弱.该结论对于实际地震资料处理、解释以及staggcred-grid high-order difference method of first-order储层预测烃类检测具有普遍的指导意义elastic wave equation. Chinese J. Gcophys. in Chinese)2000。43(6);856~864本文不足之处主要有三点「7]萤良国.复杂地表条件下地震波传播数值模找1.勘探地球(1)在三维孔隙度建模时采用的是常规阻抗信物理进展,2005,28(3);187~194息(约束稀坑脉冲反演),其纵向分辨率不够(只能分Dong L G. Numerical simulation of seismic wave propagation辨1/4波长以上的孔洞储集体),在后续工作中将尝under complex near surface conditions [J]. Progress in试使用地质统计学反演的阻抗体来约束建模以大幅Exploration Geophysics(in Chinese), 2005, 28(3):187--194提高纵向分辨能力[8奚先,姚姚,二维随机介质及波动方程正演模拟[.石油地球物理劫探,2001,36(5):546-552(2)在弹性波正演模拟时采用的是2.5维思想XiX, Yao Y. 2D random media and wavc cquation forward口前正在研制全三维算法有望更逼真的还原孔洞储modeling [J]. Oil Geophysical Prospecting in Chinese集体的真实地下情况001,35(5);546~5523)考虑到缝的各向异性更为复杂,本文尚未涉9]奚先,姚姚,二维粘弹性随机介质中的波场特征分析[刀地及,对于碳酸盐岩中这类油气运移的重要通道,将在球物埋学进展,2004,19(3):608~615今后的工作中进一步研究Xi x, Yao Y. The analysis of the wave field characteristics in2-D viscoelastic random medium LJ. Progress in Geophysics参考文献( References):hinese),2004,19(3):608~[10]奚先,姚姚,二维横各向同性弹性随机介质中的波场特征1]刘文岭.大庆宋芳屯油田芳2区块地震与地质资料综合储层J.地球物理学进展,2004,19(4):924~932地质建模研究(博土论文儿D1.北京:中国地质大学,2002Xi x, Yao Y. The wave field characteristics of 2-DLiu W I. A Study on Reservoir Geological Modeling withclo].ESeismic and Well-log Data in Fang 2 Area of DaqingGeophysics(in Chinese), 2004,19(4):924-932ongfangtun Oil Field (doctor dissertation)(in Chi[111吴永国,贺振华,黄德济.串珠状溶涧模型介质波动方程正Beijing: CUG, 2002.演与偏移[.地球物理学进展,2008,23(2);539~5444期闵小刚,等:三维孔洞储层建模及其地震波场正演模拟1297Wu Y G, He Z H, Huang d J. Wave equation forward[19]肖玉茹,何峰煜,孙义梅,等,古洞穴型碳酸盐岩储层特征modeling and migration for heads-shaped corroded cave model研究一以塔河油田奥陶系古洞穴为例匚门。石油与天然气地EJ]. Progress in Geophysics(in Chinese), 2008, 23(2): 539质,200324(1):71~80.Xiao YR, He f Y, Sun Y M, et al. Reservoir charactetistics12]股文,印兴耀,吴国忧.高特度频率域弹性波方程有限差分of paleocave carhonates-a casc study of Ordovician paleocave方法及波场模拟[」.地球物理学报,2006,49(2):561in tahe oilfield, Tarim basin UJ]. Oil Gas Geology(inChinese),2003:24(1):71-80Yinw, YinXi,WuGC. The method of finite difference of[20]姚蟋,唐文榜.深层碳酸盐岩岩溶风化壳洞缝型油气藏可检high precision elastic wave equations in the frequcncy domain测性的理论研究[门.石油地球物理勘探,2003,38(6):623and wave-field simulation [J. Chinese J, Geophys.629Chinese),2006,49(2):561~568.Yao Y, Tang W B. Theoretical study of detectable cavern[13]马贵,土尚旭,宋建勇.频率域波动方程正演中的多网格Fractured reservoir in weathered Karst of dccp carbonatite迭代箅法[门].石油地球物理勘探,2010,45(1):15[J]. Oil Geophysical Prospecting(in Chinese), 2003,38(6):Ma ZG, Wang S X. Sun J Y. Multigrid iterative algorithm in623~629,domain wave equation forward modeling [J]. Oil [21] Levander A R. Fourth-order finite difference P-SvGeophysical Prospecting(in Chinese ) 2010, 45(1): 1-5seismograms []. Geophysics, 1988, 53(11): 25-36.[14]张金海,王卫民,赵连锋,等.傅里叶有限差分法三维波动[22] Crase e. Iligh- order( space and timc) finite-difference方程正演模拟[.地球物理学报,2007,50(6):1854A, In: 60th SEG Annual1862C].1990:987~991.Zhang j H, Wang W M, Zhao L F, et aL. Modeling 3-D [23] IkelleL T, Yung SK, Daube F. 2-D random media withscalar waves using the Fourier finite-difference method.ellipsoidal autocorrelation function [J]. Geophysics, 199350(6):1854[24]奚先.随机介质模型的构造及其波场模拟(博土论文)[D][15] Qin Z, Lt武汉:中国地质mproved NPML absorbing boundary condition in elastic waveXix. Construction and scismic wave field modeling ofmodeling [J]. Applied Geophysics, 2009. 6(2): 113-121random medium model doctor dissertation ) in Chinese)[16][D].Wwave equation [J. Geophysics, 1986, 51(1): 54-61[25]吴国忱,王华忠.波场模拟中的数值频散分析与校正策略[17] Virieux J. P-Sv wave propagation in heterogeneous mediaLJ.地球物理学进展,2005,20(1):58-65velocity-strcss finite-difference methud LJ]. GeophysicsWu GC, Wang H Z. Analysis of numerical dispersion in1986,V51;889~901.wave-field simulation [J]. Progress in GreaphysiEs ( in18] Igel H, Riollet B. Mora P. Accuracy of staggered 3-D finiteChinese),2005,20(1):58~65difference grids for anisotropie wave propagation [J]. 62th [26] Cerjan C, Kosloff D, Kosloff R, et al. A nonreflectingAnn, Internat, Mtg, Soc. ExpL. Geophys, Expboundary condition for discrete acoustic and elastic- wav1992,1244~1246.equation []. Gcophysics, 1985, 50(4): 705-708.
- 2021-05-06下载
- 积分：1
|
联系客服（QQ：3324461878）
