Robust Kalman filter-based dynamic state estimation of natural gas pipeline networks

To obtain the accurate transient states of the big scale natural gas pipeline networks under the bad data and non-zero mean noises conditions, a robust Kalman filter-based dynamic state estimation method is proposed using the linearized gas pipeline transient flow equations in this paper. Firstly, the dynamic state estimation model is built. Since the gas pipeline transient flow equations are less than the states, the boundary conditions are used as supplementary constraints to predict the transient states. To increase the measurement redundancy, the zero mass flow rate constraints at the sink nodes are taken as virtual measurements. Secondly, to ensure the stability under bad data condition, the robust Kalman filter algorithm is proposed by introducing a time-varying scalar matrix to regulate the measurement error variances correctly according to the innovation vector at every time step. At last, the proposed method is applied to a 30-node gas pipeline networks in several kinds of measurement conditions. The simulation shows that the proposed robust dynamic state estimation can decrease the effects of bad data and achieve better estimating results.Comment: Accepted by Mathematical Problems in Engineerin

First Arrival Time Auto-Picking Method Based on Multi-Time Windows Energy Ratio

First arrival time auto-picking technique plays an important role in seismic exploration. It is widely used in shallow layer tomography and static correction. Conventional method that based on sliding time windows energy ratio is not stable. Here a new method based on multi-time windows energy ratio is proposed. Combining with automatic quality control and phase-domain first arrival estimation technique, our method performs perfectly on seismic records of normal S/N ratio. In the computational process of conventional sliding time windows energy ratio method, first arrivals are often determined by the maximum energy ratio of two adjacent sliding time windows. It is well known that for low S/N ratio data the conventional picking is not effective, and for high S/N ratio data weak reflections are hardly detected. The reason is that first arrival time does not correspond to the maximum energy ratio. Meanwhile conventional method sometime picks local secondary extreme of energy ratio. The new method of multi-time windows energy ratio method takes both maximum and local secondary extreme in consideration. Hence new method promotes the stability and accuracy of first arrival picking. Combined with automatic quality control and phase-domain first arrival estimation, the new method performs well in its application in the middle part of Dzungarian Basin(Northwest China)

The Improvement of 3D Traveltime Tomographic Inversion Method

As 3D high-precision seismic exploration is more and more widely used in seismic data acquisition, traveltime tomographic inversion based on first arrivals is developed from 2D to 3D. However, magnanimity data of 3D traveltime inversion brings about the problem of data storage; the absence of first arrivals with near offset reduces the precision of shallow layer; the utilization of prior information, such as small refraction and micro-logging data, can improve the precision of 3D traveltime inversion. Therefore, we make some improvements in 3D traveltime inversion method. We take compression storage for large and sparse matrix, propose virtual receivers technology, and add prior information to tomographic inversion linear equations. The application in 3D real data indicates that the improvements can effectively improve 3D traveltime tomographic inversion.Key words: 3D seismic exploration; 3D traveltime inversion method; 3D traveltime tomographic inversio

A new family-based association test via a least-squares method

To test the association between a dichotomous phenotype and genetic marker based on family data, we propose a least-squares method using the vector of phenotypes and their cross products within each family. This new approach allows covariate adjustment and is numerically much simpler to implement compared to likelihood- based methods. The new approach is asymptotically equivalent to the generalized estimating equation approach with a diagonal working covariance matrix, thus avoiding some difficulties with the working covariance matrix reported previously in the literature. When applied to the data from Collaborative Study on the Genetics of Alcoholism, this new method shows a significant association between the marker rs1037475 and alcoholism

Responses of seasonal indicators to extreme droughts in southwest China

Significant impact of extreme droughts on human society and ecosystem has occurred in many places of the world, for example, Southwest China (SWC). Considerable research concentrated on analyzing causes and effects of droughts in SWC, but few studies have examined seasonal indicators, such as variations of surface water and vegetation phenology. With the ongoing satellite missions, more and more earth observation data become available to environmental studies. Exploring the responses of seasonal indicators from satellite data to drought is helpful for the future drought forecast and management. This study analyzed the seasonal responses of surface water and vegetation phenology to drought in SWC using the multi-source data including Seasonal Water Area (SWA), Permanent Water Area (PWA), Start of Season (SOS), End of Season (EOS), Length of Season (LOS), precipitation, temperature, solar radiation, evapotranspiration, the Palmer Drought Severity Index (PDSI), the Normalized Difference Vegetation Index (NDVI), the Enhanced Vegetation Index (EVI), Gross Primary Productivity (GPP) and data from water conservancy construction. The results showed that SWA and LOS effectively revealed the development and recovery of droughts. There were two obvious drought periods from 2000 to 2017. In the first period (from August 2003 to June 2007), SWA decreased by 11.81% and LOS shortened by 5 days. They reduced by 21.04% and 9 days respectively in the second period (from September 2009 to June 2014), which indicated that there are more severe droughts in the second period. The SOS during two drought periods delayed by 3~6 days in spring, while the EOS advanced 1~3 days in autumn. All of PDSI, SWA and LOS could reflect the period of droughts in SWC, but the LOS and PDSI were very sensitive to the meteorological events, such as precipitation and temperature, while the SWA performed a more stable reaction to drought and could be a good indicator for the drought periodicity. This made it possible for using SWA in drought forecast because of the strong correlation between SWA and drought. Our results improved the understanding of seasonal responses to extreme droughts in SWC, which will be helpful to the drought monitoring and mitigation for different seasons in this ecologically fragile region