估算模式 AERSCREEN的参数敏感性研究

《环境影响评价技术导则大气环境》（HJ2.2-2018）推荐的估算模式AERSCREEN在气象和地形资料的处理以及建筑物下洗等多个方面做了改进。利用估算模式AERSCREEN,针对30 m左右高度的点源,进行了不同排放参数、不同气象条件下最大落地浓度的敏感性试验。结果表明：随着烟气出口流速的增大,地面浓度最大值逐渐减小;随着波文比的变化,地面浓度最大值没有明显的变化;随着地面粗糙度的增大,地面浓度最大值逐渐减小;随着烟气出口温度的增高,地面浓度最大值逐渐减小;当烟气温度为75℃,粗糙度达到1.3 m时,地面浓度达到最小;随着反照率的增大,地面浓度最大值逐渐减小;随着烟囱高度的增大,地面浓度最大值逐渐减小;在各种烟囱高度条件下,随着最高环境温度的增高,地面浓度最大值逐渐增大;而在各种环境温度条件下,随着烟囱高度的增高,地面浓度最大值在逐渐减小;模式中,随着最低环境温度的增高,地面浓度最大值没有变化;但随着最小风速的增大,模拟得到的地面浓度最大值会逐渐减小。     

Identification-based simplified model of large container ships using support vector machines and artificial bee colony algorithm

The 6 degrees of freedom (DOF) model with a high degree of complexity for capturing ship dynamics is generally able to track the nonlinear and coupling dynamics of ships. However, the 6 DOF model makes challenges in estimating model coefficients and designing the model-based control. Therefore, simplified ship dynamic models within allowed accuracy are essential. This paper simplified the 6 DOF nonlinear dynamic model of ships into two decoupled models including the speed model and the steering model through reasonable assumptions. Those models were tested through maneuvering simulations of a container ship with a 4 DOF dynamic model. Support vector machines (SVM) optimized by the artificial bee colony algorithm (ABC) was used to identify parameters of speed and steering models by analyzing the rudder angle, propeller shaft speed, surge and sway velocities, and yaw rate from simulated data extracted from a series of maneuvers made by the container ship. Comparisons with the first order linear and nonlinear Nomoto models show that the simplified nonlinear steering model can capture more complicated dynamics and performs better. Additionally, comparisons among three different parameter identification methods demonstrate similar identification results but the different performance involving the applicability and effectiveness. SVM optimized by ABC is relatively convenient and effective for parameter identification of ship simplified dynamic models.     

提高地震定位精度方法的研究

本文对地震定位产生误差的原因进行分析,并且对提高地震定位精度的方法进行了研究。随着沈阳数字化地震台网即将投入运行,需要将原有的ＤＯＳ环境下运行,用ＢＡＳＩＣ语言编写的地震定位软件,移植到新的系统运行环境,并且在吸收原定位系统的优点的基础上作了重大改进,进一步提高了网缘地震定位精度和可靠性,对提高地震台网监测能力有重大的意义。     

Development and application of a coupled-process parameter inversion model based on the maximum likelihood estimation method

The coupled flow-mass transport inverse problem is formulated using the maximum likelihood estimation concept. An evolutionary computational algorithm, the genetic algorithm, is applied to search for a global or near-global solution. The resulting inverse model allows for flow and transport parameter estimation, based on inversion of spatial and temporal distributions of head and concentration measurements. Numerical experiments using a subset of the three-dimensional tracer tests conducted at the Columbus, Mississippi site are presented to test the model's ability to identify a wide range of parameters and parametrization schemes. The results indicate that the model can be applied to identify zoned parameters of hydraulic conductivity, geostatistical parameters of the hydraulic conductivity field, angle of hydraulic conductivity anisotropy, solute hydrodynamic dispersivity, and sorption parameters. The identification criterion, or objective function residual, is shown to decrease significantly as the complexity of the hydraulic conductivity parametrization is increased. Predictive modeling using the estimated parameters indicated that the geostatistical hydraulic conductivity distribution scheme produced good agreement between simulated and observed heads and concentrations. The genetic algorithm, while providing apparently robust solutions, is found to be considerably less efficient computationally than a quasi-Newton algorithm.     

