0 引言

相控处理就是按照沉积微相对某些地球物理参数的控制规律，对地震反演得到的参数剖面和数据体做进一步处理，以获得沉积微相剖面和数据体。相控的目的就是把地震反演的参数剖面按一定规律转换成能反映沉积微相展布的沉积相特征剖面，当该特征剖面所表征的优势沉积微相与储层含油气岩性相对应时，就可用沉积微相定量预测岩性的分布[1-3]。从理论上讲，按照每种沉积微相与地球物理参数的关系和相控数学规律，可把每种沉积微相的参数信息从参数剖面或数据体中分别提取出来，以实现“沉积微相的分离处理”[4-5]。

在川东南地区，用地震波主特征参数和波阻抗联合进行沉积相的预测研究和沉积相的划分已取得初步成效；但由于沉积相形成的复杂性，目前该技术方法还很难对沉积相进行准确的定量预测和划分，其在实际生产处理中有局限性。相控技术目前主要有2种方法，其一是用沉积相对油气有利区块进行约束，其二是用已知井的沉积相来控制未知相带的划分，这两种方法对沉积相正确的划分都很重要[6-7]。

目前沉积相的研究基本上采用“三相”法，也就是将地震相的模式与区域地质的沉积规律进行对比,然后再结合单井相的分析结果实现目的层段的地震相向沉积相转换。该方法是勘探阶段纵向有利勘探层和有利勘探位置比较有效的方法。利用地震进行沉积相的预测，有3种预测方法：①根据利用已知井统计出各种沉积微相的类型与地层速度或波阻抗的关系，而后经过数学方法拟合得到有利油气藏的沉积相与这些参数随着深度变化的规律，直接用参数采用多次滤波来预测沉积相；②利用统计得到的优势沉积相地震波型与剖面的地震波型进行相关的分析，由于沉积相变化会引起地震波型特征的变化，当沉积相与速度、波阻抗等地层参数存在一定的关系时，地震波型与沉积微相就存在对应关系；③在井旁的地震记录中，提取出优势相所对应的地震子波，经研究发现沉积相的变化和子波的变化存在一定关联度，因此利用该优势相的地震子波对所要研究的剖面进行相关分析，也可预测这种优势相。本文在川东南官渡地区采用的相控处理方法是第三种，即子波相关沉积相预测方法[8]。

为了说服他改变主意继续留在队上，我提出的问题或者解决方法必须要赢得包某的认同，但是出国施工对比留在国内，无论从收入还是休息上都有很大的优势，很难让人不动心。他家中只有一个女儿，经济上并不困难，所以我分析他想出国并不是因为缺钱。我只有另辟一条道路，那就是替他分析一下出国的缺点，让他觉得出国在获得高薪的同时，还有不少不可避免的要面对的问题，而这些问题是他不能不考虑的因素。

1 子波相关沉积相预测方法

(1)子波振幅谱计算

由图2可知，合格的水泥材料在显微镜下可观察发现大量的黑灰色玻璃体，同时黑色的煤渣等杂质较少。由于水泥的生产中允许加入少量的混合材，故还可以观察到少量粉煤灰形态的球形体。

理论研究表明，在地震波反射系数是白噪音的情况下，地震记录功率谱与子波功率谱相等。通常将多道的地震记录平均功率谱当作是子波的功率谱。实际的子波提取中，常常采用多道的平均相关谱作为子波的振幅谱。

则有

2.要做好水质调控，培育加州鲈鱼苗的水质不能过肥，池水透明度应保持在25～35cm，养殖期间要保持清新的水质和较高的溶氧量，每10～15天需换水1次，每次换水量约占池水的1/3，如能采用微流水或增氧措施养殖，则效果更好。

rxixj(t)=Tww(t)

(1)

在式(1)中，rxixj(t)是第i道地震记录的互相关，Tww(t)是子波的自相关。频率域中，式(1)为：

|Rxixj(w)|=|Rww(w)|=

(2)

因此，应用多道平均相关谱统计来提取子波振幅谱公式是：

