Constraints on Characteristics and Distribution of Gas Hydrate and Free Gas Using Broad-Band Processing of Three-Dimensional Seismic Data

Gas hydrate drilling expeditions in the Pearl River Mouth Basin,South China Sea,have identified concentrated gas hydrates with variable thickness.Moreover,free gas and the coexistence of gas hydrate and free gas have been confirmed by logging,coring,and production tests in the foraminifera-rich silty sediments with complex bottom-simulating reflectors(BSRs).The broad-band processing is conducted on conventional three-dimensional(3D)seismic data to improve the image and detection accuracy of gas hydratebearing layers and delineate the saturation and thickness of gas hydrate-and free gas-bearing sediments.Several geophysical attributes extracted along the base of the gas hydrate stability zone are used to demonstrate the variable distribution and the controlling factors for the differential enrichment of gas hydrate.The inverted gas hydrate saturation at the production zone is over 40% with a thickness of 90 m,showing the interbedded distribution with different boundaries between gas hydrate-and free gas-bearing layers.However,the gas hydrate saturation value at the adjacent canyon is 70%,with 30-m-thick patches and linear features.The lithological and fault controls on gas hydrate and free gas distributions are demonstrated by tracing each gas hydrate-bearing layer.Moreover,the BSR depths based on broad-band reprocessed 3D seismic data not only exhibit variations due to small-scale topographic changes caused by seafloor sedimentation and erosion but also show the upward shift of BSR and the blocky distribution of the coexistence of gas hydrate and free gas in the Pearl River Mouth Basin.

国家地震烈度速报江苏子系统仪器烈度产出分析

对国家地震烈度速报与预警工程部署的烈度速报江苏子系统进行了介绍,简述了2021年12月系统部署以来的系统产出情况,同时以系统部署以来产出的震级最大的事件即2021年12月22日江苏天宁M4.2地震为例,使用f-k方法合成各观测台站上的理论地震图,将实际峰值加速度与理论峰值加速度进行比较,评估烈度速报江苏子系统产出结果的可靠性。结果表明,烈度速报系统产出的各台站峰值加速度随震中距逐渐减小,且和理论地震图的结果较为接近,一定程度上证明了烈度速报江苏子系统产出结果的可靠性。

A broad-band sub-harmonic mixer for W-band applications

In order to improve the bandwidth of the conventional sub-harmonic mixer, a broad-band, high intermediate frequency（IF） sub-harmonic mixer for W-band applications is proposed. Replacing the open and short stubs that are used in the convertional sub-harmonic mixer with a broad-band band-pass filter and a low-pass filter, respectively, a wide operating frequency band is achieved. Furthermore, without the use of the edge-coupled band-pass filter at radio frequency（RF） port, the proposed structure can be realized by common hybrid microwave integrated circuit technology at W- band. The measured results show that the proposed subharmonic mixer can operate from 80 to 107.5 GHz for RF frequency and support up to 18 GHz for the IF bandwidth. Also, the measured results show that the single-sideband conversion loss is less than 13. 7 dB over the available RF frequency band, while the minimum conversion loss is about 9 dB at an RF of 92. 5 GHz and an 1F of 3 GHz. Thus, a large operating bandwidth performance at W-band can be achieved by the orooosed mixer.

基于近场记录的强震台站场地波速结构反演

基于近场加速度记录反演场地浅部速度结构,可为补充强震台站场地信息提供一种简便的途径,有利于记录的充分利用。本文详细推导了基于单一台站记录的P波地震图法,分析了其中的反演目标函数;通过数值模拟大量近场场地上的理论地震图,分析了场地覆盖层剪切波速、厚度等特征参数及震源持时对PS波与P波延时的影响,给出了地震记录的筛选原则;通过增加P波初到时刻径向分量与竖向分量幅值比函数作为反演目标函数,对单台地震图法进行改进;最后,利用汶川地震主震及余震记录,反演了武都强震台的场地波速结构,并与表面波法和钻孔法进行对比,结果表明,改进的单台地震图法搜索能力得到了增强,反演效率和精度得到显著提高。

Numerical Calculation of Seafloor Synthetic Seismograms Caused by Low Frequency Point Sound Source

