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

    “嫦娥5號”登陸候選地Mons Rümker的光照與溫度特征分析*

    2020-06-30 12:14:10鐘振張騰張杰陳世國
    物理學(xué)報 2020年11期
    關(guān)鍵詞:風(fēng)化層熱流月球

    鐘振 張騰 張杰 陳世國

    1) (貴州師范大學(xué)物理與電子科學(xué)學(xué)院, 貴陽 550001)2) (貴州師范學(xué)院地理與資源學(xué)院, 貴陽 550018)(2020 年1 月16日收到; 2020 年3 月15日收到修改稿)

    即將開展的“嫦娥5號”探月任務(wù), 將使我國在月球上首次實現(xiàn)無人鉆井取樣. 考慮到實際探測活動與探測區(qū)域的光照和溫度有關(guān), 有必要對研究區(qū)域的光照和溫度特征進(jìn)行分析, 為此, 本文利用SPICE系統(tǒng)對“嫦娥5號”候選登陸區(qū)Mons Rümker高原的實時光照進(jìn)行計算. 發(fā)現(xiàn)求解的相對光強(qiáng)度分布與日本SELENE衛(wèi)星提供的早晨光照影像一致, 驗證了光照算法及計算程序的合理性. 以此為基礎(chǔ), 利用1維熱傳導(dǎo)模型, 對候選登陸區(qū)風(fēng)化層不同深度的溫度進(jìn)行仿真分析. 結(jié)果表明風(fēng)化層溫度在近表面區(qū)域受光照的影響較大, 隨著深度的增加, 光照影響逐漸減弱. 到達(dá)0.57 m深度時, 風(fēng)化層溫度不再變化. 為確保鉆井任務(wù)的開展, 實際鉆井作業(yè)應(yīng)考慮風(fēng)化層內(nèi)外溫度差異引起的應(yīng)力不均. 考慮“嫦娥5號”的實際鉆井深度遠(yuǎn)大于0.57 m, 應(yīng)能測量到常溫層的熱流值, 后續(xù)探月任務(wù)可考慮搭載熱流探測設(shè)備, 以促進(jìn)月球科學(xué)研究的發(fā)展.

    1 引 言

    根據(jù)我國探月工程計劃, 繼“嫦娥”1號、2號、3號和4號的成功實施, 我國即將開展探月工程第三期“回”[1-3]. 該階段計劃發(fā)射“嫦娥5號”月球車,在月表進(jìn)行軟著陸, 對月表風(fēng)化層進(jìn)行鉆井, 提取相關(guān)巖石樣品并返回地球, 屆時, 我國將成為繼蘇聯(lián)和美國之后, 全球第3個實現(xiàn)月球采樣并返回的國家[4]. 有關(guān)“嫦娥5號”登陸區(qū)的選擇, 文獻(xiàn)[5]認(rèn)為風(fēng)暴洋(oceanus procellarum)北面的Rümker區(qū)域, 曾經(jīng)擁有較長的火山活動, 形成了多重地質(zhì)單元, 具有不同的元素組成, 采集該區(qū)域的巖石樣本具有重要的地質(zhì)和科學(xué)意義. 該區(qū)域的Mons Rümker位于風(fēng)暴洋北部, 中心坐標(biāo)為(40°N, 58°W),擁有古火山活動時遺留的月溪、穹隆和熔巖流管道. 其地形起伏在1000—2000 m之間, 文獻(xiàn)[6-8]認(rèn)為它屬于火山中心, 文獻(xiàn)[9]認(rèn)為它屬于大型盾型火山. 文獻(xiàn)[10]對Mons Rümker的地形地貌,以及地質(zhì)特征進(jìn)行分析, 提出將該區(qū)域作為“嫦娥5號”登陸候選區(qū), 采集該區(qū)域的巖石樣品, 有助于理解月球火山作用特征及其演化歷史.

    文獻(xiàn)[11]對Mons Rümker高原的累積光照條件進(jìn)行研究, 發(fā)現(xiàn)3年時間擁有30%至53%的光照率, 且地基觀測系統(tǒng)的通信條件較好. 優(yōu)良的光照條件可以快速地提升月表溫度, 但由于月表風(fēng)化層具有良好的絕熱特性, 使得光照對風(fēng)化層內(nèi)部溫度的貢獻(xiàn)有限[12-15]. 盡管文獻(xiàn)[12-15]對風(fēng)化層溫度分布進(jìn)行了估算, 但太陽光照方向主要依據(jù)理論公式求得, 與現(xiàn)實情況存在一定的偏差. 另外,他們在估算光照引起的溫度分布時, 并沒有考慮地形對太陽光的遮擋, 而地形遮擋效應(yīng)對實際探測活動至關(guān)重要. 為此, 本文利用美國NASA (National Aeronautics and Space Administration)開發(fā)的SPICE(Spacecraft, Planet, Instrument, Cameramatrix, Events)系統(tǒng), 精確確定不同時刻的太陽光照方向. 結(jié)合文獻(xiàn)[16,17]提供的高分辨率月球數(shù)字地形模型DTM (digital terrain model), 估算顧及地形遮擋的實時太陽光照. 同時, 考慮月球內(nèi)部熱流, 對Mons Rümker高原風(fēng)化層的溫度分布進(jìn)行計算, 以期為“嫦娥5號”鉆井采樣活動, 以及后繼探月活動的開展提供一定程度的參考.

    2 地形模型和光照條件

    Mons Rümker高原的平均面積約為4000 km2,地形圖如圖1白色方框區(qū)域所示. 本文采用的數(shù)字地形模型, 來自LRO (lunar reconnaissance orbiter)的激光測高數(shù)據(jù)LOLA (lunar orbiter laser altimeter). LRO自2009年發(fā)射以來, 其激光測距載荷LOLA取得了高精度的月球全球地形數(shù)據(jù),促進(jìn)了月表永久陰影區(qū)及水冰的研究, 同時也確保了月球全球高精度大地參考框架的建立, 為后續(xù)載人登月提供了安全保障[18]. LRO目前仍處于在軌狀態(tài), 盡管后期拓展任務(wù)階段因調(diào)軌操作, 無法采集月球北半球地形數(shù)據(jù), 但激光測距載荷LOLA仍然在不斷地更新數(shù)據(jù), 地形格網(wǎng)數(shù)據(jù)已更新至2019年(https://pds-geosciences.wustl.edu/missions/lro/lola.htm). 本文采用的地形模型, 來自LOLA的格網(wǎng)數(shù)據(jù)LDEM_512_00N_45N_270_360 (1/512° × 1/512°), 其分辨率約為59.2 m ×59.2 m. 其中, 黑色方塊表示文獻(xiàn)[10]建議的“嫦娥5號”候選登陸點(303.34°E, 40.11°N).

    圖 1 Mons Rümker區(qū)域地形圖, 如白色方框所示, 其中黑色方塊表示文獻(xiàn)[10]建議的“嫦娥5號”登陸點(303.34°E,40.11°N)Fig. 1. Topography around Mons Rümker region, which is figured out with a white box. The black box indicates the candidate landing site of CE-5 proposed by reference [10],and this site is centered at (303.34°E, 40.11°N).

    如圖2所示, 為了判斷目標(biāo)點A是否有太陽光照射, 以及判斷A點是否被鄰區(qū)地形點B遮擋,可以通過點A的太陽高度角, 以及比較BC和DC的大小來判斷. 過點A作等效球面(圓弧型虛線所示), 與OB相交于點C, 點A和點B間的中心角為θ. 假定點A的太陽高度角為α, 它表示太陽入射光線與目標(biāo)點水平地面間的夾角. 若α ≤ 0,文獻(xiàn)[11]的研究表明, 對于高海拔區(qū)域, 太陽光有可能照射至目標(biāo)點A. 如圖1所示, 本文研究區(qū)域Mons Rümker高原, 其最大高程不超過—1.5 km,不屬于高海拔區(qū)域, 因此, 本文不考慮文獻(xiàn)[11]的高海拔特殊情況, 認(rèn)為α ≤ 0時, 太陽光不能照射至目標(biāo)點A. 當(dāng)α > 0時, 若DC < BC, 則目標(biāo)點A被鄰區(qū)地形點B遮擋, 太陽光無法到達(dá)點A; 當(dāng)DC≥BC時, 則不受點B的地形遮擋. 考慮月球表面地形高低的最大落差, 文獻(xiàn)[11]表明最大中心角θ = 8°, 通過不斷地減小中心角, 可以依次判斷距離目標(biāo)點A遠(yuǎn)近不同的地形遮擋效應(yīng). BC的大小可以根據(jù)數(shù)字地形模型求得, 依據(jù)文獻(xiàn)[11], 可得

    圖 2 太陽光照條件示意圖Fig. 2. Schematic of illumination condition.

    其中AC可以根據(jù)中心角θ和OA求得; OA可以根據(jù)數(shù)字地形模型DTM求得; α為實時太陽光高度角, 可以根據(jù)SPICE系統(tǒng)由行星歷表求得.

    3 熱傳導(dǎo)方程及邊界條件

    已知光照強(qiáng)度, 根據(jù)熱傳導(dǎo)方程, 可以估算風(fēng)化層的溫度變化. 月表風(fēng)化層的厚度, 就整個研究區(qū)域而言相對較小, 因此, 通常采用一維熱傳導(dǎo)模型來計算溫度分布[12-15]. 就溫度T和風(fēng)化層深度z而言, 一維模型的控制方程為

    其中, ρ和cp分別表示風(fēng)化層的密度和比熱容,k表示熱傳導(dǎo)系數(shù), 這些參數(shù)通常與風(fēng)化層的深度和溫度分布有關(guān). 根據(jù)文獻(xiàn)[13,14,19-21], 有

    式中相關(guān)參數(shù)的取值如表1所列.

    表 1 相關(guān)參數(shù)取值Table 1. Values of parameters used in study.

    為了求解(2)式, 需要指定兩個邊界條件, 一個在月表, 另一個在風(fēng)化層底部. 月表邊界取決于太陽光照及月表紅外輻射, 風(fēng)化層底部取決于月球內(nèi)部熱流分布. 對于深度為zb的底部邊界z = zb,有如下關(guān)系[14]:

    其中, Q表示風(fēng)化層底部的熱流值, 根據(jù)Apollo任務(wù)測量值[22], 取Q = 0.018 W·m—2.

    對于月表邊界溫度Ts, 假定太陽光照加熱率為Qs, 考慮紅外輻射后的月表溫度梯度有[13-15]:

    其中a, b為系數(shù), 參考文獻(xiàn)[13], 取a = 0.06, b = 0.25;A0表示月球表面的邦德反照率(又稱球面反照率),該參數(shù)用于描述太陽于正中天, 太陽光垂直入射時, 天體表面的反射系數(shù), 參考文獻(xiàn)[12], 本文取A0= 0.12. 太陽實時入射流Fsun和入射角β通常與月表計算點和太陽之間的實時距離r有關(guān), 該距離可以通過SPICE系統(tǒng)精確求得. 假定地球與太陽間的平均距離為r0, 對應(yīng)的太陽常數(shù)為S, 則入射流Fsun為

    其中, r0為日-地平均間距. 參考文獻(xiàn)[25], 取r0=1.49598261 × 1011m, S = 1361 W·m—2. 由(8)式和(9)式可知, 求取太陽光照加熱率, 關(guān)鍵在于入射角β的求取. 參考文獻(xiàn)[26], 入射角β與研究點的坡度λ、方位角γ、太陽天頂距θz、太陽方位角γs有關(guān), 關(guān)系為

    其中, 方位角γ可由研究點的坡向求出; 太陽天頂距θz及其方位角γs可由SPICE系統(tǒng)提供的實時太陽星下點坐標(biāo)以及文獻(xiàn)[26]的方法求出; 研究點的坡度λ可由文獻(xiàn)[17,18]提供的月球全球坡度分布圖求出.

    4 結(jié)果與分析

    為了測試算法與程序的合理性, 圖3給出了Mons Rümker高原的早晨光照分布. 該圖對應(yīng)月球地方時tm= 06:30:30, Mons Rümker高原在早晨相應(yīng)時刻, 地表對太陽光反射的相對光強(qiáng)度分布. 其中圖3(a)來自日本SELENE衛(wèi)星提供的早晨光照影像(http://darts.isas.jaxa.jp/planet/pdap/selene/), 由 影 像TCO_MAPm04_N42E 300N39E303SC和TCO_MAPm04_N42E303N39 E306SC拼接得到, 其中白色和黑色分別代表光照強(qiáng)度最強(qiáng)和最弱, 其他顏色表示光照強(qiáng)度在兩者間的變化. 圖3(a)表明太陽光自東向西, 東邊光照強(qiáng)度強(qiáng), 西邊光照弱; 迎光方向出現(xiàn)光亮區(qū)域, 而背光方向則由于地形遮擋出現(xiàn)陰暗區(qū)域. 采用與圖3(a)相同的光照時刻, 圖3(b)表示基于LOLA地形模型、(1)式和(10)式求解的相對光強(qiáng)度RII(relative intensity of illumination). 圖3(b)中等高線表示Mons Rümker高原, 相對半徑為1737.151 km的參考球面的地形高. 由圖3(b)可知: 1)東邊相對光強(qiáng)度強(qiáng), 西邊相對光強(qiáng)度弱;2)東邊迎光方向部分區(qū)域, 由于坡向正對太陽光,相對光強(qiáng)度較強(qiáng), 呈現(xiàn)出白色; 3)西邊背光方向部分區(qū)域, 由于地形遮擋, 相對光強(qiáng)度接近于零, 呈現(xiàn)出黑色. 對比圖3(a)和圖3(b), 可以發(fā)現(xiàn)除少部分小區(qū)域外, 兩者光照強(qiáng)度的分布總體上一致. 部分小區(qū)域出現(xiàn)偏差, 主要是由于光照計算采用的LOLA地形模型的分辨率, 低于日本SELENE衛(wèi)星提供的早晨光照影像所致. 本文計算的光強(qiáng)度分布, 在總體上與觀測結(jié)果一致, 說明整個算法及計算程序具有一定的合理性, 可以進(jìn)一步應(yīng)用于Mons Rümker高原風(fēng)化層溫度的估算.

    圖 3 Mons Rümker區(qū)域早晨光照圖, 對應(yīng)月球地方時tm =06:30:30 (a)日本SELENE衛(wèi)星提供的早晨光照圖; (b)本文計算的與圖3(a)相同時刻的實時光照圖Fig. 3. Morning map of illumination over Mons Rümker at the lunar local time tm = 06:30:30: (a) Japan’s SELENE morning map of illumination; (b) our estimated relative intensity of illumination at the same time of Fig. 3(a).

    利用本文計算的實時光照強(qiáng)度, 根據(jù)(2)—(7)式, 可以求出Mons Rümker高原風(fēng)化層的實時溫度分布. 為了計算風(fēng)化層不同深度處的溫度,利用文獻(xiàn)[14]提供的標(biāo)準(zhǔn)化有限差分方法對(2)式進(jìn)行數(shù)值計算. 圖4給出了Mons Rümker高原不同時刻的表面溫度, 其中圖4(a)對應(yīng)協(xié)調(diào)世界時UTC(universal time coordinated): 2020年10月28日11點30分00秒, 簡寫為UTC 2020-10-28T 11:30:00, 后文所有時間依此格式表述. 此時, 由于沒有太陽光照, 大部分區(qū)域的溫度在80 K左右. 至圖4(b)時(UTC 2020-10-29T 06:45:00),隨著太陽光的到來, 東邊溫度逐漸升高, 特別是坡向正對太陽光的區(qū)域, 溫度一度接近200 K, 西邊由于背光, 溫度維持在120 K左右. 隨著太陽高度角的增大, 類似情況也出現(xiàn)在圖4(c)中, 大部分區(qū)域的溫度升至200 K左右, 高溫區(qū)域一度接近280 K, 該圖對應(yīng)的時刻為UTC 2020-10-30T 06:45:00. 到達(dá)正午時分(如圖4(d)所示, 對應(yīng)時刻為UTC 2020-11-02T 04:45:00), 大部區(qū)域的溫度升至360 K左右, 少部分區(qū)域接近或超過400 K. 至圖4(e)時, 由于太陽高度角的下降, 光照減少, 表面溫度也逐漸下降. 此時, Mons Rümker高原處于月球地方時的下午時分, 太陽光照至西向東, 西邊迎光方向光照強(qiáng), 東邊背光方向光照弱,對應(yīng)時刻為UTC 2020-11-12T 02:45:00. 如圖4(f)所示(對應(yīng)時刻為UTC 2020-11-12T 17:00:00), 隨著太陽高度角的進(jìn)一步降低, 大部分區(qū)域的溫度降至120 K左右, 西邊迎光方向部分區(qū)域的溫度維持在200 K左右. 盡管溫度下降, 部分背光區(qū)域的溫度一度高于200 K, 這些區(qū)域恰好對應(yīng)圖4(b)—(d)中的高溫區(qū)域. 在表面光照強(qiáng)時, 這些高溫區(qū)域的熱量不斷地向風(fēng)化層底部傳遞, 至圖4(f)時, 盡管表面光照減弱, 已傳遞至風(fēng)化層內(nèi)部的熱量反向傳遞至風(fēng)化層表面, 使得表面溫度升高. 由于這些區(qū)域內(nèi)部溫度高于其他區(qū)域, 至圖4(f)時, 盡管光照減弱, 得到內(nèi)部熱量的補(bǔ)充, 表面溫度仍然高于其他區(qū)域.

    圖 4 Mons Rümker區(qū)域表面溫度分布圖 (a)對應(yīng)時刻為UTC 2020-10-28T 11:30:00; (b)對應(yīng)時刻為UTC 2020-10-29T 06:45:00; (c)對應(yīng)時刻為UTC 2020-10-30T 06:45:00; (d)對應(yīng)時刻為UTC 2020-11-02T 04:45:00; (e)對應(yīng)時刻為UTC 2020-11-12T 02:45:00; (f)對應(yīng)時刻為UTC 2020-11-12T 17:00:00Fig. 4. Surface temperature distribution with time over Mons Rümker plateau: (a) Time at UTC 2020-10-28T 11:30:00; (b) time at UTC 2020-10-29T 06:45:00; (c) time at UTC 2020-10-30T 06:45:00; (d) time at UTC 2020-11-02T 04:45:00; (e) time at UTC 2020-11-12T 02:45:00; (f) time at UTC 2020-11-12T 17:00:00.

    這種結(jié)論也可以由圖5得到佐證. 圖5表示Mons Rümker高原風(fēng)化層5 cm深度的溫度分布,對應(yīng)時刻與圖4各分圖一致. 由圖5(a)—5(c)可知, 隨著太陽光照增加, 盡管表面溫度在升高(如圖4(a)—4(c)所示), 但風(fēng)化層5 cm深度處的溫度始終維持在200 K左右, 至正午時分圖5(d)時, 底部溫度才有所升高. 這說明月球風(fēng)化層的導(dǎo)熱能力較弱, 具有一定的絕熱性能, 該結(jié)果與文獻(xiàn)[14,15]的一致. 隨著下午時分太陽光照的減弱, 盡管表面溫度下降至120 K左右(如圖4(e)和圖4(f)所示), 但5 cm深度的溫度始終維持在260 K左右,部分區(qū)域甚至達(dá)到300 K. 此時, 隨著表面溫度下降, 內(nèi)部溫度高于月表溫度, 熱量由內(nèi)向外傳遞,使得圖4(f)部分背光區(qū)域的溫度高于鄰近區(qū)域.

    為了進(jìn)一步探究風(fēng)化層的溫度變化, 圖6給出了溫度隨深度變化的剖面圖, 剖面方向沿圖1藍(lán)色直線的方向, 橫軸表示經(jīng)度, 縱軸表示深度, 對應(yīng)時刻與圖4和圖5一致. 由圖6(a)—(c)可知, 隨著太陽光照的增加, 風(fēng)化層5 cm深度以內(nèi)的溫度不斷地增加, 5 cm以下的溫度變化較弱, 這與圖5(a)—(c)的結(jié)果一致. 至正午時刻圖6(d)時,表面溫度升至最大值, 熱量向風(fēng)化層內(nèi)部傳遞, 內(nèi)部溫度不斷升高. 圖6(d)中黑色曲線表示沿剖面方向的表面地形輪廓, 其幅度表示地形的相對變化, 其大小與圖6縱軸刻度范圍無關(guān), 該圖表明溫度的變化與地形有關(guān). 至圖6(e)和圖6(f)時, 內(nèi)部儲存的熱量向表面?zhèn)鬟f, 使得表面對應(yīng)區(qū)域的溫度高于鄰近區(qū)域.

    為了進(jìn)一步研究0.2 m深度以下的溫度變化,圖7給出了兩個參考點不同深度處溫度隨時間的變化關(guān)系. 其中, 圖7(a)表示“嫦娥5號”候選登陸點(圖1所示黑色方塊)的溫度變化, 圖7(b)表示另一參考點(圖1所示黑色五角星)的溫度變化.由圖7可知, 兩個參考點的溫度在0.27 m深度的溫度最小值在239 K左右. 溫度變化幅度隨著深度的降低逐漸減小, 到達(dá)0.57 m深度時, 溫度幾乎不再變化, 保持在241.5 K左右. 結(jié)合圖4—圖7,保守估計風(fēng)化層的常溫層深度在0.6 m. 參考文獻(xiàn)[4], “嫦娥5號”的鉆井深度接近2 m, 在進(jìn)行鉆井作業(yè)及仿真分析時, 有必要考慮風(fēng)化層內(nèi)外溫度的差異; 另外, 考慮到常溫層深度在0.6 m左右,“嫦娥5號”的鉆井深度達(dá)2 m, 能探測到常溫層的熱流. 結(jié)合“嫦娥5號”的鉆井經(jīng)驗, 后期探月工程可考慮搭載探測月球內(nèi)部熱流值的載荷, 以促進(jìn)月球內(nèi)部結(jié)構(gòu)及熱演化研究的發(fā)展.

    圖 5 Mons Rümker區(qū)域底部5 cm深度的溫度分布圖 (a)對應(yīng)時刻為UTC 2020-10-28T 11:30:00; (b)對應(yīng)時刻為UTC 2020-10-29T 06:45:00; (c)對應(yīng)時刻為UTC 2020-10-30T 06:45:00; (d)對應(yīng)時刻為UTC 2020-11-02T 04:45:00; (e)對應(yīng)時刻為UTC 2020-11-12T 02:45:00; (f)對應(yīng)時刻為UTC 2020-11-12T 17:00:00Fig. 5. Subsurface temperature distribution at the depth of 5 cm over Mons Rümker plateau: (a) Time at UTC 2020-10-28T 11:30:00; (b) time at UTC 2020-10-29T 06:45:00; (c) time at UTC 2020-10-30T 06:45:00; (d) time at UTC 2020-11-02T 04:45:00;(e) time at UTC 2020-11-12T 02:45:00; (f) time at UTC 2020-11-12T 17:00:00.

    圖 6 剖面溫度沿圖1所示經(jīng)度方向的分布 (a)對應(yīng)時刻為UTC 2020-10-28T 11:30:00; (b)對應(yīng)時刻為UTC 2020-10-29T 06:45:00; (c)對應(yīng)時刻為UTC 2020-10-30T 06:45:00; (d)對應(yīng)時刻為UTC 2020-11-02T 04:45:00; (e)對應(yīng)時刻為UTC 2020-11-12T 02:45:00; (f)對應(yīng)時刻為UTC 2020-11-12T 17:00:00. 圖6(d)所示黑色曲線表示圖1藍(lán)線方向的表面地形輪廓Fig. 6. Temperature variation along the longitude direction shown in Fig.1: (a) Time at UTC 2020-10-28T 11:30:00; (b) time at UTC 2020-10-29T 06:45:00; (c) time at UTC 2020-10-30T 06:45:00; (d) time at UTC 2020-11-02T 04:45:00; (e) time at UTC 2020-11-12T 02:45:00; (f) time at UTC 2020-11-12T 17:00:00. The black carves in Fig. 6(d) represents the surface topography along the same blue line direction displayed in Fig. 1.

    5 討 論

    前文在估算Mons Rümker高原風(fēng)化層的溫度分布時, 控制風(fēng)化層密度分布的H參數(shù)取值為0.06 m, 熱流值Q選用的是Apollo 15號和Apollo 17號測量的平均值(Q = 0.018 W·m—2), 不同參數(shù)取值有可能引起溫度波動, 為此, 有必要討論相關(guān)參數(shù)的適用性. 參考文獻(xiàn)[14], 相關(guān)參數(shù)如熱傳導(dǎo)系數(shù)、風(fēng)化層密度和比熱容等參數(shù)主要受H參數(shù)的調(diào)控, 因此, 僅討論H參數(shù)對溫度的影響. 參考文獻(xiàn)[14], H參數(shù)的最小值為0.02 m, 最大值為0.09 m. 參考文獻(xiàn)[22], 熱流測量的最大值在Apollo 15號登陸點, 約為0.021 W·m—2. 文獻(xiàn)[27]利用數(shù)值方法, 估算的全球熱流值約為0.012 W·m—2, 小于Apollo任務(wù)的測量值. 為了直觀方便, 圖8和表2給出了H參數(shù)和熱流值Q取極值時, 參考點(如圖1黑色五角星所示)溫度隨深度的變化. 其中黑色實線和紅色虛線, 表示Q =0.018 W·m—2, H參數(shù)分別取0.02 m和0.09 m時,參考點溫度隨深度的變化; 綠色和藍(lán)色虛線表示H = 0.06 m, 熱流值分別取Q = 0.012 W·m—2和Q = 0.021 W·m—2時的溫度分布. 圖8(a)—8(c)分別對應(yīng)月球地方時tm= 06:30:30, tm= 12:30:30和tm= 18:30:30, 對應(yīng)參考點的早晨、正午和傍晚時刻. 1)由圖8(a)—8(c)可知, H參數(shù)和熱流值Q取值的不同, 對溫度的變化只有細(xì)微的影響, 絕對誤差不超過10 K; 2)如圖8(c)以及表2所示,H參數(shù)對溫度的影響超過熱流值Q, H參數(shù)主要對風(fēng)化層0.2 m以內(nèi)區(qū)域產(chǎn)生較大的影響, 對底層溫度的影響較弱; 3)如表2所列, 熱流值Q取值的不同, 主要對風(fēng)化層基底溫度產(chǎn)生較小的影響, 絕對誤差不超過1 K. 綜上所述, 同時考慮實際鉆進(jìn)深度超過1 m, H參數(shù)和熱流值Q的差異, 對溫度分布的影響可以忽略不計, 因此, 本文計算結(jié)果具有一定的參考價值.

    圖 8 參考點(圖1中黑色五角星)溫度隨深度的變化 (a)對應(yīng)月球地方時tm = 06:30:30; (b)對應(yīng)月球地方時tm = 12:30:30;(c)對應(yīng)月球地方時tm = 18:30:30Fig. 8. Subsurface temperature variations with depth for the point of black star in Fig. 1: (a) Temperature variations at the lunar local time tm = 06:30:30; (b) temperature variations at the lunar local time tm = 12:30:30; (c) temperature variations at the lunar local time tm = 18:30:30.

    表 2 參考點(圖1中黑色五角星)溫度(單位為K)在不同時刻隨深度的變化Table 2. Temperature (in K) variations with depth for the point of black star in Fig. 1 at various lunar local time.

    6 結(jié) 論

    針對“嫦娥5號”候選登陸區(qū)Mons Rümker高原風(fēng)化層的溫度變化, 利用SPICE系統(tǒng)對實時光照進(jìn)行計算, 在此基礎(chǔ)上, 結(jié)合1維熱傳導(dǎo)模型對風(fēng)化層的溫度進(jìn)行仿真分析. 結(jié)果表明: 1)本文計算的實時相對光照強(qiáng)度分布與日本SELENE衛(wèi)星提供的早晨光照影像總體一致, 說明本文的算法及計算程序具有一定的合理性; 2)依據(jù)求解的光照作為邊界條件, 發(fā)現(xiàn)0.5 cm深度以內(nèi)的風(fēng)化層溫度, 受表面光照的影響較大, 隨著深度的增加, 光照影響逐漸減弱; 3)風(fēng)化層底部0.57 m深度以下的溫度不隨時間變化, 保守估計常溫層的深度在0.6 m左右; 4)決定風(fēng)化層物理特征的H參數(shù), 其大小的改變并不會引起溫度的大幅波動, 對常溫層產(chǎn)生的影響較弱; 5)熱流值Q的不同, 僅對基底邊界溫度產(chǎn)生較小的影響, 實際應(yīng)用中可以不考慮熱流值Q選取的差異.

    猜你喜歡
    風(fēng)化層熱流月球
    到月球上“飆車”
    軍事文摘(2022年12期)2022-07-13 03:11:44
    巖質(zhì)邊坡風(fēng)化層的關(guān)鍵滑動面及其穩(wěn)定性研究
    近地表速度模型精度影響分析
    陪我去月球
    飛碟探索(2022年1期)2022-03-12 14:11:20
    月球上的另一個我
    軍事文摘(2021年18期)2021-12-02 01:28:54
    內(nèi)傾斜護(hù)幫結(jié)構(gòu)控釋注水漏斗熱流道注塑模具
    空調(diào)溫控器上蓋熱流道注塑模具設(shè)計
    流域風(fēng)化層雨洪調(diào)蓄生態(tài)機(jī)理
    聚合物微型零件的熱流固耦合變形特性
    中國塑料(2017年2期)2017-05-17 06:13:24
    坡度在巖石風(fēng)化層解譯中的應(yīng)用
    河北遙感(2015年2期)2015-07-18 11:11:14
    国产伦人伦偷精品视频| 精品一区二区三区视频在线观看免费 | 免费在线观看黄色视频的| 欧美日韩成人在线一区二区| 女警被强在线播放| 久久久精品国产亚洲av高清涩受| e午夜精品久久久久久久| 人成视频在线观看免费观看| 国产欧美日韩一区二区精品| 久久久久国产一级毛片高清牌| 夜夜夜夜夜久久久久| 18禁黄网站禁片午夜丰满| 成年人黄色毛片网站| 天天操日日干夜夜撸| 亚洲中文av在线| 不卡一级毛片| 高潮久久久久久久久久久不卡| 大片电影免费在线观看免费| 久久久水蜜桃国产精品网| 电影成人av| 国产单亲对白刺激| av不卡在线播放| ponron亚洲| 免费在线观看完整版高清| 丝袜美腿诱惑在线| 动漫黄色视频在线观看| 久99久视频精品免费| 中国美女看黄片| 老司机影院毛片| 亚洲国产毛片av蜜桃av| 日韩制服丝袜自拍偷拍| 日韩中文字幕欧美一区二区| 男男h啪啪无遮挡| 久久热在线av| 国产在线精品亚洲第一网站| 国产区一区二久久| 在线观看午夜福利视频| 一区二区日韩欧美中文字幕| 婷婷精品国产亚洲av在线 | 亚洲精品中文字幕一二三四区| 国产精品美女特级片免费视频播放器 | 一a级毛片在线观看| 俄罗斯特黄特色一大片| 久久久久国产一级毛片高清牌| 亚洲中文日韩欧美视频| 久久久国产欧美日韩av| 男女之事视频高清在线观看| 久久久国产一区二区| 午夜福利影视在线免费观看| 亚洲欧美色中文字幕在线| 午夜福利在线免费观看网站| 亚洲成国产人片在线观看| 天天躁日日躁夜夜躁夜夜| www.自偷自拍.com| 香蕉丝袜av| 亚洲人成伊人成综合网2020| 成人永久免费在线观看视频| av在线播放免费不卡| 亚洲七黄色美女视频| 在线av久久热| 成年版毛片免费区| 亚洲第一青青草原| 一级作爱视频免费观看| 99久久99久久久精品蜜桃| 999久久久国产精品视频| 一二三四社区在线视频社区8| 久久天躁狠狠躁夜夜2o2o| 亚洲精品国产一区二区精华液| 免费日韩欧美在线观看| 国产免费av片在线观看野外av| cao死你这个sao货| 久久ye,这里只有精品| 久久国产精品影院| 1024香蕉在线观看| 91国产中文字幕| 久久ye,这里只有精品| 水蜜桃什么品种好| 99re在线观看精品视频| 久久天堂一区二区三区四区| 露出奶头的视频| 天堂√8在线中文| 免费观看人在逋| 黄片大片在线免费观看| 两人在一起打扑克的视频| 在线av久久热| 免费人成视频x8x8入口观看| 免费人成视频x8x8入口观看| 国产国语露脸激情在线看| 女人精品久久久久毛片| 午夜成年电影在线免费观看| 亚洲精品国产色婷婷电影| 色尼玛亚洲综合影院| 777米奇影视久久| 国产精品秋霞免费鲁丝片| 在线观看免费日韩欧美大片| 亚洲人成伊人成综合网2020| 亚洲va日本ⅴa欧美va伊人久久| 国产一区有黄有色的免费视频| 亚洲国产精品sss在线观看 | 久久精品国产清高在天天线| 欧美人与性动交α欧美精品济南到| 国产片内射在线| 亚洲专区国产一区二区| 人人妻人人爽人人添夜夜欢视频| bbb黄色大片| 色综合婷婷激情| 女人高潮潮喷娇喘18禁视频| 免费人成视频x8x8入口观看| 成年女人毛片免费观看观看9 | 国产亚洲精品久久久久5区| 香蕉丝袜av| 人人妻人人澡人人爽人人夜夜| 视频区欧美日本亚洲| 精品久久久久久久毛片微露脸| 亚洲黑人精品在线| 久久精品成人免费网站| 99精品久久久久人妻精品| 老鸭窝网址在线观看| 国产午夜精品久久久久久| 欧美午夜高清在线| 国产精品.久久久| 18禁裸乳无遮挡免费网站照片 | 日本撒尿小便嘘嘘汇集6| 久久久久国内视频| 黄色视频不卡| 一区二区三区精品91| 午夜成年电影在线免费观看| 搡老熟女国产l中国老女人| 国产在视频线精品| 香蕉久久夜色| 丰满人妻熟妇乱又伦精品不卡| 久久精品国产99精品国产亚洲性色 | 一级作爱视频免费观看| 精品卡一卡二卡四卡免费| 不卡一级毛片| 成人18禁在线播放| 久久久久久久国产电影| 午夜福利欧美成人| 亚洲精品美女久久久久99蜜臀| 日本a在线网址| 精品熟女少妇八av免费久了| 麻豆国产av国片精品| 欧美丝袜亚洲另类 | 成人18禁高潮啪啪吃奶动态图| 在线十欧美十亚洲十日本专区| 国产精品影院久久| 1024香蕉在线观看| 日本黄色日本黄色录像| 天堂中文最新版在线下载| 一本一本久久a久久精品综合妖精| 久久中文看片网| 国产免费现黄频在线看| 天天操日日干夜夜撸| 女警被强在线播放| 一级a爱视频在线免费观看| 欧美黄色片欧美黄色片| 少妇的丰满在线观看| 九色亚洲精品在线播放| 成人免费观看视频高清| 男女高潮啪啪啪动态图| 日本五十路高清| 高清黄色对白视频在线免费看| 制服人妻中文乱码| 天天躁狠狠躁夜夜躁狠狠躁| 久久九九热精品免费| 久久国产亚洲av麻豆专区| 国产野战对白在线观看| 欧美日韩视频精品一区| 又黄又爽又免费观看的视频| 女性生殖器流出的白浆| 在线播放国产精品三级| videos熟女内射| 色精品久久人妻99蜜桃| 在线国产一区二区在线| 成年动漫av网址| 欧美乱妇无乱码| 在线观看免费高清a一片| 亚洲国产精品sss在线观看 | videosex国产| 亚洲午夜精品一区,二区,三区| 免费在线观看日本一区| 真人做人爱边吃奶动态| 水蜜桃什么品种好| 久久婷婷成人综合色麻豆| 成年动漫av网址| av网站免费在线观看视频| 18禁裸乳无遮挡动漫免费视频| 久久精品亚洲精品国产色婷小说| 香蕉久久夜色| 91成年电影在线观看| 高清视频免费观看一区二区| 亚洲熟女精品中文字幕| 黄片播放在线免费| 欧美不卡视频在线免费观看 | 天天操日日干夜夜撸| 亚洲av欧美aⅴ国产| √禁漫天堂资源中文www| 成年女人毛片免费观看观看9 | 曰老女人黄片| 黄网站色视频无遮挡免费观看| 69精品国产乱码久久久| 91老司机精品| 美女视频免费永久观看网站| 久久国产乱子伦精品免费另类| 精品高清国产在线一区| 欧美国产精品va在线观看不卡| 亚洲欧美日韩高清在线视频| 亚洲色图 男人天堂 中文字幕| 久久人妻av系列| 国产视频一区二区在线看| 国产欧美日韩一区二区三区在线| 国产一区二区三区综合在线观看| 十八禁人妻一区二区| 久久久久国产一级毛片高清牌| 免费观看精品视频网站| 国产精品98久久久久久宅男小说| 国产一区二区三区视频了| 午夜日韩欧美国产| 日日夜夜操网爽| 999久久久精品免费观看国产| 一二三四社区在线视频社区8| 亚洲国产精品一区二区三区在线| 亚洲国产欧美一区二区综合| 亚洲免费av在线视频| 高清毛片免费观看视频网站 | 天堂俺去俺来也www色官网| 夜夜躁狠狠躁天天躁| 九色亚洲精品在线播放| 午夜精品国产一区二区电影| a级毛片在线看网站| 色综合婷婷激情| 热99久久久久精品小说推荐| 欧美一级毛片孕妇| 女人高潮潮喷娇喘18禁视频| 国产精品美女特级片免费视频播放器 | 精品免费久久久久久久清纯 | 一个人免费在线观看的高清视频| 我的亚洲天堂| 一区二区日韩欧美中文字幕| 亚洲第一欧美日韩一区二区三区| 国产在线一区二区三区精| 日韩大码丰满熟妇| 激情在线观看视频在线高清 | 亚洲午夜精品一区,二区,三区| 亚洲熟妇熟女久久| 最近最新中文字幕大全电影3 | 亚洲人成电影免费在线| 丰满饥渴人妻一区二区三| 51午夜福利影视在线观看| 看免费av毛片| 亚洲av美国av| svipshipincom国产片| 嫩草影视91久久| 亚洲一区二区三区欧美精品| 国产精品国产高清国产av | 老司机在亚洲福利影院| 午夜亚洲福利在线播放| 亚洲国产日韩欧美精品在线观看 | 91在线观看av| 色综合站精品国产| 免费在线观看日本一区| 久久久久久久久久黄片| 国产久久久一区二区三区| 日韩欧美国产一区二区入口| 1024手机看黄色片| 日韩中文字幕欧美一区二区| 18美女黄网站色大片免费观看| 国产亚洲精品一区二区www| 亚洲精品美女久久久久99蜜臀| 亚洲不卡免费看| 国产三级在线视频| 在线免费观看不下载黄p国产 | 757午夜福利合集在线观看| 成年女人永久免费观看视频| 一级毛片女人18水好多| 99精品欧美一区二区三区四区| 校园春色视频在线观看| 五月伊人婷婷丁香| 又爽又黄无遮挡网站| 91字幕亚洲| 免费观看的影片在线观看| 手机成人av网站| 俄罗斯特黄特色一大片| 免费av观看视频| 一二三四社区在线视频社区8| 高清在线国产一区| 搡老岳熟女国产| 在线免费观看不下载黄p国产 | 在线观看免费视频日本深夜| 最近最新免费中文字幕在线| 国产乱人伦免费视频| 我要搜黄色片| 国产亚洲精品一区二区www| 国产黄a三级三级三级人| 尤物成人国产欧美一区二区三区| 久久婷婷人人爽人人干人人爱| 国产亚洲欧美98| 两个人的视频大全免费| 色精品久久人妻99蜜桃| 窝窝影院91人妻| 亚洲 欧美 日韩 在线 免费| 精品人妻一区二区三区麻豆 | 欧美一级毛片孕妇| 一级毛片女人18水好多| 午夜福利欧美成人| 桃红色精品国产亚洲av| xxx96com| 日韩免费av在线播放| 超碰av人人做人人爽久久 | 又黄又爽又免费观看的视频| 国产主播在线观看一区二区| 免费看十八禁软件| 欧美色视频一区免费| 内地一区二区视频在线| 国产精品国产高清国产av| 国产午夜福利久久久久久| 亚洲无线在线观看| 亚洲天堂国产精品一区在线| 香蕉av资源在线| 精品日产1卡2卡| 久久国产乱子伦精品免费另类| 中文字幕精品亚洲无线码一区| 久久精品亚洲精品国产色婷小说| av欧美777| 欧美黄色淫秽网站| 亚洲欧美日韩无卡精品| 亚洲一区高清亚洲精品| 日本与韩国留学比较| 国产麻豆成人av免费视频| av在线蜜桃| 亚洲五月婷婷丁香| 久久精品91蜜桃| 免费看十八禁软件| 国产亚洲精品久久久久久毛片| 国产精品女同一区二区软件 | 国产成人a区在线观看| 国产免费av片在线观看野外av| 校园春色视频在线观看| 热99在线观看视频| 国产av一区在线观看免费| eeuss影院久久| 性欧美人与动物交配| 亚洲精品一区av在线观看| 少妇人妻精品综合一区二区 | 午夜a级毛片| 亚洲av不卡在线观看| 亚洲美女黄片视频| 亚洲一区高清亚洲精品| 亚洲精品成人久久久久久| 少妇的丰满在线观看| 亚洲美女黄片视频| 国产v大片淫在线免费观看| 精品久久久久久,| 亚洲欧美激情综合另类| 国产三级在线视频| 中出人妻视频一区二区| 欧美av亚洲av综合av国产av| 精品午夜福利视频在线观看一区| 精品国产三级普通话版| 国产av一区在线观看免费| 黄色女人牲交| 最近最新中文字幕大全电影3| 香蕉av资源在线| 精品人妻一区二区三区麻豆 | 熟女少妇亚洲综合色aaa.| 啦啦啦免费观看视频1| 久久草成人影院| 搡老岳熟女国产| 国产v大片淫在线免费观看| 美女大奶头视频| 欧美一级毛片孕妇| 搡女人真爽免费视频火全软件 | 国产真实伦视频高清在线观看 | 亚洲五月婷婷丁香| 757午夜福利合集在线观看| 国产精品精品国产色婷婷| 国产高清三级在线| 少妇丰满av| 免费在线观看成人毛片| 精品久久久久久,| 成年女人毛片免费观看观看9| 丁香欧美五月| svipshipincom国产片| 久久精品夜夜夜夜夜久久蜜豆| 欧美另类亚洲清纯唯美| 99久久精品热视频| 国产 一区 欧美 日韩| 久久久久久久久久黄片| av专区在线播放| 国产亚洲欧美在线一区二区| 搡老熟女国产l中国老女人| 亚洲av一区综合| 亚洲欧美日韩东京热| 久久久久久久精品吃奶| 女同久久另类99精品国产91| 久久九九热精品免费| 麻豆久久精品国产亚洲av| 午夜精品久久久久久毛片777| 日韩免费av在线播放| 少妇人妻精品综合一区二区 | 国产精品一及| 脱女人内裤的视频| 露出奶头的视频| 欧美日韩国产亚洲二区| 免费人成视频x8x8入口观看| 熟女电影av网| 在线观看免费午夜福利视频| 午夜免费成人在线视频| 色综合欧美亚洲国产小说| 日本 av在线| 久久久久精品国产欧美久久久| 丰满的人妻完整版| 天堂网av新在线| 有码 亚洲区| 亚洲国产欧美网| 少妇裸体淫交视频免费看高清| 一个人看视频在线观看www免费 | 亚洲专区中文字幕在线| 国产亚洲欧美98| 97人妻精品一区二区三区麻豆| 熟女少妇亚洲综合色aaa.| 丰满的人妻完整版| 脱女人内裤的视频| av在线天堂中文字幕| 日本熟妇午夜| 国产高清有码在线观看视频| 好看av亚洲va欧美ⅴa在| 国产主播在线观看一区二区| 亚洲中文字幕一区二区三区有码在线看| 午夜免费观看网址| 非洲黑人性xxxx精品又粗又长| 国产老妇女一区| av欧美777| 十八禁人妻一区二区| 成人精品一区二区免费| 最新中文字幕久久久久| 欧美日韩一级在线毛片| 国产视频一区二区在线看| 免费看美女性在线毛片视频| 国内精品久久久久精免费| 丰满人妻一区二区三区视频av | 中文在线观看免费www的网站| 日韩亚洲欧美综合| 在线播放国产精品三级| 午夜a级毛片| 亚洲av电影不卡..在线观看| 亚洲欧美日韩无卡精品| 热99re8久久精品国产| 91久久精品电影网| 俄罗斯特黄特色一大片| 午夜视频国产福利| 亚洲精品粉嫩美女一区| 亚洲在线自拍视频| 婷婷六月久久综合丁香| 首页视频小说图片口味搜索| 国产乱人视频| 国产成人av教育| 国产高潮美女av| 亚洲欧美日韩高清专用| 国产视频一区二区在线看| 国产精品精品国产色婷婷| 一级毛片高清免费大全| 国产成年人精品一区二区| 国内久久婷婷六月综合欲色啪| 国产精品av视频在线免费观看| 亚洲人成网站在线播| 日本 av在线| 色尼玛亚洲综合影院| 黄色女人牲交| 国产综合懂色| 国产免费av片在线观看野外av| 12—13女人毛片做爰片一| 中亚洲国语对白在线视频| 久久精品国产综合久久久| 19禁男女啪啪无遮挡网站| 国产av麻豆久久久久久久| 亚洲18禁久久av| 色综合欧美亚洲国产小说| ponron亚洲| 久9热在线精品视频| 国产99白浆流出| av女优亚洲男人天堂| 国产免费一级a男人的天堂| 久久久国产成人免费| 日本黄大片高清| 国产淫片久久久久久久久 | 亚洲国产欧美网| 色综合站精品国产| 天堂av国产一区二区熟女人妻| 国产亚洲欧美在线一区二区| 午夜免费男女啪啪视频观看 | 欧美另类亚洲清纯唯美| av欧美777| 日本三级黄在线观看| 草草在线视频免费看| 尤物成人国产欧美一区二区三区| 在线免费观看的www视频| 淫妇啪啪啪对白视频| 亚洲av美国av| 免费人成在线观看视频色| 国产精品三级大全| av专区在线播放| 亚洲av日韩精品久久久久久密| 色综合婷婷激情| 观看美女的网站| 欧美激情久久久久久爽电影| 日本 欧美在线| 熟女电影av网| 麻豆国产av国片精品| 日韩精品青青久久久久久| 很黄的视频免费| 最好的美女福利视频网| 亚洲成人久久爱视频| 999久久久精品免费观看国产| 成人午夜高清在线视频| 中文资源天堂在线| 久久久久久大精品| 高清在线国产一区| 黄片小视频在线播放| 国产又黄又爽又无遮挡在线| 国产高潮美女av| 国内毛片毛片毛片毛片毛片| 一级毛片高清免费大全| 叶爱在线成人免费视频播放| 亚洲国产日韩欧美精品在线观看 | 精品久久久久久久久久免费视频| 哪里可以看免费的av片| 午夜视频国产福利| 欧美性猛交╳xxx乱大交人| 内射极品少妇av片p| 国产高清有码在线观看视频| 国语自产精品视频在线第100页| 欧美成人一区二区免费高清观看| 亚洲人成网站在线播放欧美日韩| 超碰av人人做人人爽久久 | 两个人的视频大全免费| 熟女人妻精品中文字幕| 国产av在哪里看| 成熟少妇高潮喷水视频| 久久久久久久久大av| 人妻丰满熟妇av一区二区三区| 国产精品三级大全| 国产精品亚洲美女久久久| 搡老妇女老女人老熟妇| 色老头精品视频在线观看| 嫩草影院入口| 人妻夜夜爽99麻豆av| 51国产日韩欧美| 丁香欧美五月| 校园春色视频在线观看| 国产精品,欧美在线| 精品国内亚洲2022精品成人| 噜噜噜噜噜久久久久久91| 国产精品99久久久久久久久| 一个人观看的视频www高清免费观看| 白带黄色成豆腐渣| 美女cb高潮喷水在线观看| 精品一区二区三区人妻视频| 99riav亚洲国产免费| 国产精品久久久久久久电影 | 91久久精品电影网| 免费在线观看成人毛片| 中文字幕人成人乱码亚洲影| 成人18禁在线播放| 亚洲av免费在线观看| 中文字幕久久专区| 国产成人福利小说| 在线播放无遮挡| 99国产综合亚洲精品| 成人永久免费在线观看视频| 亚洲国产精品成人综合色| 成人特级黄色片久久久久久久| 国产麻豆成人av免费视频| 在线视频色国产色| av黄色大香蕉| 九色成人免费人妻av| 国产伦在线观看视频一区| 国产精品久久久久久人妻精品电影| 日韩精品中文字幕看吧| 亚洲国产精品成人综合色| 亚洲电影在线观看av| 很黄的视频免费| 男人舔女人下体高潮全视频| 校园春色视频在线观看| 国产一区在线观看成人免费| 欧美国产日韩亚洲一区| 男女床上黄色一级片免费看| 亚洲中文日韩欧美视频| 欧美性感艳星| 亚洲成av人片在线播放无| 51午夜福利影视在线观看| 黄色成人免费大全| 国产亚洲精品久久久com| 国产精品综合久久久久久久免费| 亚洲欧美激情综合另类| 日韩 欧美 亚洲 中文字幕| 最近最新中文字幕大全电影3| 亚洲精品国产精品久久久不卡| 搡老熟女国产l中国老女人| 国产精品久久久久久人妻精品电影| 日本黄色视频三级网站网址| 国产成人av教育| 欧美在线黄色| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | svipshipincom国产片| 搡女人真爽免费视频火全软件 | 99精品在免费线老司机午夜| 一个人免费在线观看的高清视频| 99国产精品一区二区三区| 69人妻影院| 黄色成人免费大全| 麻豆国产av国片精品| 国产av在哪里看| 特大巨黑吊av在线直播| 免费av毛片视频|