Characteristics of the Seismic Waves from a New Active Source Based on Methane Gaseous Detonation

Active seismic sources are critical for obtaining high resolution images of the subsurface. For active imaging in urban areas, environment friendly and green seismic sources are required. In present work, we introduce a new type of green active source based on the gaseous detonation of methane and oxygen. When fired in a closed container, the chemical reaction, i.e. gaseous detonation, will produce high pressure air over 150MPa. Seismic waves are produced when high pressure air is quickly released to impact the surroundings. The first field experiment of this active source was carried out in December, 2017 in Jingdezhen, Jiangxi Province, where a series of active sources were excited to explore their potential in mine exploration. In current work, we analyzed the seismic waves recorded by near-field accelerators and a dense short-period seismic array and compared them with those from a mobile airgun source, another kind of active source by releasing high pressure air into water. The results demonstrate that it can be used for high resolution near surface imaging. Firstly, the gaseous detonation productions are harmless CO₂ and water, making it a green explosive source. Secondly, the dominant seismic frequencies are 10-80Hz and a single shot can be recorded up to 15km, making it suitable for local structure investigations. Thirdly, it can be excited in vertical wells, similar to traditional powder explosive sources. It can also act as an additional on-land active source to airgun sources, which requires a suitable water body as intermediate media to generate repeating signals. Moreover, the short duration and high frequency signature of the source signals make it safe with no damage to nearby buildings. These make it convenient to excite in urban areas. As a new explosive source, the excitation equipment and conditions, such as gas ratio, sink depth and air-releasing directions, need further investigation to improve seismic wave generation efficiency.     

The June 2013 Alberta catastrophic flooding event – part 2: fine‐scale precipitation and associated features

Data obtained from a variety of sources including the Canadian Lightning Detection Network, weather radars, weather stations and operational numerical weather model analyses were used to address the evolution of precipitation during the June 2013 southern Alberta flood. The event was linked to a mid‐level closed low pressure system to the west of the region and a surface low pressure region initially to its south. This configuration brought warm, moist unstable air into the region that led to dramatic, organized convection with an abundance of lightning and some hail. Such conditions occurred in the southern parts of the region whereas the northern parts were devoid of lightning. Initially, precipitation rates were high (extreme 15‐min rainfall rates up to 102 mm h^?1 were measured) but decreased to lower values as the precipitation shifted to long‐lived stratiform conditions. Both the convective and stratiform precipitation components were affected by the topography. Similar flooding events, such as June 2002, have occurred over this region although the 2002 event was colder and precipitation was not associated with substantial convection over southwest Alberta. Copyright © 2016 Her Majesty the Queen in Right of Canada. Hydrological Processes. © John Wiley & Sons, Ltd.     

广西龙滩库区深部孕震结构大地电磁探测研究

在广西龙滩库区布置2条北北西-南东东向大地电磁探测剖面进行了32个测点的观测,观测频段为320 Hz~1000 s.对2条剖面上测量和计算得到的各种电磁参数、深部电性结构图像以及库区地质构造、岩性和水域分布等进行分析研究.结果显示库区发育的6条断裂在电性上表现为不同程度的电性差异边界,其中拉浪-达良断裂和拉色-辉马断裂为库区主要断裂带.库区深部电性结构总体表现为纵向上自地表到20 km深度呈现高-低-高的波浪起伏的三层结构,横向上呈明显的高低电阻块状组合样式.龙滩库区4个地震丛集区的震源都位于高、低电阻接触带附近,其中3个地震丛集区的震源位于低阻特性的石炭系地层的下部,为岩溶水体诱发地震;发生最大震级的地震丛集区的震源位于上宽下窄似"铆钉状"的高阻体下部,推测是因水库蓄水后水体压力增大和库水渗透作用下在聚集高变形能的脆性高电阻体内部发生的地震.     

气体水合物密闭采样原位饱和度评估技术研究