Elastic wave on seafloor caused by low frequency noise radiated from ship is called ship seismic wave which can be used to identify ship target. In order to analyze the wave components and the propagating properties of ship seismic wave, the numerical calculation of synthetic seismograms on seafloor aroused by a low frequency point sound source is carried out using a wave number integration technique combined with inverse Fourier transform. According to the numerical example of hard seafloor, the time series of seismic wave on seafloor are mostly composed of interface waves and normal mode waves. Each normal mode wave has a well defined low cut-off frequency, while the interface wave doesn't have. The frequency dispersion of normal mode wave is obvious when frequency is lower than 100Hz, while the interface wave is dispersive only in the infra-sound frequency range. The time series of seismic wave is dominated by the interface wave when the source frequency is less than the minimal cut-off frequency of normal mode wave.

Seismogram Synthesis in Multi-layered Half-space Part Ⅰ. Theoretical Formulations

In the past two decades numerous studies were made to develop and improve the theory and practical computation techniques of synthesizing theoretical seismograms for the model of multi-layered half-space. Today, synthesizing theoretical seismograms in multi-layered half-space is an important tool for understanding the structure of the Earth’s interior as well as the seismic source process from well-recorded seismograms data. As part of a review of the state-of-the-art, in this article I shall present a systematic and self-contained theory of elastic waves in multi-layered half-space media based on the developments in the past two decades.

A novel technique for automatic seismic data processing using both integral and local feature of seismograms

A novel technique for automatic seismic data processing using both integral and local feature of seismograms was presented in this paper. Here, the term integral feature of seismograms refers to feature which may depict the shape of the whole seismograms. However, unlike some previous efforts which completely abandon the DIAL approach, i.e., signal detection, phase identifi- cation, association, and event localization, and seek to use envelope cross-correlation to detect seismic events directly, our technique keeps following the DIAL approach, but in addition to detect signals corresponding to individual seismic phases, it also detects continuous wave-trains and explores their feature for phase-type identification and signal association. More concrete ideas about how to define wave-trains and combine them with various detections, as well as how to measure and utilize their feature in the seismic data processing were expatiated in the paper. This approach has been applied to the routine data processing by us for years, and test results for a 16 days＇ period using data from the Xinjiang seismic station network were presented. The automatic processing results have fairly low false and missed event rate simultaneously, showing that the new technique has good application prospects for improvement of the automatic seismic data processing.

A cloud-based synthetic seismogram generator implemented using Windows Azure

Synthetic seismograms generated by solving the seismic wave equation using numerical methods are being widely used in seismology. For fully three-dimensional seismic structure models, the generation of these synthetic seismograms may require large amount of computing resources. Conventional high-performance computer clusters may not provide a cost-effective solution to this type of applications. The newly emerging cloud-computing platform provides an alternative solution. In this paper, we describe our implementation of a synthetic seismogram generator based on the reciprocity principle using the Windows Azure cloud application framework. Our preliminary experiment shows that our cloud-based synthetic seismogram generator provides a costeffective and numerically efficient approach for computing synthetic seismograms based on the reciprocity principle.

The strain seismograms of P- and S-waves of a local event recorded by four-gauge borehole strainmeter

At a sampling rate of 100 samples per second,the YRY-4 four-gauge borehole strainmeters（FGBS） are capable of recording transient strains caused by seismic waves such as P and S waves or strain seismograms. At such a high sampling rate, data from the YRY-4 strainmeters demonstrate fairly satisfactory self-consistency. The strain tensor seismograms demonstrate the senses of motion of P waves, that is, the type of seismic wave travels in the direction of the maximum normal strain change. The observed strain patterns of S waves significantly differ from those of P waves and should contain information about the source mechanism. Spectrum analysis shows that the strain seismograms are consistent with conventional broadband seismograms from the same site.

Research on anisotropic parameters by syn-thetic seismogram

Based on the extensive-dilatancy anisotropy theory, the method of synthetic seismogram is used to estimate the anisotropic parameters. The advantages of the method lie in that it avoids the singularity resolution and saves calculation time of computer by using the eigenvalue and eigenvector analytical expressions of Christoffel equation, at the same time, the result is tested by coherence function. The test result reveals that there exists a fine linear relation between original records and synthetic records, indicating the anisotropic parameters estimated by synthetic seismogram can reflect and describe the anisotropic characteristics of the given region medium.

The optimum source－time function for generating finite－difference synthetic seismograms

A source-time function plays an important role in generating finite-difference synthetic seismograms.In this regard,we consider the amplitude spectra of five different time functions.This investigation shows that when the high-frequency signals of the time function carry more power than the main signal,the tailing phenomena would be associated with the main signal which is not desirable in producing finite-difference synthetic seismograms.It is illustrated that the time function denoted by is the optimal source-time function among those examined in this paper for generating finite-difference synthetic seismograms.This matter is also verified by making finite-difference synthetic seismograms over an elastic earth model consisting of a layer over a half space.

