琼东南盆地天然气水合物及其成藏模式的海洋可控源电磁研究

南海琼东南盆地是天然气水合物重要远景区之一.由于盆地大部分地区海底地形平缓、地层近于水平,增加了利用地震反射剖面识别似海底反射（BSR,bottom simulating reflector）的难度,从而影响了对水合物的评价.为了进一步开展琼东南盆地水合物调查研究,本文在研究海域进行了海洋可控源电磁探测试验,将自主研发的10台接收机以500 m的间距,投放至水深约为1360 m的海底,完成了一条4.5 km剖面的电磁数据采集.通过对采集的数据进行处理与二维（2D）反演,获得了研究剖面海底的电阻率断面图.反演结果显示,研究区海底60~330 mbsf（meter bottom of seafloor）的地层中,存在多个横向不连续分布的高阻异常体,电阻率介于2~10 Ωm之间;在海底330 mbsf之下,横向上发育了电阻率为2~4 Ωm的3个高阻体.根据研究区热力学条件,本文估算了生物成因气与热成因天然气的水合物稳定带（GHSZ,gas hydrate stability zone）厚度,结合高阻体的分布特征推断了地震剖面上BSR的位置.在此基础上,对反演的电阻率断面进行解释,推断了研究区水合物的分布及游离气运移通道.研究表明,勘探区具有形成天然气水合物矿藏的地质与地球物理条件,其成藏模式可能属于"断层、裂隙输导的下生上储型",水合物的气源为生物成因气.     

南海天然气水合物远景区海洋可控源电磁探测试验

为了测试我国自主设计与研发的海洋可控源电磁仪器性能及其在水合物探测中的适用性,本文从海洋可控源电磁法基本原理出发,首先根据试验海域水合物地质特征,建立简化地电模型开展理论研究,确定海洋可控源电磁试验的技术方案;利用研发的海洋可控源电磁仪器,在南海天然气水合物远景区开展探测试验,首次获得了我国深水海域的海洋可控源电磁数据.通过对采集数据进行处理与反演,获得了试验剖面的海底电性结构模型,揭示了4号测点下方存在一个50 m厚的高阻层,其电阻率为25 Ωm、顶部埋深为181 m,为该区天然气水合物调查提供了有价值的电性参考资料.研究结果表明,自主研发的海洋可控源电磁仪器性能达到了预期的设计指标,这标志着我国海洋可控源电磁探测技术向实用化进程迈出了重要一步.     

勘查天然气水合物资源的海洋可控源电磁发射系统

海洋可控源电磁探测(MCSEM)是勘查天然气水合物资源的有效方法之一.海洋可控源电磁发射系统是实现MCSEM的重要硬件组成部分.本文以天然气水合物在海底的一维地电模型为例,研究电偶源发射频率、发射偶极侧向漂移、源偶极矩大小对电场响应的影响,进而指导发射系统研制.该系统包括:甲板升压控制单元,可将船载大功率电能升为高压并通过万米光电复合缆输送至海底的发射机;水下变压器,可将深拖缆中的大功率高压转换为低压;拖曳式大功率电磁发射机,其内部的控制电路硬件和嵌入式驱动软件可将水下变压器输出的电能逆变为大功率矩形脉冲,并通过发射偶极将脉冲发送至海水介质中;借助水动力学设计的发射机拖体,用于装载发射系统水下部件和保持拖曳过程中的平衡与稳定;甲板端上位机监控单元,利用万米光电复合缆中的光纤实现船上计算机与海底发射机的远程数据通信.2012年5月和2013年5月海洋试验的结果均表明,所研制的发射系统可作为天然气水合物资源勘探的有效激励场源.     

新型坐底式海洋可控源电磁发射系统及其海试应用

利用拖曳式海洋可控源电磁发射系统在探测埋藏较浅的天然气水合物资源时,会遇到一些新问题.如拖曳式轴向发射偶极源不能贴紧海底,发射偶极源和海底间海水的电磁衰减,使得10 Hz以上的相对高频能量难以导入至海底以下介质;发射偶极源有可能出现水平或垂向摆动;拖曳式发射时,时间窗口内叠加的数据有限;拖曳路线可能与海底构造走向平行,不利于揭示探测目标体的异常形态.新型坐底式发射系统有望解决上述问题.坐底式发射系统的硬件部分包括甲板端供电和监控单元、长距离电力和数据通信单元,以及水下发射机主体.新型坐底式发射系统的发射电极紧贴海底,没有海水层衰减,有利于相对高频的人工源电磁信号经过海底以下介质传输至接收端;通过超短基线信标更容易精确定位发射机拖体,利用姿态方位参考系统确定发射电极的供电方向,可以精确校正电偶源水平偏角的影响;增加单点供电时间,提高接收信号信噪比;提供两对电偶极源,从两个相互垂直的方向对异常体进行人工源激发.2015年海洋试验的结果表明,坐底式与拖曳式发射系统联合作业,可从多角度对海底异常体进行宽频带电磁激发,为数据采集提供丰富的场源信息.

     