海底天然气水合物采集样品的饱和度确定通常是通过密闭保压采样技术来实现,但采样器的温度压力控制非常困难.本文用热力学模型来模拟密闭体系下温度压力变化过程中甲烷在各相态中饱和度的变化过程,并在海洋水合物热力学系统基础上,完成了一项分析水合物原位饱和度的技术,其只要求采集器密闭,不需要进行严格的温度和压力控制,降低了对采样器的要求.实例分析显示,完全可以利用密闭样品在地表下气体和液体成分的测量,包括盐度的测量来估计出甲烷水合物原位饱和度.     

多层弹性半空间中的地震波（一）

为了了解地震震源和地球介质的性质,很有必要对地震波的辐射、传播和衰减问题作仔细的分析。作为一种近似,可以暂且忽略地球的曲率,把传播地震波的地球介质视为多层半空间。为简便起见,地震波的衰减问题另作考虑。这样,便需要研究多层、均匀、各向同性和完全弹性半空间中地震震源辐射的地震波传播问题。 用哈斯克尔（Haskell）矩阵法解多层介质中弹性波的传播问题是很方便的。如果     

基于累积耗损能量的饱和粉土液化特性试验研究

饱和粉土场地在强地震作用下易发生液化现象。开展饱和粉土的循环三轴试验,以循环加载的累积耗损能量为指标,探讨黏粒含量、密实度、有效围压和循环应力比等因素对粉土液化特性的影响,试验结果表明:粉土液化所需的耗损能量随黏粒含量的增加呈先减小后增大的趋势,当黏粒含量约为8%时粉土的液化耗损能量最低;液化耗损能量随粉土密实程度的增大而逐渐增加,并随初始有效围压的增长而增加,但粉土的液化耗损能量与循环应力比之间的关系不明显。     

Acoustic directional snapshot wavefield decomposition

Up–down wavefield decomposition is effectuated by a scaled addition or subtraction of the pressure and vertical particle velocity, generally on horizontal or vertical surfaces, and works well for data given on such surfaces. The method, however, is not applicable to decomposing a wavefield when it is given at one instance in time, i.e. on snapshots. Such situations occur when a wavefield is modelled with methods like finite-difference techniques, for the purpose of, for example, reverse time migration, where the entire wavefield is determined per time instance. We present an alternative decomposition method that is exact when working on snapshots of an acoustic wavefield in a homogeneous medium, but can easily be approximated to heterogeneous media, and allows the wavefield to be decomposed in arbitrary directions. Such a directional snapshot wavefield decomposition is achieved by recasting the acoustic system in terms of the time derivative of the pressure and the vertical particle velocity, as opposed to the vertical derivative in up–down decomposition for data given on a horizontal surface. As in up–down decomposition of data given at a horizontal surface, the system can be eigenvalue decomposed and the inverse of the eigenvector matrix decomposes the wavefield snapshot into fields of opposite directions, including up–down decomposition. As the vertical particle velocity can be rotated at will, this allows for decomposition of the wavefield into any spatial direction; even spatially varying directions are possible. We show the power and effectiveness of the method by synthetic examples and models of increasing complexity.     

Quantitative sediment source ascription using magnetic measurements in a reservoir-catchment system near Nijar,S.E. Spain

The magnetic properties of sediments from a dried-out reservoir and its inflowing streams were measured to provide a basis for sediment source determination. Quantitative estimates of source component proportions in the reservoir sediments were made using linear programming procedures. The study confirms that this type of research must be based on a detailed knowledge of the relationship between magnetic properties and particle size in both the sediments and the potential sources. In the case of the site studied, the Isabel II reservoir, near Nijar, S.E. Spain, the sediments generated by contemporary surface processes within the reservoir catchment fail to reflect fully the 19th century material largely responsible for the rapid infill of the reservoir. An additional source type, namely unweathered schist-derived material not represented in the fine sediments sampled within channels upstream from the reservoir, has to be included in the model before calculations of source component contributions can be made. It is inferred that one or more high magnitude, low frequency event outside the amplitude of the most recently operating processes was largely responsible for the infilling of the reservoir.     

北京及周边气溶胶区域影响与大雾相关特征的研究进展