Grid dispersion in generating finite-differences synthetic seismograms

In order to obtain a finite-difference synthetic seismogram, the finite earth model must be subdivided into cells.If the cell size becomes too large in comparison with wavelength of the source signal, waves disperse with increasing traveltime. This phenomenon is known as grid dispersion. The grid dispersion produces a variation of velocity with different frequencies; that is the higher signal frequencies travel more slowly than the lower signal frequencies. Consequently, substantial tailing of the signal arises with increasing traveltime. This may happpen if: (1) the grid spacing is large; (2) the sampling rate is large; or (3) the source wavelength is too short compared with grid size. In other words, an important parameter in generating finite-difference synthetic seismograms is the number of grid points per wavelength of the source signal. In this paper, it is shown that the frequency of a source function has great effects on grid dispersion when P-Sv coupled waves propagate through elastic models by finite-differences. The two-dimensional elastic models considered in this paper consist of : (1) a normal fault, and (2) a layer over a half-space. This study verifies that, when generating finite-difference synthetic seismograms, the grid dispersion will be reduced to a satisfactory level if the grid points per wavelength at the half-power frequency of the source signal for the lowest velocity of the medium exceeds ten.

A new method of weighted superposition for seismogram synthesis in depth domain

The extensive use of depth-imaged seismic data produced by pre-stack depth migration(PSDM)leads to the necessity to synthesize seismogram directly in depth domain.However,since seismic wavelet in depth domain is dependent on media velocities.The time-domain convolution operation directly used in depth domain does not meet the linear time-invariant condition.In this paper,we present a new method for genuine depth-domain seismic synthesis.This method constructs the velocity-dependent seismic wavelets varying adaptively with the corresponding interval velocities in the depth direction and weights them by the reflectivities,then the synthetic seismic record is obtained by the superposition of these weighted depth-varying wavelets.We applied this method to synthesize the seismic record of both a multi-layered geological model and the field data.The results show that the method can accommodate the intrinsic velocity-dependent variation characteristics of seismic events in depth domain and avoid the redundant depth-to-time and time-to-depth transformations.

Broad-Band All-Optical Wavelength Conversion of Differential Phase-Shift Keyed Signal Using an SOA-Based Nonlinear Polarization Switch

Broad-band all-optical wavelength conversion of differential phase-shift keyed （DPSK） signal is experimentally demonstrated. This scheme is composed of a one-bit delay interferometer demodulation stage followed by a semiconductor optical amplifier （SOA） based nonlinear polarization switch. A wavelength converter for the 10 G b/s DPSK signal is presented, which has a wide wavelength range of more than 30 nm. The converted signals experience small power penalties less than 1.4 dB compared with the original signal, at a bit error rate of 10-9. Additionally, the optical spectra, the measured waveforms and the open eye diagrams of the converted signals show a high quality wavelength conversion performance.

Broad-Band FMR Linewidth of Co_2MnSi Thin Films with Low Damping Factor:The Role of Two-Magnon Scattering

The low Gilbert damping factor, which is usually measured by ferromagnetic resonance, is crucial in spintronic applications. Two-magnon scattering occurs when the orthogonMity of the ferromagnetic resonance mode and other degenerate spin wave modes was broken by magnetic anisotropy, voids, second phase, surface defects, etc., which is important in analysis of ferromagnetic resonance linewidth. Direct fitting to linewidth with Gilbert damping is advisable only when the measured linewidth is a linear function of measuring frequency in a broad band measurement. We observe the nonlinear ferromagnetic resonance linewidth of Co2MnSi thin films with respect to measuring frequency in broad band measurement. Experimental data could be well fitted with the model including two-magnon scattering with no fixed parameters. The fitting results show that two-magnon scattering results in the nonlinear linewidth behavior, and the Gilbert damping factor is much smaller than reported ones, indicating that our Co2MnSi films are more suitable for the applications of spin transfer torque.

Seismogram Analysis of Indonesian Earthquakes at DAV Observation Station

