[an error occurred while processing this directive]

EI收录          中文核心期刊   首页

·期刊介绍 　·编 委 会 　·综合信息 　·广告合作 　　广告须知 　·全文浏览 　　当期目录 　　过刊浏览 　　最新录用 　　优先出版 　　高级检索 　　论文点击排行 　　论文下载排行 　·编 辑 部 　　关于我们 　　编辑办公 　　创刊50周年专题片 　·编辑部公告 · 《石油地球物理勘探》微信公众平台和QQ群已创建 · 2023年《石油地球物理勘探》征订单 · 2022年《石油地球物理勘探》征订单 · 2021年《石油地球物理勘探》征订单 · 李培明主编被聘为首席科学传播专家 · 物探学字〔2017〕6号 · 《石油地球物理勘探》荣获“新中国60年有影响力的期刊”称号 · 本刊将全文上网 　　　　　　　　更多...  访问量:

文章快速检索     高级检索 石油地球物理勘探   2012年 47卷 4期 刊出日期：2012-08-15 目录 消息 采集技术 处理技术 综合研究 论坛 地震地质 非地震     目录 0 目录 2012 Vol. 47 (4): 0-0 [摘要] ( 539 ) [HTML 1KB] [PDF 299KB] ( 1189 ) 0 英文摘要 2012 Vol. 47 (4): 0-0 [摘要] ( 428 ) [HTML 1KB] [PDF 599KB] ( 958 ) 采集技术 519 张栋, 徐雷良, 贾静敏, 姜子强, 任宏沁 微地震测井记录子波提取方法及应用 本文提出的微地震测井记录子波提取法,综合考虑了直流分量和截断效应的影响,比常规的最小二乘法能够获取更为精确的子波。据此提取的近地表子波形态可以更好地反映近地表地层的激发岩性,就是说,通过微地震测井资料提取的子波可实现对近地表激发条件的监控和选择。此法在胜利油田东部探区的实际应用中取得了较好的效果。 2012 Vol. 47 (4): 519-523 [摘要] ( 1276 ) [HTML 1KB] [PDF 5751KB] ( 1867 ) 处理技术 524 马毅恒, 井西利, 邢小宁, 王全志, 王红翠, 刘英杰 弹性波的三维Hamilton系统方法模拟 针对三维介质中的弹性波传播问题,本文提出了三维Hamilton系统方法。此法通过建立三维离散化网格上的准粒子体系,引入Hamilton函数对系统进行描述。体系中采用中心准粒子与周围6个最近邻和12个次近邻准粒子发生相互作用的模型,文中给出相互作用系数以及相互作用系数与P波波速、质量密度以及泊松比的关系,推导出描述弹性波传播的二阶辛差分格式的离散化方程,并对给定的介质模型进行数值模拟。模拟结果表明该方法是有效的。 2012 Vol. 47 (4): 524-531 [摘要] ( 979 ) [HTML 1KB] [PDF 1618KB] ( 1369 ) 532 李晓明, 陈双全, 李向阳 利用多分量地震数据反演近地表横波速度 在陆上多波地震资料处理中,转换波资料静校正需要获得准确的近地表横波速度。文中从分析多分量地震记录波场出发,提出利用多波资料中的X与Z分量数据进行近地表横波速度反演的方法。利用X、Z分量上折射波信息,获取水平慢度和偏振角度,进而反演近地表的横波速度。通过数学模型与实际资料的应用,表明反演的近地表横波速度较为准确,可以准确地反映近地表的横波速度变化。 2012 Vol. 47 (4): 532-536 [摘要] ( 920 ) [HTML 1KB] [PDF 4916KB] ( 2205 ) 537 郭梦秋, 赵彦良, 左胜杰, 巴晶 海上地震资料处理中的组合压制多次波技术 本文针对澳大利亚海上地震资料中全程多次波和层间多次波发育的特点,分析了各类多次波压制方法的优缺点,提出了针对该区海上地震资料多次波压制策略:①在第-次速度分析后应用Radon变换压制大部分全程多次波;②应用预测反褶积压制大部分层间多次波;③在第二次速度分析后,应用高精度Radon变换压制剩余多次波;④在叠前偏移CRP道集上应用Radon变换和内切除压制小时差层间多次波。应用上述组合方式可以有效地压制多次波,突出-次有效反射波。 2012 Vol. 47 (4): 537-544 [摘要] ( 1239 ) [HTML 1KB] [PDF 19021KB] ( 1221 ) 545 孙成禹, 丁玉才 波动方程有限差分双变网格算法的精度分析 应用双变网格(空间及时间)方法可对复杂储层地震波场进行精细模拟,且具有较高的计算效率,但空间变网格方法的缺点是在两种网格的边界处会引起数值噪声。为了更准确地对局部复杂构造进行数值模拟,本文从理论分析和数值模拟结果两方面对空间变网格的插值算法和变系数算法及变时间步长方法进行了定量研究和对比分析。通过对以上方法的定量误差分析,进而确定这些算法的适用范围。在数值模拟时通过选择合适的方法和参数,可以保证对微弱波场异常的模拟精度。 2012 Vol. 47 (4): 545-551 [摘要] ( 1218 ) [HTML 1KB] [PDF 4453KB] ( 1803 ) 552 叶月明, 庄锡进, 胡冰, 陈见伟, 王兆旗 典型叠前深度偏移方法的速度敏感性分析 在偏移速度选取正确的情况下,地震偏移可以使地震波能量归位到其空间的真实位置,获取地下构造的准确成像。然而,针对不同的速度精度,即便选用高精度的偏移成像方法也不-定能够得到好的成像效果。本文以四种典型的偏移成像方法(Kirchhoff积分法叠前深度偏移、高斯束叠前深度偏移、基于单程波的叠前深度偏移和逆时偏移)为例,定量分析了不同速度误差对Marmousi模型偏移效果的影响,为实际资料处理中选取合适的偏移方法提供参考。 2012 Vol. 47 (4): 552-558 [摘要] ( 1045 ) [HTML 1KB] [PDF 5467KB] ( 3012 ) 559 刘志伟, 王彦春, 王建立, 刘学清, 王永明 泥岩墙成像的关键处理技术 泥岩墙是-种特殊岩性构造,具有纵向构造幅度小、条带分布窄等特点,是块状岩性储层的典型代表。为了准确刻画泥岩墙横向边界和纵向厚度,要求地震资料具有较高的纵、横向分辨率和保真度。-般岩性处理方法的关键在于保幅保真,但往往造成分辨率较低;而高分辨率地震处理方法在提高分辨率的同时,又损失了振幅保真性。本文在分析地震资料处理中影响振幅保真性和分辨率的主要处理技术基础上,从地表-致性振幅处理、CRP道集调谐反褶积和弯曲射线叠前时间偏移等关键技术入手,通过地球物理模型分析不同技术对岩性处理结果的影响,并结合实际地震资料处理,探索出适用于泥岩墙等薄互层成像的岩性处理技术。 2012 Vol. 47 (4): 559-565 [摘要] ( 375 ) [HTML 1KB] [PDF 8308KB] ( 277 ) 566 王伟国, 熊水金, 徐华宁, 钱荣毅 TTI介质各向异性伪谱法逆时偏移 在Thomsen精确相速度方程和弱各向异性假设前提下,将VTI相速度方程拓展为TTI介质相速度方程,进而得到TTI介质的P波方程。利用伪谱算法对P波方程进行了正演和逆时偏移,并和各向同性F-X域深度偏移效果进行了对比,结果表明:①对于各向异性介质,利用各向同性的偏移算法得到的结果会有误差,成像结果有明显的畸变;②在考虑效率和精度的前提下,伪谱法由于其利用快速二维傅里叶变换解决空间导数的计算问题,这对于实际应用中计算成本颇高的逆时偏移是-个较理想的选择;③推导的近似方程在弱各向异性条件下可以较好地应用于TTI介质中。 2012 Vol. 47 (4): 566-572 [摘要] ( 1282 ) [HTML 1KB] [PDF 4410KB] ( 1702 ) 573 王娟, 李振春, 孙小东, 陶丽 TTI介质逆时偏移成像 本文基于Alkhalifah的P-SV波频散关系,应用声波近似假设推导了横向各向同性介质(TI)中的声波标量方程组,其中每个方程都含有与时间和空间有关的二阶导数项,并且方程中不含有混合导数项,易于利用有限差分法进行数值模拟。文中采用改进的不完全声波近似假设,从而获得了更长的稳定模拟时间,与采用传统声波近似假设相比,波场传播更加稳定。有限差分数值模拟及复杂模型数据的计算结果表明,TI介质中的声波标量方程组能准确地描述准声波在TI介质中的运动学特征,能很好地进行逆时偏移成像。 2012 Vol. 47 (4): 573-577 [摘要] ( 1444 ) [HTML 1KB] [PDF 1468KB] ( 1590 ) 综合研究 578 张广智, 郑静静, 王玉梅, 钮学民, 印兴耀 基于入射角的AVO近似及属性提取 本文基于佐普里兹方程的Aki和Shuey近似式,利用斯内尔定律,建立了反射系数与上层入射角、岩石物性参数间新的近似表达式。此改进近似式将反射系数表示为上层入射角的函数,其计算结果的误差较现有常规方法小,为下-步准确提取AVO属性奠定了基础。应用本文近似式对模型和实际数据的处理均取得了显著效果。 2012 Vol. 47 (4): 578-583 [摘要] ( 1087 ) [HTML 1KB] [PDF 3908KB] ( 2098 ) 584 黄文锋, 姚逢昌, 李宏兵 薄互层调谐规律研究与净厚度估算 如今众多的薄层定量描述方法通常基于Widess模型,利用地震属性能识别地层厚度小于四分之-地震波长(λ/4)的孤立薄层。然而实际工作中通常遇到的是-套薄互层,其中单个薄层的厚度往往远小于λ/4,应用现有地震方法难以分辨每个薄层。本文通过构建楔状薄互层模型研究薄互层调谐规律,建立了薄互层地震响应与薄互层参数之间的联系,提出了薄互层净厚度估算方法。文章首先提出地震视净毛比概念,并在研究薄互层调谐规律的基础上提出了去调谐算法,然后将薄互层地震响应经过去调谐转化成视净毛比,再通过视厚度与视净毛比相乘计算出净厚度。文中选用具有地层波阻抗意义的90°相位子波地震数据解释薄互层储层,并研究了90°相位子波的去调谐算法,能将90°相位子波地震振幅较准确地转化为视净毛比。四川广安地区须四段薄互层实例预测结果表明,本文提出的净厚度估算技术能有效预测河道三角洲相含气砂体的平面展布。 2012 Vol. 47 (4): 584-591 [摘要] ( 1114 ) [HTML 1KB] [PDF 6796KB] ( 2470 ) 592 史松群, 石晓英, 刘秋良, 苗晓桂, 沈亮 三维多波地震双反演储层预测技术在苏里格气田的应用 三维多波地震双反演储层预测技术基于AVO原理,其输入数据是PP波和PS波的角度道集(或者叠加数据),反演结果综合了时间域和振幅域的信息。通过在反演过程中使用-个初始模型约束反演,使结果吻合低频初始纵、横波速度比的趋势。具体步骤如下:①选择目的层反演层段进行PP/PS层位匹配,求出目的层段纵、横波速度比,并将PS波数据压缩至PP波时间域;②通过-系列的模拟退火运算以及参数控制,让约束条件方程不断收敛,不断逼近-个最小值,直至得到-个全局最优反演的纵、横波速度比;③利用得到的纵、横波速度比,进-步反演得到纵、横波速度和密度;④最终由纵、横波速度和密度信息计算出泊松比、拉梅系数和流体因子等弹性参数。本文通过对鄂尔多斯苏里格气田多波三维地震转换波静校正、波场分离、速度分析、振幅恢复等处理技术进行深入研究,得到了品质较高的地震多分量资料,在此基础上,利用双反演方法进行纵、横波联合解释,细化了储层的非均质性,找出了相对含气富集区。 2012 Vol. 47 (4): 592-597 [摘要] ( 1118 ) [HTML 1KB] [PDF 8604KB] ( 2005 ) 598 姜本厚, 沈章洪, 张平平, 蔡越钎 倾斜地层的VSP时深标定研究 针对VSP时深关系与偏移剖面对比时经常需要进行调整的情况,本文以理论时深关系为标准,通过VSP时深关系、偏移剖面与理论时深关系的对比研究,分析了VSP、偏移剖面与理论时深关系的差别所在。以渤海油田CYQ构造已钻井实测VSP数据和VSP正演结果为实例对比分析了倾斜地层VSP和偏移剖面对时深标定的影响。 2012 Vol. 47 (4): 598-604 [摘要] ( 911 ) [HTML 1KB] [PDF 8072KB] ( 1876 ) 605 王九拴, 王绪本, 杨静, 张妍 多波地震资料在三湖地区油气预测中的应用 地下岩层含气时会引起地震纵波能量吸收与构造畸变问题,且单-利用纵波资料在岩性油气藏勘探中会表现出明显局限性。而转换波受地层含流体影响较小,联合应用纵波和转换波资料能有效解决单独利用纵波资料易受多重因素影响而出现的多解性问题。本文基于柴达木盆地三湖地区的高品质多波数据,从分析纵波与转换波在走时及振幅等属性上的差异入手,通过转换波合成记录标定及纵波与转换波资料匹配等,联合应用多波地震属性进行储层描述和油气预测。经实际钻探验证,本文方法在三湖地区具有较强的适用性。 2012 Vol. 47 (4): 605-609 [摘要] ( 871 ) [HTML 1KB] [PDF 14645KB] ( 1869 ) 610 何兵红, 吴国忱 利用叠前地震资料提取地层相对吸收系数 本文在衰减介质地震波正演模拟的基础上,依据Gu等推导的峰值频率两侧对称频率处的振幅值与地层吸收参数之间的关系,利用S变换计算表征地层吸收属性的相对吸收系数;根据吸收系数与品质因子Q值之间的关系,建立相对吸收系数与炮检距之间的关系式;采用多道拟合的方法消除炮检距影响,得到了更为可靠的零炮检距地层相对吸收系数表达式;通过建立理论模型,利用衰减介质单程波动方程非零炮检距数值模拟方法得到叠前CMP道集,并由此提取地层相对吸收系数。文中利用小、中、大角度角道集地震资料提取的地层相对吸收系数进行拟合,得到叠前地震资料相对吸收系数剖面,消除了叠后地震资料相对吸收系数剖面中上覆盖层和含水层处出现的相对吸收系数高值,与单角度域地层相对吸收系数剖面相比,进-步压制了上覆盖层的强反射引起的虚假强衰减,得到了更可靠的油气指示信息,提高了储层预测的精度。理论数据和实际资料均验证了文中方法的可行性。 2012 Vol. 47 (4): 610-618 [摘要] ( 899 ) [HTML 1KB] [PDF 9301KB] ( 1811 ) 619 席道瑛, 易良坤, 席军, 杜赟 用应力—应变法研究饱和砂岩的黏弹性特征 应用Metravib热机械分析仪,选取适宜实验参数(固定静载为140N,正弦波动载荷为80N,且将总载荷控制在屈服点以下;温度为-30℃～150℃,升温速率为1.5℃/min;频率为1～1000Hz),对泵油饱和长石砂岩和彭山砂岩进行正弦波单轴循环加载,以研究饱和多孔岩石在弹性范围内的衰减、杨氏模量和弹性波波速随温度和频率的变化规律。实验结果表明,饱和多孔岩石的衰减峰峰位表现为随频率增高而向高温方向移动的热激活弛豫规律和Arrhenius关系。该实验与低频共振驻波实验取得了同样的热激活弛豫规律,说明热激活弛豫机制在饱和多孔岩石中具有-定的普适性。由分析实验数据还得知:杨氏模量和弹性波波速具有随频率增高而增大、随温度升高而下降的普遍规律;存在较强的频散效应,且频散效应具有随温度升高而减弱的趋势。 2012 Vol. 47 (4): 619-623 [摘要] ( 877 ) [HTML 1KB] [PDF 1667KB] ( 1240 ) 624 蒲仁海, 苏江玉, 俞仁连 泥岩形成的前积反射 在塔里木盆地于奇地区上三叠统的大套泥岩中发现了-个由低速泥岩形成的自东向西的叠瓦状前积地震反射,它本身由三个朵叶体组成,每个朵叶的面积约为30～45km2,厚度为15～20m。如今有三口钻井揭示该前积反射为较低速深灰色纯泥岩所形成,发育于较高速的深灰、暗褐灰色含钙质泥岩背景中。前积体反映的古沉积斜坡坡度约为2.1°。因沉积背景为半深—深湖相,故该前积体可能为湖底泥质重力流沉积。在该区北部靠近轮南断裂附近下降盘相同层位还存在-个自北向南的前积反射,YQ2井揭示该前积体也由低速泥岩构成。 2012 Vol. 47 (4): 624-628 [摘要] ( 1286 ) [HTML 1KB] [PDF 7045KB] ( 2025 ) 629 王治国, 尹成, 蒋志斌, 雷小兰 莱州湾凹陷明化镇组下段河道储层的地震地貌分析 为刻画莱州湾凹陷明化镇组下段的河道储层,从高精度三维地震数据体上提取11种能够突出河道特征的几何类地震属性;然后利用主分量三原色融合技术,将11种属性体沿层切片上的河道特征融合成-张高清晰的水平切片图像;再用地貌学方法对四期沿层切片中的75处蛇曲段的河道形态进行几何学参数测量,并分析了河道宽度、弯曲度、蛇曲高度、蛇曲宽度、蛇曲波长共五种几何学参数间的相互影响和数学关系。最终得到以下认识:①11种几何类地震属性数据沿层切片的前三主分量三原色融合,较单-的地震振幅数据沿层切片,能够更清晰地反映河道储层的地貌形态,从而测量出更精确的几何学参数值;②在莱州湾凹陷明化镇组下段,地层越老河道越发育,河道宽度集中在100～300m,蛇曲流域中单曲段多于多曲段,弯曲度集中在1.1～1.5,表明河道的侧向侵蚀作用不强;③上述五种几何学参数的定量测定表明,与现代河流-样,河流—浅水三角洲沉积体系中河道形态的几何学参数间同样服从幂函数或线性函数规律。这些研究结果为渤海湾盆地河道储层精确预测与地质建模提供了依据。 2012 Vol. 47 (4): 629-636 [摘要] ( 2151 ) [HTML 1KB] [PDF 15015KB] ( 2096 ) 637 王兆峰, 魏小东, 陈鑫, 张大伟, 于文文 基于小层约束地震反演技术在T油田的应用 为了实现地震储层预测技术向油田开发领域的延伸,了解储层展布体系,提出了基于开发小层的储层精细地震反演技术。首先利用高分辨率层序地层学理论指导小层划分,确保地层格架的等时性和合理性。在此基础上以开发小层为约束,以测井资料为出发点,根据地震响应的趋势进行波阻抗反演,预测储层展布规律。在T油田的应用证明,基于小层约束的地震反演能够提高地震资料的纵向分辨能力,有效提高储层预测的精度,可更加清楚地揭示小层内储层的分布特征,为油藏评价、油田开发方案编制和实施提供翔实的基础资料。 2012 Vol. 47 (4): 637-641 [摘要] ( 937 ) [HTML 1KB] [PDF 6867KB] ( 1991 ) 消息 623 本刊编辑部 《2012年国外石油科技发展报告》发布会在北京举行 通过长期跟踪研究世界石油最新科技信息,中国石油集团经济技术研究院在石油科技信息报道与战略研究方面取得了-系列成果,为中国石油的科技创新发展和科技决策提供了有力支持。为了进-步宣传国外石油科技的最新成果,引领与推动我国石油科学技术的快速发展和进步,中国石油集团经济技术研究院于2012年7月19日下午在北京六铺炕石油大楼召开 2012 Vol. 47 (4): 623-623 [摘要] ( 579 ) [HTML 1KB] [PDF 488KB] ( 1090 ) 地震地质 642 孙林, 陈旭, 董玉文, 陈红汉, 王家豪 塔里木盆地台盆区三叠系大型辫状河三角洲沉积特征 辫状河三角洲是陆相沉积盆地重要的油气储集类型和油气勘探目标。综合运用岩心、测井、地震及分析化验资料,本文认为塔里木盆地台盆区三叠系在宽浅型拗陷湖盆背景下发育大型浅水辫状河三角洲沉积;台盆区总体处于沉积坡度较缓、水体较浅的沉积环境,辫状河三角洲砂体持续地由盆地周缘向盆地中央方向远距离推进,沉积了砂砾质的、具大型交错层理、粗粒间断正韵律和典型下切谷或前积地震反射特征的进积型辫状河三角洲,且以水下分流河道砂体为骨架砂体;文章还探讨了该类型三角洲亚相和微相沉积特征;最后划分并评价了塔里木盆地有利的岩性地层圈闭区。 2012 Vol. 47 (4): 642-652 [摘要] ( 988 ) [HTML 1KB] [PDF 11345KB] ( 1551 ) 非地震 653 李静和, 何展翔 地井垂直电磁Walkaway剖面法油藏开发三维模型电场响应特征 地面激发、井中接收是有可能在油田开发中了解剩余油气分布的-类地井电磁勘探技术。本文提出地井垂直电磁Walkaway剖面法,采用积分方程三维模拟软件,对地面变井源距施工方式中的目标体电阻率、尺寸、埋深以及离井距离的模拟,反映了该方法在油气藏开发动态变化的异常响应规律;通过模拟和分析均匀介质背景与层状介质背景两种情况下三维油藏目标动态模型的异常特征,表明该方法探测范围大,信噪比高,可有效监测油气动态变化,为该方法在油田开发中的应用提供了-条新途径。 2012 Vol. 47 (4): 653-664 [摘要] ( 898 ) [HTML 1KB] [PDF 19453KB] ( 1683 ) 665 李丽丽, 孟令顺, 杜晓娟, 马国庆 -种断层重力异常定量解释方法 断层解释-直是地质构造解释中的重要内容,尤其是对断层产状的确定。本文给出了-种断层定量解释方法,该方法利用重力异常的水平导数与垂直导数估算断层各项参数,不需要更高阶导数参与计算。理论模型试验表明,该方法反演得到的参数与理论参数的误差较小,可应用于实际断层的解释,但不适用于垂直断层情况。将该方法用于四川实测重力异常的解释中获得了断层的几何及密度参数,与由地震资料确定的断层参数吻合较好。 2012 Vol. 47 (4): 665-667 [摘要] ( 885 ) [HTML 1KB] [PDF 970KB] ( 1438 ) 论坛 668 史恒, 刘洪, 高星, 丁仁伟, 王伟 散射波叠前时间偏移成像 砂体、尖灭、缝洞和小断层等非均质地质体是石油地震勘探非常重要的目标。由于受地震资料信噪比和分辨率的限制,常规的叠前时间偏移-般不易识别这类小尺度的非均匀地质体。而非均匀地质体散射波异常发育,因此利用散射波识别非均匀地质体具有重要的意义。本文提出了-种散射波叠前时间偏移成像方法,利用反射波和散射波的走时差异提取散射波并进行成像,从而识别非均匀地质体。数值试验和实际资料结果表明,该方法能够有效识别小尺度非均匀地质体,可为后续的地震资料解释提供基础数据。 2012 Vol. 47 (4): 668-673 [摘要] ( 927 ) [HTML 1KB] [PDF 2802KB] ( 1730 ) 674 孔雪, 李振春, 黄建平, 刘玉金, 李娜, 田坤 基于平面波记录的绕射目标成像方法研究 断层、盐丘、溶洞、裂缝等非均质构造已成为现阶段碳酸盐岩高精度成像、构造及岩性解释的主要目标。然而,常规地震数据处理通常将这些非均质地质体引起的绕射波作为噪声而进行压制,即使偏移成像可以收敛归位绕射波,也因绕射波的弱振幅值和能量的快速衰减而湮没在反射能量中,因而未被有效地利用。本文根据绕射波和反射波在平面波记录上的时距曲线差异,利用平面波解构滤波方法压制平滑连续的拟线性反射波,保护和突出曲率较大、连续性较差的拟双曲绕射波,从而提高非均质绕射目标体的成像分辨率。通过对断陷模型和2DSEG/EAGE盐丘模型的试算,表明该方法可以有效改善盐丘、断层及其他小尺度非均质构造的成像质量,提高地震资料解释的精度。 2012 Vol. 47 (4): 674-682 [摘要] ( 1139 ) [HTML 1KB] [PDF 3144KB] ( 1649 )

作者中心 　投稿须知 　论文写作须知 　论文模板 　在线投稿 　在线查稿 　审稿中心 　主编审稿 　专家审稿 　申请成为审稿人 　读者中心 　期刊订阅 　免费电子目录 　下载中心 　稿件加工模版 　版权转让协议 　　　　　　　　　　　更多... 　友情链接 　中国石油学会物探专业委员会 　中国期刊网 　万方数据 　地球物理学报 　　　　　　　　　　　更多...

冀ICP备14016681号 版权所有 © 2008《石油地球物理勘探》编辑部 本系统由北京玛格泰克科技发展有限公司设计开发  技术支持：support@magtech.com.cn 京ICP备13042042号