针对2011年12月初北京及华北持续近一周的严重大雾天气这一热点事件,从城市群大雾过程气溶胶区域影响的视角,基于"973"项目"北京及周边地区大气-水-土环境污染机理与调控原理"的研究工作,就北京及周边地区大雾天气与大气气溶胶区域影响的关系等方面进行了讨论.研究表明,北京城市大雾前低空SO₂和NO₂浓度存在"积聚"与"突增"现象.北京及周边地区冬季雾日数和气溶胶光学厚度则呈正相关,并具有"同位相"的年际变化趋势.研究同时发现北京及其南部周边的冬季气溶胶高值区呈南北向带状分布,其与北京周边居民户数高值区有所吻合,反映了冬季北京城市气溶胶颗粒物的远距离影响源区及大尺度输送效应.统计分析指出,冬季北京气溶胶颗粒物PM₁₀、PM_2.5主要影响成分是SO₂和NO_X,且有关研究也表明,电厂、采暖和工业面源是SO₂的三大本地排放源,而机动车、电厂、工业为NO_X的三大本地排放源,上述大气PM₁₀、PM_2.5主成分污染源亦与雾水样本化学分析结果相吻合,即冬季由于燃煤在生活能源中的比例较大,北京雾水中硫元素和碳元素的含量都较高.因此,北京冬季大雾不仅与北京城区气溶胶及其污染排放影响存在相关关系,而且与北京周边天津、河北、山东等地气溶胶及大气污染物的远距离输送和气溶胶区域影响效应有着重要的联系.因此,北京雾霾天气及相关大气污染的治理工作首先要着眼于局地污染物的减排,但同时如何做好区域大气污染的协同治理也是不容忽视的问题.     

Coupled infiltration and filtration behaviours of concrete porous pavement for stormwater management

Cementitious porous pavement (CPP) is a structural low‐impact development material for rainfall–runoff management. Both infiltration and filtration are critical functions for CPP stormwater quality and quantity control. In this study, three groups of CPP specimens exposed to rainfall–runoff for 4 years and experienced with different maintenance intervals (6, 12 and 48 months, respectively) were used to examine CPP infiltration and filtration performance. Particle mass strained on CPP surface, saturated infiltration rate I_f, temporal infiltration rate I(t), suspended sediment concentration (SSC) and turbidity (τ) were measured to evaluate the process of filtration/infiltration. I(t), SSC and τ were examined less than 50 mg/l of the suspended particle loading. It was found that the CPP surface cleaning methods used in the past 4 years, namely, high pressure wash followed by vacuuming with one atmosphere (100 kPa), were effective, and a 12‐month maintenance interval was verified suitable to maintain the pore structure an acceptable infiltration rate for stormwater management. It was also found that CPP infiltration and filtration process affect each other, and the two properties are coupled in urban stormwater quality and quantity control. On the basis of the experimental measurements, the temporal infiltration rate of the cleaned CPP under a certain particle loading could be simulated by a first‐order nonlinear rational model, and effluent turbidity–SSC relationship was found following a power law. Copyright © 2012 John Wiley & Sons, Ltd.     

The Earth surface deformations caused by air pressure variations

The relationship between the tilts and strains of rocks and air pressure variations is investigated using deformation measurements from the Protvino observatory. Deformations due to air pressure are large in the observation area. The deformation response to the observed air pressure variations is compared with calculated values. In addition, the effective elasticity of the rock is determined from measured deformations in the underground galleries of the observatory after creating a load on the Earth surface. The observed tilts are considerably larger than the calculated values, for local as well as atmospheric loading, and the rock shear modulus determined from the tilt data is anomalously small compared to the mean shear modulus value for the crust. The observed horizontal strain data indicate no anomaly in the rock shear modulus value. The cavity effects at the observatory cannot explain this phenomenon. It may be caused by the rocks high fracturing. It is supposed, that the influence of cracks on the rock horizontal deformations is small if the cracks are closed and have the vertical direction. In case of large wavelength the large tilts correspond to the great displacements. These displacements may worsen the accuracy of geodetic measurements, especially when measuring vertical crustal movements.     

The June 2013 Alberta Catastrophic Flooding Event: Part 1—Climatological aspects and hydrometeorological features

