• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    Water level updating model for flow calculation of river networks

    2014-03-06 06:21:54XiaolingWUXiaohuaXIANGLiLIChuanhaiWANG
    Water Science and Engineering 2014年1期

    Xiao-ling WU*, Xiao-hua XIANG, Li LI, Chuan-hai WANG

    1. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, P. R. China

    2. State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, P. R. China

    3. Zhenjiang Engineering Survey and Design Institute Co., Ltd., Zhenjiang 212003, P. R. China

    The nodal water level in Eq. (9) can be used to update the water level and discharge of the river network connecting to the same node. With the help of Eq. (8), the correction value can be extended from observation stations to any other river cross-sections.

    It is difficult to update the river network state for a complex river network with fewer gauge

    Water level updating model for flow calculation of river networks

    Xiao-ling WU*1,2, Xiao-hua XIANG1,2, Li LI3, Chuan-hai WANG1,2

    1. College of Hydrology and Water Resources, Hohai University, Nanjing 210098, P. R. China

    2. State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering, Hohai University, Nanjing 210098, P. R. China

    3. Zhenjiang Engineering Survey and Design Institute Co., Ltd., Zhenjiang 212003, P. R. China

    Complex water movement and insufficient observation stations are the unfavorable factors in improving the accuracy of flow calculation of river networks. A water level updating model for river networks was set up based on a three-step method at key nodes, and model correction values were collected from gauge stations. To improve the accuracy of water level and discharge forecasts for the entire network, the discrete coefficients of the Saint-Venant equations for river sections were regarded as the media carrying the correction values from observation locations to other cross-sections of the river network system. To examine the applicability, the updating model was applied to flow calculation of an ideal river network and the Chengtong section of the Yangtze River. Comparison of the forecast results with the observed data demonstrates that this updating model can improve the forecast accuracy in both ideal and real river networks.

    plain river network; cyclic looped channel network; water level updating model; hydrodynamic model; error correction

    1 Introduction

    Meant to improve the initial stage of flood routing, real-time updating is a technique combining field measurements with a flood forecasting model (Mu and Zhang 2007; Li and Weng 2001). Coupled with hydrological and hydrodynamic models, the real-time updating technique has been widely applied to solving problems in single-channel and dendritic channel networks (Crissman et al. 1993; Neal et al. 2007; Romanowicz et al. 2006; Wang and Bai 2008). There are many difficulties associated with flood inundation and updating in cyclic looped channel networks, especially in downstream areas. One involves defining the boundary condition and flow direction, which are complicated due to the reciprocating water flow motion under the tidal influence (Hsu et al. 2003). As the updating model is based onmeasurements, the other difficulty is the measurement of both water quality and quantity. Compared with discharge observations, water level observations are commonly conducted in downstream areas of river networks, but few of them are performed over long periods of time (Madsen and Skotner 2005). Fewer and simpler measurements load a heavy burden on the development of an updating model for cyclic looped channel networks (Kachroo 1992; Kachroo and Liang 1992), let alone complex space states.

    To improve the accuracy of water level and discharge forecasts, a dynamic updating model based on the real-time water level observation was developed. With the developed updating model, the correction value was transmitted from a gauge station to other cross-sections nearby (Lai 2009). The applicability of the water level updating model was verified by its application to both ideal and real cyclic looped channel networks.

    2 Flow simulation for cyclic looped channel networks

    According to the connecting structure, river networks can be classified into two categories: dendritic channel networks and cyclic looped channel networks (Wu et al. 2013). Cyclic looped channel networks consist of junctions and branches with looping structures inside, which are further divided into inside and outside channels, with their flow states determined by the boundary conditions and the compatibility conditions of internal junctions (Cheng et al. 2006).

    2.1 Internal river flow calculation

    The basic continuity and momentum equations describing one-dimensional gradually varied flow are as follows:

    where q is the lateral inflow per unit channel length, Q is the discharge, A is the area of a cross-section, B is the channel width, Z is the water level, n is Manning’s roughness coefficient, R is the hydraulic radius, g is the gravitational acceleration, u is the mean flow velocity at the cross-section, and ω is the momentum correction coefficient.

    The numerical solution of Eq. (1) can be obtained when the initial and boundary conditions are available. For a river channel with several cross-sections, the first and last cross-sections were defined as N1and N2. The four-point implicit finite difference scheme (Wang and Bai 2008) was used to discretize Eq. (1) between the ith and the (i+1)th cross-sections with the weight coefficient equal to 1.0. In the discretization processes, the ith cross-section can be organized into two equations. With the three-step method, the discharge at each cross-section can be described as a linear function of water levels at the correspondinglocation and the first and last cross-sections of the river channel:

    where αi, βi, ξi, θi, ηi, and γiare the discrete coefficients of the ith cross-section, ZN1and ZN2are the water level at the first and last cross-sections, Ziis the water level at the ith cross-section, and Qiis the discharge at the ith cross-section.

    Through simultaneous solution of Eq. (2), the water level and discharge of the ith cross-section are as follows:

    2.2 Updating model for river network

    For a river network in a plain area, gauge stations are always located at junctions and cross-sections of river channels. The gauge stations provide the tidal level data, which is the basal information of the updating model, and the updating model should be switched according to different locations of gauge stations. In general, it is difficult to update the calculated value directly using the observed information from gauge stations located at cross-sections of a river channel. The purpose of this study was to develop a real-time water level updating model especially for the condition with gauge stations located inside the river channel.

    There are three steps to updating the river network state, including the water level and discharge of the network. The first step is to analyze the model error at a gauge station and obtain the corresponding correction value, the second step is to update the coefficients α and θ of the cross-section where the gauge station is located, and the third step is to extend the model correction value to all the cross-sections and update all the coefficients of the river network. With the new coefficients, the water level and discharge of all the cross-sections can

    be updated.

    First, the model error at a gauge station, located at the ith cross-section, is collected (Wu et al. 2013):

    where Y?i,kand Yi,kare the observed and calculated values at the ith cross-section at the kth time step, respectively; N is the number of samples; and ΔYiis defined as the model error at the ith cross-section at the tth time step.

    It is assumed that there is a gauge station located at the ith cross-section with the water level correction value of ΔZi. For the coefficient recursive calculation, Eq. (2) becomes

    where αi′ = αi+Δαiand θi′ = θi+Δ θi, with Δαi=βiΔZiand Δθi=ηiΔZi.

    As the recurrence goes on, the correction value of coefficient α at the jth cross-section can be Δαj= λjΔ Zi, where the influence coefficient

    For the updating coefficient θ at the sth cross-sections, the correction value Δθscan be Δθs= μsΔ Zi, where the influence coefficient

    where Cj, Ej, Fj, Gj, and Wjare the discrete coefficients (Wu et al. 2013). With the coefficients α and θ updated by the correction value of the water level at the ith cross-section, the water level and discharge of the ith cross-section defined by Eq. (3) can be replaced as follows:

    2.3 Nodal water level equation

    Based on the linear relationship between the river network discharge and nodal water level, the nodal water level equation can be set up according to the principle of water balance at key nodes:

    where,ijQ is the discharge from the jth river channel connecting to the ith node (m3/s), with inflow defined as positive, and outflow as negative; Aiis the joining area of the ith node (m2), for nodes without regulation and storage capacity, Ai= 0; Ziis the water level at the ith node (m); t is time (s); and m is the total number of river channels connecting to the ith node. If the joining area of the ith node is much smaller, withthe node is considered the absence of the regulation and storage capacity.

    The nodal water level in Eq. (9) can be used to update the water level and discharge of the river network connecting to the same node. With the help of Eq. (8), the correction value can be extended from observation stations to any other river cross-sections.

    3 Application examples

    It is difficult to update the river network state for a complex river network with fewer gauge

    stations in a plain area. Many factors confuse the flow calculation for real river networks, such as typhoons, tides, climate, and human activities. To evaluate the performance of the presented updating model, an ideal river network and a real river network were used as case studies.

    3.1 Application to ideal river network

    An artificial river network consists of cyclic looped channels based on real data was adopted in this study, as shown in Fig. 1.

    The lengths of inside and outside channels were 10 km and 40 km, respectively. For the implicit discretization, the time step could be set from 5 to 60 min. Accounting for the simulation accuracy and computation burden, we used a 10-min time step in this study. The space interval was specified as 1 km.

    Three inner branches are marked as B1, B2, and B3, with the nodes P1, P2, and P3connecting with one another; three outer branches are marked as B4, B5, and B6, with the boundary nodes P4, P5, and P6assumed to be the Jiangyin, Yanglin, and Qinglonggang stations, respectively, near the Yangtze Estuary. The flow at the three gauge stations is significantly influenced by tides. In this study, the boundary conditions were specified with the observed water level from the Jiangyin, Yanglin, and Qinglonggang stations. The hydrological observation data of flood stage from 0:00 on June 28 to 0:00 on August 10, 1998 were selected for evaluating the efficiency of the water level updating model. The flow direction is presented with arrows in Fig. 1. If the computational result of a discharge is negative, the flow direction is opposite to the arrow direction. The transverse profile of the cross-section of the river channels is displayed in Fig. 2. In the river network, Manning’s roughness coefficient, mainly affecting the dynamic routing model, is about 0.02.

    Fig. 1 Structure of artificial river network

    Fig. 2 Profile of river channel cross-section

    In this case, the high-resolution results of the dynamic routing model were taken as measurements in the updating scheme. For this purpose, the originally calculated situation should be changed to shorten the time step and space interval from 10 min and 1 km to 1 min and 50 m, respectively.

    In order to test the efficiency of the updating model in forecasting the river network state without external interference, the middle cross-section H of B1was selected as themeasurement location (Fig. 1), and the calculated results of the water level updating model based on single water level correction values were compared with those obtained without adopting the water level updating model. The calculated results of the water level at inner nodes and discharge at cross-sections with and without adopting the water level updating model during the flood period from 2:00 on June 30 to 14:20 on July 16, 1998 are shown in Tables 1 and 2.

    Table 1 Simulated results of water level at inner nodes with and without adopting updating model

    Table 2 Simulated results of discharge at cross-sections with and without adopting updating model

    From Tables 1 and 2, it can be seen that with water level updating, the values of RMSE of water level are less than 0.016 m, and those of discharge are less than 203.82 m3/s. The values of RMSE of water level with and without water level updating are almost the same at P3, the distance from which to the measurement point is larger than those from P1and P2. For the cross-sections SE2and SF3, the values of RMSE of discharge with water level updating are close to those without water level updating, and the largest value of RMSE (203.82 m3/s) occurs at SF3, indicating that the updating efficiency decreases with the increase of the distance from the cross-section to the measurement point.

    In respect to the water level and discharge characteristics in the flood period, the results of the water level updating model with single water level correction values are much better than those without adopting the updating model. Thus, the water level updating model basedon single water level correction values is effective for water level and discharge forecasts.

    3.2 Application to real river network

    The Chengtong River network is located at the longitude from 120°15′E to 120°57′E and latitude from 31°45′N to 32°05′N from Jiangyin to Xuliujing. The model was applied to the Chengtong River network downstream of the Yangtze River, which consists of river channels with circular structures. Because the study area is beside the Yangtze Estuary, the water flow movement in this cyclic looped channel network is affected by both inflow from upstream and the tidal current from downstream. In the summer and fall seasons, the torrential rain falling in the river basin causes flood discharge over a 1- to 3-h period.

    The four-point implicit finite difference approximation used in discretization of the dynamic wave equations is an unconditionally stable scheme. Hsu et al. (2006) demonstrated that the influences of the time step on dynamic routing are not significant. Accounting for the simulation accuracy and computation burden, we used a 15-min time step in this study.

    In the network, there are only three permanent water level stations, the Jiangyin (JY), Tianshenggang (TSG), and Xuliujing (XLJ) stations, providing a series of long-term observations. The other observation stations, including the Heshanggang (HSG), Rugaogang (RGG), Taiziwei (TZW), Jiulonggang (JLG), Wuganhe (WGH), Yingchuangang (YCG), Qiganhe (QGH), and Huifengmatou (HFMT) stations, shown in Fig. 3, are all temporary ones without long-term water level observations. The roughness coefficients in flood routing were classified into those of the main channel and flood plain by weight coefficients according to the real-time water level. The model was calibrated and verified with the measured water level from August 26, 2004 at 14:00 to September 18, 2004 at 11:50. The values of Manning’s roughness coefficient were calibrated against the observed flood events in August, 2004. The roughness coefficient of the main channel is around 0.025, and the value of the flood plain at both sides of the main channel is around 0.04. The upstream and downstream boundary conditions were specified based on the observed water level at the JY and XLJ stations, respectively.

    Fig. 3 Distribution of observation stations along Chengtong River section

    All measurement information from observation stations was incorporated into the water level updating model. Fig. 4 exhibits the forecasted hydrographs from 15:00 on September 6 to 15:00 on September 7, 2004, at nine observation stations based on 1- to 3-h flood routing using the water level updating model and the results of observation and obtained by the dynamic routing model. Fig. 4 reveals that the hydrographs of forecast results with water level updating are closer to those of observed data, showing a better prediction as compared with those obtained with the original dynamic routing model without water level updating, especially in a short lead-time forecast.

    Fig. 4 Comparison of hydrographs at nine cross-sections obtained using water level updating model with a lead time from 1 to 3 h with those of observations and obtained by dynamic routing model

    Detailed comparisons of the efficiency coefficients (E) and RMSE between the forecast results and observations at nine cross-sections for water level forecast with a lead time from 1 to 3 h are list in Table 3. Table 3 shows that the value of RMSE with water level updating increases with the lead time (T), while the value of E maintains a decreasing trend, demonstrating that the forecast accuracy of the water level updating model decreases with the increase of the lead time, and the effect of real-time updating gradually fades away with the increase of the forecast time.

    Table 3 Model evaluation results of E and RMSE of water level forecast from September 6 to September 7, 2004

    4 Conclusions

    Based on the dynamic wave theory, a water level updating model was developed and applied in cyclic looped channel networks. With the discrete coefficients, the correction values obtained from observation stations can be extended from measurement points to other cross-sections nearby. The model was applied to ideal and real river networks, and some conclusions are as follows:

    (1) The water level updating model shows a good performance in water level and discharge forecasts; its forecast accuracy is higher than the original dynamic routing model without water level updating, especially in a short lead-time forecast.

    (2) The updating efficiency of the model in water level and discharge forecasts decreases with the increase of distance between the forecast position and the measurement point.

    (3) Future research should be focused on systematic analysis of Manning’s roughness coefficient and its effect on flow forecasting.

    Cheng, W. H., Wang, C. H., and Zhu, Y. 2006. Taihu Basin Model, 10, 71-77. Nanjing: Hohai University Press. (in Chinese)

    Crissman, R. D., Chiu, C. L., Yu, W. Z., Mizumura, K., and Corbu, I. 1993. Uncertainties in flow modeling and forecasting for Niagara River. Journal of Hydrualic Engineering, 119(11), 1231-1250. [doi:10.1061/ (ASCE)0733-9429(1993)119:11(1231)].

    Hsu, M. H., Fu, J. C., and Liu, W. C. 2003. Flood routing with real-time stage correction method for flash flood forecasting in the Tanshui River, Taiwan. Journal of Hydrology, 283(1-4), 267-280. [doi:10.1016/S0022-1694(03)00274-9]

    Hsu, M. H., Fu, J. C., and Liu, W. C. 2006. Dynamic routing model with real-time roughness updating forflood forecasting. Journal of Hydraulic Engineering, 132(6), 605-619. [doi:10.1061/(ASCE) 0733-9429(2006)132:6(605)]

    Kachroo, R. K. 1992. River flow forecasting, part 1: A discussion of the principles. Journal of Hydrology, 133(1-2), 1-15. [doi:10.1016/0022-1694(92)90146-M]

    Kachroo, R. K., and Liang, G. C. 1992. River flow forecasting, part 2: Algebraic development of linear modelling techniques. Journal of Hydrology, 133(1-2), 17-40. [doi:10.1016/0022-1694(92)90147-N]

    Lai, X. J. 2009. Real-updating multivariate analysis for unsteady flows with ensemble Kalman filter. Advances in Water Sciences, 20(2), 241-248. (in Chinese)

    Li, Z. J., and Weng, M. H. 2001. Real-time flood forecasting modeling of 1D unsteady channel flow and kalman filter. Journal of Hydrodynamics, 13(1), 64-69.

    Madsen, H., and Skotner, C. 2005. Adaptive state updating in real-time river flow forecasting: A combined filtering and error forecasting procedure. Journal of Hydrology, 308(1-4), 302-312. [doi:10.1016/ j.jhydrol.2004.10.030]

    Mu, J. B., and Zhang, X. F. 2007. Real-time flood forecasting method with 1-D unsteady flow model. Journal of Hydrodynamics, 19(2), 150-154. [doi:10.1016/S1001-6058(07)60041-9]

    Neal, J. C., Atkinson, P. M., and Hutton, C. W. 2007. Flood inundation model updating using an ensemble Kalman filter and spatially distributed measurements. Journal of Hydrology, 336(3-4), 401-415. [doi: 10.1016/j.jhydrol.2007.01.012]

    Romanowicz, R. J., Young, P. C., and Beven, K. J. 2006. Data assimilation and adaptive forecasting of water levels in the river Severn catchment, United Kingdom. Water Resources Research, 42(6), 1-12. [doi:10.1029/2005WR004373]

    Wang, C. H., and Bai, Y. L. 2008. Algorithm for real time correction of stream flow concentration based on Kalman filter. Journal of Hydrologic Engineering, 13(5), 290-296. [doi:10.1061/(ASCE)1084-0699 (2008)13:5(290)]

    Wu, X. L., Xiang, X. H., Wang, C. H., Chen, X., Xu, C. Y., and Yu, Z. B. 2013. Coupled hydraulic and Kalman filter model for real-time correction of flood forecast in the Three Gorges interzone of Yangtze River, China. Journal of Hydrologic Engineering, 18(11), 1416-1425. [doi:10.1061/(ASCE)HE. 1943-5584. 0000473]

    (Edited by Ye SHI)

    ——

    This work was supported by the Major Program of the National Natural Science Foundation of China (Grant No. 51190091), the National Natural Science Foundation of China (Grant No. 51009045), and the Open Research Fund Program of the State Key Laboratory of Water Resources and Hydropower Engineering Science of Wuhan University (Grant No. 2012B094).

    *Corresponding author (e-mail: freebird7237@163.com)

    Received Jul. 15, 2012; accepted Dec. 7, 2012

    亚洲精品成人久久久久久| 99热这里只有精品一区| 欧美性猛交黑人性爽| 国产一区亚洲一区在线观看| 男的添女的下面高潮视频| 99久久久亚洲精品蜜臀av| 永久网站在线| 老师上课跳d突然被开到最大视频| 高清在线视频一区二区三区 | 91aial.com中文字幕在线观看| 亚洲欧美精品自产自拍| 亚洲成人av在线免费| 两性午夜刺激爽爽歪歪视频在线观看| 国产精品不卡视频一区二区| 欧洲精品卡2卡3卡4卡5卡区| 亚洲av不卡在线观看| 美女高潮的动态| 国产69精品久久久久777片| 欧美3d第一页| 精品日产1卡2卡| 99久国产av精品| 嫩草影院新地址| 久久精品影院6| 亚洲精品色激情综合| 国产一区二区亚洲精品在线观看| АⅤ资源中文在线天堂| 日本在线视频免费播放| 午夜精品在线福利| 一级毛片我不卡| 欧美成人一区二区免费高清观看| 六月丁香七月| 精品免费久久久久久久清纯| 欧美日韩乱码在线| 91久久精品国产一区二区成人| 91久久精品国产一区二区成人| 国产精品久久视频播放| 成人无遮挡网站| 在线观看午夜福利视频| 日本av手机在线免费观看| 久久人人爽人人爽人人片va| 精品久久国产蜜桃| 青春草亚洲视频在线观看| 亚洲av一区综合| 午夜亚洲福利在线播放| 亚洲七黄色美女视频| 午夜老司机福利剧场| 国产精品嫩草影院av在线观看| 一本久久精品| 亚洲欧美中文字幕日韩二区| 精品欧美国产一区二区三| 国产成人aa在线观看| 热99re8久久精品国产| 欧美色视频一区免费| 国产中年淑女户外野战色| 少妇熟女欧美另类| 久久99热6这里只有精品| 18禁在线无遮挡免费观看视频| 午夜福利在线观看吧| 欧美xxxx黑人xx丫x性爽| 国产片特级美女逼逼视频| 午夜爱爱视频在线播放| 免费在线观看成人毛片| 精品人妻熟女av久视频| 一本久久中文字幕| 国产精华一区二区三区| 综合色丁香网| 我要搜黄色片| 卡戴珊不雅视频在线播放| 国产精品久久视频播放| 久久久色成人| 久久人人爽人人片av| 午夜免费激情av| 午夜福利视频1000在线观看| 最近的中文字幕免费完整| 亚洲精品456在线播放app| 成人亚洲精品av一区二区| 又爽又黄a免费视频| 十八禁国产超污无遮挡网站| 午夜精品一区二区三区免费看| 老女人水多毛片| 啦啦啦观看免费观看视频高清| 国产伦理片在线播放av一区 | 老熟妇乱子伦视频在线观看| 99视频精品全部免费 在线| 中文字幕av成人在线电影| 最好的美女福利视频网| 18禁黄网站禁片免费观看直播| 亚洲中文字幕日韩| 在线a可以看的网站| 悠悠久久av| 国产精品无大码| 亚洲欧美成人综合另类久久久 | 国产午夜精品久久久久久一区二区三区| 国产不卡一卡二| 日韩三级伦理在线观看| 国产一区二区在线av高清观看| 菩萨蛮人人尽说江南好唐韦庄 | 国产高清不卡午夜福利| 国产精品综合久久久久久久免费| 日本五十路高清| 国产高清激情床上av| 国产高清视频在线观看网站| 久久精品国产自在天天线| 亚洲av男天堂| 一级毛片电影观看 | 97人妻精品一区二区三区麻豆| 日本与韩国留学比较| 国内揄拍国产精品人妻在线| АⅤ资源中文在线天堂| 成人午夜高清在线视频| 最后的刺客免费高清国语| 美女 人体艺术 gogo| 一级av片app| 亚洲欧美日韩高清在线视频| 国产淫片久久久久久久久| 国产69精品久久久久777片| 中文字幕人妻熟人妻熟丝袜美| 精品久久久久久久人妻蜜臀av| 久久午夜福利片| 亚洲精品影视一区二区三区av| 一进一出抽搐动态| 男人狂女人下面高潮的视频| 中文字幕人妻熟人妻熟丝袜美| 日韩人妻高清精品专区| 亚洲欧美日韩无卡精品| 国产成年人精品一区二区| 3wmmmm亚洲av在线观看| 久久久久久伊人网av| 我要看日韩黄色一级片| 午夜福利在线观看吧| 午夜激情欧美在线| 高清毛片免费看| 日韩视频在线欧美| 在线国产一区二区在线| 国产亚洲精品久久久久久毛片| 2022亚洲国产成人精品| 99久久精品热视频| 国产精品一区二区三区四区久久| 深爱激情五月婷婷| 午夜免费激情av| 男人狂女人下面高潮的视频| 亚洲在线观看片| 久久久久久久久久成人| 亚洲经典国产精华液单| 日韩av不卡免费在线播放| 99热网站在线观看| 一级二级三级毛片免费看| 男人和女人高潮做爰伦理| 国产单亲对白刺激| 亚洲欧美日韩高清在线视频| 伦理电影大哥的女人| 日韩精品青青久久久久久| 99久久久亚洲精品蜜臀av| 国产三级在线视频| 亚洲国产欧洲综合997久久,| av在线播放精品| 最近的中文字幕免费完整| 久久精品国产亚洲av天美| 少妇的逼水好多| 欧美日本视频| 亚洲欧美精品自产自拍| 亚洲成人中文字幕在线播放| 亚洲欧洲日产国产| 亚洲精品成人久久久久久| 99视频精品全部免费 在线| 白带黄色成豆腐渣| 亚洲色图av天堂| 边亲边吃奶的免费视频| 国产精品久久久久久精品电影| 日韩强制内射视频| 亚洲欧美日韩高清专用| 国产精品蜜桃在线观看 | 欧美最新免费一区二区三区| 国产精品无大码| 国产高清视频在线观看网站| 欧美一区二区亚洲| 老司机影院成人| 久久精品国产鲁丝片午夜精品| 日日摸夜夜添夜夜添av毛片| 蜜臀久久99精品久久宅男| 热99在线观看视频| 舔av片在线| 欧美色视频一区免费| 成人二区视频| 国产伦精品一区二区三区四那| 人体艺术视频欧美日本| 亚洲精品日韩av片在线观看| 一级毛片电影观看 | 免费av观看视频| 日本一本二区三区精品| 日韩欧美 国产精品| 一区福利在线观看| 中国国产av一级| 中国美女看黄片| 成人亚洲精品av一区二区| 中文资源天堂在线| 日本黄色视频三级网站网址| 国产三级在线视频| 日日撸夜夜添| 91久久精品国产一区二区成人| 国产黄色视频一区二区在线观看 | 中出人妻视频一区二区| 午夜免费男女啪啪视频观看| 日本黄大片高清| 亚洲综合色惰| 精品免费久久久久久久清纯| 51国产日韩欧美| 亚洲av成人精品一区久久| 夫妻性生交免费视频一级片| 亚洲成人中文字幕在线播放| 亚洲无线在线观看| 欧美色视频一区免费| 国产大屁股一区二区在线视频| 少妇熟女aⅴ在线视频| 久久精品91蜜桃| 国产三级在线视频| 中文字幕久久专区| 亚洲丝袜综合中文字幕| av.在线天堂| eeuss影院久久| 精品免费久久久久久久清纯| 日韩亚洲欧美综合| 一级av片app| 色5月婷婷丁香| 欧美一区二区精品小视频在线| 亚洲一级一片aⅴ在线观看| 亚洲自偷自拍三级| 春色校园在线视频观看| 国产av一区在线观看免费| 九草在线视频观看| 深夜a级毛片| 国产精品精品国产色婷婷| 国产麻豆成人av免费视频| 亚洲图色成人| 搞女人的毛片| 亚洲欧美日韩无卡精品| 国产视频首页在线观看| 国产精品久久久久久精品电影小说 | 中国美白少妇内射xxxbb| 日产精品乱码卡一卡2卡三| 亚洲最大成人av| 免费在线观看成人毛片| 亚洲av中文字字幕乱码综合| 边亲边吃奶的免费视频| 99在线人妻在线中文字幕| 卡戴珊不雅视频在线播放| 亚洲成a人片在线一区二区| 在线国产一区二区在线| 22中文网久久字幕| 亚洲内射少妇av| 国产淫片久久久久久久久| 国产精品久久电影中文字幕| 久久亚洲国产成人精品v| av在线观看视频网站免费| 国产视频首页在线观看| 久久人人爽人人片av| 好男人在线观看高清免费视频| 中文字幕制服av| 亚洲最大成人中文| 97人妻精品一区二区三区麻豆| 色综合站精品国产| 夜夜夜夜夜久久久久| 亚洲国产日韩欧美精品在线观看| 亚洲精品久久久久久婷婷小说 | 麻豆成人av视频| 免费看a级黄色片| 一级毛片久久久久久久久女| 性欧美人与动物交配| 不卡视频在线观看欧美| 自拍偷自拍亚洲精品老妇| 久久久精品欧美日韩精品| 熟女人妻精品中文字幕| 又黄又爽又刺激的免费视频.| 一级黄色大片毛片| 国产私拍福利视频在线观看| 午夜福利在线观看吧| 国产又黄又爽又无遮挡在线| 国产成人午夜福利电影在线观看| 一区二区三区高清视频在线| 天堂√8在线中文| 久久久国产成人免费| 啦啦啦韩国在线观看视频| 22中文网久久字幕| 久久人人爽人人爽人人片va| 亚洲精品国产av成人精品| 老女人水多毛片| 一区二区三区高清视频在线| 变态另类丝袜制服| 亚洲欧美中文字幕日韩二区| 久久这里有精品视频免费| 国产三级中文精品| 嫩草影院精品99| 婷婷色av中文字幕| 成人综合一区亚洲| 婷婷精品国产亚洲av| 九九久久精品国产亚洲av麻豆| 少妇猛男粗大的猛烈进出视频 | 啦啦啦观看免费观看视频高清| 在线a可以看的网站| 国产一区二区在线观看日韩| 六月丁香七月| 国产高清三级在线| 午夜福利高清视频| 美女xxoo啪啪120秒动态图| 久久久国产成人精品二区| 欧美变态另类bdsm刘玥| 国产亚洲5aaaaa淫片| 观看美女的网站| 亚洲人成网站在线播| 变态另类成人亚洲欧美熟女| 亚洲av免费在线观看| 国产成人aa在线观看| 亚洲无线观看免费| 国产精品无大码| 亚洲国产精品成人综合色| 日本一本二区三区精品| 日韩一本色道免费dvd| 免费观看在线日韩| www.色视频.com| 亚洲av免费高清在线观看| 日韩一区二区视频免费看| 亚洲成人精品中文字幕电影| 日韩成人av中文字幕在线观看| 国产久久久一区二区三区| 国产伦理片在线播放av一区 | 日本免费a在线| 99久国产av精品| .国产精品久久| 日韩三级伦理在线观看| 在线免费观看不下载黄p国产| 亚洲性久久影院| 日韩制服骚丝袜av| videossex国产| 国产在线精品亚洲第一网站| 男插女下体视频免费在线播放| 亚洲va在线va天堂va国产| 欧美丝袜亚洲另类| 免费观看精品视频网站| 国产精品综合久久久久久久免费| 一边亲一边摸免费视频| 人妻少妇偷人精品九色| 成人高潮视频无遮挡免费网站| 久久99热这里只有精品18| 午夜精品在线福利| 免费av观看视频| 少妇被粗大猛烈的视频| 18禁裸乳无遮挡免费网站照片| 国产色婷婷99| 日本五十路高清| 免费搜索国产男女视频| 日韩视频在线欧美| 国产午夜福利久久久久久| 国产av不卡久久| 美女高潮的动态| 一个人看的www免费观看视频| av又黄又爽大尺度在线免费看 | 美女xxoo啪啪120秒动态图| 日本三级黄在线观看| 99视频精品全部免费 在线| 久久久久久久久久黄片| 午夜爱爱视频在线播放| 国产伦在线观看视频一区| 国产成人91sexporn| 少妇猛男粗大的猛烈进出视频 | 国产av麻豆久久久久久久| 亚洲真实伦在线观看| 99国产精品一区二区蜜桃av| 亚洲人成网站高清观看| 99国产极品粉嫩在线观看| 亚洲国产精品成人综合色| 亚洲人成网站高清观看| 亚洲av中文av极速乱| av视频在线观看入口| av在线老鸭窝| 亚洲精品自拍成人| 91久久精品国产一区二区成人| 免费人成在线观看视频色| 成人三级黄色视频| 日韩高清综合在线| a级一级毛片免费在线观看| 久久人人精品亚洲av| 一边摸一边抽搐一进一小说| 少妇的逼水好多| 亚洲激情五月婷婷啪啪| 日日摸夜夜添夜夜爱| av黄色大香蕉| 夜夜看夜夜爽夜夜摸| 青春草亚洲视频在线观看| 久久精品国产亚洲av天美| 精品一区二区三区人妻视频| 一级毛片电影观看 | a级毛片a级免费在线| 晚上一个人看的免费电影| 日本黄色视频三级网站网址| 国产精品精品国产色婷婷| 午夜福利成人在线免费观看| 国产极品精品免费视频能看的| 国产黄色小视频在线观看| 国产精品麻豆人妻色哟哟久久 | 欧美成人精品欧美一级黄| 亚洲中文字幕一区二区三区有码在线看| 美女xxoo啪啪120秒动态图| 亚洲天堂国产精品一区在线| 日日啪夜夜撸| 能在线免费看毛片的网站| 欧美高清成人免费视频www| 中文欧美无线码| 禁无遮挡网站| 亚洲成av人片在线播放无| 九九久久精品国产亚洲av麻豆| 最近2019中文字幕mv第一页| 国内精品一区二区在线观看| 国产黄a三级三级三级人| 91精品国产九色| 三级国产精品欧美在线观看| 青春草亚洲视频在线观看| 蜜桃久久精品国产亚洲av| 亚洲精品乱码久久久久久按摩| 六月丁香七月| 国产成人午夜福利电影在线观看| 草草在线视频免费看| 51国产日韩欧美| 波野结衣二区三区在线| 青青草视频在线视频观看| 精品久久久噜噜| 草草在线视频免费看| 99久久精品国产国产毛片| 特级一级黄色大片| 亚洲av成人av| 国产精品蜜桃在线观看 | 久久精品国产清高在天天线| 国产麻豆成人av免费视频| 久久久久久久久久成人| 波多野结衣高清作品| 欧美极品一区二区三区四区| 日本免费一区二区三区高清不卡| 99久久无色码亚洲精品果冻| 欧美变态另类bdsm刘玥| av免费在线看不卡| 老司机福利观看| 精品人妻偷拍中文字幕| 天天躁日日操中文字幕| 国产精品不卡视频一区二区| 美女黄网站色视频| 内地一区二区视频在线| 日本-黄色视频高清免费观看| www.色视频.com| 欧美丝袜亚洲另类| 永久网站在线| 男女视频在线观看网站免费| 欧美日本视频| 日日撸夜夜添| 99久国产av精品| 大型黄色视频在线免费观看| 亚洲精品国产av成人精品| 久久精品人妻少妇| 国产不卡一卡二| 热99re8久久精品国产| 亚洲中文字幕日韩| 日日干狠狠操夜夜爽| 中文字幕av在线有码专区| 99热只有精品国产| 91久久精品国产一区二区成人| 国产精品野战在线观看| 99热网站在线观看| 免费看av在线观看网站| 人体艺术视频欧美日本| 日本av手机在线免费观看| 国产伦精品一区二区三区四那| 国产精品一及| 欧美极品一区二区三区四区| 亚洲精品久久久久久婷婷小说 | 深夜精品福利| 精品日产1卡2卡| 日韩欧美 国产精品| 尾随美女入室| 中文精品一卡2卡3卡4更新| 蜜桃久久精品国产亚洲av| 人人妻人人澡欧美一区二区| 国产淫片久久久久久久久| 亚洲成人av在线免费| 美女xxoo啪啪120秒动态图| 日韩大尺度精品在线看网址| 99热网站在线观看| 午夜精品国产一区二区电影 | 免费无遮挡裸体视频| 国产高清视频在线观看网站| 九九爱精品视频在线观看| 91aial.com中文字幕在线观看| 成人av在线播放网站| 亚洲国产高清在线一区二区三| 日韩av在线大香蕉| 亚洲,欧美,日韩| 亚洲最大成人av| 日韩一区二区视频免费看| 色吧在线观看| 成年av动漫网址| 亚洲国产精品sss在线观看| 麻豆精品久久久久久蜜桃| 日韩人妻高清精品专区| 国产亚洲精品久久久com| 国产精品一区二区在线观看99 | 99热这里只有精品一区| 蜜臀久久99精品久久宅男| 国产精品一区www在线观看| 欧美日韩一区二区视频在线观看视频在线 | 亚洲性久久影院| 国产一区二区亚洲精品在线观看| 欧美成人免费av一区二区三区| 99国产极品粉嫩在线观看| 国产老妇伦熟女老妇高清| 精品不卡国产一区二区三区| 亚洲av熟女| 久久久久久久久久成人| 青春草亚洲视频在线观看| 欧美最黄视频在线播放免费| 国产私拍福利视频在线观看| 九九爱精品视频在线观看| 婷婷色综合大香蕉| 亚洲精品456在线播放app| 国产 一区精品| 国产单亲对白刺激| 亚洲精品自拍成人| 美女cb高潮喷水在线观看| 在线天堂最新版资源| 成人无遮挡网站| 国产精品1区2区在线观看.| 成人午夜高清在线视频| 一个人看视频在线观看www免费| 深爱激情五月婷婷| 国产午夜精品一二区理论片| 波多野结衣高清作品| 亚洲av成人精品一区久久| 长腿黑丝高跟| 久久久久久久亚洲中文字幕| 日韩欧美一区二区三区在线观看| 精品99又大又爽又粗少妇毛片| 亚洲精品日韩av片在线观看| 国产精品伦人一区二区| www日本黄色视频网| 成人鲁丝片一二三区免费| 自拍偷自拍亚洲精品老妇| 麻豆精品久久久久久蜜桃| 又粗又爽又猛毛片免费看| 乱系列少妇在线播放| 天天躁夜夜躁狠狠久久av| 国产人妻一区二区三区在| 欧美精品一区二区大全| 亚洲成人久久爱视频| 中文字幕av成人在线电影| 三级毛片av免费| 日韩精品青青久久久久久| 99热全是精品| 99国产极品粉嫩在线观看| 日韩一区二区三区影片| 国产午夜精品久久久久久一区二区三区| 少妇熟女欧美另类| 国产亚洲5aaaaa淫片| 干丝袜人妻中文字幕| 18禁黄网站禁片免费观看直播| 女的被弄到高潮叫床怎么办| 国产精品av视频在线免费观看| 别揉我奶头 嗯啊视频| av免费在线看不卡| 国产色爽女视频免费观看| 日韩三级伦理在线观看| 久久久精品大字幕| 中文精品一卡2卡3卡4更新| 99久久精品一区二区三区| 国国产精品蜜臀av免费| 国产精品99久久久久久久久| 岛国在线免费视频观看| 少妇高潮的动态图| 最近视频中文字幕2019在线8| 国产午夜精品论理片| 亚洲av第一区精品v没综合| www.av在线官网国产| 99riav亚洲国产免费| 寂寞人妻少妇视频99o| 一进一出抽搐动态| 亚洲真实伦在线观看| 欧美成人免费av一区二区三区| 男女下面进入的视频免费午夜| 别揉我奶头 嗯啊视频| 美女黄网站色视频| 99热精品在线国产| 国国产精品蜜臀av免费| 久久这里有精品视频免费| 丰满乱子伦码专区| 尾随美女入室| 两个人视频免费观看高清| 亚洲成人av在线免费| 久久九九热精品免费| 男的添女的下面高潮视频| 在线免费观看不下载黄p国产| 国产极品天堂在线| 欧美3d第一页| 日韩中字成人| 高清在线视频一区二区三区 | 九九在线视频观看精品| 精品欧美国产一区二区三| 大型黄色视频在线免费观看| 一级黄片播放器| 美女大奶头视频| 久久精品国产亚洲av天美| 99热只有精品国产| 小蜜桃在线观看免费完整版高清| 美女xxoo啪啪120秒动态图| 岛国毛片在线播放| 久久午夜福利片| 全区人妻精品视频| 久99久视频精品免费| 亚洲真实伦在线观看| 亚洲av不卡在线观看| 国产精品久久久久久精品电影|