The S-wave velocity across the earth structure under Indonesia for Indonesia earthquakes has been investigated through seismogram analysis, simultaneously in the time domain and three Cartesian components. The data were recorded at DAV observational station, the Philippines. The main data set is the seismogram comparison between the measured and synthetic seismogram, instead of travel time data, as commonly used in other seismological research. The synthetic seismogram is calculated using the GEMINI method, which is equivalent to Mode Summation. The above seismogram comparison shows that the global earth mantle of PREMAN gives a deviating synthetic seismogram and has earlier arrival times than those of the measurement. The gradient of βh in the upper mantle layers is altered into a positive, rather than negative slope as stated in the PREMAN model, and negative corrections are imposed to the zero order of the polynomiars coefficients in all earth mantle layers. The excellent fitting, as well as travel time or waveform, is obtained from the surface waves of Love and Rayleigh, surface wave to the S and SS mantle waves as well as the core reflected waves. This result expresses that part of the earth mantle, due to a collision between India and Asia tectonic released zones, has a negative anomaly in S-wave velocity and vertical anisotropy in all of the earth mantle layers.

Seismogram Analysis of Earthquake C060394F,South Java: S-wave Velocity Structure

The S wave velocity structure between the hypocenter of C060394F earthquake,South Java and a series of observatory stations located in Australia and South-East Asia have been investigated through seismogram analysis in the time domain and the three Cartesian components.The synthetic seismogram is constructed from the PREMAN global earth model.Seismogram comparison between the measured and synthetic seismograms shows large discrepancies.A correction to the S wave velocity structure is needed to solve these discrepancies.Seismogram analysis and fitting at the observatory stations located in Australia show a positive anomaly,whereas the stations in South-East Asia show a strong negative anomaly.Vertical anisotropy occurred in the upper mantle layer and earth mantle layers beneath it.The seismogram analysis in the time domain and the three components simultaneously provides a better picture of the earth model,compared to the analysis of the seismogram with arrival time,travel time difference and dispersion analysis.

Seismogram Analysis and Fitting of Three Earthquakes,West Papua,Recorded at Stations ENH and KMI,China

The S wave velocity structure of the earth below Eastern Southeast Asia has been investigated by analyzing the seismogram from surface wave to multiple depth waves in the time domain and three Cartesian components simultaneously.The wave passes across the front area of subduction zone between the Philippine plate and the Asian plate.The main data are waveform comparisons,instead of the arrival times.The synthetic seismogram is calculated using the GEMINI method.The synthetic seismogram constructed by PREMAN global earth model deviates greatly from the measured one.To solve this problem,corrections are needed for the β speed structure.Corrections cover the gradient change of β_h,which turns from negative to positive in upper mantle layers as in the PREMAN,change of earth crust depth and change of zero order coefficients of β velocity function in all earth mantle layers.So,the fitting is obtained,as well as the arrival time or the waveform of Love and Rayleigh surface waves,the S wave and the repetitive depth waves ScS_2 and ScS_3.This result reveals that the Southeast Asia,being stretched due to tectonic release,has a mantle in some parts with negative anomaly of S wave velocity and vertical anisotropy in all earth mantle layers.

基于改进的介观裂缝岩石物理模型的频散AVO特征分析

地下储层中存在介观尺度的裂缝是导致地震波的速度频散和能量衰减的一个重要原因。从Galvin介观裂缝模型出发,重构其低频极限和高频极限的表达形式,将频散关系施加到裂缝模型上构建依赖频率的裂缝柔度参数,建立了一种改进的介观尺度的裂缝型岩石物理模型。基于该改进模型,利用附加柔度的组合性质,可以容易地构建更复杂的模型并进行依赖频率的地震响应特征分析。在传统的反射系数公式中引入频率参数,分析了具有不同裂缝长度、背景渗透率、流体粘度等频散敏感参数的反射界面上的频散AVO及地震记录响应特征。通过数值模拟得出频散敏感参数主要影响发生频散的频段,即影响岩石的特征频率;当模型的特征频率与地震波频段主频相近时,频散现象最突出,当特征频率大于地震波频段主频一个数量级时,可以忽略频散作用,此时采用传统的Gassmann方程即可较准确地分析地震剖面的特征。基于改进模型得到的认识为地震频散特征的实际应用提供了岩石物理模型的构建方面的指导。

河北历史地震模拟波形图纸数字化方法及初步研究

利用HSD软件对河北台网部分模拟图纸进行波形数字化,提取到9个地震27条SAC格式波形文件,利用软件对数字化波形进行震相到时、振幅、震级偏差等多方面校验。结果表明,所有可拾取的震相到时总体平均绝对偏差为0.31 s;振幅总体平均绝对偏差为0.04μm,与原振幅值基本相当;震级平均绝对偏差为0.08,整体偏差较小,精度较高。最终校验结果表明,数字化后的波形可较好地应用于数字化软件,为后续波形数字化以及历史数据的数字化应用提供较好的研究基础。