基于修正SWAT模型的岩溶地区非点源污染模拟初探——以横港河流域为例

岩溶流域含水系统的主要特征之一是连通地表的落水洞等垂直管道将近水平的地下暗河联系起来,降水及其形成的地表径流可以通过这些管道迅速地灌入地下河系,从而改变了水及其所携带的非点源污染物质在垂直与水平方向的传输速度与数量,使岩溶流域内地表-地下之间的物质交换与传输过程变得比较复杂;应用广泛的SWAT模型在模拟岩溶地区的水文、水质时会存在一些不足与局限.为此,本文针对岩溶水系统特征,引入落水洞、伏流、暗河的水文过程以及主要营养盐的输移过程,修正SWAT模型原有的水文循环过程及相关算法,改变其只适用于松散均匀介质流域非点源污染模拟的单一特征,并研究建立适合于岩溶流域的非点源污染模型和相应的模拟方法.选取横港河流域岩溶地区作为非点源污染的对象,应用修正后的模型通过控制性的模拟方法和敏感性性分析,定量评估落水洞、伏流、暗河等岩溶特征对氮、磷等主要非点源污染物质输移的影响及其带来的时空效应,并进一步探讨落水洞、伏流、暗河等对地表-地下水文与营养盐的交互作用及转换机理.结果表明,岩溶特征对流域的氮、磷负荷有增加作用,其中总磷的增加明显大于总氮的增加,总磷和总氮的增量分别为0.86%和2.12%;植被岩溶指数的增加会导致流域可溶性磷、有机磷的产出量增加,有机氮、地表产流中硝酸氮和沉积磷的产出量则居其次,落水洞改变了降雨的产流方式,增加了落水洞所在流域的有机磷和有机氮的产出,其增量变化在0~0.7和0~0.3 kg/hm²之间.     

Effects of aperture on Marchenko focussing functions and their radiation behaviour at depth

A focussing function is a specially constructed field that focusses on to a purely downgoing pulse at a specified subsurface position upon injection into the medium. Such focussing functions are key ingredients in the Marchenko method and in its applications such as retrieving Green's functions, redatuming, imaging with multiples and synthesizing the response of virtual sources/receiver arrays at depth. In this study, we show how the focussing function and its corresponding focussed response at a specified subsurface position are heavily influenced by the aperture of the source/receiver array at the surface. We describe such effects by considering focussing functions in the context of time-domain imaging, offering explicit connections between time processing and Marchenko focussing. In particular, we show that the focussed response radiates in the direction perpendicular to the line drawn from the centre of the surface data array aperture to the focussed position in the time-imaging domain, that is, in time-migration coordinates. The corresponding direction in the Cartesian domain follows from the sum (superposition) of the time-domain direction and the directional change due to time-to-depth conversion. Therefore, the result from this study provides a better understanding of focussing functions and has implications in applications such as the construction of amplitude-preserving redatuming and imaging, where the directional dependence of the focussed response plays a key role in controlling amplitude distortions.     

Non-point source pollution estimation in the Pingqiao River Basin,China, using a spatial hydrograph-separation approach

For efficient and targeted management, this study demonstrates a recently developed non-point source (NPS) pollution model for a year-long estimation in the Pingqiao River Basin (22.3 km²) in China. This simple but physically reasonable model estimates NPS export in terms of land use by reflecting spatial hydrological features and source runoff measurements under different land-use types. The NPS export was separately analysed by a distributed hydrological model, a spatial hydrograph-separation technique, and an empirical water quality sub-model. Simulation results suggest that 57 890 kg of total nitrogen (TN) and 1148 kg of total phosphorus (TP) were delivered. The results, validated with observed stream concentrations, show relative errors of 23.3% for TN and 47.4% for TP. Countermeasures for urban areas (5.3% of total area) were prioritized because of the high contribution rate to TN (14.1%) and TP (26.2%) which is caused by the high degree of runoff (8.5%) and pollution source.     

隐伏正断层错动致地表破裂变形特征的研究