In June 2013, excessive rainfall associated with an intense weather system triggered severe flooding in southern Alberta, which became the costliest natural disaster in Canadian history. This article provides an overview of the climatological aspects and large‐scale hydrometeorological features associated with the flooding event based upon information from a variety of sources, including satellite data, upper air soundings, surface observations and operational model analyses. The results show that multiple factors combined to create this unusually severe event. The event was characterized by a slow‐moving upper level low pressure system west of Alberta, blocked by an upper level ridge, while an associated well‐organized surface low pressure system kept southern Alberta, especially the eastern slopes of the Rocky Mountains, in continuous precipitation for up to two days. Results from air parcel trajectory analysis show that a significant amount of the moisture originated from the central Great Plains, transported into Alberta by a southeasterly low level jet. The event was first dominated by significant thunderstorm activity, and then evolved into continuous precipitation supported by the synoptic‐scale low pressure system. Both the thunderstorm activity and upslope winds associated with the low pressure system produced large rainfall amounts. A comparison with previous similar events occurring in the same region suggests that the synoptic‐scale features associated with the 2013 rainfall event were not particularly intense; however, its storm environment was the most convectively unstable. The system also exhibited a relatively high freezing level, which resulted in rain, rather than snow, mainly falling over the still snow‐covered mountainous areas. Melting associated with this rain‐on‐snow scenario likely contributed to downstream flooding. Furthermore, above‐normal snowfall in the preceding spring helped to maintain snow in the high‐elevation areas, which facilitated the rain‐on‐snow event. Copyright © 2016 John Wiley & Sons, Ltd.     

Spatial structure and scale feature of the atmospheric pollution source impact of city agglomeration

The spatial structure and multi-scale feature of the atmospheric pollution influence domain of Beijing and its peripheral areas (a rapidly developed city agglomeration) is dissected and analyzed in this paper on the basis of the atmospheric pollution dynamic-chemical process observation data of the urban building ensemble boundary layer of the Beijing City Air Pollution Observation Experiment (BECAPEX) in winter (February) and summer (August) 2003, and relevant meteorological elements and satellite retrieval aerosol optical depth (AOD), etc. comprehensive data with the dynamic-statistical integrated analysis of "point-surface" spatial structure. Results show that there existed significant difference in the contribution of winter/summer different pollution emission sources to the component character of atmospheric pollution, and the principal component analysis (PCA) results of statistical model also indicate that SO2 and NOX dominated in the component structure of winter aerosol particle; instead, CO and NOX dominated in summer. Surface layer atmospheric dynamic and thermal structures and various pollutant species at the upper boundary of building ensembles at urban different observational sites of Beijing in winter and summer showed an "in-phase" variation and its spatial scale feature of "influence domain". The power spectrum analysis (PSA) shows that the period spectrum of winter/summer particle concentration accorded with those of atmospheric wind field: the longer period was dominative in winter, but the shorter period in summer, revealing the impact of the seasonal scale feature of winter/summer atmospheric general circulation on the period of atmospheric pollution variations. It is found that from analyzing urban area thermal heterogeneity that the multiscale effect of Beijing region urban heat island (UHI) was associated with the heterogeneous expansion of tall buildings area. In urban atmospheric dynamical and thermal characteristic spatial structures, the turbulent scale feature of the urban boundary layer (UBL) of architectural complexes had important impact on the multi-scale feature of urban atmospheric pollution. The comprehensive analyses of the variational analysis field of Moderate Resolution Imaging Spectroradiometer (MODIS) AOD-surface PM10 under the condition of clear sky and the correlation resultant wind vector field for pollution source-tracing suggest that the emission sources for winter Beijing atmospheric pollution aerosols particle might be remotely traced to the south peripheral greater-scale spatial range of Hebei, Shandong, Tianjin, etc., and the spatial distribution of the high value area of AOD was associated with that of the high value area of resident family number (heating surface source). The backward trajectory feature of winter/ summer air particles exhibits analogous multi-scale feature, and depicts the difference in the scale feature of the pollution sources spatial distribution in different seasons. The peripheral source trajectory paths of urban atmospheric pollution (UAP) mainly come from the fixed industrial surface source or heating surface source in the outskirt of Beijing, and the diffusion and transport distance of peripheral sources in winter is larger than one in summer. The above conclusions depict the multi-scale spatial influence domain and seasonal features caused by UAP source influence and atmospheric dynamical structure. The high value area of the winter Total Ozone Mapping Spectrometer (TOMS) AOD lay in the Beijing region and its south peripheral area, an S-N zonal pattern, which reflects the dynamical effect of peripheral topographic pattern on the diffusion of regional scale atmospheric pollution sources. Study suggests that the extent of winter atmospheric pollution within the "valley" megarelief in Beijing and periphery was close related with the pollution emission sources of the south peripheral area; and the significant "anti-phase" variation feature of winter AOD and sunshine duration in Beijing and its peripheral areas, and regional scale correlation of low cloud cover, fog days, and aerosols reflects the local climatic effect of aerosol influence in this region. Besides, analysis of the impacts of atmospheric dry/wet deposition distributions within a valley-scale on the regional water body of Miyun reservoir also reveals the possible influence of the multi-scale spatial structure of summer water, soil and atmospheric pollution sources on the water quality of Miyun reservoir.     

