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

    濱?;?、核電廠溫排放數(shù)值模擬研究進(jìn)展

    2010-08-15 00:51:44張長寬陶建峰
    水利水電科技進(jìn)展 2010年3期
    關(guān)鍵詞:排放口濱海浮力

    張長寬,姚 靜,陶建峰

    (1.河海大學(xué)水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098;2.河海大學(xué)港口海岸與近海工程學(xué)院,江蘇南京 210098)

    濱?;?、核電廠溫排放數(shù)值模擬研究進(jìn)展

    張長寬1,2,姚 靜1,2,陶建峰1,2

    (1.河海大學(xué)水文水資源與水利工程科學(xué)國家重點(diǎn)實(shí)驗(yàn)室,江蘇南京 210098;2.河海大學(xué)港口海岸與近海工程學(xué)院,江蘇南京 210098)

    回顧近年來國內(nèi)外濱海火、核電廠溫排放數(shù)值模擬研究的進(jìn)展,著重討論現(xiàn)有溫排放數(shù)值計(jì)算模型的發(fā)展以及存在的問題,主要包括溫排放系統(tǒng)模式的選取和模型邊界控制等。認(rèn)為三維數(shù)值模型在一定程度上代表溫排放數(shù)值模擬的發(fā)展方向,并應(yīng)考慮斜壓效應(yīng),且將水動(dòng)力場與物質(zhì)擴(kuò)散場耦合求解。

    水環(huán)境;溫排放;熱污染;數(shù)值模擬

    隨著國民經(jīng)濟(jì)的發(fā)展,我國的火、核電力建設(shè)發(fā)展迅速。根據(jù)工藝流程的要求,火、核電廠必須連續(xù)不斷大量供應(yīng)機(jī)組冷卻水,同時(shí)向環(huán)境水體(湖泊、水庫、河流、海洋等天然水域)排放等量的超溫?zé)崴?即溫排放,超溫?zé)崴Q之為溫排水[1]。溫排水排入受納水體后,經(jīng)過摻混和擴(kuò)散,使得一定范圍內(nèi)的受納水體水溫高于自然水溫2~14℃,改變了原有水體的物理、化學(xué)性質(zhì)以及生態(tài)環(huán)境,使水生生態(tài)系統(tǒng)平衡失調(diào),水生生物的生存受到威脅,水產(chǎn)資源受到破壞。另一方面,就火、核電廠自身而言,從天然水體中抽取的冷卻水溫度直接關(guān)系到機(jī)組的效率,炎熱季節(jié)冷卻水溫度每升高2℃,機(jī)組效率降低1%,當(dāng)水溫超過一定限度時(shí),還會(huì)形成水循環(huán)短路,影響發(fā)電機(jī)組的安全。因此,對(duì)火、核電廠溫排水在受納水域中的時(shí)空分布規(guī)律進(jìn)行深入了解非常必要,進(jìn)行溫排水?dāng)U散研究不僅具有重要的學(xué)術(shù)價(jià)值,而且對(duì)防止受納水體熱污染、生態(tài)環(huán)境保護(hù)、海洋環(huán)境預(yù)測以及電廠選址和方案比選等都具有重要的實(shí)際意義。

    影響濱?;?、核電廠溫排水時(shí)空分布的因素包括:溫排水流量、流速、初始溫度以及排放口附近海域動(dòng)力條件等。為方便起見,根據(jù)溫排水的物理特性和運(yùn)動(dòng)機(jī)理,人們從環(huán)境水力學(xué)角度將受納水域分為近區(qū)和遠(yuǎn)區(qū)。溫排水排放口附近,受出流初始動(dòng)量的影響,溫排水與環(huán)境水體之間發(fā)生卷吸和劇烈摻混,同時(shí)受浮力作用,溫排水上浮。一般把受紊動(dòng)及浮力雙重控制的區(qū)域稱為近區(qū)。近區(qū)之外,隨著溫排水沿程擴(kuò)散和向大氣散熱,初始動(dòng)量和浮力衰減殆盡,以對(duì)流、擴(kuò)散作用為主,不受動(dòng)量、浮力影響,此區(qū)域稱為遠(yuǎn)區(qū)。在濱海水域,水下地形變化劇烈,水體流動(dòng)主要表現(xiàn)為非恒定性、非均勻性、斜壓性及其三維特征,其運(yùn)動(dòng)尺度在時(shí)空上相差甚大,決定了水動(dòng)力特征的區(qū)域異常可變性和隨機(jī)性。電廠溫排水使近岸天然水體中增加了溫度變量以及由溫差引起的浮力、傳質(zhì)、傳熱等各種效應(yīng)。由于有廣闊的受納條件,濱海電廠溫排水多采用深取淺排的直接排放方式,深層高鹽冷卻水經(jīng)過電廠機(jī)組后,高溫高鹽溫排水直接排放到海洋表層,并與周圍較冷海水混合,呈現(xiàn)出復(fù)雜的溫差密度分層和紊流流動(dòng)特性。溫差引起的浮力效應(yīng)、天然水體的復(fù)雜動(dòng)力條件以及密度層結(jié)效應(yīng)成為影響溫排水排放口近區(qū)流場的重要因素,這也使得濱海水域溫排水運(yùn)動(dòng)規(guī)律變得更為復(fù)雜。因此,適用于水動(dòng)力條件相對(duì)簡單的受納水體近、遠(yuǎn)區(qū)劃分概念已不再適用于濱海水域。為敘述方便,下文中將溫排水排放口附近溫升影響水域統(tǒng)稱為溫排水排放口近區(qū)。

    溫排水以及由此引起的水質(zhì)、生態(tài)環(huán)境變化研究已成為世界各國海洋環(huán)境學(xué)家十分關(guān)注的重要課題。國內(nèi)外海洋環(huán)境學(xué)家們從20世紀(jì)70年代開始就有關(guān)電廠溫排水對(duì)附近海域的影響進(jìn)行調(diào)查和研究。1984年聯(lián)合國海洋污染專家組(GESAMP)[2]對(duì)世界各國科學(xué)家前十幾年在該領(lǐng)域取得的研究成果進(jìn)行總結(jié)。王麗霞等[3-4]從確定溫排水與環(huán)境水溫的變化關(guān)系出發(fā),詳細(xì)分析了多種影響海洋水溫變化的因素,并對(duì)近年來熱擴(kuò)散預(yù)測方法的研究現(xiàn)狀做了回顧。蔣爽等[5]給出了近年來海水熱擴(kuò)散研究進(jìn)展的主要文獻(xiàn)目錄。20世紀(jì)中后期,電廠溫排水研究大多以數(shù)學(xué)模型解析解[6-9]和物理模型試驗(yàn)[10]為主。隨著電子計(jì)算機(jī)技術(shù)和數(shù)值計(jì)算方法的發(fā)展,數(shù)值模型計(jì)算已成為電廠溫排水研究的重要手段。但因?yàn)闇嘏欧帕W(xué)機(jī)理的復(fù)雜性以及濱海動(dòng)力條件的復(fù)雜性,濱海溫排放數(shù)值模型仍處在發(fā)展過程中,筆者主要針對(duì)濱?;?、核電廠溫排放數(shù)值模型的發(fā)展進(jìn)行回顧和評(píng)述。

    1 溫排放數(shù)值模型

    描述溫排水運(yùn)動(dòng)一般規(guī)律的理論為時(shí)均化的N-S方程和溫度輸移擴(kuò)散方程,其控制方程有2種,即二維模型方程和三維模型方程。電廠溫排水?dāng)?shù)值模擬研究起始于20世紀(jì)70年代中期,數(shù)值模型具有投資小、周期短且不存在物理模型中模型相似問題等優(yōu)越性。多年來人們系統(tǒng)地開展了溫排水輸移擴(kuò)散的數(shù)學(xué)模型研究,提出了多種較為實(shí)用、有效的理論模型和計(jì)算方法。但針對(duì)自然條件復(fù)雜的濱海水域,其模式選取、模型邊界控制等問題一直成為能否精確反映溫排水排放口近區(qū)溫升分布規(guī)律的關(guān)鍵。

    1.1 紊流模式的選取

    紊流模式能詮釋很多水動(dòng)力現(xiàn)象,諸如射流、異重流、浮力流和回流等。當(dāng)電廠溫排水的受納水域?yàn)樗畮?、湖泊、直流河段等天然水域時(shí),環(huán)境水體的動(dòng)力條件單一,溫排水主要受其初始動(dòng)量和浮力控制,在排放口附近形成紊流狀態(tài)。隨著動(dòng)量和浮力沿程衰減殆盡,溫排水隨環(huán)境水體運(yùn)移。在濱海水域,遠(yuǎn)離排放口不受溫排水初始動(dòng)量影響的水流仍表現(xiàn)為復(fù)雜的紊流狀態(tài)。因此,紊流模式的正確選取對(duì)數(shù)值模擬精度相當(dāng)重要。模擬溫排水運(yùn)動(dòng)的紊流模式主要有k-ε雙方程模式、雷諾應(yīng)力方程模式(DSM)和代數(shù)應(yīng)力模式(ASM)。

    1.1.1k-ε雙方程模式

    k-ε雙方程模式采用紊流動(dòng)能k和能量耗散率ε雙方程對(duì)時(shí)均化的N-S方程進(jìn)行封閉模擬,多年來已得到廣泛應(yīng)用。McGuirk等[11]較早用深度平均的k-ε雙方程模式模擬溫水側(cè)排問題,成功模擬了溫排水在明渠中的流場和溫升場,并能較好地處理回流區(qū),而后又發(fā)展了三維模式[12],用來模擬靜水中的溫排水。陽昌陸等[13]根據(jù)k-ε模型建立完全三維計(jì)算,超越了傳統(tǒng)的準(zhǔn)三維近似,成功地預(yù)測了重疊式取排水口的三維湍流流動(dòng)。

    在工程實(shí)踐中,溫差引起的浮力作用不可忽視,浮力有助于形成穩(wěn)定的分層,對(duì)紊動(dòng)有抑制和削弱作用,運(yùn)用紊流模式時(shí)應(yīng)考慮浮力影響。由于浮力使紊流導(dǎo)熱減少的程度比使紊流黏性減少的程度更甚,Rodi[14]分別在k和ε方程中添加浮力源項(xiàng)以修正k-ε方程,對(duì)紊流Prandtl數(shù)采用Munk-Anderson公式[15]進(jìn)行進(jìn)一步的修正。倪浩清等[16]在動(dòng)量方程和k方程及ε方程中均引入浮力項(xiàng),同時(shí)假設(shè)溫躍層中熱交換系數(shù)與速度和溫度的垂向梯度有關(guān),提出了一個(gè)半經(jīng)驗(yàn)公式對(duì)紊流Prandtl數(shù)加以修正,建立了浮力修正的k-ε模式,成功模擬了溫差異重流的溫躍層現(xiàn)象以及突擴(kuò)湍浮力回流問題,但經(jīng)多次檢驗(yàn)計(jì)算,ε方程中計(jì)及密度項(xiàng)的作用并不明顯。

    浮力修正使得標(biāo)準(zhǔn)k-ε模式無法反映浮力作用這一問題得以解決,但針對(duì)參數(shù)的修正關(guān)系式往往具有一定的經(jīng)驗(yàn)性,缺乏通用性,對(duì)某些各向異性的紊動(dòng)特性依然無能為力。

    1.1.2 雷諾應(yīng)力方程模式

    雷諾應(yīng)力方程模式從各向異性的前提出發(fā),直接封閉求解雷諾應(yīng)力的輸運(yùn)方程,計(jì)算應(yīng)力分量。Pathak等[17]分別用雷諾應(yīng)力方程模式和標(biāo)準(zhǔn)的k-ε模式模擬垂直于橫向來流的熱射流,結(jié)果表明雷諾應(yīng)力方程模式的模擬效果要優(yōu)于k-ε模式的模擬效果。范麗麗等[18]采用考慮浮力的雷諾應(yīng)力紊流模式對(duì)不同長寬比的排放口進(jìn)行三維模擬,得到擴(kuò)散寬度和排放口的長寬比相關(guān)關(guān)系。從實(shí)際應(yīng)用角度來講,雷諾應(yīng)力方程模式雖然精度高,但十分復(fù)雜,計(jì)算量過于龐大,在工程實(shí)踐中難以普及。

    1.1.3 代數(shù)應(yīng)力模式

    代數(shù)應(yīng)力模式是雷諾應(yīng)力方程模式的一種簡化,保留了紊流各向異性的特征,但計(jì)算量比雷諾應(yīng)力方程模式大幅度減少,而精度比k-ε模式高,結(jié)合了前者的通用性和后者的經(jīng)濟(jì)性。倪浩清等[19]基于Rodi[14]提出的代數(shù)應(yīng)力模式,通過簡化,將部分源項(xiàng)人為地置入到擴(kuò)散項(xiàng)中進(jìn)行計(jì)算,建立了簡化的紊流代數(shù)應(yīng)力模式,較好地反映了溫排水各向異性的特性。江春波等[20]基于代數(shù)應(yīng)力模式,改進(jìn)了求解流場的變量分割解法和求解對(duì)流擴(kuò)散方程的流線迎風(fēng)有限元格式,模擬了明渠中的二次流及溫躍層。金海生等[21]用一個(gè)簡化的應(yīng)力代數(shù)模式,成功地模擬了同時(shí)存在溫度和鹽度梯度的密度分層流動(dòng),并用水槽試驗(yàn)結(jié)果進(jìn)行驗(yàn)證,結(jié)果吻合良好。雖然代數(shù)應(yīng)力模式仍然存在收斂速度受限[22]的問題,但能夠較好地反映密度躍層這類各向異性的紊流流動(dòng)現(xiàn)象,特別是在分層流的流速、密度梯度不大時(shí),可給出精確的結(jié)果。因此,代數(shù)應(yīng)力模式是當(dāng)前較為實(shí)用的模式。

    1.2 時(shí)均模型的選取

    如前所述,描述溫排水運(yùn)動(dòng)一般規(guī)律的數(shù)值模型有二維模型和三維模型。在寬淺型水域,水體分層不明顯,二維模型基本上能夠反映出水流的水平運(yùn)動(dòng)規(guī)律和溫升分布,已在許多工程實(shí)際問題中取得一定效果[23-24]。河海大學(xué)[23]在大亞灣核電站的遠(yuǎn)區(qū)模擬中,為同時(shí)反映大洋潮波控制下的近區(qū)水力特性,首次通過嵌套耦合技術(shù)從大水域推算到局部水域。

    溫排水與受納水體的摻混及傳熱過程是一個(gè)三維過程。在溫排水排放口附近水域,溫差引起的浮力效應(yīng)是影響溫排水排放口近區(qū)流場和溫升場的重要因素,濱海水域復(fù)雜的海底地形、潮汐動(dòng)力條件以及明顯的鹽度層結(jié)效應(yīng)亦使得水域的溫排水在環(huán)境水體中的擴(kuò)散過程異常復(fù)雜,二維模型難以反映出水流和溫排水的時(shí)空分布規(guī)律。要了解溫排水垂向水流特性和溫度分布,必須進(jìn)行三維模擬。黃平[25]建立了汕頭港三維溫排放對(duì)流擴(kuò)散模型,用特征差分法求解,并用水槽試驗(yàn)結(jié)果進(jìn)行驗(yàn)證。王麗霞等[26]根據(jù)一階湍封閉理論建立了三維熱擴(kuò)散預(yù)測模型,模型中引入了計(jì)算網(wǎng)格無法分辨的次網(wǎng)格能量密度,將均勻流體的湍黏性與局部能量聯(lián)系起來,使湍流黏性系數(shù)的模擬較一般擴(kuò)散模型更加精確。因此,三維數(shù)值模型在一定程度代表了溫排放數(shù)值模擬的發(fā)展方向。

    現(xiàn)有的濱海三維溫排放數(shù)值模型中,控制方程大多采用Blumberg等[27]提出的基于靜水壓強(qiáng)假定和Boussinesq假定的三維水流動(dòng)力學(xué)和熱力學(xué)原始方程組。模型將動(dòng)力學(xué)和熱力學(xué)方程組耦合求解,并考慮海水密度變化所引起的斜壓項(xiàng)。實(shí)際工程應(yīng)用中,人們?cè)趯?duì)流場及溫升場進(jìn)行模擬時(shí),大多先模擬水動(dòng)力場,然后在水動(dòng)力場模擬的基礎(chǔ)上模擬溫升場,并且忽略由密度差異引起的密度梯度和斜壓效應(yīng),即正壓模型。這對(duì)水庫、湖泊等水動(dòng)力場相對(duì)較弱、密度均勻的環(huán)境水體來說,其計(jì)算精度是可以接受的。但是,淡鹽水交匯的濱海水域,水動(dòng)力場非常復(fù)雜,鹽度在水平方向和垂向上的分布不均勻,使得環(huán)境水體產(chǎn)生斜壓效應(yīng),而溫排水和環(huán)境水體的溫差亦會(huì)產(chǎn)生密度分層。研究表明,0.4%鹽度差的等溫度海水所具有的密度差,大體上相當(dāng)于10℃溫差的等鹽度海水所產(chǎn)生的密度差。由此可知,鹽度差引起的密度差大幅度超過溫差引起的密度差。因而,對(duì)于深取淺排的電廠取、排水口布置方案而言,溫排水的垂向分布狀態(tài)受鹽度的垂向梯度影響很大,導(dǎo)致溫排水不一定只存在于水體上層,也有可能在水體下層[28]。因此,在三維溫排水模型中,考慮斜壓效應(yīng)尤為重要,且需要將水動(dòng)力場和物質(zhì)擴(kuò)散場耦合求解。

    在實(shí)際應(yīng)用中,三維溫排水模型的許多計(jì)算參數(shù)的選取都與垂向坐標(biāo)系的選取直接相關(guān)。目前國內(nèi)外三維模型主要采用等平面(z)坐標(biāo)和地形擬合(σ)坐標(biāo),如國際上較為流行的 TRIM3D,POM,EFDC,DELFT-3D,FVCOM等。近年來,國內(nèi)已有多家科研單位引進(jìn)了這些模型和商業(yè)軟件進(jìn)行溫排水?dāng)?shù)值模擬研究,取得了一些成果[29-32]。但這些模型并不能很好地適用于地形變化劇烈、密度分層明顯的水域。陶建峰等[33]對(duì)三維水流數(shù)值模型中的多種垂向坐標(biāo)模式進(jìn)行了評(píng)述,同時(shí)指出,σ坐標(biāo)系下的三維模型在處理水平斜壓梯度力時(shí)會(huì)引入較大計(jì)算誤差,容易引起局部流態(tài)的失真。針對(duì)這一問題,陶建峰等[34]在垂向引入雙σ坐標(biāo)變換,對(duì)地形擬合坐標(biāo)系下的水平斜壓梯度力的計(jì)算誤差改進(jìn)進(jìn)行了嘗試。

    1.3 模型邊界的處理

    溫排放數(shù)值模型的邊界處理主要分為控制邊界處理和自由面散熱邊界處理2種。

    1.3.1 控制邊界處理

    控制邊界主要涉及2個(gè)方面,一是遠(yuǎn)區(qū)模型的外邊界控制,二是近區(qū)、遠(yuǎn)區(qū)模型的耦合邊界。近區(qū)、遠(yuǎn)區(qū)模型在運(yùn)用上簡單方便,但忽視了兩區(qū)之間的相互作用,近區(qū)模型往往忽略了遠(yuǎn)區(qū)回流的影響,而遠(yuǎn)區(qū)模型也無法精確反映排放口近區(qū)的特性。因此耦合模型和全場模型應(yīng)運(yùn)而生。

    耦合模型是指在近區(qū)、遠(yuǎn)區(qū)采用不同的模型,通過一定的方法將兩者連接、耦合。早期耦合模型往往在近區(qū)采用射流或浮射流積分模型,遠(yuǎn)區(qū)則采用淺水環(huán)流對(duì)流擴(kuò)散模型。Kaufman等[35]建立的耦合模型在近區(qū)采用S-H模型,遠(yuǎn)區(qū)采用二維有限元數(shù)學(xué)模型CAFE和DISPER,但網(wǎng)格尺寸和時(shí)間步長都較小,比較耗機(jī)時(shí),成本較高。Adams等[36]在其基礎(chǔ)上用潮汐海灣分析模型TEA代替CAFE模型計(jì)算二維水流,用歐拉-拉格朗日分析模型 ELA代替DISPER模型計(jì)算溫度場,成功地用于2座電廠的溫排放模擬。近區(qū)、遠(yuǎn)區(qū)模型在范圍劃定上存在一定的困難,如何銜接也具有一定的難度,因而銜接方法的選取就變得十分關(guān)鍵。Zhang等[37]針對(duì)近、遠(yuǎn)區(qū)銜接方法做了很好地闡述。Zhang[38]用粒子追蹤法耦合近區(qū)RSB模型和遠(yuǎn)區(qū)ECOMSI模型,將粒子密度轉(zhuǎn)化為濃度,并考慮兩區(qū)之間的相互作用。Suh[39]用高斯puff算法耦合近區(qū)CORMIX模型和遠(yuǎn)區(qū)TEA模型以及ELA模型。Maderich等[40]用DESA法在近區(qū)拉格朗日積分模型JETLAG和遠(yuǎn)區(qū)THREETOX模型之間建立雙向水動(dòng)力連接。紀(jì)平等[41]采用湍射流模型計(jì)算近、遠(yuǎn)區(qū)交界上的速度值,同時(shí)將其作為遠(yuǎn)區(qū)二維淺水環(huán)流模型的邊界進(jìn)行耦合計(jì)算。

    全場模型避免了分區(qū)模型帶來的劃分及銜接問題,使流場和溫升場的預(yù)測結(jié)果更為合理。20世紀(jì)90年代初,倪浩清等[42-43]通過引入考慮速度、溫度垂向分布不均勻性的流散系數(shù),建立深度平均的k-ε紊流全場模式,不僅模擬了大水域冷卻池,而且對(duì)室內(nèi)冷卻池的模擬效果也較好。梁書秀等[44-45]提出了深度平均的代數(shù)應(yīng)力全場模式,假設(shè)流體不可壓縮,且不考慮浮力、表面風(fēng)應(yīng)力及柯氏力作用,對(duì)明渠和潮汐水域的側(cè)排問題進(jìn)行模擬,與深度平均的k-ε紊流模式的計(jì)算結(jié)果對(duì)比,該模式優(yōu)勢顯著。Yu等[46]建立了深度平均的~k-~w模式,用3個(gè)溫排放算例進(jìn)行檢驗(yàn),并與深度平均的k-ε紊流模式比較,結(jié)果表明當(dāng)排放口尺寸相對(duì)較小或排放速度相對(duì)較小時(shí),前者要優(yōu)于后者。

    1.3.2 自由面散熱邊界處理

    溫排水排入到濱海水域后,通過近區(qū)摻混稀釋、遠(yuǎn)區(qū)對(duì)流擴(kuò)散,其廢熱直接進(jìn)入到天然海水中,通過水溫升高的自由水面蒸發(fā)、對(duì)流和輻射3種機(jī)制散失于大氣中。因此,水面散熱系數(shù)的選取直接關(guān)系到水面散熱能力和環(huán)境水體受熱容量的計(jì)算精度。

    水面散熱機(jī)制十分復(fù)雜,陳惠泉等[47]、毛世民等[48]根據(jù)氣象可控專用風(fēng)洞系統(tǒng)試驗(yàn)研究成果和全國各地代表性水域(以湖泊和水庫為主)現(xiàn)場實(shí)測數(shù)據(jù),通過分析氣溫、水面風(fēng)速、水溫及水汽溫差等各因子對(duì)水面蒸發(fā)系數(shù)的相關(guān)關(guān)系,提出了充分考慮溫差效應(yīng)的全國通用水面蒸發(fā)散熱系數(shù)公式,統(tǒng)一了自由、強(qiáng)迫對(duì)流的散熱機(jī)理,已推薦成為《電力工程水文技術(shù)規(guī)程》及《工業(yè)循環(huán)水冷卻設(shè)計(jì)規(guī)范》等規(guī)范公式。國外學(xué)者將水體表面熱通量考慮為短波輻射、長波輻射、感熱通量和潛熱通量之和,并給出了相應(yīng)的計(jì)算公式[49]。目前大多數(shù)溫排放研究工作中,水面蒸發(fā)散熱系數(shù)通常采用規(guī)范公式計(jì)算甚至用常數(shù)代替,這對(duì)于湖泊、水庫等水域的溫排放數(shù)值模擬,其計(jì)算精度是可以接受的,但對(duì)于河口海岸水域,其海洋性氣候條件與湖泊、水庫等的氣候條件差異較大,水面上時(shí)有波浪發(fā)生,波浪條件下的散熱問題變得十分復(fù)雜,在此條件下規(guī)范公式的通用性存在不少值得推敲之處,需接受實(shí)測資料的進(jìn)一步檢驗(yàn)。

    2 結(jié) 語

    迄今為止,溫排放數(shù)值模型的研究為海洋水環(huán)境的實(shí)用性研究提供了許多寶貴的經(jīng)驗(yàn)和方法,但溫排水三維數(shù)值模擬方法研究仍處于探索階段,濱海電廠溫排水的擴(kuò)散規(guī)律已引起國內(nèi)外學(xué)者的極大重視。濱海水域岸線曲折,海底地形變化劇烈,水動(dòng)力條件復(fù)雜且密度層結(jié)效應(yīng)明顯,加之溫排水與環(huán)境水域溫差引起的斜壓效應(yīng),使得溫排水排放口近區(qū)的局部流態(tài)和溫升分布變得非常復(fù)雜。我國濱?;?、核電廠建設(shè)發(fā)展迅速,針對(duì)現(xiàn)有溫排放數(shù)值模式的不足,嘗試建立濱海水域溫度、鹽度大梯度水體高精度三維斜壓溫排水?dāng)?shù)值模型,探求更加符合實(shí)際的紊流模型和有效的計(jì)算方法,也許是今后的研究方向。但由于需要同時(shí)考慮海水斜壓效應(yīng)和溫差引起的斜壓效應(yīng),其理論推導(dǎo)和數(shù)值算法實(shí)現(xiàn)也較湖泊、水庫、河流等溫排放模型更為復(fù)雜。因此三維溫排放數(shù)值模式的發(fā)展將是一個(gè)十分漫長的過程。

    [1]金嵐.水域熱影響概論[M].北京:高等教育出版社,1993.

    [2]GESAMP.Thermal discharges in the marine environment[R].Roma:UHEP Regional Seas Reports andStudies,1984:1-44.

    [3]王麗霞,孫英蘭,田暉.熱擴(kuò)散預(yù)測方法研究概況:影響海洋水溫的因素[J].海洋科學(xué),1997,5:24-25.

    [4]王麗霞,孫英蘭,田暉.熱擴(kuò)散預(yù)測方法研究概況:熱擴(kuò)散的研究現(xiàn)狀[J].海洋科學(xué),1997,6:16-17.

    [5]蔣爽,端木琳,王樹剛.海水熱擴(kuò)散研究進(jìn)展與新問題分析[J].能源環(huán)境保護(hù),2006,20(5):5-9.

    [6]張書農(nóng),孔佩文.溫差異重流運(yùn)行長度的研究[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,1987,15(2):1-11.

    [7]蔣春風(fēng).感潮河段電廠溫排水熱影響范圍的特征分析[J].河海大學(xué)學(xué)報(bào):自然科學(xué)版,1990,18(6):114-121.

    [8]JIRKA G H,DONEKER R L,HINTON S W.User's manual for CORMIX:a hydrodynamic mixing zone model and decision support system for pollutant discharges into surface waters[M].Washington D.C.:U.S.EPA,2006.

    [9]STOLZENBACH K D,HARLEMAN D R F.An analytical and experimental investigation of surface discharges of heated water[R].Cambridge:MassachusettsInstitute of Technology,1971.

    [10]陳惠泉,許玉麟,賀益英.火、核電廠冷卻水試驗(yàn)研究50年的進(jìn)展和體驗(yàn)[J].中國水利水電科學(xué)研究院學(xué)報(bào),2008,6(4):288-298.

    [11]McGUIRK J J,RODI W.A depth-averaged mathematical model for the near field of the side discharge into openchannel flow[J].Journal of Fluid Mechanics,1978,86(4):761-781.

    [12]McGUIRK J J,RODI W.Mathematical modeling of threedimensional heated surface jets[J].Journal of Fluid Mechanics,1979,95(4):609-633.

    [13]陽昌陸,倪浩清.完全三維湍流數(shù)值模擬及其在冷卻水工程中的應(yīng)用[J].水利學(xué)報(bào),1993(5):12-22.

    [14]RODI W.Turbulent model for environmental problem[C]//Prediction Methods for Turbulent Flows.Von Karman:Hemisphere Publishing Coorperation,1980:260-349.

    [15]MUNK W H,ANDER SON E R.Notes on a theory of the thermocline[J].Journal of Marine Research,1948,7:276-295.

    [16]倪浩清,賀益英,王能家,等.明渠溫差異重流中湍浮力回流的數(shù)學(xué)模型[J].中國科學(xué):A輯,1987,2(4):437-448.

    [17]PATHAK M,DEWAN A,DASS A K.Computational prediction ofaslightlyheated turbulent rectangularjet discharged into a narrow channel crossflow using two different turbulence models[J].International Journal of Heat and Mass Transfer,2006,49:3914-3928.

    [18]范麗麗,沙海飛,楊海燕,等.熱污水表面?zhèn)扰诺慕鼌^(qū)特性[J].江蘇大學(xué)學(xué)報(bào):自然科學(xué)版,2007,28(2):151-155.

    [19]倪浩清,王能家,周力行.應(yīng)力代數(shù)模型在各向異性湍浮力回流的應(yīng)用[J].力學(xué)學(xué)報(bào),1989,21(1):26-34.

    [20]江春波,許協(xié)慶.各向異性湍浮力流的三維有限元計(jì)算[J].水利學(xué)報(bào),1991(7):8-18.

    [21]金海生,張書農(nóng).同時(shí)存在溫度和鹽度的分層流的數(shù)值模擬[J].水動(dòng)力學(xué)研究與進(jìn)展:A輯,1991,6(4):72-78.

    [22]張?jiān)?楊永全,吳持恭.混合k-ε模型研究[J].水動(dòng)力學(xué)研究與進(jìn)展:A輯,1992,7(增刊):571-576.

    [23]河海大學(xué).大亞灣核電站排放的廢熱水對(duì)遠(yuǎn)區(qū)的影響[R].南京:河海大學(xué),1988.

    [24]中國水利水電科學(xué)研究院水力學(xué)所.秦山核電站冷卻水排放數(shù)學(xué)模型研究[R].北京:中國水利水電科學(xué)研究院,1990.

    [25]黃平.汕頭港水域溫排水熱擴(kuò)散的三維數(shù)值模擬[J].海洋環(huán)境科學(xué),1996,15(1):59-65.

    [26]王麗霞,孫英蘭,鄭連遠(yuǎn).三維熱擴(kuò)散預(yù)測模型[J].青島海洋大學(xué)學(xué)報(bào),1998,28(1):29-35.

    [27]BLUMBER G A F,MELLOR G L.A description of a threedimensionalcoastal ocean circulation modelin threedimensional coastal ocean models[C]//HEAPS N.Coastal and EstuarineSciences.Washington D.C.:American Geophysical Union,1987:1-16.

    [28]李瑞生,劉長貴.鹽淡水交匯水域電廠溫排水及航道取水的試驗(yàn)研究[J].水動(dòng)力學(xué)研究與進(jìn)展:A輯,1996,11(6):661-671.

    [29]汪一航,魏澤勛,王永剛,等.潮汐潮流三維數(shù)值模擬在莊河電廠溫排水問題中的應(yīng)用[J].海洋通報(bào),2006,25(1):8-15.

    [30]周巧菊.大亞灣海域溫排水三維數(shù)值模擬[J].海洋湖沼通報(bào),2007(4):37-46.

    [31]朱軍政.強(qiáng)潮海灣溫排水三維數(shù)值模擬[J].水力發(fā)電學(xué)報(bào),2007,26(4):56-60.

    [32]曹穎,朱軍政.基于FVCOM模式的溫排水三維數(shù)值模擬研究[J].水動(dòng)力學(xué)研究與進(jìn)展:A輯,2009,24(4):432-439.

    [33]陶建峰,張長寬.河口海岸三維水流數(shù)值模型中幾種垂向坐標(biāo)模式研究述評(píng)[J].海洋工程,2007,25(1):133-142.

    [34]陶建峰,張長寬.河口海岸三維雙σ坐標(biāo)斜壓水流數(shù)值模型研究[J].海洋工程,2008,26(1):71-76.

    [35]KAUFMAN J T,ADAMS E E.Coupled near and far field thermal plume analysis using finite element techniques[R].Cambridge:Massachusetts Institute of Technology Energy Laboratory,1981.

    [36]ADAMS E E,COSLER D J.Predicting circulation and dispersion nearcoastal power plants:applications using models TEA and ELA[R].Cambridge:Massachusetts Institute of Technology Energy Laboratory,1987.

    [37]ZHANG Xue-yong,ADAMS E E.Prediction of near field plume characteristics using far field circulation model[J].Journal of Hydraulic Engineering,1999,125(3):233-241.

    [38]ZHANG Xue-yong.Ocean outfall modeling-interfacing near and far field models with particle tracking method[D].Massachusetts:Massachusetts Institute of Technology,1995.

    [39]SUH S W.A hybrid near-field/far-field thermal discharge model for coastal areas[J].Marine Pollution Bulletin,2001,43(7/8/9/10/11/12):225-233.

    [40]MADERICH V,HELING R,BEZHENARR,etal.Developmentand application of3D numericalmodel THREETOX to the prediction of cooling water transport and mixing in the inland and coastal waters[J].Hydrological Processes,2008,22:1000-1013.

    [41]紀(jì)平,秦曉,袁玨,等.冷卻池近區(qū)湍射流與遠(yuǎn)區(qū)環(huán)流耦合計(jì)算模型[J].水利水電技術(shù),2006,37(9):23-26.

    [42]倪浩清,龍中林,周力行.深度平均的紊流全場水環(huán)境新模型及其在大水域冷卻池中的應(yīng)用[J].水利學(xué)報(bào),1991(9):7-13.

    [43]倪浩清,沈永明,陳惠泉.深度平均的k-ε紊流全場模型及其驗(yàn)證[J].水利學(xué)報(bào),1994(11):8-17.

    [44]梁書秀,沈永明,孫昭晨.深度平均的應(yīng)力-通量代數(shù)全場模型及其驗(yàn)證[J].海洋環(huán)境科學(xué),1999,18(3):33-38.

    [45]梁書秀,沈永明,孫昭晨.潮汐水域中污水側(cè)向排放的數(shù)值模擬[J].海洋學(xué)報(bào),2000,22(2),113-119.

    [46]YU Li-ren,RIGHETTO A M.Depth-averaged turbulence~k-~wmodel and applications[J].Advancesin Engineering Software,2001,32(5):375-394.

    [47]陳惠泉,毛世民.水面蒸發(fā)系數(shù)全國通用公式的驗(yàn)證[J].水科學(xué)進(jìn)展,1995,6(2):116-120.

    [48]毛世民,陳惠泉.超溫水體蒸發(fā)系數(shù)通用公式的研究與應(yīng)用[J].水利水電科技進(jìn)展,2007,27(6):85-89.

    [49]SIMPSON J J,DICKEY T D.Alternative parameterizations of downward irradiance and their dynamical significance[J].Journal of Physical Oceanography,1981,11(6):876-882.

    Advances in numerical simulation of thermal discharge of heat or nuclear power plants in coastal areas

    ZHANG Changkuan1,2,YAO Jing1,2,TAO Jian-feng1,2(1.State Key Laboratory of Hydrology-Water Resources and Hydraulic Engineering,Hohai University,Nanjing210098,China;2.College of Harbour,Coastal and Offshore Engineering,Hohai University,Nanjing210098,China)

    The advances in numerical simulation of thermal discharge of worldwide heat or nuclear power plants in coastal areas in recent yearswere reviewed.The development situations and some basic issues of the existing numerical models for the thermal discharge were highlighted,mainly including selection of turbulent modes of the thermal discharge and control of model boundaries.A conclusion is drawn that 3D numerical models are regarded to represent the development tendency of those for the thermal discharge to a certain extent.The oblique pressure effects should be considered.The hydrodynamic fields and the diffusion fields of matters should be coupled for the solution.

    water environment;thermal discharge;thermal pollution;numerical simulation

    TK124

    A

    1006-7647(2010)03-0084-05

    10.3880/j.issn.1006-7647.2010.03.021

    國家自然科學(xué)基金(50779016);高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金(20090094120005)

    張長寬(1954—),江蘇興化人,教授,從事海岸動(dòng)力學(xué)、近海水環(huán)境及物質(zhì)輸運(yùn)研究。E-mail:ckzhang@hhu.edu.cn

    2009-10-21 編輯:駱超)

    猜你喜歡
    排放口濱海浮力
    水下航行器熱射流多孔排放口結(jié)構(gòu)仿真分析
    遮流罩排放口對(duì)水下航行器熱射流溫度場影響的數(shù)值仿真
    “浮力”知識(shí)鞏固
    我們一起來“制服”浮力
    浮力大小由誰定
    煤制油項(xiàng)目有組織排放口VOCs排放源強(qiáng)分析
    化工管理(2021年19期)2021-08-02 01:25:38
    濱海白首烏
    濱海頂層公寓
    岳濱海 藏石欣賞
    寶藏(2018年12期)2019-01-29 01:50:56
    第一類污染物廢水總排放口“超標(biāo)”應(yīng)如何處理
    精品一品国产午夜福利视频| 国产欧美日韩一区二区三| 一a级毛片在线观看| 国产黄a三级三级三级人| 十八禁人妻一区二区| 十分钟在线观看高清视频www| 黄色片一级片一级黄色片| 黄色a级毛片大全视频| 国产精品1区2区在线观看.| 97人妻精品一区二区三区麻豆 | 亚洲 欧美 日韩 在线 免费| 久久国产亚洲av麻豆专区| 一个人免费在线观看的高清视频| 在线观看午夜福利视频| 一区在线观看完整版| 99在线人妻在线中文字幕| 免费少妇av软件| 黄色 视频免费看| 黄网站色视频无遮挡免费观看| 黄色 视频免费看| 国产av在哪里看| 亚洲片人在线观看| 国产一区二区在线av高清观看| 99国产综合亚洲精品| а√天堂www在线а√下载| 日本 欧美在线| 国产欧美日韩一区二区精品| 亚洲成a人片在线一区二区| 色综合欧美亚洲国产小说| 最新美女视频免费是黄的| 免费一级毛片在线播放高清视频 | 国产成人系列免费观看| 亚洲无线在线观看| 亚洲成a人片在线一区二区| 一级毛片高清免费大全| 亚洲伊人色综图| 琪琪午夜伦伦电影理论片6080| 欧美乱妇无乱码| 久久人人爽av亚洲精品天堂| 每晚都被弄得嗷嗷叫到高潮| or卡值多少钱| 色综合欧美亚洲国产小说| 欧美 亚洲 国产 日韩一| 涩涩av久久男人的天堂| 高清黄色对白视频在线免费看| 极品人妻少妇av视频| 真人做人爱边吃奶动态| 黑人巨大精品欧美一区二区蜜桃| 中文字幕另类日韩欧美亚洲嫩草| 欧美黄色淫秽网站| 99国产精品一区二区三区| 亚洲成人精品中文字幕电影| 啪啪无遮挡十八禁网站| √禁漫天堂资源中文www| 怎么达到女性高潮| 国产片内射在线| 午夜视频精品福利| 少妇 在线观看| www日本在线高清视频| 国产亚洲av高清不卡| 美女大奶头视频| 欧美黄色片欧美黄色片| 国产精品av久久久久免费| 久久九九热精品免费| 色尼玛亚洲综合影院| 国产1区2区3区精品| 动漫黄色视频在线观看| 视频区欧美日本亚洲| 国产精品秋霞免费鲁丝片| 精品久久久久久久人妻蜜臀av | 久久香蕉国产精品| 中文字幕色久视频| 欧美日韩精品网址| 黄色女人牲交| 少妇裸体淫交视频免费看高清 | 9色porny在线观看| 麻豆久久精品国产亚洲av| 天堂√8在线中文| 在线观看免费午夜福利视频| 美国免费a级毛片| 久久久国产欧美日韩av| 亚洲国产欧美一区二区综合| 久久人妻av系列| 身体一侧抽搐| 天天添夜夜摸| 美女扒开内裤让男人捅视频| 国产精品1区2区在线观看.| 丰满的人妻完整版| 99精品久久久久人妻精品| 日日摸夜夜添夜夜添小说| 亚洲 国产 在线| 国产精品av久久久久免费| 中文字幕色久视频| 一区二区三区激情视频| 美女高潮到喷水免费观看| 久久久精品国产亚洲av高清涩受| 久久精品亚洲熟妇少妇任你| 色综合亚洲欧美另类图片| 久久久久国内视频| 99精品在免费线老司机午夜| 亚洲美女黄片视频| 高清黄色对白视频在线免费看| 欧美国产日韩亚洲一区| 国产欧美日韩一区二区三区在线| 天堂√8在线中文| 国产精品av久久久久免费| 激情在线观看视频在线高清| 一边摸一边抽搐一进一出视频| 97人妻天天添夜夜摸| 色播在线永久视频| 亚洲一区二区三区色噜噜| avwww免费| 99在线人妻在线中文字幕| 最近最新免费中文字幕在线| 国产精品亚洲av一区麻豆| 一级毛片精品| 国产亚洲av高清不卡| 国产主播在线观看一区二区| av有码第一页| 久久性视频一级片| 少妇的丰满在线观看| 国内精品久久久久久久电影| √禁漫天堂资源中文www| 久久久久久久久中文| 国产亚洲精品久久久久久毛片| 亚洲自拍偷在线| 国产午夜精品久久久久久| 1024视频免费在线观看| 十分钟在线观看高清视频www| 中文字幕久久专区| 黄色成人免费大全| 免费少妇av软件| 久久久精品欧美日韩精品| 欧美色视频一区免费| 成人18禁在线播放| 757午夜福利合集在线观看| 波多野结衣高清无吗| 老汉色av国产亚洲站长工具| 久久精品成人免费网站| 精品卡一卡二卡四卡免费| 亚洲国产毛片av蜜桃av| 国产色视频综合| 色综合婷婷激情| 亚洲国产看品久久| 青草久久国产| 国产亚洲精品综合一区在线观看 | 又黄又爽又免费观看的视频| 亚洲av美国av| 精品久久蜜臀av无| 不卡一级毛片| 亚洲九九香蕉| 啦啦啦观看免费观看视频高清 | 亚洲国产中文字幕在线视频| 纯流量卡能插随身wifi吗| 久久人妻福利社区极品人妻图片| 人人澡人人妻人| 嫩草影视91久久| 亚洲精品中文字幕一二三四区| 亚洲精品国产色婷婷电影| 亚洲成人久久性| 久久人妻福利社区极品人妻图片| 黄色a级毛片大全视频| 亚洲少妇的诱惑av| 久久九九热精品免费| 国产亚洲精品av在线| 国产高清激情床上av| 欧美色视频一区免费| 熟妇人妻久久中文字幕3abv| 深夜精品福利| 国产单亲对白刺激| ponron亚洲| 亚洲av第一区精品v没综合| 久久久久精品国产欧美久久久| 少妇粗大呻吟视频| 国产精品 国内视频| 成人欧美大片| 免费av毛片视频| 久99久视频精品免费| 欧美 亚洲 国产 日韩一| 久久中文看片网| 99久久精品国产亚洲精品| 久热爱精品视频在线9| 欧美激情极品国产一区二区三区| 欧美中文综合在线视频| 一区二区三区高清视频在线| 亚洲黑人精品在线| 丝袜美足系列| 久久天堂一区二区三区四区| 国内精品久久久久久久电影| 黄网站色视频无遮挡免费观看| 久久久精品欧美日韩精品| www.熟女人妻精品国产| 伊人久久大香线蕉亚洲五| 高清黄色对白视频在线免费看| 夜夜躁狠狠躁天天躁| 欧美av亚洲av综合av国产av| 好男人在线观看高清免费视频 | 一a级毛片在线观看| 欧美色视频一区免费| 久久精品亚洲熟妇少妇任你| 国产三级黄色录像| 亚洲男人的天堂狠狠| 午夜精品国产一区二区电影| 在线观看午夜福利视频| 一本大道久久a久久精品| 亚洲第一青青草原| 精品午夜福利视频在线观看一区| 日韩欧美国产一区二区入口| 老鸭窝网址在线观看| 一级毛片精品| 婷婷六月久久综合丁香| www国产在线视频色| 精品国产超薄肉色丝袜足j| 91av网站免费观看| 99国产精品99久久久久| 99精品欧美一区二区三区四区| 欧美中文日本在线观看视频| 999久久久国产精品视频| 一个人观看的视频www高清免费观看 | 亚洲国产中文字幕在线视频| 久久人人97超碰香蕉20202| 国产精品久久久久久人妻精品电影| 日韩大尺度精品在线看网址 | 美女扒开内裤让男人捅视频| 18禁观看日本| 后天国语完整版免费观看| 欧美精品啪啪一区二区三区| 露出奶头的视频| 中文字幕最新亚洲高清| 日韩欧美国产一区二区入口| www.999成人在线观看| 男女下面插进去视频免费观看| 中文字幕色久视频| 大型av网站在线播放| 国产成人精品在线电影| 久久九九热精品免费| 精品久久久久久久久久免费视频| 亚洲专区国产一区二区| 亚洲欧美日韩无卡精品| 久久久精品欧美日韩精品| e午夜精品久久久久久久| 亚洲专区字幕在线| 啦啦啦免费观看视频1| 亚洲精品中文字幕一二三四区| 欧美日韩亚洲国产一区二区在线观看| 曰老女人黄片| 午夜免费观看网址| 国产色视频综合| 日日摸夜夜添夜夜添小说| 嫩草影视91久久| 国内毛片毛片毛片毛片毛片| 久久久久国产一级毛片高清牌| 久久国产精品影院| 亚洲av美国av| 久久午夜亚洲精品久久| 国产精品免费视频内射| 国产成人精品久久二区二区91| 一级毛片高清免费大全| 99在线人妻在线中文字幕| 日本三级黄在线观看| 涩涩av久久男人的天堂| 久久这里只有精品19| 亚洲av五月六月丁香网| 制服丝袜大香蕉在线| 成人18禁高潮啪啪吃奶动态图| 亚洲成人精品中文字幕电影| 午夜精品久久久久久毛片777| 一二三四在线观看免费中文在| 丝袜美腿诱惑在线| 看黄色毛片网站| 男女床上黄色一级片免费看| 午夜成年电影在线免费观看| 精品国产一区二区三区四区第35| 夜夜夜夜夜久久久久| 欧美一区二区精品小视频在线| 视频在线观看一区二区三区| 9热在线视频观看99| 国产精品爽爽va在线观看网站 | 亚洲一卡2卡3卡4卡5卡精品中文| av电影中文网址| 亚洲片人在线观看| 久久性视频一级片| 一级a爱视频在线免费观看| 久久中文字幕一级| 亚洲国产日韩欧美精品在线观看 | bbb黄色大片| 亚洲 欧美 日韩 在线 免费| 日韩高清综合在线| 久久国产精品影院| 久久天躁狠狠躁夜夜2o2o| 无限看片的www在线观看| 女性生殖器流出的白浆| 中文字幕av电影在线播放| 久久天堂一区二区三区四区| 亚洲国产中文字幕在线视频| 男女下面进入的视频免费午夜 | 91大片在线观看| 亚洲第一av免费看| 精品欧美一区二区三区在线| 天天躁狠狠躁夜夜躁狠狠躁| 夜夜夜夜夜久久久久| 天天躁夜夜躁狠狠躁躁| 可以在线观看的亚洲视频| 夜夜看夜夜爽夜夜摸| 母亲3免费完整高清在线观看| av电影中文网址| videosex国产| 成人精品一区二区免费| 成人免费观看视频高清| 一区二区三区国产精品乱码| 女性生殖器流出的白浆| 中文字幕高清在线视频| 久久久久久久午夜电影| 精品高清国产在线一区| 国产精品爽爽va在线观看网站 | 欧美日韩瑟瑟在线播放| 精品免费久久久久久久清纯| 欧美 亚洲 国产 日韩一| 黄色视频,在线免费观看| 看片在线看免费视频| 国产精品二区激情视频| 高潮久久久久久久久久久不卡| 91麻豆精品激情在线观看国产| 一卡2卡三卡四卡精品乱码亚洲| 天天躁狠狠躁夜夜躁狠狠躁| 可以免费在线观看a视频的电影网站| 久久欧美精品欧美久久欧美| 欧美色欧美亚洲另类二区 | 亚洲男人的天堂狠狠| 搡老岳熟女国产| 欧美丝袜亚洲另类 | 国产亚洲欧美精品永久| 国产一区在线观看成人免费| 亚洲精品国产一区二区精华液| 亚洲欧美日韩无卡精品| 精品午夜福利视频在线观看一区| 一a级毛片在线观看| 大码成人一级视频| 精品国产一区二区三区四区第35| 国产精品久久电影中文字幕| 午夜福利一区二区在线看| 人妻久久中文字幕网| 国产在线观看jvid| 欧美另类亚洲清纯唯美| 一级a爱片免费观看的视频| 亚洲成a人片在线一区二区| 乱人伦中国视频| 国产极品粉嫩免费观看在线| 欧美+亚洲+日韩+国产| 别揉我奶头~嗯~啊~动态视频| 亚洲中文字幕一区二区三区有码在线看 | 看片在线看免费视频| 在线天堂中文资源库| 久久青草综合色| 国产精品一区二区三区四区久久 | 欧美日本视频| 午夜亚洲福利在线播放| 搞女人的毛片| 69av精品久久久久久| 一边摸一边抽搐一进一出视频| 亚洲av熟女| 亚洲最大成人中文| av超薄肉色丝袜交足视频| 一边摸一边抽搐一进一出视频| 久久久久久久精品吃奶| 在线观看66精品国产| 欧美成人一区二区免费高清观看 | 久久久精品国产亚洲av高清涩受| 免费看a级黄色片| 一级作爱视频免费观看| 777久久人妻少妇嫩草av网站| 视频区欧美日本亚洲| 欧美国产日韩亚洲一区| 自线自在国产av| 人人妻人人澡人人看| 亚洲全国av大片| 热99re8久久精品国产| 韩国精品一区二区三区| e午夜精品久久久久久久| 午夜影院日韩av| 国产单亲对白刺激| 亚洲精品久久成人aⅴ小说| 欧洲精品卡2卡3卡4卡5卡区| 国产精品免费一区二区三区在线| 久久午夜综合久久蜜桃| www.999成人在线观看| 久久青草综合色| 亚洲精品粉嫩美女一区| 中文字幕最新亚洲高清| 亚洲成人精品中文字幕电影| 欧美乱妇无乱码| 岛国在线观看网站| 99国产精品免费福利视频| 热99re8久久精品国产| 国语自产精品视频在线第100页| 琪琪午夜伦伦电影理论片6080| 欧美国产日韩亚洲一区| 国产午夜福利久久久久久| 亚洲一码二码三码区别大吗| 91精品国产国语对白视频| av视频免费观看在线观看| 国产区一区二久久| 99久久久亚洲精品蜜臀av| 色播亚洲综合网| 欧美乱码精品一区二区三区| 岛国视频午夜一区免费看| 久久伊人香网站| 天天躁狠狠躁夜夜躁狠狠躁| 精品一区二区三区四区五区乱码| 麻豆久久精品国产亚洲av| 久久精品成人免费网站| 色婷婷久久久亚洲欧美| 天天添夜夜摸| 窝窝影院91人妻| 欧美黄色片欧美黄色片| 日韩欧美在线二视频| 97人妻精品一区二区三区麻豆 | 淫秽高清视频在线观看| 欧美大码av| 在线观看一区二区三区| 久久久精品国产亚洲av高清涩受| 久久精品91蜜桃| 嫁个100分男人电影在线观看| 欧美绝顶高潮抽搐喷水| 国产午夜精品久久久久久| 制服人妻中文乱码| 男女下面进入的视频免费午夜 | 男人操女人黄网站| 亚洲中文字幕日韩| 久久精品91无色码中文字幕| 不卡av一区二区三区| 午夜福利免费观看在线| 日本一区二区免费在线视频| 亚洲va日本ⅴa欧美va伊人久久| 国产一区二区在线av高清观看| 久久人妻av系列| 精品人妻在线不人妻| 老司机福利观看| 99国产极品粉嫩在线观看| 精品熟女少妇八av免费久了| 国产av在哪里看| 人妻久久中文字幕网| 欧美午夜高清在线| 人妻丰满熟妇av一区二区三区| 久久久久久大精品| www日本在线高清视频| 一边摸一边抽搐一进一小说| 国产精品综合久久久久久久免费 | 亚洲成人久久性| 亚洲欧美激情综合另类| 国产不卡一卡二| 久久热在线av| 一个人免费在线观看的高清视频| 美女国产高潮福利片在线看| 如日韩欧美国产精品一区二区三区| 好男人电影高清在线观看| 最新在线观看一区二区三区| 国产精品久久久人人做人人爽| 怎么达到女性高潮| 亚洲性夜色夜夜综合| 色婷婷久久久亚洲欧美| 91国产中文字幕| 97碰自拍视频| 亚洲国产中文字幕在线视频| 老汉色av国产亚洲站长工具| 欧美乱色亚洲激情| 亚洲精品久久国产高清桃花| 自线自在国产av| 婷婷丁香在线五月| 免费在线观看完整版高清| 久久国产亚洲av麻豆专区| 午夜亚洲福利在线播放| 亚洲熟女毛片儿| 无人区码免费观看不卡| 午夜a级毛片| 国产精品一区二区免费欧美| 露出奶头的视频| 手机成人av网站| 久久九九热精品免费| av天堂在线播放| 精品久久久久久成人av| 亚洲熟妇熟女久久| 一本大道久久a久久精品| 久久精品成人免费网站| 伊人久久大香线蕉亚洲五| 精品日产1卡2卡| 男女床上黄色一级片免费看| 成人永久免费在线观看视频| 高清毛片免费观看视频网站| 久久欧美精品欧美久久欧美| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品国产精品久久久不卡| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜免费成人在线视频| 日本vs欧美在线观看视频| 久久久久亚洲av毛片大全| 一进一出抽搐动态| 久久久久久免费高清国产稀缺| 日韩av在线大香蕉| 中文字幕精品免费在线观看视频| av有码第一页| 怎么达到女性高潮| 日韩欧美一区二区三区在线观看| 性欧美人与动物交配| 日韩欧美国产一区二区入口| 日韩有码中文字幕| 日本五十路高清| 亚洲第一电影网av| 久久久国产精品麻豆| 亚洲av五月六月丁香网| 可以在线观看的亚洲视频| 午夜福利欧美成人| 欧美最黄视频在线播放免费| 亚洲七黄色美女视频| 久久久久久免费高清国产稀缺| 国产亚洲欧美98| 国产精品久久久人人做人人爽| 一区二区三区国产精品乱码| 亚洲五月天丁香| 身体一侧抽搐| 色综合欧美亚洲国产小说| 免费高清视频大片| 日韩精品中文字幕看吧| 欧美日韩一级在线毛片| 99re在线观看精品视频| 亚洲精品av麻豆狂野| 久久久久久久精品吃奶| 日韩免费av在线播放| 美女扒开内裤让男人捅视频| 国产亚洲欧美在线一区二区| 国产精品99久久99久久久不卡| 国产精品美女特级片免费视频播放器 | 婷婷六月久久综合丁香| 国产精品一区二区在线不卡| 久久久精品欧美日韩精品| 又紧又爽又黄一区二区| x7x7x7水蜜桃| 午夜免费观看网址| 国产亚洲精品一区二区www| www.www免费av| 男女午夜视频在线观看| 91老司机精品| 男女之事视频高清在线观看| 丝袜在线中文字幕| 亚洲av电影不卡..在线观看| 国产精品九九99| 欧美不卡视频在线免费观看 | 亚洲成av人片免费观看| 一级a爱片免费观看的视频| 脱女人内裤的视频| 久久亚洲真实| 人人妻,人人澡人人爽秒播| 正在播放国产对白刺激| 国内毛片毛片毛片毛片毛片| 国产一区二区三区视频了| 色综合欧美亚洲国产小说| 国产黄a三级三级三级人| 亚洲色图av天堂| 国产精品秋霞免费鲁丝片| 欧美午夜高清在线| 无限看片的www在线观看| 好男人在线观看高清免费视频 | 亚洲五月天丁香| 婷婷精品国产亚洲av在线| 久久午夜综合久久蜜桃| 搡老妇女老女人老熟妇| 老司机靠b影院| 欧美另类亚洲清纯唯美| 淫秽高清视频在线观看| 一区二区三区激情视频| 自拍欧美九色日韩亚洲蝌蚪91| 国产一区二区三区综合在线观看| 老司机在亚洲福利影院| 国内精品久久久久久久电影| 亚洲av美国av| 久久久久国产精品人妻aⅴ院| 国产极品粉嫩免费观看在线| 免费av毛片视频| 无限看片的www在线观看| 久久久国产成人免费| 男女做爰动态图高潮gif福利片 | 女生性感内裤真人,穿戴方法视频| 亚洲一区二区三区不卡视频| 老鸭窝网址在线观看| 国产一级毛片七仙女欲春2 | 女人被躁到高潮嗷嗷叫费观| 国产极品粉嫩免费观看在线| av天堂在线播放| 女人精品久久久久毛片| 亚洲欧美日韩无卡精品| 亚洲国产精品久久男人天堂| 亚洲av美国av| 亚洲av成人一区二区三| 日本a在线网址| 黑人巨大精品欧美一区二区mp4| 99香蕉大伊视频| 午夜影院日韩av| www.熟女人妻精品国产| 免费高清视频大片| 91九色精品人成在线观看| 满18在线观看网站| 亚洲成人免费电影在线观看| 在线观看免费午夜福利视频| 久久国产乱子伦精品免费另类| 久久伊人香网站| 黄色视频不卡| 久久精品亚洲熟妇少妇任你| 国产成人影院久久av| 老熟妇仑乱视频hdxx| 国产单亲对白刺激| 1024香蕉在线观看| 在线国产一区二区在线| 性色av乱码一区二区三区2| 少妇 在线观看| 人人妻人人澡人人看|