双偏振雷达差分传播相移的五种滤波方法对比分析

为了对双线偏振雷达的差分传播相移ΦDP进行滤波降噪处理,提高双线偏振雷达识别粒子、估测降水的能力。通过对C波段双线偏振雷达差分传播相移分别采用滑动平均、中值滤波、FIR滤波、卡尔曼滤波和小波分析5种不同方法地介绍,利用中国气象科学研究院灾害天气国家重点实验室移动式C波段双线偏振雷达(POLC)2013年外场试验观测资料,对比分析了经5种方法处理后的差分传播相移ΦDP径向数据,及其分别在降水估测中应用的效果。结果表明:采用FIR滤波方法、卡尔曼滤波方法和小波分析方法对ΦDP径向数据进行处理的效果比传统的滑动平均和中值滤波方法好;其中小波分析处理的ΦDP径向数据波动最小,再经过变距离法拟合的KDP数据的负值能够得到最有效地抑制,最后应用到R(KDP)降水估测中的精度最高。     

Precise, automatic and fast method for vanishing point detection

A new automated approach for vanishing point detection in images of urban scenes is described. This method is based on the theorem of Thales. The main contribution of this paper is the automatic and simultaneous detection of all vanishing points of the image, achieved by converting this problem into the detection of circles in a complex cloud of points, in which each point corresponds to a segment and is associated with an uncertainty. This extraction of circles uses a RANSAC method, modified to improve its speed by using accumulation techniques (Hough transform or otherwise). This robust estimation is then refined by least squares error propagation using the individual variances of each segment. The algorithm is robust, its accuracy is optimised and it is entirely automatic. Its successful operation has been tested on a large number of images of varied urban scenes.     

局部角度域波传播步进算法研究

本文从非均匀介质中波动方程出发,提出了基于一般标架的相空间（局部角度域）波传播的步进算法. 该方法在构造单程波的步进算法时,在选择标架或正交基等方面有更大的自由度. 我们以不随频率及深度变化且具有变尺度特性的Gabor_Daubechies紧标架为例,给出了单程波传输算子的具体形式及相应的波场步进算法;详细讨论了基于Gabor_Daubechies标架的传输算子的高频渐近展开问题,得出了在高频、小传输步长条件下传输算子的近似解析表达式,并给出使用条件. 通过模型算例,比较了精确传输算子与高频近似传输算子用于非均匀介质中波传播的结果,说明在一定条件下由两者得出的波场几乎是相同的.     

基于分数阶时间导数的黏弹性衰减VTI介质中平面波传播

目前在地震勘探频带范围内通常假设品质因子Q与频率无关,且呈衰减各向同性.事实上,相比较速度各向异性,介质的衰减各向异性同样不可忽视.本文将衰减各向异性和速度各向异性二者与常Q模型相结合,建立了黏弹性衰减VTI介质模型,并基于分数阶时间导数理论,给出了对应的本构关系和波动方程.利用均匀平面波分析和Poynting定理,推...     

Modeling 3D thermal fracture propagation by transient cooling using virtual multidimensional internal bonds

Thermal fracturing can play an important role in development of unconventional petroleum and geothermal resources. Thermal fractures can result from the nonlinear deformation of the rock in response to thermal stress related to cold water injection as well as heating. Before the rock reaches the final failure stage, material softening and bulk modulus degradation can cause changes in the thermo‐mechanical properties of the solid. In order to capture this aspect of the rock fracture, a virtual multidimensional internal bond‐based thermo‐mechanical model is derived to track elastic, softening, and the failure stages of the rock in response to the temporal changes of its temperature field. The variations in thermo‐mechanical properties of the rock are derived from a nonlinear constitutive model. To represent the thermo‐mechanical behavior of pre‐existing fractures, the element partition method is employed. Using the model, numerical simulation of 3D thermal fracture propagation in brittle rock is carried out. Results of numerical simulations provide evidence of model verification and illustrate nonlinear thermal response and fracture development in rock under uniform cooling. In addition, fracture coalescence in a cluster of fractures under thermal stress is illustrated, and the process of thermal fracturing from a wellbore is captured. Results underscore the importance of thermal stress in reservoir stimulation and show the effectiveness of the model to predict 3D thermal fracturing. Copyright © 2016 John Wiley & Sons, Ltd.     

Variability of temperature,salinity and tidally‐averaged density in the middle estuary of the St Lawrence

Abstract

