The effect of expansion ratio on the critical Reynolds number in single fracture flow with sudden expansion

Flow in a single fracture (SF) is an important research subject in groundwater hydrology, hydraulic engineering, radioactive nuclear waste repository and geotechnical engineering. An abruptly changing aperture is a unique type of SF. This study discusses the relation between the values of the critical Reynolds number (Re_c) for the onset of symmetry breaking of flow and the expansion ratio (E) of SF, which is defined as the ratio between the outlet (D) and inlet (d) apertures. This study also investigates the effect of inlet aperture d on Re_c for flow in an SF with abruptly changing apertures (SF‐ACA) using the finite volume method. Earlier numerical and experimental results showed that flow is symmetric in respect to the central plane of the SF‐ACA at small Reynolds number (Re) but becomes asymmetric when Re is sufficiently large. Our simulations show that the value of Re_c decreases with the increasing E, and the relationship between the logarithm of Re_c and E can be described accurately using either a quadratic polynomial function or a logarithmic function. However, the relationship of Re_c and d for a given E value is vague, and Re_c becomes even less sensitive to d when E increases. This study also reveals that the hydraulic gradient (J) and flow velocity (v) follow a super‐linear relationship that can be fitted almost perfectly by the Forchheimer equation. The inertial component (J_i) of J increases monotonically with Re, whereas the viscous component (J_v) of J decreases monotonically with Re. The Re value corresponding to equal inertial and viscous components of J (named as the transitional point Re) decreases when E increases, and such a transitional point Re should be closely related to the critical Reynolds number Re_c, although a rigorous theoretical proof is not yet available. Copyright © 2015 John Wiley & Sons, Ltd.     

Assessing the effectiveness of imperviousness on stormwater runoff in micro urban catchments by model simulation

Imperviousness, considered as a critical indicator of the hydrologic impacts of urbanization, has gained increasing attention both in the research field and in practice. However, the effectiveness of imperviousness on rainfall–runoff dynamics has not been fully determined in a fine spatiotemporal scale. In this study, 69 drainage subareas <1 ha of a typical residential catchment in Beijing were selected to evaluate the hydrologic impacts of imperviousness, under a typical storm event with a 3‐year return period. Two metrics, total impervious area (TIA) and effective impervious area (EIA), were identified to represent the impervious characteristics of the selected subareas. Three runoff variables, total runoff depth (TR), peak runoff depth (PR), and lag time (LT), were simulated by using a validated hydrologic model. Regression analyses were developed to explore the quantitative associations between imperviousness and runoff variables. Then, three scenarios were established to test the applicability of the results in considering the different infiltration conditions. Our results showed that runoff variables are significantly related to imperviousness. However, the hydrologic performances of TIA and EIA were scale dependent. Specifically, with finer spatial scale and the condition heavy rainfall, TIA rather than EIA was found to contribute more to TR and PR. EIA tended to have a greater impact on LT and showed a negative relationship. Moreover, the relative significance of TIA and EIA was maintained under the different infiltration conditions. These findings may provide potential implications for landscape and drainage design in urban areas, which help to mitigate the runoff risk. Copyright © 2015 John Wiley & Sons, Ltd.     

Frequency analysis of nonstationary annual maximum flood series using the time‐varying two‐component mixture distributions

The most popular practice for analysing nonstationarity of flood series is to use a fixed single‐type probability distribution incorporated with the time‐varying moments. However, the type of probability distribution could be both complex because of distinct flood populations and time‐varying under changing environments. To allow the investigation of this complex nature, the time‐varying two‐component mixture distributions (TTMD) method is proposed in this study by considering the time variations of not only the moments of its component distributions but also the weighting coefficients. Having identified the existence of mixed flood populations based on circular statistics, the proposed TTMD was applied to model the annual maximum flood series of two stations in the Weihe River basin, with the model parameters calibrated by the meta‐heuristic maximum likelihood method. The performance of TTMD was evaluated by different diagnostic plots and indexes and compared with stationary single‐type distributions, stationary mixture distributions and time‐varying single‐type distributions. The results highlighted the advantages of TTMD with physically‐based covariates for both stations. Besides, the optimal TTMD models were considered to be capable of settling the issue of nonstationarity and capturing the mixed flood populations satisfactorily. Copyright © 2016 John Wiley & Sons, Ltd.     

L1范数与IQR统计量组合的GNSS坐标序列粗差探测算法

