关于JCX 3型井磁仪器及其资料整理

经参加全国危机矿山示范项目的实际工作,认为目前重庆地质仪器厂生产的新型JCX 3井中三分量磁测仪探头倾斜面向北时X读数正值最大,倾斜面向南时X读数负值最大,X分量指向井斜方向。两者在直角坐标系中虽然同样符合左手法则,但新老仪器的井下磁测传感器x、y坐标相差90°。因而,对于当前广泛在用的数十台JCX 3型仪器的测量结果,如果仍然套用建立在垂向系统基础上的矢量计算公式,就会导致应用和解释错误,因此有必要共同探讨。     

广义非线性动态处理模型及其解算方法

本文针对当今各国高层领导和科学家十分关注并大力倡导的“数字地球”、“数字国家”、“数字城市”、“数字矿山”等科学工程构建中遇到的大量的多源、多维、多类型、多时态、多精度并具有非线性特征的联合数据处理问题的特点,建立了一个广义非线性动态联合数据处理模型及其相应的广义非线性最小二乘模型。针对该模型规模大、维数高的特点,借鉴多变量函数寻优的“变量轮换法”或“因素交替法”的思想、结合无记忆牛顿法,建立了一个解算算法,该算法将大规模的优化问题分解为两个较低规模的优化问题进行解算,降低了问题的规模,借助无记忆牛顿法,减少了存储量,特别适合大规模问题的解算。     

等高线悖论与广义等高线定义的研究

根据等高线是地图上可视化的线并与地图比例尺密切相关的条件,推出了等高线悖论,即等高线是地面上一定范围内高程相等或不相等的点的连线在水平面上的投影。提出了点元概念,推导了等高线是地面上的等高点元集在地图平面上的同胚或拓朴映射的广义等高线定义。本定义包容了传统的等高线定义并可解释制图综合规律。因此,广义等高线定义更适用于等高线可视化的地图应用环境。     

Validation of ground‐motion simulations for historical events using MDoF systems

The study presented in this paper addresses the issue of engineering validation of Graves and Pitarka's (2010) hybrid broadband ground motion simulation methodology with respect to some well‐recorded historical events and considering the response of multiple degrees of freedom (MDoF) systems. Herein, validation encompasses detailed assessment of how similar is, for a given event, the seismic response due to comparable hybrid broadband simulated records and real records. In the first part of this study, in order to investigate the dynamic response of a wide range of buildings, MDoF structures are modeled as elastic continuum systems consisting of a combination of a flexural cantilever beam coupled with a shear cantilever beam. A number of such continuum systems are selected including the following: (1) 16 oscillation periods between 0.1 and 6 s; (2) three shear to flexural deformation ratios to represent respectively shear‐wall structures, dual systems, and moment‐resisting frames; and (3) two stiffness distributions along the height of the systems, that is, uniform and linear. Demand spectra in terms of generalized maximum interstory drift ratio (IDR) and peak floor acceleration (PFA) are derived using simulations and actual recordings for four historical earthquakes, namely, the 1979 M_w 6.5 Imperial Valley earthquake, 1989 M_w 6.8 Loma Prieta earthquake, 1992 M_w 7.2 Landers earthquake, and 1994 M_w 6.7 Northridge earthquake. In the second part, for two nonlinear case study structures, the IDR and PFA distributions over the height and their statistics, are obtained and compared for both recorded and simulated time histories. These structures are steel moment frames designed for high seismic hazard, 20‐story high‐rise and 6‐story low‐rise buildings. The results from this study highlight the similarities and differences between simulated and real records in terms of median and intra‐event standard deviation of logs of seismic demands for MDoF building systems. This general agreement, in a broad range of moderate and long periods, may provide confidence in the use of the simulation methodology for engineering applications, whereas the discrepancies, statistically significant only at short periods, may help in addressing improvements in generation of synthetic records. Copyright © 2013 John Wiley & Sons, Ltd.     

Trimmed L-moments (1,0) for the generalized Pareto distribution

Abstract

Statistical analysis of extremes is often used for predicting the higher return-period events. In this paper, the trimmed L-moments with one smallest value trimmed—TL-moments (1,0)—are introduced as an alternative way to estimate floods for high return periods. The TL-moments (1,0) have an ability to reduce the undesirable influence that a small value in the statistical sample might have on a large return period. The main objective of this study is to derive the TL-moments (1,0) for the generalized Pareto (GPA) distribution. The performance of the TL-moments (1,0) was compared with L-moments through Monte Carlo simulation based on the streamflow data of northern Peninsular Malaysia. The result shows that, for some cases, the use of TL-moments (1,0) is a better option as compared to L-moments in modelling those series.