Surveys in the Middle Estuary of the St Lawrence have yielded a data base consisting of more than 15,000 T‐S pairs distributed over 62 13‐h profiling stations. Although the T‐S curves at each station are remarkably linear, the variability of the slopes and intercepts of the lines is considerable. The means and standard deviations of the temperature and salinity at each individual station are not explicable in terms of linear combinations of the parameters for location in the Estuary, the upstream water properties, the phase of the spring‐neap cycle and the tidal energies.

It is shown that the tidally‐averaged density structure is separable into horizontal and vertical components and that its vertical variation over the whole Estuary may be explained by any one of three different functional forms. However, its horizontal variation is not explicable in terms of linear combinations of the parameters mentioned in the paragraph above.

Plots of the horizontal variations in temperature, salinity or density may only be meaningful if the data are collected synoptically, and even then cannot be considered to be accurate over time‐scales longer than one tidal cycle.     

Water waves generated by instantaneous disturbances at the bed of a sloping beach

The two-dimensional problem of the generation of water waves due to instantaneous disturbances prescribed at the bed of a beach sloping at an arbitrary angle is studied here. It is formulated in terms of an initial-boundary-value problem for the velocity potential describing the motion in the fluid region assuming the linear theory. Using the Laplace transform in time and the Mellin transform in distance, the problem is reduced to solving a difference equation whose method of solution is of considerable importance in the literature. The form of the free surface is obtained in terms of a multiple infinite integral that is evaluated by the method of steepest-descent. For some prescribed forms of the disturbance at the bed of the beach, the free surface is depicted in a number of figures for different beach angles. It is observed that as the beach angle decreases, the maximum wave height increases, which is plausible.     

Shear wave propagation along infinite slopes: A theoretically based numerical study

In this paper, the seismic response of ‘infinitely’ long slopes is numerically analysed via the formulation of a 1D analytical/numerical model, in which the soil mechanical behaviour is assumed to be elasto‐perfectly viscoplastic and simple shear (SS) kinematical constraints are imposed. In order to simplify the problem, a theoretically based procedure to set up a fully 1D shear constitutive model is defined, within which the mechanical response of a multiaxial relationship is condensed. The use of a 1D shear constitutive model is aimed at reducing the number of unknowns and, therefore, the computational costs. In particular, the case of the Mohr–Coulomb yield criterion is considered, while an enhanced Taylor–Galerkin finite element algorithm is employed to simulate the seismic wave propagation within the soil stratum. The proposed ‘condensation’/calibration procedure captures both the ‘pseudo’‐hardening pre‐failure behaviour and the influence of dilation on the occurrence of strain‐localization, which characterize, under SS conditions, the static response of virgin perfectly plastic soils. The effectiveness of the conceived method is shown with reference to freshly deposited deposits, while, in the case of highly overconsolidated strata, some difficulties arise because of the brittle behaviour induced both by unloading and non‐associativeness. Copyright © 2011 John Wiley & Sons, Ltd.     

Theoretical analyses of acidization dissolution front instability in fluid‐saturated carbonate rocks

This paper presents an instability theory that can be used to understand the fundamental behavior of an acidization dissolution front when it propagates in fluid‐saturated carbonate rocks. The proposed theory includes two fundamental concepts, namely the intrinsic time and length of an acidization dissolution system, and a theoretical criterion that involves the comparison of the Zhao number and its critical value of the acidization dissolution system. The intrinsic time is used to determine the time scale at which the acidization dissolution front is formed, while the intrinsic length is used to determine the length scale at which the instability of the acidization dissolution front can be initiated. Under the assumption that the acidization dissolution reaction is a fast process, the critical Zhao number, which is used to assess the instability likelihood of an acidization dissolution front propagating in fluid‐saturated carbonate rocks, has been derived in a strictly mathematical manner. Based on the proposed instability theory of a propagating acidization dissolution front, it has been theoretically recognized that: (i) the increase of the mineral dissolution ratio can stabilize the acidization dissolution front in fluid‐saturated carbonate rocks; (ii) the increase of the final porosity of the carbonate rock can destabilize the acidization dissolution front, while the increase of the initial porosity can stabilize the acidization dissolution front in fluid‐saturated carbonate rocks; (iii) the increase of the mineral dissolution ratio can cause an increase in the dimensionless propagation speed of the acidization dissolution front; (iv) the increase of the initial porosity can enable the acidization dissolution front to propagate faster, while the increase of the final porosity can enable the acidization dissolution front to propagate slower in the acidization dissolution system. Copyright © 2012 John Wiley & Sons, Ltd.