• <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)如何處理
    99riav亚洲国产免费| 十八禁网站免费在线| 国产一区在线观看成人免费| 亚洲熟妇熟女久久| 国产亚洲欧美在线一区二区| 99国产综合亚洲精品| 亚洲黑人精品在线| 国产精华一区二区三区| 他把我摸到了高潮在线观看| 波多野结衣av一区二区av| 少妇粗大呻吟视频| av天堂在线播放| 国产视频一区二区在线看| 在线观看66精品国产| 国产精品亚洲一级av第二区| 在线观看66精品国产| 欧美日韩黄片免| 久久人人精品亚洲av| 精品一品国产午夜福利视频| 一区二区三区高清视频在线| 国产精品一区二区精品视频观看| 九色国产91popny在线| 桃色一区二区三区在线观看| 香蕉丝袜av| 男人舔女人下体高潮全视频| 精品一区二区三区av网在线观看| 欧美另类亚洲清纯唯美| 国产精品亚洲av一区麻豆| 午夜福利视频1000在线观看 | 免费一级毛片在线播放高清视频 | 中文字幕高清在线视频| 亚洲在线自拍视频| 91在线观看av| 亚洲色图av天堂| 亚洲精品中文字幕在线视频| 日韩 欧美 亚洲 中文字幕| 日韩三级视频一区二区三区| 色综合站精品国产| 最近最新中文字幕大全免费视频| 午夜亚洲福利在线播放| 国产精品久久久久久人妻精品电影| 91av网站免费观看| videosex国产| 精品久久久久久,| 亚洲电影在线观看av| 亚洲色图综合在线观看| 国产aⅴ精品一区二区三区波| 女人精品久久久久毛片| 日韩有码中文字幕| 欧美日韩福利视频一区二区| 亚洲色图 男人天堂 中文字幕| 亚洲性夜色夜夜综合| 日韩成人在线观看一区二区三区| 亚洲国产毛片av蜜桃av| 嫩草影视91久久| 国产精品影院久久| 色在线成人网| 97人妻精品一区二区三区麻豆 | 国产成+人综合+亚洲专区| 亚洲欧美激情在线| 国产私拍福利视频在线观看| 最近最新中文字幕大全免费视频| 老司机午夜十八禁免费视频| 99精品在免费线老司机午夜| 黄网站色视频无遮挡免费观看| 国产又色又爽无遮挡免费看| 露出奶头的视频| 99在线视频只有这里精品首页| svipshipincom国产片| 又黄又粗又硬又大视频| 成人三级黄色视频| 国产伦人伦偷精品视频| 国产国语露脸激情在线看| 91老司机精品| 日本黄色视频三级网站网址| 国产男靠女视频免费网站| 99精品欧美一区二区三区四区| 一个人观看的视频www高清免费观看 | 99国产精品99久久久久| 纯流量卡能插随身wifi吗| 欧美老熟妇乱子伦牲交| 日本a在线网址| 香蕉国产在线看| 国产高清videossex| 免费无遮挡裸体视频| 91成人精品电影| 国产精品一区二区三区四区久久 | 久久久久亚洲av毛片大全| 天堂√8在线中文| 91av网站免费观看| 久久久国产成人精品二区| 9191精品国产免费久久| 自线自在国产av| 欧美黑人欧美精品刺激| 欧美成狂野欧美在线观看| 国产精品永久免费网站| 久热爱精品视频在线9| avwww免费| 搡老妇女老女人老熟妇| 色综合欧美亚洲国产小说| 香蕉久久夜色| 国产真人三级小视频在线观看| 精品久久久久久久毛片微露脸| 麻豆国产av国片精品| 最新美女视频免费是黄的| 欧美乱码精品一区二区三区| 美女大奶头视频| 亚洲无线在线观看| 黑人操中国人逼视频| 老熟妇乱子伦视频在线观看| 黄色女人牲交| 一个人观看的视频www高清免费观看 | 啦啦啦 在线观看视频| 亚洲五月婷婷丁香| 国产激情欧美一区二区| 日韩精品青青久久久久久| 精品一品国产午夜福利视频| cao死你这个sao货| 桃色一区二区三区在线观看| 美国免费a级毛片| 亚洲成人久久性| 99国产精品免费福利视频| 男人的好看免费观看在线视频 | 最近最新中文字幕大全免费视频| 人妻久久中文字幕网| 91九色精品人成在线观看| 日韩精品青青久久久久久| 嫩草影院精品99| 午夜老司机福利片| 满18在线观看网站| 美女扒开内裤让男人捅视频| 熟妇人妻久久中文字幕3abv| 宅男免费午夜| 美女 人体艺术 gogo| 国产精品一区二区免费欧美| 天天躁夜夜躁狠狠躁躁| 亚洲人成77777在线视频| 操美女的视频在线观看| 免费无遮挡裸体视频| 欧美激情 高清一区二区三区| 波多野结衣一区麻豆| 精品国产乱子伦一区二区三区| 此物有八面人人有两片| 久久亚洲真实| 国产99久久九九免费精品| 国产欧美日韩一区二区三区在线| 51午夜福利影视在线观看| 欧美中文日本在线观看视频| 精品国产亚洲在线| 国产精品精品国产色婷婷| 性少妇av在线| 麻豆成人av在线观看| 99re在线观看精品视频| 久久九九热精品免费| 国产精品野战在线观看| 国产av在哪里看| 大陆偷拍与自拍| 成年女人毛片免费观看观看9| 国产精品98久久久久久宅男小说| 久久香蕉国产精品| 久久久久国产一级毛片高清牌| 成人国产综合亚洲| 国产一卡二卡三卡精品| 欧美丝袜亚洲另类 | 免费高清视频大片| 成年女人毛片免费观看观看9| 在线观看66精品国产| 长腿黑丝高跟| 成人三级做爰电影| 亚洲全国av大片| 高清黄色对白视频在线免费看| 级片在线观看| 亚洲美女黄片视频| tocl精华| 久久精品亚洲熟妇少妇任你| 纯流量卡能插随身wifi吗| 韩国精品一区二区三区| 免费在线观看亚洲国产| 一卡2卡三卡四卡精品乱码亚洲| 又黄又粗又硬又大视频| 97人妻精品一区二区三区麻豆 | 免费在线观看视频国产中文字幕亚洲| 国产一区二区三区视频了| 亚洲视频免费观看视频| 日本a在线网址| 亚洲av电影不卡..在线观看| 亚洲午夜精品一区,二区,三区| 亚洲 国产 在线| 午夜福利,免费看| av网站免费在线观看视频| 亚洲五月婷婷丁香| 操美女的视频在线观看| 国产成人av激情在线播放| 一边摸一边抽搐一进一出视频| 亚洲精品在线美女| 村上凉子中文字幕在线| 1024香蕉在线观看| 亚洲精品美女久久av网站| 亚洲国产精品成人综合色| 欧美色视频一区免费| 99久久精品国产亚洲精品| 精品高清国产在线一区| 嫩草影视91久久| 精品一品国产午夜福利视频| 日本a在线网址| 三级毛片av免费| 欧美精品啪啪一区二区三区| videosex国产| 国产精品一区二区免费欧美| 在线观看免费视频网站a站| 欧美人与性动交α欧美精品济南到| 精品国产乱子伦一区二区三区| 中文亚洲av片在线观看爽| 午夜福利一区二区在线看| 88av欧美| 69精品国产乱码久久久| 亚洲第一欧美日韩一区二区三区| 亚洲第一电影网av| 夜夜躁狠狠躁天天躁| 国产精品98久久久久久宅男小说| 久久久国产欧美日韩av| 久久久久久久午夜电影| 此物有八面人人有两片| 99热只有精品国产| 大香蕉久久成人网| 丰满人妻熟妇乱又伦精品不卡| 成人国产综合亚洲| 国产av一区二区精品久久| 精品久久久久久久毛片微露脸| 一区在线观看完整版| 亚洲第一欧美日韩一区二区三区| 热re99久久国产66热| 九色亚洲精品在线播放| 动漫黄色视频在线观看| 日本精品一区二区三区蜜桃| 亚洲五月婷婷丁香| 久久青草综合色| 久久久久亚洲av毛片大全| 丝袜美腿诱惑在线| 女人被狂操c到高潮| 亚洲国产日韩欧美精品在线观看 | 日韩国内少妇激情av| 国产精品香港三级国产av潘金莲| 电影成人av| 日本在线视频免费播放| 女同久久另类99精品国产91| 亚洲,欧美精品.| 无人区码免费观看不卡| 国产成人精品久久二区二区免费| 精品福利观看| 免费在线观看完整版高清| 欧美日韩福利视频一区二区| 91九色精品人成在线观看| 乱人伦中国视频| 女性生殖器流出的白浆| 国产亚洲av嫩草精品影院| 成人18禁在线播放| 精品福利观看| 在线视频色国产色| 亚洲精品美女久久av网站| 久久精品91无色码中文字幕| 女人精品久久久久毛片| 久久久久久久久中文| 国产午夜福利久久久久久| 国产私拍福利视频在线观看| 中国美女看黄片| 欧美成人免费av一区二区三区| netflix在线观看网站| 欧美一区二区精品小视频在线| 99热只有精品国产| 精品国产一区二区久久| 日韩大尺度精品在线看网址 | av片东京热男人的天堂| 午夜两性在线视频| 国产亚洲欧美在线一区二区| 国产97色在线日韩免费| 丝袜美足系列| 亚洲欧美激情在线| 久久午夜综合久久蜜桃| 国产精品1区2区在线观看.| 精品人妻1区二区| 国产精品二区激情视频| 午夜精品在线福利| 亚洲天堂国产精品一区在线| 亚洲专区字幕在线| 最好的美女福利视频网| 国产精品爽爽va在线观看网站 | 国产蜜桃级精品一区二区三区| 日韩欧美免费精品| 女人精品久久久久毛片| 成人永久免费在线观看视频| 亚洲成人精品中文字幕电影| 他把我摸到了高潮在线观看| 嫁个100分男人电影在线观看| 免费一级毛片在线播放高清视频 | 国产精品美女特级片免费视频播放器 | 国产亚洲欧美精品永久| 亚洲中文字幕一区二区三区有码在线看 | 国内久久婷婷六月综合欲色啪| 嫩草影视91久久| 黄频高清免费视频| 成人18禁在线播放| 纯流量卡能插随身wifi吗| 久久精品人人爽人人爽视色| 热99re8久久精品国产| av电影中文网址| av中文乱码字幕在线| 琪琪午夜伦伦电影理论片6080| 欧美色视频一区免费| 久久国产亚洲av麻豆专区| 亚洲少妇的诱惑av| 亚洲av第一区精品v没综合| 国内久久婷婷六月综合欲色啪| www日本在线高清视频| 国产精品久久久久久精品电影 | 国产麻豆69| 一本大道久久a久久精品| 国产不卡一卡二| 男人的好看免费观看在线视频 | 亚洲黑人精品在线| 一级毛片女人18水好多| 久热这里只有精品99| 国产亚洲av嫩草精品影院| 亚洲国产日韩欧美精品在线观看 | 亚洲成人久久性| 国产一区二区激情短视频| 老司机午夜福利在线观看视频| 悠悠久久av| 成人亚洲精品一区在线观看| 大码成人一级视频| 中文亚洲av片在线观看爽| 欧美日本中文国产一区发布| 亚洲欧美激情综合另类| 亚洲自偷自拍图片 自拍| 日韩欧美一区二区三区在线观看| 99精品久久久久人妻精品| 亚洲久久久国产精品| 久久这里只有精品19| 国产乱人伦免费视频| 国产黄a三级三级三级人| 亚洲成a人片在线一区二区| 亚洲欧美精品综合一区二区三区| 99在线人妻在线中文字幕| 操美女的视频在线观看| 精品不卡国产一区二区三区| 国产午夜精品久久久久久| 桃色一区二区三区在线观看| 一边摸一边做爽爽视频免费| АⅤ资源中文在线天堂| 午夜福利成人在线免费观看| 亚洲成人免费电影在线观看| 成人亚洲精品一区在线观看| 色播在线永久视频| 亚洲九九香蕉| 亚洲一区二区三区不卡视频| 91av网站免费观看| 欧美日韩瑟瑟在线播放| 夜夜爽天天搞| 亚洲五月色婷婷综合| 久久天躁狠狠躁夜夜2o2o| 亚洲av片天天在线观看| 叶爱在线成人免费视频播放| 免费一级毛片在线播放高清视频 | 9191精品国产免费久久| 99国产综合亚洲精品| 欧美日本视频| 午夜福利免费观看在线| 99国产精品99久久久久| 亚洲性夜色夜夜综合| 亚洲午夜理论影院| 91av网站免费观看| 色老头精品视频在线观看| 俄罗斯特黄特色一大片| 99国产精品一区二区三区| 国产精品爽爽va在线观看网站 | 亚洲最大成人中文| 免费搜索国产男女视频| 免费无遮挡裸体视频| 国产精品久久久久久人妻精品电影| 亚洲一区中文字幕在线| 久久久久久久久久久久大奶| 欧美日韩瑟瑟在线播放| 日韩av在线大香蕉| 在线观看免费视频日本深夜| 日韩高清综合在线| 久久婷婷成人综合色麻豆| 69精品国产乱码久久久| 日韩欧美国产一区二区入口| 亚洲一区高清亚洲精品| 久久狼人影院| 日韩欧美国产在线观看| 天堂动漫精品| 一级毛片高清免费大全| 国内毛片毛片毛片毛片毛片| 国产精品国产高清国产av| 琪琪午夜伦伦电影理论片6080| 一级毛片高清免费大全| 丝袜美腿诱惑在线| 美女国产高潮福利片在线看| 一边摸一边抽搐一进一小说| 亚洲精品中文字幕在线视频| 99久久国产精品久久久| 无遮挡黄片免费观看| 久久久久国产一级毛片高清牌| 午夜福利欧美成人| 伦理电影免费视频| 十八禁网站免费在线| 窝窝影院91人妻| 久久久国产精品麻豆| 久久精品国产综合久久久| 视频区欧美日本亚洲| 黑人操中国人逼视频| 免费观看人在逋| 午夜日韩欧美国产| 久久精品国产亚洲av高清一级| 美女国产高潮福利片在线看| 亚洲人成电影观看| 午夜视频精品福利| 日本一区二区免费在线视频| 亚洲成人精品中文字幕电影| 熟女少妇亚洲综合色aaa.| 国产精品久久久av美女十八| 成人18禁在线播放| 91精品三级在线观看| 看片在线看免费视频| 国产欧美日韩精品亚洲av| 精品久久久久久久人妻蜜臀av | 看免费av毛片| 成在线人永久免费视频| 色老头精品视频在线观看| 麻豆av在线久日| 精品国产超薄肉色丝袜足j| 别揉我奶头~嗯~啊~动态视频| 成人国产一区最新在线观看| 久久草成人影院| 精品一品国产午夜福利视频| 国产在线精品亚洲第一网站| 午夜福利免费观看在线| www日本在线高清视频| 色精品久久人妻99蜜桃| 国产精品99久久99久久久不卡| 黄片大片在线免费观看| 999久久久国产精品视频| 精品久久久久久久人妻蜜臀av | 两性夫妻黄色片| 国产精品乱码一区二三区的特点 | 精品高清国产在线一区| 免费在线观看完整版高清| 欧美丝袜亚洲另类 | 精品国产美女av久久久久小说| 桃色一区二区三区在线观看| 久久香蕉精品热| 日韩高清综合在线| 亚洲电影在线观看av| 大香蕉久久成人网| 在线观看www视频免费| 国产精品精品国产色婷婷| 9191精品国产免费久久| 视频区欧美日本亚洲| 中文字幕av电影在线播放| 国产蜜桃级精品一区二区三区| 老司机福利观看| 18禁观看日本| 久久久精品国产亚洲av高清涩受| 成人国产一区最新在线观看| 久久久国产成人免费| 久热爱精品视频在线9| 色播亚洲综合网| 激情在线观看视频在线高清| 一夜夜www| 亚洲午夜精品一区,二区,三区| 高清黄色对白视频在线免费看| 国产亚洲精品一区二区www| 男人舔女人的私密视频| 老鸭窝网址在线观看| 久久精品国产综合久久久| 亚洲精品粉嫩美女一区| 午夜久久久久精精品| 日韩 欧美 亚洲 中文字幕| 在线播放国产精品三级| 麻豆久久精品国产亚洲av| 男人的好看免费观看在线视频 | 男人的好看免费观看在线视频 | 99精品久久久久人妻精品| 日本黄色视频三级网站网址| 1024视频免费在线观看| 欧美国产精品va在线观看不卡| 国产精品av久久久久免费| 热re99久久国产66热| 欧美日本中文国产一区发布| 亚洲国产欧美网| 亚洲黑人精品在线| 12—13女人毛片做爰片一| 两人在一起打扑克的视频| 午夜福利欧美成人| 两个人看的免费小视频| 久久中文字幕人妻熟女| 国产视频一区二区在线看| 一进一出抽搐动态| 男女下面进入的视频免费午夜 | 色综合婷婷激情| 满18在线观看网站| 精品国产超薄肉色丝袜足j| aaaaa片日本免费| 欧美日本亚洲视频在线播放| 久久久国产成人精品二区| 久久国产乱子伦精品免费另类| 黄色丝袜av网址大全| 欧美日韩中文字幕国产精品一区二区三区 | 天天躁夜夜躁狠狠躁躁| 校园春色视频在线观看| 中文字幕人妻丝袜一区二区| 欧美乱色亚洲激情| a在线观看视频网站| 久久欧美精品欧美久久欧美| 丰满人妻熟妇乱又伦精品不卡| 国产三级黄色录像| 精品福利观看| 18禁黄网站禁片午夜丰满| 一a级毛片在线观看| 人妻丰满熟妇av一区二区三区| 国产亚洲精品久久久久5区| 精品不卡国产一区二区三区| 性少妇av在线| 宅男免费午夜| 一本综合久久免费| 禁无遮挡网站| 亚洲自拍偷在线| 中文字幕最新亚洲高清| 日韩视频一区二区在线观看| 一级黄色大片毛片| 午夜老司机福利片| 丰满的人妻完整版| 久久精品国产综合久久久| 又大又爽又粗| 91成人精品电影| 亚洲情色 制服丝袜| 黑丝袜美女国产一区| 亚洲精品一卡2卡三卡4卡5卡| 国产国语露脸激情在线看| 久久精品aⅴ一区二区三区四区| 黄色毛片三级朝国网站| 亚洲黑人精品在线| 97碰自拍视频| 级片在线观看| 久久中文字幕人妻熟女| 激情视频va一区二区三区| e午夜精品久久久久久久| 在线永久观看黄色视频| 中文字幕av电影在线播放| 夜夜爽天天搞| 久久久久亚洲av毛片大全| 亚洲天堂国产精品一区在线| 精品国产乱码久久久久久男人| 动漫黄色视频在线观看| 久久久久国产一级毛片高清牌| 长腿黑丝高跟| 老司机在亚洲福利影院| 精品一区二区三区四区五区乱码| 久久婷婷人人爽人人干人人爱 | 母亲3免费完整高清在线观看| 免费在线观看影片大全网站| 欧美日韩亚洲国产一区二区在线观看| 国产精品日韩av在线免费观看 | 中文字幕人妻熟女乱码| 久久亚洲精品不卡| 看黄色毛片网站| 国产成人影院久久av| 亚洲成人精品中文字幕电影| 一级毛片女人18水好多| 夜夜爽天天搞| 精品人妻1区二区| xxx96com| 男女下面插进去视频免费观看| 18禁观看日本| 午夜久久久久精精品| 欧美激情高清一区二区三区| 老司机靠b影院| 亚洲男人的天堂狠狠| 久久久久久大精品| 日韩中文字幕欧美一区二区| 91麻豆av在线| videosex国产| 黄网站色视频无遮挡免费观看| 老司机午夜福利在线观看视频| 亚洲熟妇熟女久久| 欧美激情久久久久久爽电影 | 久久国产精品男人的天堂亚洲| 中文字幕久久专区| av超薄肉色丝袜交足视频| 国产男靠女视频免费网站| 我的亚洲天堂| 男女下面插进去视频免费观看| 亚洲av成人不卡在线观看播放网| 法律面前人人平等表现在哪些方面| 久久久国产精品麻豆| 成人特级黄色片久久久久久久| 99久久综合精品五月天人人| 啦啦啦 在线观看视频| 久久国产精品影院| 极品人妻少妇av视频| 中文字幕人妻丝袜一区二区| 欧美日韩精品网址| 十分钟在线观看高清视频www| 变态另类丝袜制服| 一级黄色大片毛片| 久久草成人影院| 久久天堂一区二区三区四区| 国产激情欧美一区二区| 老汉色av国产亚洲站长工具| 日韩av在线大香蕉| 国产亚洲精品久久久久久毛片| 真人一进一出gif抽搐免费| 麻豆久久精品国产亚洲av| 真人一进一出gif抽搐免费|