Research Note: Internal/external seismic source wavefield separation and cancellation

A compact representation is obtained for the separation of a scalar wavefield on a closed surface into parts due to internal and external sources. The formula assumes that the total field and its gradients are known on the surface, as is the exact Green function of the medium. The derivation involves four rather straightforward applications of Green’s theorem or the representation theorem, though it is a remarkable result in that waves from either source that traverse the boundary many times are appropriately separated. The intermediate results at the four steps of the derivation also shed light on the possibility of acoustic shielding from unwanted sources without knowledge of the Green function for the medium.     

Combining two seismic experiments to attenuate free-surface multiples in OBC data*

The current inverse scattering solution used for multiple attenuation of marine seismic reflection data assumes that sources and receivers are located in the water. To adapt this solution to the ocean-bottom cable (OBC) experiment where receivers are located on the sea-floor, we have proposed combining the conventional marine surface seismic reflection data (streamer data) with OBC data. The streamer data add to the OBC data some of the wave paths needed for multiple attenuation. This combination has allowed us to develop a multiple attenuation method for OBC data which does not require any knowledge of the subsurface and which takes into account all free-surface multiples, including receiver ghosts. A non-linear synthetic data example consisting of pressure and particle velocity fields is used to illustrate the procedure.     

基于萤火虫算法的雷瑞波非线性反演(英文)

雷瑞波具有强振幅、低频和低速的特点,在反射地震勘探中通常是需要被压制的强噪声。本文研究如何利用雷瑞波获取近地表地层的横波速度和地下结构,选取萤火虫优化算法进行面波的反演,萤火虫优化算法是一种新的粒子群算法理论,具有稳定、快捷、全局搜索等特点。针对萤火虫优化算法优缺点进行了讨论和改进,通过对理论模型和野外数据的测试应用,将提取的瑞利面波频散曲线反演得到横波速度信息。结果表明萤火虫优化算法能实现面波非线性反演,并具有分辨率高、抗干扰能力强等优点和实际使用前景。     

Passive seismic data primary estimation and noise removal via focal‐denoising closed‐loop surface‐related multiple elimination based on 3D L1‐norm sparse inversion

Passive seismic has recently attracted a great deal of attention because non‐artificial source is used in subsurface imaging. The utilization of passive source is low cost compared with artificial‐source exploration. In general, constructing virtual shot gathers by using cross‐correlation is a preliminary step in passive seismic data processing, which provides the basis for applying conventional seismic processing methods. However, the subsurface structure is not uniformly illuminated by passive sources, which leads to that the ray path of passive seismic does not fit the hyperbolic hypothesis. Thereby, travel time is incorrect in the virtual shot gathers. Besides, the cross‐correlation results are contaminated by incoherent noise since the passive sources are always natural. Such noise is kinematically similar to seismic events and challenging to be attenuated, which will inevitably reduce the accuracy in the subsequent process. Although primary estimation for transient‐source seismic data has already been proposed, it is not feasible to noise‐source seismic data due to the incoherent noise. To overcome the above problems, we proposed to combine focal transform and local similarity into a highly integrated operator and then added it into the closed‐loop surface‐related multiple elimination based on the 3D L1‐norm sparse inversion framework. Results proved that the method was capable of reliably estimating noise‐free primaries and correcting travel time at far offsets for a foresaid virtual shot gathers in a simultaneous closed‐loop inversion manner.     