通过建立三维计算模型,对隐伏正断层在均匀错动、倾斜错动和翘倾错动方式下地表土体的应力路径、破裂和变形特征进行了研究。根据地表破裂临界值,分析了工程建设“避让带”的宽度和起始位置的变化特征。根据行业规范,提出工程建设“关注带”的确定方法,分析了“关注带”的宽度和起始位置的变化特征,得到以下主要结论:①在断层错动过程中,位于两侧的地表土体应力路径变化明显不同,下盘一侧和上盘一侧分别以三轴拉伸和三轴压缩为主;②地表强变形带与地表破裂带的分布并不一致,需要综合考虑等效塑性应变和总位移比2个指标来评价同震地表错动对建筑物的影响;③当隐伏断层错动的垂直位移达到3m时,工程建设“避让带”的宽度在10—90m范围内变化,受上覆土体厚度和断层倾角的影响最大,而工程建设“关注带”的宽度在150—400m范围内变化,受上覆土体的性质影响最大。     

THE ACTIVITY OF WESTERN LISHAN FAULT SINCE THE LATE PLEISTOCENE

Along the northern piedmont of Mt. Lishan, the characteristics and locations of the active normal Lishan fault in west of Huaqing Pool provide important evidences for determining the seismotectonic environment, seismic stability evaluation of engineering in the eastern Weihe Basin. After reviewing the results from high-density resistivity method, seismic profile data, geological drillhole section and trenching in west of the Huaqing Pool, it is found that the strike of western normal Lishan Fault changes from EW direction at the eastern part to the direction of N60°W, and the fault consists of two branches, dipping NE with a high dip angle of~75°. The artificial shallow seismic profile data reveals that the attitude of strata near Lishan Fault mainly dips to south, which is presumed to be related to the southward tilt movement of Mt. Lishan since the Cenozoic. The section of geological drillhole reveals that since the late middle Pleistocene, the displacement of the paleo-soil layer S₂ is about 10m. And the maximum displacement of western Lishan Fault recorded in the paleo-soil layer S₁ reaches 7.8m since the late Pleistocene. In addition, evidences from trench profile show that the western Lishan Fault was active at least 3 times since Malan loess deposition with ¹⁴ C dating age(32 170±530)Cal a BP. The multiple activities of the Lishan Fault result in a total displacement about 3.0m in the Malan loess layer L₁. The latest activity of the western Lishan Fault produced a displacement of about 0.9m in the early Holocene loess layer L₀((8 630±20)Cal a BP)and caused obvious tensile cracks in the Holocene dark leoss layer S₀((4 390±20)Cal a BP). Briefly, we have obtained a vertical movement rate of about 0.11~0.19mm/a since the Holocene((8 630±20)Cal a BP)in the western extension of the Lishan Fault, the recurrence interval of earthquakes on the fault is about(10.7±0.5)ka, and the co-seismic surface rupture in a single event is inferred to be about 0.9m.     

GEODETIC AND TELESEISMIC CONSTRAINTS ON SLIP DISTRIBUTION OF 2015 MW6.4 PISHAN EARTHQUAKE

On July 3^rd, 2015, a M_W6.4 earthquake occurred on Pishan County, Xinjiang, located in the front of western Kunlun thrust belt, which is the largest earthquake(M_W6.0~7.0)in the past 40 years in this region. In this study, we collected both the near-filed geodetic coseismic deformation observations including 4 GPS sites and one high-resolution ALOS-2 InSAR imagery, and far-field teleseismic P waveforms from 25 stations provided by IRIS/USGS, to invert the fault parameters(strike and dip)and coseismic rupture model of 2015 M_W6.4 Pishan earthquake. Using the finite fault theory, a non-linear simulated annealing algorithm was employed to resolve our joint inversion problem. The strike (120°~130°) and dip angle(35°~40°)of optimal models are different from that of some previous studies, and the dip change is strongly constrained by combined data than that of strike. In fixing the geometric parameters of optimal fault model, we also considered data weight(5)(geodetic data/teleseismic P waveforms)and constrained weight from moment and smooth factor(2.5). Clearly, our results indicate that the slip distribution mainly concentrates in the depth range from 9 to 16km and a length range of 20km along the strike direction, which is similar to the spatial distribution of the relocated aftershocks. The maximum slip is~95cm. The seismic moment release is 5.45×10¹⁸N·m, corresponding to M_W6.42. Compared with the single data set, geodetic data or teleseismic waveform, our joint inversion model could simultaneously constrain the seismic moment and slip distribution well, thus avoiding effectively a lower-resolution rupture distribution determined by teleseismic-only inversion and a bias released moment estimated by the geodetic-only inversion. Importantly, we should consider both the near-field geodetic data and far-field teleseismic data in retrieving the rupture model for accurately describing the seismogenic structure of active fault in western Kunlun region.     