海洋天然气水合物可控源电磁法的响应特征:一维正演模拟

We discuss the feasibility of using controlled-source electromagnetic （CSEM） in the frequency domain for prospecting marine gas hydrates. Based on the Ocean Drilling Program （ODP） Leg 164 log data, we have established several 1-D resistivity models which have different gas hydrate concentrations. Meanwhile, we analyzed the electromagnetic response of marine gas hydrates in the frequency domain based on these models. We also studied the relationship between electrical field magnitude or phase and parameters such as receiver-transmitter distance and frequency. Our numerical modeling results provide us with a quantitative reference for exploration and resource evaluation of marine gas hydrates.     

海洋可控源电磁法发射源姿态影响研究

海洋电磁理论是建立在水平发射偶极的基础上.然而,由于海底洋流作用海洋可控源发射偶极可能产生水平摆动、倾斜和水平旋转,造成测量数据与理想状态的水平电偶极子的电磁响应有较大的偏差.为了研究发射源姿态变化对场的影响,本文提出利用欧拉旋转将发射偶极由源坐标系转换到地球坐标系,进而分析发射源姿态变化造成海洋电磁观测数据的误差分布...     

海洋可控源电磁法发射源姿态影响研究

海洋电磁理论是建立在水平发射偶极的基础上.然而,由于海底洋流作用海洋可控源发射偶极可能产生水平摆动、倾斜和水平旋转,造成测量数据与理想状态的水平电偶极子的电磁响应有较大的偏差.为了研究发射源姿态变化对场的影响,本文提出利用欧拉旋转将发射偶极由源坐标系转换到地球坐标系,进而分析发射源姿态变化造成海洋电磁观测数据的误差分布特征.研究发现,对于发射源水平摆动和倾斜造成的海洋电磁观测数据的畸变,其分布规律取决于接收机所在区域.以发射频率1 Hz为例,对于由发射源水平摆动或倾斜导致的地球坐标系中的两个等效电偶极子分量产生的水平电场Ex,在旁线区域,二者的相位差约为360°,接收机接收到的电场信号约为两个偶极子场之和;而在同线区域,二者的相位差接近于180°,接收机接收到的信号为两者之差.在两个区域的交界处,发射源姿态对电磁信号影响最大.相比之下,发射源水平摆动只改变发射源的位置,不会改变源的性质,其影响主要发生在收发距较小的区域.本文的研究成果有助于海洋电磁系统设计时确定发射源姿态变化的参数指标,同时也为未来实际海洋电磁观测时分析和校正发射源姿态变化造成的观测数据误差提供理论基础.     

海底可控源电磁测量电路的Linux驱动程序

国外的水合物探测试验证明,可控源电磁法是探测天然气水合物的有效方法之一.开展海底可控源电磁法工作的重点之一在于研制高精度、高可靠、低功耗观测海底可控源电磁信号的设备.为可靠地观测海底可控源电磁场宽频带、大动态范围信号,设计了以Linux与ARM9相结合的测量电路,并开发了Linux2.4内核下的设备驱动程序.本文分别论述了采集硬件原理,Linux字符设备驱动模型,SSC和GPS对钟驱动实例.     

考虑海底地形的三维频率域可控源电磁响应有限体积法模拟

为了提高海洋可控源电磁的地形响应模拟精度,探讨了三维交错网格剖分有限体积法求解频率域电磁场解的数值方法.在似稳场近似条件下,推导了有限差分法和有限体积法低频电磁场控制方程的离散表达式.通过对离散表达式的分析比较表明,有限体积法较之有限差分具有更高的计算精度,同时也具有与有限差分相当的计算效率.对含地形界面网格单元电导率采用加权平均处理,通过与二维有限元程序的计算结果对比,验证了该方法可以有效提高有限体积法对地形变化的模拟精度.应用有限体积法计算了一个三维带地形储层模型的可控源电磁响应,分析表明地形变化对电场分量影响明显,磁场分量对地形变化不敏感.     

On the determination of marine magnetotelluric receiver orientation using land magnetic observatory data