塔巴庙低渗致密砂岩渗透率有效应力定律实验研究

为了从实验角度探索低渗砂岩是否存在很小的ESCK值及重新认识低渗砂岩ESCK的变化规律,用两种修正的析因设计方案开展了塔巴庙低渗致密砂岩渗透率有效应力定律实验研究.一种方案包含了3个循环,每个循环是在孔隙流体压力不变,通过加载和卸载围压完成的;另一种方案包含4个循环,每个循环是在围压不变,通过降低和增加孔隙流体压力完成的.采用稳态法采集每个测点的渗透率值,并选择合适的经验模型拟合实验数据.为了使模型更好地拟合实验数据,本文采用最大似然函数法确定的转换系数转换实验得到的渗透率数据,使拟合得到的经验模型计算的渗透率值与实验值偏差的联合概率密度趋于极大值且残余平方和最小.拟合得到的经验模型可以用渗透率-围压-内压响应面直观地表示,再用Bernabe的ESCK计算式将这一响应面转换成渗透率有效应力系数ESCK-围压-内压响应面.ESCK响应面的响应特征表明,ESCK随围压和孔隙流体压力的变化而变化,随围压的增加而降低,随内压的增加而增加,反之亦然;实验研究还表明,ESCK的变化范围在0.0~1.02之间,这一变化范围和以往实验的结果存在巨大的差异,为此,文中分析了产生这一差异原因,同时提出一新的机理模型解释了实验获取的ESCK的变化规律是低渗岩石中微裂缝和孔隙变形共同作用的结果.     

HFT events: Shallow moonquakes?

A few large distant seismic events of distinctly high signal frequency, designated HFT (high-frequency teleseismic) events, are observed yearly by the Apollo lunar seismic network. Their sources are located on or near the surface of the moon, leaving a large gap in seismic activity between the zones of HFT sources and deep moonquakes. No strong regularities are found in either their spatial or temporal distributions. Several working hypotheses for the identity of these sources have been advanced, but many characteristics of the events seem to favor a hypothesis that they are shallow moonquakes. Simultaneous observations of other lunar phenomena may eventually enable the determination of their true identity.     

Spacecraft potential control aboard Equator-S as a test for Cluster-II

The payload of Equator-S was complemented by the potential control device (PCD) to stabilise the electric potential of the spacecraft with respect to the ambient plasma. Low potentials are essential for accurate measurements of the thermal plasma. The design of PCD is inherited from instruments for Geotail and Cluster and utilises liquid metal ion sources generating a beam of indium ions at several keV. The set-up of the instrument and its interaction with the plasma instruments on board is presented. When the instrument was switched on during commissioning, unexpectedly high ignition and operating voltages of some ion emitters were observed. An extensive investigation was initiated and the results, which lead to an improved design for Cluster-II, are summarised. The cause of the abnormal behaviour could be linked to surface contamination of some emitters, which will be monitored and cured by on-board procedures in future. The mission operations on Equator-S were not at all affected, because of the high redundancy built into the instrument so that a sufficient number of perfectly operating emitters were available and were turned on routinely throughout the mission. Observations of the effect of spacecraft potential control on the plasma remained limited to just one event on January 8, 1998, which is analysed in detail. It is concluded that the ion beam lead to the predicted improvement of the particle measurements even outside the low density regions of the magnetosphere where the effect of spacecraft potential control would have been much more pronounced, and that the similar instruments for the four Cluster-II spacecraft to be launched in 2000 will be very important to ensure accurate plasma data from this mission.