1.2 文字图像二值化 经过灰度处理的彩色图像还须经过二值化处理将文字与背景进一步分离开。所谓二值化，就是将彩色值图像信号转化成只有黑（0）和白（1）的二值图像信号。二值化效果的好坏，会直接影响灰度文本图像的识别率。二值化处理的具体经验有很多〔11-12〕，方法大致可以分为局部阈值二值化和整体阈值二值化。本文参考基于k中心点聚类的图像二值化方法〔13〕，采用局部阈值二值化，将方块苗文图像分割为互不相交的若干个w×w的小窗口，在每个窗口中求得所辖像素的灰度平均值，经添加惩罚项后作为初步阈值进行二值化处理，继而对各小方块进行二值化变换，各个阈值的计算公式如下：

(3)

求互相关函数的时候，应该对该函数作指数衰减，为：

主特征参数的提取克服了傅氏分析长时窗及瞬时分析刻画局部信息的不足，主特征参数计算可提取子波的主极值，能够更好地反映地质体岩性、物性等引起的差异。

(4)

式(4)中，Txixi+1(τ)是作指数衰减后的互相关函数，τ为相关的变量，β为衰减的系数。对函数作指数衰减的主要目的是为了使地层的反射系数自相关更加逼近δ(t)函数。

(2)子波相位谱计算

在通常子波提取的方法中，一般假设子波为最小相位或者是零相位。若子波为零相位，则：

φw(ω)=0

(5)

爱国主义始终是中华民族团结统一、自强不息的精神源泉。弘扬爱国主义精神，是实现中国梦的内在要求，是赢得具有许多新的历史特点的伟大斗争的现实要求，是应对经济全球化时代各种挑战的必然要求，不仅如此，还有助于维护世界和平发展、构建人类命运共同体。

(6)

式(6)中，φmin(ω)是最小相位子波的相位，|B(ω)|是子波的振幅谱，P是柯西(Cauchy)主值。利用式(6)可求得最小相位子波相位谱，再利用相位值和振幅谱就能够求得最小相位的子波实序列。

然而，实际的子波通常都是混合相位的。对于这种情况，采用了展开相位计算来计算地震道的连续相位谱。展开相位谱的值可表示为：

arg[x(ejω)]=Arg[X(ejω)]+2πle

(7)

等式成立，则所求展开相位为：

因此可直接利用多道平均方法去求子波振幅谱，经过傅氏变换后可得到零相位子波。假设子波满足最小相位，那么根据最小相位实序列因果性，该子波的相位谱和该子波对数振幅谱互为Hilbert变换对：

图3为最小相位子波和零相位子波不同噪音下的反演结果，可见噪音对反演结果的影响不大。

(8)

在得到地震道的相位谱连续曲线后，对多道相位谱曲线进行求和取平均即可得到子波的相位谱。N道地震记录的平均相位谱可以表示为：

(9)

可以证明：

式(9)中，φxi(ω)是第i道地震记录的展开相位谱，φHi(ω)是第i道地震记录中的反射系数的相位谱成分，φwi(ω)是子波相位谱。

(10)

可以证明：

另外，院区内的各间病房均配备了独立卫生间和阳台。卫生间靠外墙设计，并设有窗户，这既能采集自然光线，又能通风透气，保持病房内空气清新，有效防止异味滋生。病房里的阳台使用大落地玻璃门，自然光能直接透射进病房里，患者可以按需要在阳台进行舒展活动，呼吸新鲜空气，配合周边的花槽，令患者充分享受大自然的阳光与绿化，保持良好的心态，有助于患者病情的控制及治疗。

(11)

显然，地震道数越多，统计特性就会越好，所求得的子波相位值精度也会越高。

课题主要采用行动研究的方法，以11级机制1班学生为研究对象进行研究。在此基础上实施规划、展开调研、收集数据、开展行动研究，并最终形成生产性综合实习课程标准。再对其效果进行新一轮的验证以进一步检测前期成果，进而得出研究结论，撰写相应的研究报告。

最后，将所求得的子波振幅谱与相位谱通过反傅氏变换后就可得到子波的序列，它是在没有任何相位特征假设条件下从多道地震记录中统计提取出来的子波。

3.HBO治疗方法：HBO治疗采用面罩吸氧法，压力为0.20～0.25 MPa (2.0～2.5 ATA)，1次/d，吸氧时间为1 h，以5～10 d为1个疗程，每个疗程之间休息1～2 d，共治疗60次。治疗期间发生中耳调压不良或感冒者，给予呋麻滴鼻液，仍不能缓解者可以酌情暂停HBO治疗3～5 d；若HBO治疗中止时间超过10 d仍无法继续再次行HBO治疗者，剔除出组。