钢板阻尼墙基于Bouc-Wen模型的模拟和工程应用研究

钢板阻尼墙抗震性能好,设计自由度大,在工程中使用越来越多。为了更高精度地模拟钢板阻尼墙的力学性能,提出基于Bouc-Wen模型的模拟分析方法。给出一种基于Bouc-Wen模型钢板阻尼墙参数辨识方法,根据试验数据和仿真结果的验证对比,文中提出的模拟方法可以较好地反映钢板阻尼墙的的力学性能曲线。该模拟方法可得到较好的仿真效果,对实际工程设计分析钢板阻尼墙对结构的抗震性能影响有高精度、数据化的帮助,可为进一步工程应用和设计分析提供参考。     

An efficiency-improved genetic algorithm and its application on multimodal functions and a 2D common reflection surface stacking problem

Although Genetic Algorithms have found many successful applications in the field of exploration geophysics, the convergence speed remains a big challenge as Genetic Algorithms usually require a huge amount of fitness function evaluations. In this paper, we propose an efficiency-improved Genetic Algorithm, which has both a good global search capability and a good local search capability, and is also capable of robustly handling the premature convergence challenge commonly seen in linear and directed non-linear optimization methods. In our new genetic algorithm, the global search capability is performed via a modified island model, while the local search capability is provided by a novel self-adaptive differential evolution fine tuning scheme. Premature convergence is dealt with via a local exhaustive search method. We first demonstrate the much improved convergence speed of this efficiency-improved Genetic Algorithm over that of our previously proposed advanced Genetic Algorithm on several multimodal functions. We further demonstrate the effectiveness of our efficiency-improved Genetic Algorithm by applying it to a two-dimensional common reflection surface stacking problem, which is a highly nonlinear geophysical optimization problem, to obtain very encouraging results.     

黏滑断层隧道减错措施参数对减错效果的影响分析

为提高黏滑断层隧道的结构安全性和稳定性,对黏滑断层隧道设置不同缝宽减错缝、不同刚度减错层的减错效果进行研究。研究结果表明:断层黏滑错动对上盘隧道的影响远大于下盘;减错缝对上盘部分隧道结构的减错效果优于下盘,其中上盘减错效果最大为24.50%,下盘减错效果最大为9.26%;减错层对下盘部分隧道结构的减错效果略优于上盘,其中下盘减错效果最大为105.32%,上盘减错效果最大为78.07%;随着减错缝宽度的增加,隧道上盘减错效果变好,下盘缝宽10—15cm减错效果最好;随着减错层弹性模量的增加,隧道上下盘减错效果降低,当减错层弹性模量增加到一定程度（约100MPa）,减错效果趋于稳定。研究成果可为黏滑断层隧道的减错结构设计及施工提供参考。     

Reflection-constrained 2D and 3D non-hyperbolic moveout analysis using particle swarm optimization

The objective of moveout parameter inversion is to derive sets of parameter models that can be used for moveout correction and stacking at each common midpoint location to increase the signal-to-noise ratio of the data and to provide insights into the kinematic characteristics of the data amongst other things. In this paper, we introduce a data-driven user-constrained optimization scheme that utilizes manual picks at a point on each reflector within a common midpoint gather to constrain the search space in which an optimization procedure can search for the optimal parameter sets at each reflection. The picks are used to create boundary curves which can be derived approximately via an optimization technique or analytically via the derivation of an analytical bounds function. In this paper, we derive analytical forms of bounds functions for four different moveout cases. These are normal moveout, non-hyperbolic moveout, azimuthally dependent normal moveout and azimuthally dependent non-hyperbolic moveout. The optimization procedure utilized here to search for the optimal moveout parameters is the particle swarm optimization technique. However, any metaheuristic optimization procedure could be modified to account for the constraints introduced in this paper. The technique is tested on two-layer synthetic models based on three of the four moveout cases discussed in this paper. It is also applied to an elastic forward modelled synthetic model called the HESS model, and finally to real 2D land data from Alaska. The resultant stacks show a marked improvement in the signal-to-noise ratio compared to the raw stacks. The results for the normal moveout, non-hyperbolic moveout and azimuthally dependent normal moveout tests suggest that the method is viable for said models. Results demonstrate that our method offers potential as an alternative to conventional parameter picking and inversion schemes, particularly for some cases where the number of parameters in the moveout approximation is 2 or greater.     