Marine magnetotelluric measurements using “free‐fall’’ instruments without effective compasses suffer from the problem of unknown orientation of the receivers at the seafloor. While past works indicate that marine magnetotelluric orientation of the instruments can be estimated by reference to land deployments of known orientation using the transfer tensor method, there is limited published information on how this is implemented in practice. We document this method and propose a set of new time‐ and frequency‐domain approaches to solve this orientation problem of the seafloor receivers. We test these methodologies in onshore and offshore magnetotelluric data whose orientations are well known and apply these techniques to marine magnetotelluric data with unknown orientation. For the controlled tests, both time‐ and frequency‐domain approaches produce overall comparable results. To investigate the effects of the subsurface structure distribution on the orientation process, a dimensionality analysis of a controlled dataset is carried out. In subsequent analysis using the available disoriented marine magnetotelluric data from offshore Brazil and from the Vassouras magnetic observatory on the mainland for remote referencing, frequency‐domain methods yield approximate orientation angles among themselves with low standard deviation each. Time‐domain results are consistent for most cases but differ from frequency‐domain results for some situations.     

基于改进的接收点插值算法的频率域海洋可控源电磁法2.5维正演

在海洋可控源电磁法勘探中,接收站常置于海底.在进行海洋电磁场模拟时,由于海水和海底介质存在显著电性差异,这给海底接收点处场值的求取带来困难.本文提出一种新的接收点插值算法,该算法考虑到海底电场法向分量不连续性问题,用法向电流分量进行插值以准确求取海底任意接收点处电磁场值.本文利用交错网格有限差分法实现了二维介质中频率域海洋可控源法（CSEM）正演.对构造走向做傅里叶变换,将三维电磁模拟问题转换为波数域2.5维问题,即三维场源激励下针对二维地电模型的电磁模拟问题.使用交错网格有限差分法,基于一次场/二次场分离方法导出波数域二次电场离散形式,并进一步求得波数域电磁场.采用本文提出的改进的插值算法可求得海底任意接收点处波数域电磁场,采用傅里叶逆变换对波数域电磁场进行积分可得到接收点处空间域电磁场.模型算例表明,与常规的线性插值和严格插值算法相比,本文提出的改进的插值算法具有更高的精度.

     

A numerical comparison of time and frequency‐domain marine electromagnetic methods for hydrocarbon exploration in shallow water

In shallow water the frequency domain controlled source electromagnetic method is subject to airwave saturation that strongly limits the sensitivity to resistive hydrocarbon targets at depth. It has been suggested that time‐domain CSEM may offer an improved sensitivity and resolution of these deep targets in the presence of the airwave. In order to examine and test these claims, this work presents a side‐by‐side investigation of both methods with a main focus on practical considerations, and how these effect the resolution of a hydrocarbon reservoir. Synthetic noisy data for both time‐domain and frequency domain methods are simulated using a realistic frequency dependent noise model and frequency dependent scaling for representative source waveforms. The synthetic data studied here include the frequency domain response from a compact broadband waveform, the time‐domain step‐response from a low‐frequency square wave and the time‐domain impulse response obtained from pseudo‐random binary sequences. These data are used in a systematic resolution study of each method as a function of water‐depth, relative noise and stacking length. The results indicate that the broadband frequency domain data have the best resolution for a given stacking time, whereas the time‐domain data require prohibitively longer stacking times to achieve similar resolution.     

海底地震仪的信号传递效果和噪声水平

单球式海底地震仪（以下简称OBS）由于其成本低、操作简便的优点在天然地震研究、人工地震探测中获得了广泛应用.本文首先分析多型进口和国产OBS在台湾海峡西部采集的地震数据,发现同一台OBS上的垂直向速度检波器（Z分量）的信噪比常显著低于压力检波器（H分量）,由于这两种检波器记录的都是海底的垂向振动信号,推测速度检波器的低信噪比更多的与仪器特性有关.然后从信号传递和噪声水平两方面分析影响速度检波器信噪比的因素：为检测速度检波器与OBS壳体的耦合效果,对某型宽频带OBS和陆上地震仪进行了同址同步观测试验,发现OBS的整机灵敏度有较大的差异;为分析速度检波器的水底噪声特征,以H分量记录作为基准,对比分析了同一台仪器不同站位的Z分量噪声水平,发现速度检波器在浅海区受到较大的次生干扰.本文指出OBS的内部耦合和水流次生干扰是至今尚未引起大家重视而又严重影响资料品质和多波探测成效的两个关键问题,这一研究结果对于改进OBS结构设计和制造工艺,以及OBS数据多分量处理方法研究有重要的参考意义.     

西南印度洋中脊三维地震探测中炮点与海底地震仪的位置校正