2 子波提取实验

图1为最小相位子波提取实验，第1道是反射系数序列，第2道是理论子波，第4道是反演子波，第3道和第5道分别是理论子波和反演子波的合成记录。可以看到反演子波与理论子波的误差很小，理论子波和反演子波的合成记录误差同样很小。

本文以油石比5.4%为最佳油石比，采用浸水马歇尔试验和冻融劈裂试验进行AC—10型常温彩色沥青混合料的水稳定性验证。由试验结果得出，AC—10型常温彩色沥青混合料的水稳定性能验证结果：残留稳定度为99.8%，冻融劈裂强度比为79.2%，满足规范的要求。

图1 最小相位子波反演实验 Fig.1 Inversion test of the least phase wavelet

图2 零相位子波反演实验 Fig.2 Inversion test of zero phase wavelet

图2为用零相位子波反演实验，第1道是反射系数序列，第2道是理论子波，第4道是反演子波，第3道和第5道分别是理论子波和反演子波的合成记录。可以看到反演子波与理论子波的误差很小，理论子波和反演子波的合成记录误差同样很小。

Arg[x(ejω1)]+2πle(ω1)

图3 最小相位子波(a)和零相位子波(b)不同噪音下的反演结果 Fig.3 Inversion of the least phase wavelet (a) and zero phase wavelet(b) under different noise condition 1.加10%噪音；2.加20%噪音；3.加30%噪音；4.加40%噪音

3 主特征参数剖面相特征及井旁子波提取

Txixi+1(τ)=Txixi+1(τ)·exp(-βτ)2

图4为373测线二维主振幅剖面，图中官8井(Guan8井)中有3个边滩沉积相，分别对应时间剖面中的1 043 ms、1 077 ms和1 108 ms，且时间厚度分别是12 ms、10 ms和14 ms。3个边滩沉积相分别对应主振幅剖面1 000～1 100 ms时间段3个强能量值。图5为从官渡地区373测线中提取的滩相子波。

图6为386测线二维主振幅剖面，图中官16井(Guan16井)中有4个边滩沉积相，分别对应时间剖面中的1 086 ms、1 093 ms、1 098 ms和1 112 ms，时间厚度分别是7 ms、2 ms、3 ms和6 ms。4个边滩沉积相分别对应主振幅剖面1 000～1 200 ms时间段中部的4个主振幅，可看到井中边滩沉积相的主振幅呈中弱能量，横向的变化较大。图7为从官渡地区386测线中提取的滩相子波。

4 实际应用

官渡构造带位于四川盆地川东南赤水地区东部，构造带总体呈NNW向展布，系川东高陡褶皱带

图4 373测线过官8井主振幅剖面(横向为cdp号，纵向为时间ms) Fig.4 Section of the principal amplitude of well Guan 8 through which the line 373 passes

图5 373测线中提取的滩相子波 Fig.5 Wavelet of floodland extracted from line 373

图6 386测线过官16井主振幅剖面(横向为cdp号，纵向为时间ms) Fig.6 Section of the principal amplitude of well Guan 16 through which the line 386 passes

图7 386测线中提取的滩相子波 Fig.7 Wavelet of floodland extracted from line 386

图8 4个井统计的滩相平均子波 Fig.8 Average wavelet of floodland of 4 wells

向南西方向帚状延伸的一部分。通过前期勘探及油气成藏研究，官渡地区具有良好的成油地质条件，烃源岩发育，陆相、海相烃源岩兼有。从目前的钻井分析，有官3、官8、官9、官10等多口井见油气显示，尤其在陆相地层中已获工业油气流。本区储层较为发育，有海相、陆相多套储层，陆相主要有下沙溪庙组和须家河组碎屑岩储层，海相主要有雷口坡组、嘉陵江组孔隙-裂缝型储层以及茅口组岩溶缝洞型储层。根据官渡地区须家河组有利油气藏的沉积相主要分布于边滩中，本文利用子波相关方法预测沉积相，预测控制边滩砂体的分布，这对砂体十分发育的以河流相沉积为主的工区有重要的研究意义和实用价值[6]。图8是官渡地区须家河组4个过井剖面地震记录中统计得到的滩相平均子波，应用该优势相的地震子波对所要研究的剖面进行相关处理，就可预测这种优势相。