Citation Ahmad, U.N., Shabri, A. & Zakaria, Z.A. (2011) Trimmed L-moments (1,0) for the generalized Pareto distribution. Hydrol.Sci. J. 56(6), 1053–1060.     

Dynamic model of mesoscale eddies. Eddy parameterization for coarse resolution ocean circulation models

In the framework of the eddy dynamic model developed in two previous papers (Dubovikov, M.S., Dynamical model of mesoscale eddies, Geophys. Astophys. Fluid Dyn., 2003, 97, 311–358; Canuto, V.M. and Dubovikov, M.S., Modeling mesoscale eddies, Ocean Modelling, 2004, 8, 1–30 referred as I–II), we compute the contribution of unresolved mesoscale eddies to the large-scale dynamic equations of the ocean. In isopycnal coordinates, in addition to the bolus velocity discussed in I–II, the mesoscale contribution to the large scale momentum equation is derived. Its form is quite different from the traditional down-gradient parameterization. The model solutions in isopycnal coordinates are transformed to level coordinates to parameterize the eddy contributions to the corresponding large scale density and momentum equations. In the former, the contributions due to the eddy induced velocity and to the residual density flux across mean isopycnals (so called Σ-term) are derived, both contributions being shown to be of the same order. As for the large scale momentum equation, as well as in isopycnal coordinates, the eddy contribution has a form which is quite different from the down-gradient expression.     

Simulating hourly rainfall occurrence within an equatorial rainforest,Borneo Island / Simulation de l'occurrence de pluie horaire au sein de la forêt équatoriale,Ile de Bornéo

Abstract

The ability to simulate characteristics of the diurnal cycle of rainfall occurrence, and its evolution over the seasons is important to the forecasting of hydrological impacts resulting from land-use and climate changes within the humid tropics. This stochastic modelling study uses a generalized linear model (GLM) solution to second-order Markov chain models, as these discrete models are better at describing binary occurrence processes on an hourly time-scale than continuous-time approaches such as stochastic state-space models. We show that transition probabilities derived by the Markov chain method need to be time-varying rather than stationary to simulate the evolution of the diurnal cycle of rainfall occurrence over a Southeast Asian monsoon sequence. The conceptual and pragmatic links between discrete diurnal processes and continuous processes occurring over seasonal periods are thereby simulated within the same model.     

Statistics of extremes and estimation of extreme rainfall: II. Empirical investigation of long rainfall records / Statistiques de valeurs extrêmes et estimation de précipitations extrêmes: II. Recherche empirique sur de longues séries de précipitations

Abstract

Abstract In the first part of this study, theoretical analyses showed that the Gumbel distribution is quite unlikely to apply to hydrological extremes and that the extreme value distribution of type II (EV2) is a more consistent choice. Based on these theoretical analyses, an extensive empirical investigation is performed using a collection of 169 of the longest available rainfall records worldwide, each having 100–154 years of data. This verifies the theoretical results. In addition, it shows that the shape parameter of the EV2 distribution is constant for all examined geographical zones (Europe and North America), with value κ = 0.15. This simplifies the fitting and the general mathematical handling of the distribution, which become as simple as those of the Gumbel distribution.     

地方坐标系到CGCS 2000的转换在电网GIS中的应用

随着国家电网公司空间信息服务平台的建设和2000国家大地坐标系(CGCS 2000)的启用,迫切需要将各省不同坐标系的空间数据统一转换到CGCS 2000,从而实现电网空间信息服务平台的纵向贯通及空间数据资源的整合和共享。本文根据电力GIS空间数据资源的特点及行业情况,阐述了电网GIS中现有成果从独立坐标系向CGCS 2000转换的实现,并提出了对转换成果的验证方法。     

Polarimetric synthetic aperture radar change detection for specific land cover types

This paper presents a supervised polarimetric synthetic aperture radar (PolSAR) change detection method applied to specific land cover types. For each pixel of a PolSAR image, its target scattering vector can be modeled as having a complex multivariate normal distribution. Based on this assumption, the joint distribution of two corresponding vectors in a pair of PolSAR images is derived. Then, a generalized likelihood ratio test statistic for the equality of two likelihood functions of such joint distribution is considered and a maximum likelihood distance measure for specific land cover types is presented. Subsequently, the Kittler and Illingworth minimum error threshold segmentation method is applied to extract the specific changed areas. Experiments on two repeat-pass Radarsat-2 fully polarimetric images of Suzhou, China, demonstrate that the proposed change detection method gives a good performance in determining the specific changed areas in PolSAR images, especially the areas that have changed to water.