在海上实施三维地震探测过程中,人工震源枪阵中心与船上GPS的距离及地震探测作业中的船行方向造成炮点实际位置与预设位置有一定偏差;自由落体投放的OBS由于海流的影响会偏离原定设计位置(投放点),因此,炮点与海底地震仪(OBS)的位置校正是三维地震结构研究中的基本环节.本文利用艏向信息校正了炮点位置;采用蒙特卡洛和最小二乘法方法对海底地震仪的位置进行了校正,并探讨了直达水波曲线特征.结果表明 OBS位置一般偏离设计点1 km左右,其误差范围在20 m以内,校正后的OBS记录剖面展示了真实的记录情况.该研究结果为下一步西南印度洋的三维层析成像研究提供了坚实数据基础,同时为今后南海的三维深部地壳结构探测提供经验与借鉴.     

The oxygen isotopic composition and temperature of Southern Ocean bottom waters during the last glacial maximum

We provide two new determinations of the oxygen isotopic composition of seawater during the last glacial maximum (LGM). High-resolution oxygen isotopic measurements were made on interstitial waters from Ocean Drilling Program (ODP) Sites 1168 and 1170 in the southeast Indian Ocean sector of the Southern Ocean. We use a diffusion-advection numerical model to calculate the glacial-interglacial change in bottom-water δ¹⁸O_sw from the pore water δ¹⁸O profiles; the first such determinations from this part of the oceans. Statistical analyses of the model runs indicate that Circumpolar Deep Water (CDW) δ¹⁸O_sw changed by 1.0-1.1±0.15‰ since the last glacial maximum (LGM). Our results are consistent with a previous calculation from a South Atlantic Southern Ocean location (ODP Site 1093) also situated within CDW. The new values determined in this study, together with previous estimates, are converging on a global average Δδ¹⁸O_sw of 1.0-1.1‰.Using the calculated bottom-water δ¹⁸O_sw, we have extracted the temperature component from the benthic foraminiferal δ¹⁸O record at Sites 1168 and 1170. Since the LGM, bottom waters at these two sites warmed by 2.6 and 1.9°C, respectively. The absolute temperature estimates for the LGM (−0.5°C [Θ=−0.6°C] at Site 1168 and −0.2°C [Θ=−0.4°C] at Site 1170) are slightly warmer than those reported from previous studies using the same technique, but are consistent with more homogenous deep-ocean temperatures during the LGM relative to the modern.     

The size, mass, and composition of plastic debris in the western North Atlantic Ocean

This study reports the first inventory of physical properties of individual plastic debris in the North Atlantic. We analyzed 748 samples for size, mass, and material composition collected from surface net tows on 11 expeditions from Cape Cod, Massachusetts to the Caribbean Sea between 1991 and 2007. Particles were mostly fragments less than 10 mm in size with nearly all lighter than 0.05 g. Material densities ranged from 0.808 to 1.24 g ml⁻¹, with about half between 0.97 and 1.04 g ml⁻¹, a range not typically found in virgin plastics. Elemental analysis suggests that samples in this density range are consistent with polypropylene and polyethylene whose densities have increased, likely due to biofouling. Pelagic densities varied considerably from that of beach plastic debris, suggesting that plastic particles are modified during their residence at sea. These analyses provide clues in understanding particle fate and potential debris sources, and address ecological implications of pelagic plastic debris.     

小孔径台阵垂直分量P波接收函数研究及应用

尝试用垂直分量P波接收函数技术,分离莫霍面多次反射透射震相PpPmp,以此确定台阵下方的地壳厚度．利用兰州小孔径台阵记录的16次深远地震资料,分离了PpPmp震相,根据其到时计算台阵下方的地壳厚度约(51.6plusmn;1.8) km,与前人研究结果一致．      

Habitat associations of floating debris and marine birds in the North East Pacific Ocean at coarse and meso spatial scales

While many surface foraging seabirds ingest plastic, the spatial overlap of these far-ranging predators with debris aggregations at-sea is poorly understood. We surveyed concurrent distributions of marine birds and debris along a 4400 km cruise track within a debris accumulation area in the North East Pacific Ocean using line and strip transect methods. Analysis of debris and bird distributions revealed associations with oceanographic and weather variables at two spatial scales: daily surveys and hourly transects. Hourly bird abundance (densities; 0-9 birds km⁻²) was higher in lower wind and shallower water. Hourly debris abundance (densities; 0-15,222 pieces km⁻²) was higher in lower wind, higher sea-level atmospheric pressure and deeper water. These results suggest that debris and seabird abundance and community structure are influenced by similar environmental processes, but in opposing ways, with only three far-ranging seabird species (Black-footed Albatross, Cook’s Petrel and Red-tailed Tropicbird) overlapping with high debris concentrations over meso-scales.