Multiwell study of seismic attenuation at the CO2CRC Otway project geosequestration site: Comparison of amplitude decay,centroid frequency shift and 1D waveform inversion methods

At the CO2CRC Otway geosequestration site, the abundance of borehole seismic and logging data provides a unique opportunity to compare techniques of Q (measure of attenuation) estimation and validate their reliability. Specifically, we test conventional time-domain amplitude decay and spectral-domain centroid frequency shift methods versus the 1D waveform inversion constrained by well logs on a set of zero-offset vertical seismic profiles. The amplitude decay and centroid frequency shift methods of Q estimation assume that a seismic pulse propagates in a homogeneous medium and ignore the interference of the propagating wave with short-period multiples. The waveform inversion explicitly models multiple scattering and interference on a stack of thin layers using high-resolution data from sonic and density logs. This allows for stable Q estimation in small depth windows (in this study, 150 m), and separation of the frequency-dependent layer-induced scattering from intrinsic absorption. Besides, the inversion takes into account band-limited nature of seismic data, and thus, it is less dependent on the operating frequency bandwidth than on the other methods. However, all considered methods of Q estimation are unreliable in the intervals where subsurface significantly deviates from 1D geometry. At the Otway site, the attenuation estimates are distorted by sub-vertical faults close to the boreholes. Analysis of repeated vertical seismic profiles reveals that 15 kt injection of the CO₂-rich fluid into a thin saline aquifer at 1.5 km depth does not induce detectable absorption of P-waves at generated frequencies 5–150 Hz, most likely because the CO₂ plume in the monitoring well is thin, <15 m. At the Otway research site, strong attenuation Q ≈ 30–50 is observed only in shaly formations (Skull Creek Mudstone, Belfast Mudstone). Layer-induced scattering attenuation is negligible except for a few intervals, namely 500–650 m from the surface, and near the injection interval, at around 1400–1550 m, where Q_scat ≈ 50–65.     

On possible fluid detection in equivalent transversely isotropic media

We consider a transversely isotropic medium that is long-wave equivalent to a stack of thin, parallel, isotropic layers and is obtained using the Backus average. In such media, we analyse the relations among anisotropy parameters; Thomsen parameters, ε and δ, and a new parameter ϕ. We discuss the last parameter and show its essential properties; it equals 0 in the case of isotropy of equivalent medium and/or constant Lamé coefficient λ in layers. The second property occurs to make ϕ sensitive to variations of λ in thin-bedded sequences. According to Gassmann, in isotropic media the variation of fluid content affects only the Lamé coefficient λ, not μ; thus, the sensitivity to changes of λ is an essential property in the context of possible detection of fluids. We show algebraically and numerically that ϕ is more sensitive to these variations than ε or δ. Nevertheless, each of these parameters is dependent on the changes of μ; to understand this influence, we exhibit comprehensive tables that illustrate the behaviour of anisotropy parameters with respect to specific variations of λ and μ. The changes of μ in layers can be presented by the Thomsen parameter γ that depends on them solely. Hence, knowing the values of elasticity coefficients of equivalent transversely isotropic medium, we may compute ϕ and γ, and based on the aforementioned tables, we predict the expected variation of λ; in this way, we propose a new method of possible fluid detection. Also, we show that the prior approach of possible detection of fluids, proposed by Berryman et al., may be unreliable in specific cases. To establish our results, we use the Monte Carlo method; for the range and chosen variations of Lamé coefficients λ and μ – relevant to sandstones – we generate these coefficients in thin layers and, after the averaging process, we obtain an equivalent transversely isotropic medium. We repeat that process numerous times to get many equivalent transversely isotropic media, and – for each of them  – we compute their anisotropy parameters. We illustrate ϕ, ε and δ in the form of cross-plots that are relevant to the chosen variations of λ and μ. Additionally, we present a table with the computed ranges of anisotropy parameters that correspond to different variations of Lamé coefficients.