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

    The coupling between hydrodynamic and purification efficiencies of ecological porous spur-dike in field drainage ditch *

    2018-07-06 10:01:34LeiRao饒磊PeifangWang王沛芳QingsongDai戴青松ChaoWang王超
    關(guān)鍵詞:王超青松

    Lei Rao (饒磊), Pei-fang Wang (王沛芳), Qing-song Dai (戴青松), Chao Wang (王超)

    1. College of Mechanics and Materials, Hohai University, Nanjing 211100, China

    2. Key Laboratory of Integrated Regulation and Resource Development on Shallow Lakes, Ministry of Education,Hohai University, Nanjing 210098, China

    3. College of Environment, Hohai University, Nanjing 210098, China

    This work was supported by the National Science Funds for Creative Research Groups of China:Ecological effects of water, sand and contaminants in rivers on water resources of high development (Grant No. 51421006). Our group engages researches closely related with the complex background conditions of“high development of water resources”, especially with the core problem of “transfer processes of water,sand and contaminants, interactive mechanisms and ecological effects”. The researches include the major impacts of the intensive development of water resources on hydrological and sedimentation changes in the middle and lower reaches of the Yangtze River,and the establishment of a hydrological simulation model, applicable to the Yangtze River Basin, as well as the migration laws of pollutants and colloids in the bed subsurface flow, the influential mechanism of the subsurface flow and their own settlement, deposition,adsorption and release behaviors. The mechanism of the flow variation in the migration of organic matter,phosphorus and heavy metals in the “sediment-wateraquatic organism” system is revealed and the biological response and the ecological effects are elucidated.The correlation between the aquatic organism community structure and the variations of water, sediment and pollutants is analyzed, along with the toxic effects of organic pollutants on the plankton and their physiological responses under the influence of water and sediment. The mechanism of the riverine microbe growth and the degradation characteristics and the mechanism of the river hydrodynamic succession process of the microbial community are also studied.

    This paper is mainly related to the research work of “the mechanism of flow variation in the migration of organic matter, phosphorus and heavy metals in the“sediment-water-aquatic organism” system”.

    Introduction

    The field drainage ditch system is a complex water network between farmlands and rivers, which is also the major migratory pathway of various pollutants from farmlands to downstream rivers. The field drainage ditch usually has a small cross-sectional area and a long path length in a widespread distribution and with a low flow velocity. Thus, it is the best region to block and remove the escaping pollutant from the farmland[1]. Exposed in the nature environment for a long time, the field drainage ditch is easy to be destroyed under harsh service conditions. So, most of field drainage ditches are built by cement for its structure stability and durability[2]. Compared with the nature ditch, the cemented slope and bed of this kind of artificial ditch are smooth and waterproof. It not only decreases the specific surface area and biological capacity on its surface, but also brings about a high flow velocity and a short hydraulic retention time(HRT) in the ditch[3]. These negative factors weaken the self-purification capacity of the artificial ditch greatly. So, it is a research target for a long time to enhance the purification capacity of the cemented field drainage ditch system[4-5].

    In recent years, the manmade porous material is gradually used as the road pavement brick to replace the traditional solid waterproof brick in building a sponge city[6-7]. The porous brick has a complexed 3-D porous framework structure, with most holes opened and connected to each other. The porous brick has a great specific surface area with excellent adsorption capacity, which is beneficial to improve the microorganism movement. Meanwhile, the surrounding water can flow into the brick, which increases the effective contact interface area between the polluted water and the biological film greatly. Practice indicates that the purification capacity of the ecological porous brick is much larger than that of the traditional solid brick[8]. Recently, the pollution water purifying technology is being developed by using the ecological porous material, especially the pollution removal technology in the artificial field drainage ditch[9-10].

    The amount of the microorganism attached on the material surface is an important factor for the pollution removal rate (PRR) of the ecological porous brick, which is largely decided by the effective porosity rate and the specific surface area of the porous brick. Meanwhile, the inner flow field of the porous brick and the exchanging water quantity between the porous brick and the surrounding water are closely related with the PRR too[11]. Therefore, a complex coupling process is involved, including the macro flow field in the surrounding water, the micro flow field in the ecological porous brick and the degradation of the biological film. It is well known that the exchanging water quantity and the HRT are the two major parameters to decide the PRR of the porous material. A larger anount of exchanging water can make more pollutant passing through the porous brick region, which is good for the pollutant adsorption. However, the HRT decreases with the increase of the exchanging water amount at the same time. Too short HRT is disadvantageous to the biochemical reaction process. So far, there are many experimental and simulation studies of the flow field in the porous material, but the coupling studies of the hydrodynamic condition and the purification efficiency in the ecological porous materials are still few[12].

    In this paper, a great number of ecological porous spur-dikes are arranged in an experiment channel to simulate the real field drainage ditch. Based on a series of hydrodynamics and degradation experiments,the Ammonia-Nitrogen (NH3-N) degradation evaluation model is coupled with the hydrodynamic model to calculate the PRR of the ecological porous spurdike. The variations of the total pollutant removal rate(TPRR) in the channel are discussed for different porosities ( =P15.73%, 21.84% and 26.25%), gap distances between spur-dikes ( =D0.06 m, 0.09 m and 0.12 m) and inlet flow velocities (=v0.005 m/s,0.020 m/s, 0.030 m/s, 0.040 m/s and 0.050 m/s).

    1. Experimentation design

    1.1 Experimental equipment

    A rectangular section circulation water channel is used to simulate the real field drainage ditch. It is fixed on a steel frame base with a dip angle of 1°along the length direction. The water flow in the channel is driven by a circulating water pump, and the flow rate in the channel can be controlled by the inverter motor driver and the accessory valve system.To eliminate the inlet velocity fluctuation, the water is pumped into an inlet reservoir first and then flows into the water channel through an inlet weir. At the outlet side of the channel, a moveable outlet weir gate is used to control the water level of the channel (of the height of 0.115 m). Twenty ecological porous spurdikes are arranged along the both side of the channel alternatively in the middle region with the gap distance of “D”. The middle section of the channel is selected as the computation domain. The sketch of the experimental channel and the dimension of the computation domain are shown in Fig. 1.

    1.2 Experimental materials

    Fig. 1 The sketches of experimental channel and computational domain (m)

    The dimension of each spur-dike is 0.20 m×0.06 m×0.12 m, which is a combination of two pieces of porous brick. The porous bricks are made of gravel,cement and some kinds of additional material.Different size of gravel can achieve various porosities and specific surface areas of porous bricks. All porous bricks are dealkalized in pure water for 20 d firstly,and then immersed in nitrobacterium and nutrient solution at 10 L/min circular flow rate and under 25°C℃-30°C ambient temperature condition. The nitrobacterium biological film covers on the internal and external surfaces of porous bricks after 20 d. To study the influence of porosity on the inner flow field and the PRR of the ecological spur-dike, three styles of porous brick with different porosities ( =P15.75%,21.84% and 26.25%) are used in the experiment (Fig.2).

    Fig. 2 The ecological porous bricks with different porosities

    1.3 Physical experiments

    To study the influence of the pollutant degradation efficiency of the ecological porous spur-dike under different hydrodynamics conditions, a series of experiments are performed with different inlet flow velocities (v= 0.005 m/s, 0.020 m/s, 0.030 m/s, 0.040 m/s and 0.050 m/s), gap distances ( =D0.06 m, 0.09 m and 0.12 m) and brick porosities ( =P15.75%, 21.84%and 26.25%). For describing the flow velocity distributions in the experimental channel accurately,the flow velocities at 21 sampling points are measured by Acoustic Doppler Velocimeter (ADV), and the water level of each sampling point is also recorded. To deal with the fluctuation during the flow velocity and water level measurement, a statistical average method is used in these cases.

    2. Numerical models

    In this study, each spur-dike is assumed to be a continuous porous medium with possible flows in the region. In addition, to simplify the calculation, the porosity of the porous brick keeps unchanged in the calculation process, and the separation and the regrowth of the biological film on the brick surface are neglected here.

    2.1 Flow field governing equations

    The continuity equation is as follows[13]

    wherexu,yuandzurepresent the superficial velocity components inx,yandzdirections,which can be represented as

    where the subscriptirepresents the components inx,yandzdirections of a respected variable,iuandivare the superficial velocity and physical velocity components inx,yandzdirections,γis the porosity of the media, which is equal to the porosity of the spur-dike within the spur-dike region,and 1 in other regions.

    For describing the water flow in the porous spurdike, an additional source item is required in the Navier-Stokes equations to describe the drag resistance of the porosities of the spur-dike to the water flow. The momentum equations are as follows[14-15]:

    whereρis the density of the water,tis the time,pis the time-averaged pressure,gx,gyandgzare the acceleration of the gravity component,μLis the molecular viscosity coefficient, andμtis the turbulence eddy viscous coefficient of the water.Sx,SyandSzare the additional source item compo-

    nents of the drag resistance of the porous spur-dike,described by the momentum variation due to the viscosity and the inertia of the fluid in the continuous porous media[16-17].

    whereiSare the additional source item components inx,yandzdirections (xS,ySandzS),ais the viscosity drag coefficient in the porous spurdike, and2Cis the inertia drag coefficient. Based on the Darcy’s law, these two parameters can be determined by experiments.

    The porous spur-dikes can produce a turbulent flow in the channel. Here, the -kεequations are used[18-19]:

    wherekis the turbulence kinetic energy,εis the dissipation rate of the turbulence kinetic energy.kσ,are the Prandtl numbers corresponding tokandε, respectively.1Cε,2Cεare empirical constants,Gkis the turbulent kinetic energy produced by the gradient of the average velocity, described as follows:

    Based on thekandε, the turbulence eddy viscous coefficient of the watertμ can be expressed as follows

    whereCμis an empirical constant.

    To determine the water level distribution in the channel, the volume of fluid (VOF) model is used to simulate the motion of the interface between the air and the water. The continuous equation of the volume fraction of the water phase or the air phase can be written as follows

    whereFis the volume function for the water phase and the air phase. Based on theFfunction, the free surface of the water flow can be constructed and traced.

    2.2 Degradation governing equations

    The content of Ammonia-Nitrogen (NH3-N) is an important water quality indicator, which is often used to evaluate the water pollutant state and its selfcleaning capacity. The nitrobacterium of the biological film can absorb and degrade NH3-N in the water effectively. To evaluate the purification efficiency of the ecological porous spur-dike, the NH3-N degradation rate equations of the porous filler are as follows[20]:

    where the subscriptjis the serial number of the subregions in a piece of spur-dike,jηis the NH3-N degradation rate of each subregion,1jC-is the initial NH3-N concentration of the inflow water,jCis the final NH3-N concentration of the outflow water,Kis the total degradation coefficient,asis the total interface area between the water and the porous spur-dike per unit area,Mis the total area of a piece of spur-dike,Ais the sectional area of the channel,usis the average flow velocity in the channel, PRR is the pollutant removal rate of a piece of spur-dike,nis the number of subregions in a piece of spur-dike,Qjis the water flow rate of each subregion in the spur-dike, the subscriptbis the serial number of the spur-dike in the channel,Nis the total number of spur-dikes in the channel, PRRbis the pollutant removal rate of the piece of spur-dike with label “b”,TPRR is the total pollutant removal rate of all spur-dikes in the channel,Ncycleis the water cycle numbers in an experiment period,Qinis the inlet flow rate of the channel,Tis the experimental period,Vchannelis the volume of the channel.

    2.3 Boundary conditions

    The average inlet flow velocity is set as the inlet boundary condition, measured by the ADV. The values ofkandεof the inlet flowcan be calculated using the following empirical equations[13]

    where0his the hydraulic head of the inlet side. In addition, the outlet of the water channel is taken as the outflow boundary, and an atmosphere pressure is set on the interface between the water and the air.Furthermore, the two side walls and the bottom of the water channel are taken as no-slip wall boundaries in the calculation.

    3. Results of simulation and experimentand analyses

    3.1 The flow velocity distribution

    Under the influence of the ecological porous spur-dikes, the flow pattern in the channel is complex,especially in the spur-dikes region. The flow velocity distribution in the channel at the inlet velocity of 0.03 m/s is shown in Fig. 3.

    Figure 3 shows that the water flows around the spur-dikes and an “S” shaped pattern appears. The flow direction changes quickly due to the staggered arrangement of spur-dikes. Meanwhile, because the flow through the area is blocked by the spur-dikes, the flow velocity increases sharply in this region. For clearlydescribing the characteristics ofthe flow pattern in the region of spur-dikes, the velocity and water level distributions along the channel obtained by measuring and simulating are shown in Fig. 4,respectively.

    Fig. 3 The flow field in the channel

    Fig. 4 The flow velocity and water level distributions in channel

    In Fig. 4, the channel is separated into six regions from the inlet to the outlet, marked as “Z1” through“Z6”. The “Z1” is the starting section of the water flow. The flow field in this section is smooth and uniform, and the water level is approximately stabilized. In the “Z2” region, the flow direction changes sharply due to the blocking by the first spurdike. The water flow is compressed into a narrow region between the spur-dike and the channel wall.The local flow velocity is, therefore, increased quickly in this region. The “Z3” is the front section of the spur-dike region. The water flows around the spur-dikes and the velocity increases significantly.Meanwhile, the water level arises due to the quickly increased flow resistance in this region. The “Z4” is the section of the stabilized flow pattern in the middle of the channel, and the water level keeps steady there too. The main stream of the water flow passes through the winding path quickly around the spur-dikes.Meanwhile, a small part of water passes through the porous spur-dikes and then mixes with the main stream repeatedly. The “Z5” is the wake flow influenced region of the spur-dikes. The flow velocity increases slightly in this region due to the decreased water level. The “Z6” is the flow structure recovering region. The water level decreases slightly due to the influence of the tail-weir. The flow pressure releases quickly when the water flows out from the end piece of spur-dike. Then, the flow velocity decreases sharply and its distribution tends to be uniform gradually.

    The degradation efficiency of the biological film is greatly influenced by the micro flow field inside the spur-dike, and it is also influenced by the macro flow field of the surrounding water. To see the relationship between the macro and micro flow fields clearly, the local flow field and the relative pressure field between two spur-dikes are showed in Fig. 5.

    Fig. 5 The flow velocity and the relative pressure distribution between two adjoining spur-dikes

    Based on the characteristics of the flow velocity and the relative pressure distribution between two adjoining spur-dikes, eight regions are identified in Fig. 5, marked as “R1” through “R8”. The “R1” is the spur-dike region. In this region, the water flows into the spur-dike body and the velocity decreases quickly due to the drag effect of the porous medium.Meanwhile, the pressure in the upper left area is higher than in other areas significantly. Therefore, the flow rate in the upper left area is the highest in this spur-dike. The “R2” and “R3” are two slim regions beside the spur-dike boundary, and the purified water flows into these regions from two boundaries after passing through the spur-dike. Here, they are named as the lateral outflow region and the forward outflow region, respectively. In the two regions, although the flow direction is the same, the velocity distribution is highly non-uniform. The “R4” and “R5” are the lateral and forward mixing flow regions between the outflow regions (“R2” and “R3”) and the main flow regions(“R6”, “R7” and “R8”) . In the two regions, the pressure gradient is small, but the velocity distribution is strongly non-uniform. The purified water from the outflow regions and the non-purified water from the main flow regions mix here. There are several eddies generated in these regions, which can enhance the pollutant convection diffusion in the water. The “R6”is the main inflow region. Most water flows into the gap of the spur-dikes through this region. The velocity is much higher than that in the “R4” region. The “R7”is the high pressure permeating region. Due to the blocking by the next spur-dike, the flow direction changes sharply, and the pressure near the spur-dike surface is higher than those in other areas. The larger pressure gradient between the outside and the inside of the spur-dike surface can make more water permeating into the porous spur-dike. The “R8” is the tangential outflow region. In this region, the major flow direction is parallel to the spur-dike surface.Meanwhile, the pressure gradient between this region and the spur-dike body makes a small part of water permeating into the porous spur-dike. Furthermore,the larger velocity difference can enhance the mass convection between “R5” and “R8”, “R4” and “R6”regions, as is good to improve the mixing and dilution effect of the residual pollutant in the water.

    3.2 The coupling analysis of flow field and PRR in spur-dike

    For studying the relationship between the velocity distribution and the local PRR, the micro fluid field inside the porous spur-dike is obtained by the discrete element method. The spur-dike body is divided into 3 layers and 24 subregions, which are named as “A1/B1/C1” through “A8/B8/C8”, respectively. The flow velocity distributions and the subregion locations are shown in Fig. 6.

    The water flows into the spur-dike from the layer“A” first, and the velocity distribution in this layer is quite non-uniform. Figure 6 shows that the velocity in the bordering subregions (“A1”, “A2”, “A5” and “A6”)near the channel wall is significantly higher than that in the center subregions (“A3”, “A4”, “A7” and “A8”).And then, the water goes on passing through the layer“B”. With the dragging effect of porosity, the velocity in the layer “B” decreases quickly and its distribution trends to be uniform. Finally, the water passes through the layer “C” and flows out from the back surface of the spur-dike. Based on the degradation governing equations and the flow field in the spur-dike, the PRR values in each subregion within a single cycle period are calculated, respectively. In this study, a 28 h circulation degradation experiment is executed and the water cycle numbers are considered into the TPRR value of each spur-dike. Figure 7 shows the statistics of the PRR distribution in one piece of spur-dike within 28 hours under the condition ofv=0.03m/s ,D= 0.09 m andP= 21.84%.

    Fig. 6 The flow velocity distribution and the subregions in a spur-dike

    Fig. 7 The PRR distributions of each subregion

    The HRT and the water flow rate in each subregion are two important influencing factors directly related with the PRR of each subregion. Figure 7 shows that the contribution rate in the PRR of each subregion is quite different. The PRR of the layers A,B and C are 0.975%, 1.370% and 1.505%, respectively, and the PRR of this piece of spur-dike reaches 3.85% within 28 h. It means that more pollutant is absorbed on the middle layer and behind the layer of spur-dike. In this experiment, the water circulates in the channel many times within an experimental period.The water cycle numbers are considered in the calculation of the TPRR in an experimental period.Increasing cycle numbers can lead to a higher purification efficiency.

    4. Discussions

    In practice, the pollution purifying capacity of the field drainage ditch is influenced by various factors. The influences of tthe porosity of the spurdike and the gap distance between spur-dikes on the TPRR are discussed, respectively, under different hydrodynamic conditions.

    4.1 The influence of porosity and flow veoci TPRR

    In a contrast with the traditional impermeable material, the abundant inner pores of the ecological porous material are the major pollutant removing place. The porosity of the porous spur-dike is an important parameter, which not only influences the special surface area and the biological capacity, but also influences the exchanging water quantity and the inner flow field significantly. For studying the influence of the porosity on the TPRR under different hydrodynamic conditions, three kinds of spur-dike with different porosities are considered in the experiment and the simulation. The TPRR-time variation curves for different porosities are shown in Fig. 8, and the TPRR values under conditions of three different porosities and various inlet flow velocities are shown in Fig. 9, respectively.

    Figure 8 shows that for the spur-dike with a porosity of 21.84% , a significantly larger TPRR value is obtained than under the conditions of the porosities of 15.73%, 26.25%. Based on the previous macro and micro flow field analysis, it is shown that, with a high porosity spur-dike, smaller drag resistance is achieved than with a low porosity under the same hydrodynamic condition. Hence the exchanging water quantity and the biological capacity of a higher porosity spur-dike (P=26.25%)are larger than those withP=21.84%,P= 15.73%. However, too high porosity also brings about shorter HRT as compared with that under the low porosity conditions, which are not good for improving the purifying capacity of the porous spur-dike. However, although with a lo w porosity spur-dike, a longer HRT can be obtained as compared with one of a high porosity, the exchanging water quantity and the reaction area decrease greatly under the same hydrodynamic conditions. It also limits the purifying efficiency improvement. So, the porosity of the spur-dike should be limited in a reasonable range to achieve a higher TPRR.

    Fig. 8 The TPRR-time variation curves under different porosities ( D =0.09m, v =0.030m/s )

    Fig. 9 The TPRR for different porosities and inlet flow velocities (D =0.09m)

    Moreover, the flowing water condition is more beneficial to improve the purifying capacity of the ecological porous spur-dike than under a hydrostatic condition. The inner flow field and the exchanging water quantity of each spur-dike are influenced by the macro flow field of the channel. Figure 9 shows that the TPRR variations under three porosity conditions are different with the increase of the flow velocity. For a low porosity spur-dike (P=15.73%), the highest TPRR is achieved at a high flow velocity (v=0.040 m/s). However, the highest TPRR values ofP=21.84%,P=26.25%are achieved atv=0.030 m/s,v=0.020 m/s , respectively. Meanwhile,the TPRR will be limited within a lower level under a slow velocity condition (v=0.010 m/s)in the three cases of porosities. It indicates that the flow velocity condition of the field drainage ditch should be considered seriously, and the porosity of the spur-dike should be chosen within a reasonable range under the real flow velocity condition for achieving a high TPRR.

    4.2 The influence of the gap distance and flow velocity on TPRR

    As an important construction technological parameter, the gap distance between spur-dikes not only influences the discharge capacity of the field drainage ditch, but also influences the pollutant purification capacity of the spur-dikes significantly. It is hard to achieve a high TPRR under conditions of too short or too large gap distances. For studying the influence of the gap distance on the TPRR, three gap distances are considered here under the condition of the same delivery volume of porous brick. The TPRR-time variation curves for different gap distances are shown in Fig. 10, and the TPRR values under conditions of three different gap distances and various flow velocities are shown in Fig. 11, respectively.

    Fig. 10 The TPRR-time curves for different gap distances ( P =21.84%, v =0.030m/s )

    Figure 10 shows that the TPRR value forD=0.09 m is significantly higher than those forD=0.06 m,D= 0.12 m. Under the condition of short gap distance (D=0.06 m), the flow velocity decreases in the region of spur-dikes due to a high flow resistance. Although with a slow flow velocity, the HRT can be prolonged in the spur-dike, the low exchanging water quantity also greatly limits the TPRR improvemente. Moreover, a large gap distance(D=0.12 m)can bring about a lower flow resistance than a short gap distance can do. Most water flows around the spur-dikes with a high velocity. It is also hard for the TPRR to reach a satisfactory level due to the low exchanging water quantity and the short HRT in the spur-dike. Hence, the gap distance should be controlled within a reasonable range to reach a balance between the exchanging water quantity and the HRT.

    Fig. 11 The TPRR for different gap distances and inlet flow velocities (P = 21.84%)

    Furthermore, Fig. 11 shows that the TPRR takes various values with different inlet flow velocities under the same condition of the gap distance. Due to the low exchanging water quantity under the condition of low inlet flow velocity (v=0.010 m/s), the TPRR values are limited at a low level (TPRR<50%)under the condition of the three gap distances. While the inlet flow velocity increases to 0.02 Nm/s, the TPRR for the three gap distances increases sharply and reaches more satisfactory values (TPRR>70%).With the inlet flow velocity increasing from 0.03 m/s to 0.05 m/s, the TPRR under the condition of the three gap distances decreases from 70.62% to 54.74%(D=0.06 m), from 78.18% to 68.43% (D=0.09 m)and from 66.77% to 60.11% (D=0.12 m), respectively. Hence, the gap distance should be chosen carefully based on the actual flow velocity condition of the ditch. A reasonable gap distance of the spurdike can bring about a satisfactory purifying efficiency with the same delivery quantity of ecological porous materials.

    5. Conclusions

    Based on the coupling analysis of the flow field and the PRR of the ecological porous spur-dike under various hydrodynamic conditions, some results are obtained, as summarized as follows:

    (1) The ecological porous spur-dike can effectively improve the pollutant removal capacity of the cemented field drainage ditch. Based on the coupling method and the discrete analysis of the inner flow field and the PRR of the ecological porous spur-dike,the TPRR in the field drainage ditch can be evaluated accurately.

    (2) The porosity of the porous spur-dike is an important parameter, which not only influences the special surface area and the biological capacity, but also influences the quantity of the exchanging water and the inner flow field significantly. In practice, the porosity of the spur-dike should be limited in a reasonable range to achieve a high TPRR.

    (3) The gap distance between spur-dikes can influence the discharge capacity and the TPRR of the field drainage ditch greatly. Too short or too large gap distances are not good to achieve a balance between the exchanging water quantity and the HRT in the spur-dike. It should be controlled within a reasonable range to improve the TPRR.

    (4) The porosity and the gap distance of the spur-dike should be chosen carefully based on the actual flow velocity condition of the field drainage ditch. Reasonable parameter matching of the porosity and the gap distance by the flow velocity of the ditch can bring about a satisfactory purifying efficiency with a amall delivery quantity of ecological porous materials.

    [1] Zhang Z., Robin B., Jay B. et al. Nutrient dynamics within drainage ditches under recent, medium, and long-term drainage in the Black soil zone of southeastern Saskatchewan [J].Geoderma, 2017, 289(11): 66-71.

    [2] Tao Y., Wang S., Xu D. et al. Field and numerical experiment of an improved subsurface drainage system in Huaibei plain [J].Agricultural Water Management, 2017,194(9): 24-32.

    [3] Farzana A., John S. G., Nieber J. L. Field infiltration measurements in grassed roadside drainage ditches:Spatial and temporal variability [J].Journal of Hydrology,2015, 530(10): 604-611.

    [4] Zhang S., Liu F., Xiao R. et al. Effects of vegetation on ammonium removal and nitrous oxide emissions from pilot-scale drainage ditches [J].Aquatic Botany, 2016,130(1): 37-44.

    [5] Mathieu N. K., Zhu B., Fidele S. et al. Long-term impact of primary domestic sewage on metal/loid accumulation in drainage ditch sediments, plants and water: Implications for phytoremediation and restoration [J]. Science of the Total Environment, 2017, 581-582(1): 773-781.

    [6] Shao W., Zhang H., Liu J. et al. Data integration and its application in the sponge city construction of China [C].12th International Conference on Hydroinformatics,Incheon, Korea, 2016, 779-786.

    [7] Wang Y., Sun M., Song B. Public perceptions of and willingness to pay for sponge city initiativesin China [J].Resources, Conservation and Recycling, 2017, 122(2):11-20.

    [8] Yuan X., Tang Y., Li Y. et al. Environmental and economic impacts assessment of concrete pavement brick and permeable brick production process-A case study in China[J].Journal of Cleaner Production, 2017,171(10):198-208.

    [9] Zhang S., Li M., Cheng J. et al. Response and environmental assessment of two Chinese conventional carps to water quality regulation in recirculating aquaculture [J].Aquacultural Engineering, 2016, 74(9): 180-185.

    [10] Chen S., Chu Z. Purification efficiency of nitrogen and phosphorus in Ottelia acuminata on four kinds of simulated sewage [J].Ecological Engineering, 2016, 93(5):159-165.

    [11] Long Y., Bing Y., Zhang Z. et al. Influence of plantation on microbial community in porous concrete treating polluted surface water [J].International Biodeterioration and Biodegradation, 2017, 117(11): 8-13.

    [12] Rao L., Qian J., Ao Y. H. Influence of artificial ecological floating beds on river hydraulic characteristics [J].Journal of Hydrodynamics, 2014, 26(3): 474-481.

    [13] Si Salah S. A., Filali E. G., Djellouli S. Numerical investigation of Reynolds number and scaling effects in microchannels flows [J].Journal of Hydrodynamics, 2017,29(4): 647-658.

    [14] Xu Z. H., Ma G. W., Li S. C. A graph theoretic pipe network method for water flow simulation in a porous medium: GPNM [J].International Journal of Heat and Fluid Flow, 2014, 45(2): 81-97.

    [15] Alves B., Barletta A., Hirata S. et al. Effects of viscous dissipation on the convective instability of viscoelastic mixed convection flows in porous media [J].International Journal of Heat and Mass Transfer, 2014, 70(3): 586-598.

    [16] Zhan J. M., Hu W. Q., Cai W. H. et al. Numerical simulation of flow through circular array of cylinders using porous media approach with non-constant local intertial resistance coefficient [J].Journal of Hydrodynamics, 2017,29(1): 168-171.

    [17] Ku X., Lin J., Annaland M. V. S. et al. Numerical simulation of the accumulation of heavy particles in a circular bounded vortex flow [J].International Journal of Multiphase Flow, 2016, 87(8): 80-89.

    [18] Aghaie M., Gheisari F., Zolfaghari A. Investigation of flapper valve malfunction in Tehran research reactor with porous media model [J].Progress in Nuclear Energy,2018, 105(5): 221-235.

    [19] Tang M., Lu S., Zhan H. et al. The effect of a microscale fracture on dynamic capillary pressure of two phase flow in porous media [J].Advances in Water Resources, 2018,113(5): 272-284.

    [20] Tian W. J. The direct biological strengthen purification mechanism and experimental study on weak pollution river [D]. Doctoral Thesis, Nanjing, China: Hohai University, 2005(in Chinese).

    猜你喜歡
    王超青松
    王超美術(shù)作品
    Experimental Study on Grinding Force of Electrostatic Coated Grinding Wheel
    豐收歌
    青松、紅帽與圣誕老人
    峨眉武術(shù)傳承人王超的堅(jiān)守
    林青松
    寶藏(2018年6期)2018-07-10 02:26:44
    惟有骨氣映青松
    寶藏(2017年7期)2017-08-09 08:15:19
    延伸小游戲
    Three-Dimensional Planning of Arrival and Departure Route Network Based on Improved Ant-Colony Algorithm
    Asymptotic Behavior of the Drift Coefficient Estimator of Stochastic Differential Equations Driven by Small Noises
    免费在线观看成人毛片| 亚洲成人av在线免费| 亚洲一区二区三区欧美精品| 一边亲一边摸免费视频| a 毛片基地| 高清在线视频一区二区三区| 色94色欧美一区二区| 精华霜和精华液先用哪个| 午夜久久久在线观看| 大香蕉97超碰在线| 乱码一卡2卡4卡精品| 日韩,欧美,国产一区二区三区| freevideosex欧美| 国国产精品蜜臀av免费| 国产片特级美女逼逼视频| 国产亚洲av片在线观看秒播厂| 久久鲁丝午夜福利片| 国产精品偷伦视频观看了| 精品一区在线观看国产| 97精品久久久久久久久久精品| 亚洲三级黄色毛片| 伊人亚洲综合成人网| 日日撸夜夜添| 一本一本综合久久| 亚洲丝袜综合中文字幕| 久久人妻熟女aⅴ| 精品久久久久久电影网| 春色校园在线视频观看| 国产日韩欧美在线精品| 中文字幕人妻丝袜制服| 国产国拍精品亚洲av在线观看| 婷婷色综合www| 国产日韩欧美亚洲二区| 美女脱内裤让男人舔精品视频| 国内揄拍国产精品人妻在线| 一级毛片我不卡| 久热久热在线精品观看| 国产精品一区二区在线观看99| 国产欧美日韩综合在线一区二区 | 插阴视频在线观看视频| 在线观看一区二区三区激情| 三级国产精品欧美在线观看| 97在线视频观看| 国产精品久久久久久精品电影小说| 国产精品麻豆人妻色哟哟久久| 最近中文字幕2019免费版| 国产精品秋霞免费鲁丝片| 黄色毛片三级朝国网站 | 国产精品人妻久久久影院| 精品久久久精品久久久| 国产淫语在线视频| 交换朋友夫妻互换小说| 久久ye,这里只有精品| 少妇人妻精品综合一区二区| 久久av网站| 日本与韩国留学比较| 成人18禁高潮啪啪吃奶动态图 | 午夜精品国产一区二区电影| 99精国产麻豆久久婷婷| 一级毛片 在线播放| 国模一区二区三区四区视频| 精品久久久噜噜| 午夜福利在线观看免费完整高清在| av免费在线看不卡| 我要看黄色一级片免费的| 一个人免费看片子| 久久免费观看电影| 色视频在线一区二区三区| 男人狂女人下面高潮的视频| 国产白丝娇喘喷水9色精品| 最新中文字幕久久久久| 性色avwww在线观看| 涩涩av久久男人的天堂| 亚洲av欧美aⅴ国产| 国产亚洲午夜精品一区二区久久| 国产乱来视频区| 久久亚洲国产成人精品v| 蜜臀久久99精品久久宅男| 久久影院123| 韩国高清视频一区二区三区| 看免费成人av毛片| 十分钟在线观看高清视频www | 国产精品一区二区三区四区免费观看| 色吧在线观看| 最近2019中文字幕mv第一页| 免费大片黄手机在线观看| av福利片在线| 春色校园在线视频观看| 97超碰精品成人国产| 亚洲精品,欧美精品| 校园人妻丝袜中文字幕| 国产日韩欧美亚洲二区| 中文天堂在线官网| 99精国产麻豆久久婷婷| 欧美精品国产亚洲| 麻豆乱淫一区二区| 不卡视频在线观看欧美| 免费播放大片免费观看视频在线观看| 99久久综合免费| 精品人妻一区二区三区麻豆| 久久久亚洲精品成人影院| 午夜激情久久久久久久| 热re99久久精品国产66热6| 精品少妇黑人巨大在线播放| 亚洲伊人久久精品综合| 国产 精品1| 亚洲精品一二三| 日本av免费视频播放| av在线播放精品| 香蕉精品网在线| 三级经典国产精品| 免费久久久久久久精品成人欧美视频 | 成人午夜精彩视频在线观看| 日日啪夜夜撸| 99久久综合免费| 亚洲国产毛片av蜜桃av| 中文天堂在线官网| 视频中文字幕在线观看| 极品少妇高潮喷水抽搐| 男女啪啪激烈高潮av片| av卡一久久| 精品亚洲成国产av| 久久99一区二区三区| 成人特级av手机在线观看| 免费看av在线观看网站| 一级毛片黄色毛片免费观看视频| 欧美另类一区| 国产无遮挡羞羞视频在线观看| 国产在线免费精品| 一级片'在线观看视频| 亚洲精品日本国产第一区| 午夜视频国产福利| 观看av在线不卡| 大片免费播放器 马上看| 亚洲第一av免费看| 如日韩欧美国产精品一区二区三区 | 国产亚洲午夜精品一区二区久久| 特大巨黑吊av在线直播| 久久人妻熟女aⅴ| 激情五月婷婷亚洲| 国产欧美日韩一区二区三区在线 | 中国三级夫妇交换| 免费在线观看成人毛片| 久久久国产精品麻豆| 亚州av有码| 最近的中文字幕免费完整| 99re6热这里在线精品视频| 亚洲伊人久久精品综合| 国产乱来视频区| 国产白丝娇喘喷水9色精品| 国产一区亚洲一区在线观看| 国产精品久久久久久精品古装| 成人亚洲精品一区在线观看| av一本久久久久| 如何舔出高潮| videos熟女内射| av有码第一页| 久久人人爽人人爽人人片va| 国产亚洲一区二区精品| av.在线天堂| 国国产精品蜜臀av免费| 亚洲欧美日韩另类电影网站| 国产男女超爽视频在线观看| 亚洲av电影在线观看一区二区三区| 久久久久久久大尺度免费视频| 欧美日韩一区二区视频在线观看视频在线| 观看美女的网站| 精品亚洲成a人片在线观看| 日韩av在线免费看完整版不卡| 最近中文字幕高清免费大全6| 极品少妇高潮喷水抽搐| 亚洲精品久久久久久婷婷小说| 卡戴珊不雅视频在线播放| 夜夜爽夜夜爽视频| 中文欧美无线码| 最后的刺客免费高清国语| 99久久精品一区二区三区| 亚洲天堂av无毛| 一区二区三区乱码不卡18| 赤兔流量卡办理| 少妇丰满av| 一本一本综合久久| 内射极品少妇av片p| 七月丁香在线播放| 丝袜在线中文字幕| a级毛片在线看网站| 免费观看性生交大片5| av在线老鸭窝| 3wmmmm亚洲av在线观看| 日韩 亚洲 欧美在线| 黄色配什么色好看| 国产国拍精品亚洲av在线观看| 丝袜在线中文字幕| 成人毛片60女人毛片免费| 视频中文字幕在线观看| 国产伦精品一区二区三区四那| 日韩中字成人| 国产精品99久久99久久久不卡 | 国内揄拍国产精品人妻在线| 国产乱来视频区| 大陆偷拍与自拍| 成人特级av手机在线观看| 黑人巨大精品欧美一区二区蜜桃 | 全区人妻精品视频| 亚洲av电影在线观看一区二区三区| 国产乱人偷精品视频| 亚洲精品亚洲一区二区| 91精品伊人久久大香线蕉| 在线观看免费视频网站a站| av福利片在线| 午夜激情福利司机影院| 亚洲第一区二区三区不卡| 欧美性感艳星| 两个人免费观看高清视频 | 黄色欧美视频在线观看| 成人黄色视频免费在线看| 欧美 亚洲 国产 日韩一| 亚洲无线观看免费| 国国产精品蜜臀av免费| 亚洲国产毛片av蜜桃av| 99热国产这里只有精品6| 乱系列少妇在线播放| 成人亚洲精品一区在线观看| 久久久久久久精品精品| 欧美变态另类bdsm刘玥| 国产成人免费观看mmmm| 国产精品久久久久久久久免| 蜜臀久久99精品久久宅男| 国产在视频线精品| 国产精品.久久久| 激情五月婷婷亚洲| 日本91视频免费播放| 久久99精品国语久久久| 国产在线男女| 久久久久久伊人网av| 人妻系列 视频| 全区人妻精品视频| 99久久综合免费| 久久亚洲国产成人精品v| 免费不卡的大黄色大毛片视频在线观看| 欧美3d第一页| 男人舔奶头视频| 婷婷色综合www| 日韩,欧美,国产一区二区三区| 免费久久久久久久精品成人欧美视频 | 国产精品国产三级国产av玫瑰| 另类精品久久| 大陆偷拍与自拍| 成人18禁高潮啪啪吃奶动态图 | 人妻一区二区av| xxx大片免费视频| 国产欧美日韩一区二区三区在线 | 欧美另类一区| 99热这里只有是精品50| 国产成人免费观看mmmm| 九九爱精品视频在线观看| 欧美亚洲 丝袜 人妻 在线| 亚洲欧美清纯卡通| 夫妻性生交免费视频一级片| 日日摸夜夜添夜夜爱| 国产精品国产av在线观看| 一级毛片 在线播放| 极品教师在线视频| 国产成人91sexporn| 狂野欧美白嫩少妇大欣赏| 一级毛片黄色毛片免费观看视频| 国精品久久久久久国模美| 精品人妻一区二区三区麻豆| 人妻夜夜爽99麻豆av| 日本午夜av视频| 亚洲av不卡在线观看| 亚洲久久久国产精品| a级毛片在线看网站| 一二三四中文在线观看免费高清| 精品亚洲成a人片在线观看| 午夜福利在线观看免费完整高清在| 亚洲内射少妇av| 亚洲国产欧美日韩在线播放 | 亚洲图色成人| 91午夜精品亚洲一区二区三区| 少妇裸体淫交视频免费看高清| 极品教师在线视频| 日本-黄色视频高清免费观看| 18禁在线播放成人免费| 欧美精品高潮呻吟av久久| 国内少妇人妻偷人精品xxx网站| 91精品国产国语对白视频| 中文字幕人妻熟人妻熟丝袜美| 午夜福利在线观看免费完整高清在| 成人影院久久| 久久精品国产亚洲av天美| 你懂的网址亚洲精品在线观看| 九色成人免费人妻av| 日韩一区二区三区影片| 岛国毛片在线播放| 18禁裸乳无遮挡动漫免费视频| 狂野欧美激情性bbbbbb| 国产成人精品福利久久| 免费久久久久久久精品成人欧美视频 | 午夜91福利影院| 乱人伦中国视频| 熟女电影av网| 国产成人精品婷婷| 国产乱来视频区| 国产成人精品福利久久| 亚洲国产精品专区欧美| videossex国产| 成人午夜精彩视频在线观看| 只有这里有精品99| 日本爱情动作片www.在线观看| 高清在线视频一区二区三区| 免费av不卡在线播放| 丰满人妻一区二区三区视频av| 男人爽女人下面视频在线观看| a级片在线免费高清观看视频| av免费观看日本| 亚洲国产最新在线播放| 国产免费又黄又爽又色| 日产精品乱码卡一卡2卡三| 久久午夜福利片| www.av在线官网国产| 美女xxoo啪啪120秒动态图| 精品一品国产午夜福利视频| 国产精品人妻久久久影院| 91成人精品电影| 亚洲av成人精品一二三区| 91精品国产九色| 国产精品成人在线| 日本黄色片子视频| 在线观看免费日韩欧美大片 | 亚洲国产精品999| 777米奇影视久久| 乱码一卡2卡4卡精品| 国产一区二区在线观看av| 中文字幕亚洲精品专区| 国产精品国产三级专区第一集| 亚洲婷婷狠狠爱综合网| 精品一品国产午夜福利视频| 一区二区三区精品91| 久久ye,这里只有精品| 亚洲人成网站在线播| 日韩欧美 国产精品| 一区二区三区四区激情视频| tube8黄色片| 91aial.com中文字幕在线观看| 男人和女人高潮做爰伦理| 99re6热这里在线精品视频| 免费av中文字幕在线| 国产精品蜜桃在线观看| 免费观看av网站的网址| 国产一区亚洲一区在线观看| a级毛色黄片| 99视频精品全部免费 在线| 精品午夜福利在线看| 丰满饥渴人妻一区二区三| 久久女婷五月综合色啪小说| 亚洲精品色激情综合| 免费看av在线观看网站| 国产精品一区二区在线观看99| 大片电影免费在线观看免费| 热re99久久精品国产66热6| 最近中文字幕2019免费版| 国产视频首页在线观看| 纯流量卡能插随身wifi吗| 两个人的视频大全免费| 91在线精品国自产拍蜜月| 国产一区有黄有色的免费视频| 街头女战士在线观看网站| 国产精品久久久久成人av| 精品少妇黑人巨大在线播放| 狂野欧美激情性xxxx在线观看| 久久久欧美国产精品| 成人特级av手机在线观看| 国产精品国产av在线观看| 国产视频首页在线观看| 伦精品一区二区三区| 色视频在线一区二区三区| 不卡视频在线观看欧美| 热re99久久国产66热| 国产白丝娇喘喷水9色精品| 毛片一级片免费看久久久久| 日日啪夜夜爽| 草草在线视频免费看| 内地一区二区视频在线| 日日撸夜夜添| 成人综合一区亚洲| 亚洲精品国产成人久久av| 18+在线观看网站| 国产视频首页在线观看| 亚洲在久久综合| 啦啦啦中文免费视频观看日本| 黄色一级大片看看| 六月丁香七月| 高清午夜精品一区二区三区| 曰老女人黄片| 最新的欧美精品一区二区| 国产亚洲一区二区精品| 26uuu在线亚洲综合色| 成人二区视频| 少妇人妻 视频| 国产乱人偷精品视频| 丰满人妻一区二区三区视频av| 国产日韩欧美在线精品| 国国产精品蜜臀av免费| 伦精品一区二区三区| 啦啦啦中文免费视频观看日本| 女人久久www免费人成看片| 美女福利国产在线| 全区人妻精品视频| 黄色视频在线播放观看不卡| 免费人妻精品一区二区三区视频| 纯流量卡能插随身wifi吗| 国产伦精品一区二区三区视频9| 精品人妻一区二区三区麻豆| 日韩av免费高清视频| 久久久精品免费免费高清| 成人18禁高潮啪啪吃奶动态图 | 内射极品少妇av片p| 啦啦啦视频在线资源免费观看| av福利片在线观看| 春色校园在线视频观看| a级毛片在线看网站| 国产精品一二三区在线看| 2022亚洲国产成人精品| 极品教师在线视频| 七月丁香在线播放| 日韩欧美一区视频在线观看 | 三级经典国产精品| 国内精品宾馆在线| 久久久午夜欧美精品| 超碰97精品在线观看| 日韩精品有码人妻一区| 国产成人精品一,二区| 国产精品人妻久久久影院| 99久久精品一区二区三区| 亚洲国产精品一区三区| 五月开心婷婷网| videossex国产| 亚洲av日韩在线播放| 少妇人妻精品综合一区二区| 日韩中字成人| 又黄又爽又刺激的免费视频.| 久久久国产一区二区| 91精品一卡2卡3卡4卡| 黄色怎么调成土黄色| 国产精品一区二区在线观看99| 亚洲四区av| 91午夜精品亚洲一区二区三区| 亚洲成人一二三区av| 国内揄拍国产精品人妻在线| 亚洲第一av免费看| 亚洲国产精品一区三区| 黄色日韩在线| 国产淫片久久久久久久久| 只有这里有精品99| 99久久人妻综合| 中文精品一卡2卡3卡4更新| 国产精品福利在线免费观看| 国产精品一区www在线观看| 亚洲精品,欧美精品| 晚上一个人看的免费电影| 久久久久久久久大av| 日韩精品免费视频一区二区三区 | 国产极品天堂在线| 国产黄色视频一区二区在线观看| 午夜福利网站1000一区二区三区| 午夜91福利影院| 精品少妇黑人巨大在线播放| 亚洲精品国产色婷婷电影| 国产色婷婷99| 亚洲一级一片aⅴ在线观看| av线在线观看网站| 婷婷色综合大香蕉| 亚洲av综合色区一区| 汤姆久久久久久久影院中文字幕| 曰老女人黄片| 久久99一区二区三区| 美女cb高潮喷水在线观看| 大陆偷拍与自拍| 亚洲成人手机| 女性生殖器流出的白浆| 国产av国产精品国产| 国产在线视频一区二区| 视频中文字幕在线观看| 色5月婷婷丁香| 亚洲av电影在线观看一区二区三区| a级一级毛片免费在线观看| 高清av免费在线| 在线观看一区二区三区激情| 我的女老师完整版在线观看| 亚洲国产精品国产精品| 草草在线视频免费看| 精品99又大又爽又粗少妇毛片| 成人黄色视频免费在线看| 丰满饥渴人妻一区二区三| av在线观看视频网站免费| 国产中年淑女户外野战色| 91午夜精品亚洲一区二区三区| 黄色视频在线播放观看不卡| 亚洲av欧美aⅴ国产| 日韩 亚洲 欧美在线| 两个人的视频大全免费| 老司机影院成人| 人妻夜夜爽99麻豆av| 亚洲精品日本国产第一区| 高清黄色对白视频在线免费看 | 亚洲国产精品一区二区三区在线| 最近手机中文字幕大全| av.在线天堂| 老司机影院成人| 精品熟女少妇av免费看| videossex国产| 欧美最新免费一区二区三区| 热re99久久国产66热| 熟女av电影| 久久午夜综合久久蜜桃| 一个人看视频在线观看www免费| 久久97久久精品| 女人久久www免费人成看片| 王馨瑶露胸无遮挡在线观看| 人体艺术视频欧美日本| 国产免费视频播放在线视频| .国产精品久久| 久久午夜福利片| 国产高清不卡午夜福利| 美女福利国产在线| 女性生殖器流出的白浆| 精品少妇内射三级| 国产男女超爽视频在线观看| 卡戴珊不雅视频在线播放| 欧美日韩综合久久久久久| 婷婷色综合www| 亚洲精品成人av观看孕妇| 性高湖久久久久久久久免费观看| 观看美女的网站| 青春草视频在线免费观看| 街头女战士在线观看网站| 亚州av有码| 久久久久久人妻| 少妇被粗大猛烈的视频| 妹子高潮喷水视频| 国产黄频视频在线观看| 亚洲精品第二区| 亚洲一区二区三区欧美精品| 下体分泌物呈黄色| 婷婷色麻豆天堂久久| 丰满饥渴人妻一区二区三| 久久国产乱子免费精品| 91精品国产九色| 亚洲欧洲国产日韩| 新久久久久国产一级毛片| 各种免费的搞黄视频| 日韩视频在线欧美| 国产精品欧美亚洲77777| 熟女av电影| 日韩成人av中文字幕在线观看| 一个人看视频在线观看www免费| av在线播放精品| 久久久久久久久大av| 久久久国产欧美日韩av| 大片电影免费在线观看免费| 久久狼人影院| 伦理电影免费视频| 内射极品少妇av片p| 91精品伊人久久大香线蕉| 国产av码专区亚洲av| 男人爽女人下面视频在线观看| 国产成人精品无人区| 水蜜桃什么品种好| 国产精品偷伦视频观看了| 成年av动漫网址| 亚洲精品乱久久久久久| 国产伦精品一区二区三区四那| 人人澡人人妻人| .国产精品久久| 国产淫片久久久久久久久| 欧美精品一区二区免费开放| 国产乱人偷精品视频| 免费av中文字幕在线| 欧美日本中文国产一区发布| 日韩免费高清中文字幕av| 男女边摸边吃奶| 特大巨黑吊av在线直播| 国产精品国产三级国产av玫瑰| 一级毛片aaaaaa免费看小| 美女xxoo啪啪120秒动态图| 欧美最新免费一区二区三区| 欧美 亚洲 国产 日韩一| 国产极品天堂在线| 在线亚洲精品国产二区图片欧美 | 日本wwww免费看| 久久婷婷青草| 国产高清有码在线观看视频| 夫妻性生交免费视频一级片| 男人爽女人下面视频在线观看| 自拍偷自拍亚洲精品老妇| 高清在线视频一区二区三区| 欧美国产精品一级二级三级 | 美女福利国产在线| 人妻夜夜爽99麻豆av| 性色av一级| 99久久精品热视频| h视频一区二区三区| 国产精品女同一区二区软件| 国产男女超爽视频在线观看| 国产精品熟女久久久久浪| 亚洲精品乱码久久久久久按摩| 亚洲成人一二三区av| 搡女人真爽免费视频火全软件| 亚洲av成人精品一区久久| 欧美少妇被猛烈插入视频| 欧美激情国产日韩精品一区| 午夜老司机福利剧场| 高清黄色对白视频在线免费看 | 日韩av不卡免费在线播放|