2017年3月27日漾濞5.1级地震前后气枪震源观测走时数据的变化

以2017年3月27日漾濞5.1级地震为例,根据区域特性和信噪比要求,选取数据较为完整的6个台站记录的2017年1月1日～6月6日期间的宾川地震信号发射台气枪震源波形资料,采用互相关检测技术提取6个台站各自稳定震相的走时数据,并对漾濞5.1级地震前后走时数据的变化情况进行分析。结果表明,漾濞5.1级地震前后6个台站各自稳定震相存在较为明显的走时变化,且短期内走时变化具有较好的同步性,相关台站异常幅度大小和异常出现时间存在细小差异。地震发生前,6个台站走时低值异常过程明显,以YUL台最为显著。地震发生前后走时变化形态特征为双“V”型,漾濞5.1级地震发生在第1个“V”型末端。地震发生后,不同方位相关台站受地震的影响程度不同,走时波动大小存在差异。     

基于数学建模的潜在地震破裂面源检测方法研究

大多数地震破裂面源检测方法都是通过简化地震震源,将地震震源表示成线源或者点源,无法有效描述地震带地震破裂面源产状和大小,不适用地震震级较大的情况下地震危险性检测。因此提出基于数学建模的潜在地震破裂面源检测方法,在地震震级较大时仍能检测出地震危险性概率。选取适宜的地震基岩水平峰值加速度衰减关系,分析地震震级、破裂长度、破裂宽度相互关系,确定地震引起的潜在地震破裂面源大小,计算给定地震动小于在场点处产生地震动的概率,将该概率同地震动加速度衰减关系结合,得到地震动年超越概率,分析地震危险性。经过实验检测发现,所提方法检测出的年超越概率与峰值加速度、最大震级有关,该概率能精准表示地震带地震破裂面源产状和大小,说明该方法检测地震危险性是合理的。     

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.     

Attenuation Characteristics of Earthquake Ground Motion for Large Volume Airgun

In order to further deepen the understanding of seismic wave propagation characteristics induced by the large volume airgun source, experimental data from multiple fixed excitation points in Fujian Province were used to obtain the equivalent single excitation high signal-to-noise ratio velocity and displacement records through linear stacking and simulation techniques. Then the peak displacements of different epicentral distances were used to calculate the equivalent magnitude of the airgun source excitation at different fixed excitation points so as to establish the attenuation relationship between equivalent magnitude,epicenter distance and velocity peak. Our results show that:① Within 270 km of epicentral distance,for the large-volume airgun's single shot,the peak velocity range is about 700-4 nm/s,and the peak displacement range is about 200. 0-0. 2 nm;② The equivalent magnitude of the P-wave from the airgun source with a total capacity of 8,000 in 3 is 0. 181-0. 760,and the equivalent magnitude of the S-wave is 0. 294-0. 832. By contrast,the equivalent magnitude of the P-wave from the airgun source with a total capacity of 12,000 in 3 is 0. 533-0. 896,and the equivalent magnitude of the S-wave is 0. 611-0. 946. The S-wave energy is greater than the P-wave energy, and the excitation efficiency varies greatly with different excitation environment;③ The peak velocity increases with the equivalent magnitude,and decreases with the epicentral distance. The vertical component of the P-wave peak velocity is the largest among those three components,while the S-wave has the smallest vertical component and similar horizontal components. Hence,our research can provide an important basis for the quantitative judgment of the seismic wave propagation distance using the airgun and the design of the observation system in deep exploration or monitoring with airgun.