图9为过官16井(Guan16)NE386测线须家河组剖面经过子波沉积相控处理后得到的微相剖面，图中红、黄色表示的是优势滩相(边滩)区域，红色表示的是优势滩相(边滩)特征的一级，黄色表示的是优势滩相(边滩)特征的二级。从剖面可以看到边滩相较薄，且边滩相的横向走向连续性较好，符合沉积规律。图中黄色线圈定的部位就是官16井发现气层的位置，可以看到在该位置有一层不错的边滩相。

图10为须家河组四段经过子波沉积相控处理后得到的微相平面切片展布规律图，边滩相主要分布于须家河组四段，图中的红色与黄色区域分别表示的是一级和二级边滩相，可看出在须家河四段边滩相的分布区域。

5 结论

(1)该方法只能用于同一个沉积环境的地区，且随着地层深度的不同，各沉积相对应的地震相也会不同，因此，对不同的地层进行沉积相预测时，必须提取该层位的地震子波。

(2)在同一沉积环境区域，如果有多个井，可以对不同井同一沉积层位提取出来的地震子波做相关处理，把相关性好的地震子波求和取平均值作为沉积相预测的子波。

(3)在利用子波相关沉积相预测方法时，准确提取到优势滩相是关键，但由于钻井原始资料的缺失及中下部地层的钻井少等原因，要想较准确地得到不同目的层的滩相子波较困难。

图9 NE386须家河组子波沉积相预测剖面 Fig.9 Predicted section of wavelet sedimentary facies of Xujiahezi formation from NE line 386

图10 须家河组四段三维沿层沉积相切片图 Fig.10 3D slicing diagram along layer of the 4th member of Xujiahezi formation

(4)子波相关沉积相预测只是众多储层预测方法之一，在该地区实际的储层预测工作中，要综合沉积相、地震相、相控处理、储层特征、裂缝预测、地震属性和波组抗反演结果，有利于对川东南地区有利储集层区块进行划分及评价。

9年前，奈飞内部曾经流传出一份有关它企业文化的PPT，遭到疯传，桑德伯格称它为“硅谷最重要的文件”。硅谷人都相信，这份文件中一定蕴藏着奈飞成功的魔法。今年，这份PPT的升级版本——由奈飞前任首席人才官帕蒂·麦考德撰写的《奈飞文化手册》颠覆重启，精英们正在追寻的答案和真相就在这里。

窦性心动过缓可发生于健康人群,也可由窦房结的结构功能改变或迷走神经兴奋的生理或病理因素引起[8]。有研究表明窦性心动过缓对心律失常或其他心血管疾病也有一定影响[9]，可导致突发性心脏性猝死[10] ，临床医生常常需要判断窦缓的临床意义，识别出高危人群加以干预治疗。

(5)论文研究的技术方法已形成可应用于实际生产的储层预测软件，可处理显示任意数据格式的三维数据体和二维数据体。在本次沉积相预测方法中，井旁采用反演方法取得优势相子波，无井剖面采用展开相位谱方法提取子波，反演可以得到正确的井旁子波。实验研究中取得的反演子波是稳定的，无井剖面展开相位谱方法提取的子波符合地震波的传播规律。

参考文献：

[1] 高林, 杨勤勇, 地震属性技术的新进展[J]. 石油物探, 2004, 43(增刊): 10-16.

[2] 王开燕, 徐清彦, 张桂芳, 等. 地震属性分析技术综述[J]. 地球物理学进展, 2013, 28(2): 0815-0823.

[3] 张延玲, 杨长春, 贾曙光. 地震属性技术的研究和应用[J].地球物理学进展, 2005, 20(4): 1129-1133.

[4] Barnes A E. Seismic attributes past, present, and future[J]. Expanded Abstracts of 69th Annual Internet SEG Mtg.1999: 892-895.

[5] Chopra S, Marfurt K J. 地震属性的历史回顾:上[J]. 赵剑敏, 译. 油气地球物理, 2006, 4(3): 60-65.

[6] 施泽进, 王长城. 官渡三维区相控分析及储层预测[R]. 成都: 成都理工大学, 2007.

[7] 施泽进, 王长城. 川东南地区重点层位储层沉积相研究[R]. 成都: 成都理工大学, 2009.

[8] 夏庆龙, 赵志超, 赵宪生. 渤海浅部储层沉积微相与地球物理参数关系的研究[J]. 天然气工业, 2004, 24(5): 51-53.