GNSS坐标时间序列中不可避免地含有粗差,未剔除的粗差将会导致参数估计有偏。因此,粗差探测与剔除是GNSS坐标序列分析中一项重要的数据预处理工作。针对GNSS坐标时间序列特点,提出了一种将L1范数(L1-norm)估计与四分位距统计量IQR(interquartile range)组合的移动开窗粗差探测算法,称之为L1_Mod IQR。该方法的主要思想是,首先利用L1范数估计得到较"真实"的残差,然后再对残差采用IQR统计量进行粗差探测。将L1_Mod IQR法与"3σ"法、基于最小二乘的τ检验法等粗差探测算法进行了模拟计算与对比,验证了该算法的有效性。进一步采用L1_Mod IQR算法对中国区域10个IGS站的高程时间序列进行了分析,结果表明中国区域IGS站高程序列的粗差剔除率最小为0.1%,最大为2.6%。并且以WUHN站为例与SOPAC提供的结果进行了对比,结果表明SOPAC提供的"Clean"数据仍含有大量的粗差,而L1_Mod IQR算法能够有效地剔除粗差。     

GNSS相应坐标系之间的关系及其影响

不同的GNSS采用的坐标系定义几乎相近,但参考椭球及其坐标实现不同,这将影响多GNSS融合导航定位效果。根据各GNSS坐标系所采用参考椭球的基本常数,计算比较了不同坐标系参考椭球参数的差异;导出了相应的正常重力公式,比较了这些正常重力公式确定的正常重力值差异;最后分别从坐标系统的定义与实现两个方面分析了其对定位结果的影响。结果表明:1)GPS(BDS)与Galileo和GLONASS所使用的参考椭球引起正常重力差约为0.15和0.30 mgal;2)GPS与BDS,Galileo及GLONASS所使用参考椭球引起纬度分量最大差异约为0.1 mm,3 cm和3 cm,高程分量约为0.1 mm,0.5 m和1 m;3)各GNSS所使用坐标框架间转换参数引起的坐标变化达到厘米级。     

川北旺苍唐家河剖面仙女洞组灰泥丘沉积特征及造丘环境分析

目前川北地区早寒武世微生物藻与古杯参与的生物礁(丘)发育特征存在较大争议。通过实测剖面与镜下分析,川北旺苍县鼓城乡唐家河剖面仙女洞组发育三期灰泥丘。对该剖面灰泥丘沉积特征进行精细研究认为:第一期灰泥丘发育在滑塌角砾灰岩之上,造丘生物主要为枝状和房室状附枝藻,主要岩性为砾屑灰岩、藻格架凝块灰岩和含云生物碎屑灰岩,发育丘基、丘核、丘盖和丘翼,后期演化为钙质砂泥岩结束生长。第二期生物丘发育在鲕粒滩间生屑灰岩之上,造丘生物主要为肾状藻,次为附枝藻,主要岩性为生物碎屑灰岩、藻格架凝块灰岩、漂浮砾岩和亮晶鲕粒灰岩,纵向上发育丘基、丘核、丘坪、丘盖和丘翼,最终再度沉积鲕粒灰岩。第三期灰泥丘发育在藻叠层灰岩之上,造丘生物为古杯和微生物藻,岩性主要为古杯灰岩、生物碎屑泥晶灰岩、藻格架凝块灰岩和藻叠层灰岩,发育丘基、丘核和丘盖。进一步分析三期灰泥丘的成丘环境认为第一期灰泥丘属于台缘斜坡环境内的斜坡灰泥丘低能沉积,经历三个演化阶段;第二期属于台地边缘灰泥丘,经历四期演化阶段;第三期为台缘高能鲕粒滩向陆一侧能量相对较低的台缘古杯灰泥丘沉积,经历三个演化阶段。剖面中灰泥丘多与鲕粒滩共存,构成丘滩复合体,为后期研究储集性能和油气富集关系具有重要的油气地质意义。     

内蒙古达茂旗塔林宫地区找矿前景分析

内蒙古达茂旗塔林宫地区位于华北板块北缘与西伯利亚板块南缘的结合带,天山―兴蒙造山系之温都尔庙俯冲增生杂岩带(内蒙古自治区地质矿产局. 1991),白乃庙―哈达庙铜、金、萤石成矿带上(邵和明,张履桥. 2001) 。区域上矿产资源较为丰富,金属矿产主要有铁、稀土、金、铜等,非金属矿产主要有萤石、石灰岩、石英岩等,比较重要的矿床有白云鄂博钒钛稀土矿、查干哈达庙铜矿、乌花敖包金矿、赛乌苏金矿、白乃庙铜金矿、石宝铁矿等。