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

    中低緯度電離層偶發(fā)E層電波傳播建模?

    2017-08-09 07:34:48郝書吉張文超張雅彬楊巨濤馬廣林
    物理學(xué)報 2017年11期
    關(guān)鍵詞:電子密度電波電離層

    郝書吉 張文超張雅彬 楊巨濤 馬廣林

    (中國電波傳播研究所,電波環(huán)境特性及模化技術(shù)重點實驗室,青島 266107)

    中低緯度電離層偶發(fā)E層電波傳播建模?

    郝書吉 張文超?張雅彬 楊巨濤 馬廣林

    (中國電波傳播研究所,電波環(huán)境特性及?;夹g(shù)重點實驗室,青島 266107)

    (2017年1月4日收到;2017年3月8日收到修改稿)

    基于電離層偶發(fā)E層(sporadic-E,Es)對電波傳播的多條鏈路實測數(shù)據(jù),認(rèn)為電離層Es對入射電波的作用是反射和散射兩種機制,且反射/散射的比例常數(shù)隨著電離層Es臨界頻率的變化而變化,進而建立了包含反射、散射以及反射與散射共同作用的三段式電離層Es電波傳播模型,并與國際電聯(lián)(International Telecommunications Union,ITU)給出的Es層傳播的電波場強預(yù)測模型進行了對比,驗證了本文所建模型的正確性.該模型尤其適用于中低緯地區(qū)甚高頻信號經(jīng)Es層的傳播研究.

    電離層Es,電波傳播,ITU模型

    1 引言

    偶發(fā)E層(sporadic-E,Es)是電離層內(nèi)一種短暫而不規(guī)則的強電離薄層,其高度范圍為90—140 km,厚度從數(shù)百米至1 km不等,水平尺度從幾十千米到幾百千米不等[1].電離層Es電子密度非常高,可達(dá)到常規(guī)E層電子密度的100倍[2].利用Es層可以實現(xiàn)單跳2000 km的超短波超視距通信,這對于建立超短波超視距應(yīng)急通信保障系統(tǒng)具有重要意義.

    研究人員從20世紀(jì)40年代就開始關(guān)注電離層Es對入射電波的影響及其信道特性,主要基于電離層探測儀,通過甚高頻(very high frequency,VHF)后向和斜入射鏈路對不同緯度的Es層進行觀測和統(tǒng)計研究.但由于觀測手段和統(tǒng)計方法的不同,結(jié)論也不盡相同.早在1964年,日本學(xué)者Kobayashi[3]通過長1100 km的VHF鏈路觀測,初步認(rèn)為電離層Es的反射損耗與電波入射角(θ)和電離層Es臨界頻率(foEs)的關(guān)系為Γ=25[f cos(θ)/foEs]2.67,其中θ為高度105 km的鏡反射入射角;1972年,Bramley[4]認(rèn)為電離層Es中VHF信號傳播的反射系數(shù)與(foEs/f)j成正比,其中j的值在7—30之間;1972年,Althouse和Davis[5]研究了Es層的散射特性,并對其后向散射進行了評估;1978年,日本的Miya和Shimizu[6]通過Tokyo-Wakkanai鏈路觀測,得到經(jīng)電離層Es的電波傳播經(jīng)驗公式;1980年,Kerblay和Makarenko[7]研究了Es層的大尺度水平不均勻性對反射系數(shù)的影響,并給出了中間點位置的反射系數(shù);1999年,國際電聯(lián)(International Telecommunications Union,ITU)在ITU-R p.1240建議書中給出了Es的臨界頻率發(fā)生概率分布圖,可以預(yù)測世界各地Es的發(fā)生;2007年,Sherstyukov和Akchurin[8]對中緯度Es層單跳傳播方式建模,用來預(yù)測Es層傳播的可用頻率分布;2012年,ITU發(fā)布了計算經(jīng)過Es層電波場強的建議書,在建議書中給出了經(jīng)過Es傳播的電離層Es衰減經(jīng)驗公式,其衰減值是入射電波頻率與電離層臨界頻率比值的二次方成正比;2015年,Sherstyukov等[9,10]通過對高頻(HF)和VHF頻段的無線電波測量數(shù)據(jù)的統(tǒng)計分析,也獲得了電離層Es的反射損耗.

    上述研究主要從反射角度并基于大量實驗觀測數(shù)據(jù)研究電離層Es對入射電波的作用,從電離層Es對電波的反射系數(shù)和傳播損耗方面進行了分析和建模研究.本文基于多條鏈路實驗數(shù)據(jù),認(rèn)為電離層Es對入射電波的作用很難用單一的反射或散射作用機制進行解釋,結(jié)合中低緯地區(qū)電離層Es電子密度分布特性及其中各向異性結(jié)構(gòu)的尺度關(guān)系,判斷電離層Es對入射電波的作用可能是散射作用和反射作用兩者機制都存在,并且反射作用的比例常數(shù)隨著電離層Es臨界頻率的變化而變化.據(jù)此,根據(jù)不同鏈路上的測量數(shù)據(jù),建立了包含反射、散射以及反射與散射共同作用的三段式電離層Es電波傳播模型,并與ITU模型進行了對比,驗證了本文所建模型的正確性.

    2 電離層Es電波傳播理論

    理論和實驗研究表明:電離層Es電子密度非常高,會對入射的VHF波段信號產(chǎn)生反射作用,同時Es層內(nèi)部存在大量電子密度不規(guī)則體,會對入射的信號產(chǎn)生散射作用.

    2.1 反射理論

    理論和試驗研究表明[2,11]:電離層Es層電子密度梯度在垂直方向上相對比較陡峭,據(jù)此特點,可假定Es層是由水平方向上為有限尺度的、垂直方向上具有不同介電常數(shù)的無數(shù)個薄介電層組成,針對每一個薄介電層,其VHF信號經(jīng)過Es層的反射示意圖見圖1.其中d為發(fā)射點和接收點的距離;θ為掠射角;a和b為兩個方向上的水平尺度,(λ是入射波波長).若相鄰層的介電常數(shù)的增量為dε,且滿足1?θ2?dε,則薄層間的反射系數(shù)為dρ=dε/(4θ2).

    圖1 薄層反射幾何示意Fig.1.Geometry of a thin layer model.

    目前,普遍認(rèn)為“風(fēng)剪切”理論是中緯度Es層形成的物理機制.設(shè)Es層電子密度在“風(fēng)剪切”的作用下其垂直分布滿足如下分布[12]:

    其中,Nmax是電離層Es最大電子密度;?h是Es層半厚度;z′=z?h,h為Es的中心高度,|z′|6?h;n=1,2,3···.通常Es層介電常數(shù)ε可寫為[13]

    其中fN是Es層最大等離子體頻率,fN=√是自由空間介電常數(shù),f是入射波頻率,w=2πf;e是電子電量,me是電子質(zhì)量.根據(jù)反射理論,總反射系數(shù)為

    其中λ為入射波波長;L=4πsinθ?h/λ.根據(jù)雷達(dá)方程,則經(jīng)過Es層反射的VHF信號接收功率為

    其中Pt為發(fā)射功率,Gt和Gr分別為發(fā)射和接收天線在指向Es高度方向的增益,D為傳播路徑.

    2.2 散射理論

    與對流層或常規(guī)F層不規(guī)則體相比[14,15],電離層Es不規(guī)則體尺度體要大得多[16,17],因此可運用其電子密度自相關(guān)函數(shù)來分析其散射現(xiàn)象[18].一般Es層電子密度自相關(guān)函數(shù)為

    其中:C(0)是大氣折射系數(shù)的時間方差,在Es中C(0)=(λ/λN)4(?N/N)/4;Γ為伽馬函數(shù),γ為Es不規(guī)則體尺度,l為一般湍流的尺度,p為常數(shù)因子,Km(p)為二類高階修正貝塞爾函數(shù).由于電離層Es電子密度通常比一般散射體的要大,用高階修正貝塞爾函數(shù)可更好地反映其散射特性[19].

    根據(jù)散射理論[20],電離層Es不規(guī)則體散射截面積為

    其中q=2klsinθ;χ為散射方向與入射電場方向的夾角;L1,L2是Es層不規(guī)則體在水平方向上平行、正交于水平漂移的不同尺度,L3是Es不規(guī)則體在垂直方向的尺度.

    根據(jù)(5)式,Es層不規(guī)則體尺度可表示為γ=pl,p由Es層電子密度自相關(guān)函數(shù)和修正貝塞爾函數(shù)的階數(shù)決定,通常電子密度自相關(guān)函數(shù)[21]取值為0.5,按照中低緯地區(qū)的經(jīng)驗取值,修正貝塞爾函數(shù)的階數(shù)m取5[12],由(5)式可確定出p=3.5,γ=3.5l.則電離層Es不規(guī)則體散射截面進一步表示為

    根據(jù)雷達(dá)方程,則經(jīng)電離層Es散射的VHF信號接收功率為

    3 電離層Es電波傳播模型

    為了深入研究中緯度地區(qū)電離層Es層VHF傳播特性,選擇在夏季Es高發(fā)的時段開展了電離層Es層VHF信號測量試驗,實驗共進行了20余天,有效實驗數(shù)據(jù)約100 h.為覆蓋不同的區(qū)域和距離,試驗采用南北和東西兩種鏈路布局.其中南北布局為昆明-西安的鏈路1,距離為1065 km;東西布局有兩條鏈路,其中德宏-懷化的鏈路2,距離為1240 km,德宏-郴州的鏈路3,距離為1590 km.試驗使用了多個測量頻率,并在每條鏈路中重點使用了一個頻率,分別為50.5,45.5,31 MHz.測量鏈路布局如圖2所示.

    實驗期間典型VHF信號接收功率以及對應(yīng)的頻高圖如圖3所示.由圖3可以看出:在無Es期間VHF接收信號大多淹沒在噪聲中,且信號起伏較大,不穩(wěn)定;隨著電離層Es由弱變強,接收信號逐漸增強.因此判定正常E區(qū)以及F區(qū)不規(guī)則結(jié)構(gòu)相對突發(fā)Es層的散射能力很弱,本文有效數(shù)據(jù)要求:1)明顯出現(xiàn)Es且持續(xù)至少1 min時段;2)信號至少高于背景噪聲10 dB.同時參照信道時延功率譜將較強的F層不均勻體散射的數(shù)據(jù)剔除.經(jīng)上述數(shù)據(jù)篩選后,三條測量鏈路的VHF信號平均接收功率隨著Es層臨界頻率與入射頻率的比值foEs/f的變化如圖4所示.

    由圖4可以看出:在三條測量鏈路中VHF信號平均接收功率隨著foEs/f的增長而上升.但VHF信號接收功率低端增長緩慢,中間段快速上升,高端增長又變緩.利用第2節(jié)給出的反射和散射理論對鏈路1進行仿真,根據(jù)電離層Es的厚度通常在幾百米到1 km,其不規(guī)則體水平方向尺度在幾十米到幾百米,垂直方向尺度幾十米,選擇反射作用機制中的n=5,電離層Es半厚度?h取50和500 m,散射作用機制中垂直尺度L3取30 m,水平方向尺度L1取50和300 m,仿真結(jié)果如圖5所示.

    圖2 實驗測量鏈路布局Fig.2.Link-arrangement of test and measurement.

    圖3 (網(wǎng)刊彩色)實驗期間典型測量結(jié)果(a)VHF接收信號功率;(b)鏈路中點垂測儀觀測結(jié)果Fig.3.(color online)Typical measurement results during the test:(a)Received power of VHF signal;(b)results of ionosonode in midpoint of measurement-link.

    圖5 VHF信號接收功率仿真Fig.5.Simulation of the received power of VHF signal.

    由圖5可以看出,在foEs/f高端,VHF接收信號平均功率與反射作用機制相符,在低端與散射機制相符.而在中間過渡區(qū),VHF信號接收功率快速上升,其斜率和foEs/f較高和較低區(qū)域明顯不同,即在該區(qū)域已經(jīng)不能用單一的反射作用或散射作用機制解釋,據(jù)此可初步判斷在該區(qū)域可能是反射與散射的共同作用機制,并且隨著foEs/f的上升,反射作用機制逐漸增強.

    下面通過實測數(shù)據(jù)并結(jié)合電離層Es對無線電波的反射和散射作用機制,研究并建立電離層Es電波傳播模型.首先定義實測數(shù)據(jù)與仿真結(jié)果(基于反射作用機制或散射作用機制)之間的平均誤差ˉδ為

    其中,N是在考察區(qū)域內(nèi)樣本點的個數(shù);Pr,P′r是實測接收功率和仿真接收功率.為均衡不同鏈路中的誤差,定義本次測量試驗中三條鏈路實測數(shù)據(jù)與仿真結(jié)果之間的總平均誤差ˉδ總為

    3.1 反射作用分析

    根據(jù)(3)式,對于入射到電離層Es單一頻率電波,在收發(fā)兩點距離、Es層高度確定后,影響反射信號強度的主要是Es層的半厚度?h和電離層Es的電子密度分布指數(shù)n.當(dāng)n相同時,?h越大,電波衰減越大,接收功率越小.事實上,相同電離層Es反射系數(shù)具有多組n和?h的組合,本文按照中低緯地區(qū)的經(jīng)驗[22]取值n=5,通過實驗數(shù)據(jù)確定等效半厚度?h.根據(jù)前面論述,經(jīng)Es傳播的VHF接收信號在foEs/f高端主要是Es反射作用,設(shè)這一區(qū)間為foEs/f>fr,fr為完全可用反射作用機制解釋的門限值.圖6給出了不同等效半厚度下,三條鏈路經(jīng)電離層Es反射的實測結(jié)果與基于反射作用機制仿真結(jié)果的平均誤差δˉ隨fr的變化,表1給出了三條鏈路中均滿足實測數(shù)據(jù)與仿真結(jié)果δˉ<2%時的不同鏈路門限值fr和總平均誤差δˉ總.

    通常δˉ<2%在實際測量中已經(jīng)是比較苛刻的要求,由圖6可以看出滿足此要求的條件為:鏈路1中90 m

    圖6 (網(wǎng)刊彩色)不同等效半厚度下平均誤差ˉδ隨fr的變化(a)鏈路1;(b)鏈路2;(c)鏈路3Fig.6.(color online)Average errorˉδ under di ff erent half-thickness changed with fr:(a)Link 1;(b)Link 2;(c)Link 3.

    表1 實測數(shù)據(jù)與基于反射作用仿真結(jié)果誤差ˉδ均小于2%的門限值與平均誤差Table 1.The threshold and total average error under the average errorˉδ between measured power and simulated power based on re fl ecting process less than 2%.

    如果在整個電離層Es持續(xù)時間內(nèi)?h變化不大,上述結(jié)果表明:某一頻率的VHF電波入射Es層,當(dāng)foEs/f達(dá)到某個門限值時,接收端接收到的信號主要來自于Es層對VHF電波的反射.當(dāng)foEs/f大于此門限值時,可用電離層Es反射作用機制進行解釋,此時Es反射系數(shù)等效表示為

    當(dāng)foEs/f低于該值時,需尋求其他機制去解釋.

    3.2 散射作用分析

    根據(jù)(10)式,電離層Es中不規(guī)則體水平尺度相同時,垂直方向上尺度越大,衰減越大;而當(dāng)垂直方向上尺度相同時,水平尺度越大,衰減越小.根據(jù)前面的論述,經(jīng)Es傳播的VHF接收信號在foEs/f低端主要是Es散射作用,設(shè)這一區(qū)間為foEs/f

    根據(jù)3.1節(jié)中δˉ<2%的誤差要求,由圖7可以看出滿足此要求的條件為:鏈路1中160 m

    表2 實測數(shù)據(jù)與基于散射作用仿真結(jié)果誤差ˉδ均小于2%的門限值與總平均誤差Table 2.The threshold and total average error under the average errorˉδ between measured power and simulated power based on scattering process less than 2%.

    圖7 (網(wǎng)刊彩色)不同水平尺度下平均誤差ˉδ隨fs的變化(a)鏈路1;(b)鏈路2;(c)鏈路3Fig.7.(color online)Average errorˉδ under di ff erent scale changed with fs:(a)Link 1;(b)Link 2;(c)Link 3.

    如果在整個電離層Es持續(xù)時間內(nèi)L1,L2,L3變化不大,上述結(jié)果表明:當(dāng)某一頻率的VHF電波入射Es層,當(dāng)foEs/f低于某個門限值時,接收端接收到的信號主要來自于Es層對VHF電波的散射,此時Es散射截面積可以等效表示為

    當(dāng)foEs/f大于該門限值時,需尋求其他機制去解釋.

    3.3 電離層Es電波傳播模型

    綜合分析3.1節(jié)和3.2節(jié)仿真結(jié)果,發(fā)現(xiàn)在鏈路1中當(dāng)foEs/f在0.11—0.14區(qū)間時,鏈路2中foEs/f在0.14—0.21區(qū)間時,鏈路3中foEs/f在0.13—0.25區(qū)間時,利用反射作用機制和散射作用機制都不能有效解釋實測結(jié)果.根據(jù)前面的分析,電離層Es對入射的電波既有反射作用機制,又有散射作用機制,兩種作用隨電離層的背景條件而發(fā)生變化.設(shè)反射作用機制與散射作用機制比例常數(shù)為kr,ks,且kr+ks=1,則此時接收端VHF信號強度Pr可表示為

    其中,Pt是VHF信號發(fā)射功率;Hr,Hs為對應(yīng)Es層對VHF信號反射作用和散射作用機制的傳輸參數(shù).

    根據(jù)上述分析,此區(qū)間即為fs

    圖8 (網(wǎng)刊彩色)過渡區(qū)域內(nèi)反射作用所占比重的擬合結(jié)果Fig.8.(color online)The fi tting result of coefficient re fl ection process in transition section.

    由圖8可以看出,在此區(qū)間內(nèi)反射作用所占比重隨著f′的上升而逐漸增大,并且反射作用所占比重更符合指數(shù)分布,

    確立過渡區(qū)域反射作用和散射作用所占比重后,整個Es存在期間內(nèi)電離層Es對電波的作用隨foEs/f變化的作用機制即確定.根據(jù)上面的分析,將(11),(12)和(14)式代入(13)式,建立由散射作用、散射與反射共同作用、反射作用三段式電離層Es電波傳播模型,表示為

    其中:Pt為發(fā)射功率;Gt和Gr分別為發(fā)射和接收天線在指向Es高度方向的增益;D為傳播路徑;λ是入射波波長;θ為掠射角,L=460πsinθ/λ.kr是反射信號所占比重,符合如下分布:

    其中:fr是完全反射作用機制解釋的門限值,fs是完全散射作用機制解釋的門限值.在本次測量中,鏈路1中fr=0.14,fs=0.11;鏈路2中fr=0.21,fs=0.14;鏈路3中fr=0.25,fs=0.13.根據(jù)該模型的(15)式和(16)式對三條鏈路進行模擬,結(jié)果如圖9所示.

    由圖9可以看出,本文給出的三段式動態(tài)模型很好地解釋了實驗期間接收端VHF信號幅度的變化:foEs/f>fr區(qū)間,VHF信號傳播符合電離層Es反射作用機制;foEs/f

    圖9 (網(wǎng)刊彩色)基于本文Es電波傳播模型的仿真結(jié)果和實測數(shù)據(jù)Fig.9.(color online)The result of simulated on the proposed model of radio wave propagation under Es and the result of measurement.

    3.4 模型對比

    根據(jù)本文所建立的模型,對三條鏈路中電離層Es對入射電波的衰減進行仿真,并與ITU模型比較.設(shè)Es高度110 km,鏈路1中f=50.5 MHz,fr=0.14,fs=0.11;鏈路2中f=45.5 MHz,fr=0.21,fs=0.14;鏈路3中f=31.5 MHz,fr=0.25,fs=0.13.仿真結(jié)果如圖10所示.

    利用本文模型計算的VHF信號衰減值與ITU模型結(jié)果一致,特別是在foEs/f高端時,一致程度更好;但在foEs/f取較小的值時,利用本文模型計算的VHF信號衰減值要小很多.分析其原因,本文所建模型是針對中低緯地區(qū)電離層Es對電波作用的實測數(shù)據(jù),電離層Es特性電離層背景關(guān)系比較密切,不同緯度地區(qū)電離層Es特性有差別,而ITU模型是基于中高緯地區(qū)的實測數(shù)據(jù)而建立,其應(yīng)用在中低緯緯度誤差會增大;其次,ITU模型預(yù)測時采用實驗數(shù)據(jù)統(tǒng)計后的經(jīng)驗值,與實測結(jié)果也存在差異,因此這種差異會導(dǎo)致誤差;最后,ITU預(yù)測模型誤差會隨著foEs/f減小而增大,在特別是對于foEs/f取較小值時,電離層Es對入射電波主要是散射作用機制,而散射作用遠(yuǎn)小于反射作用,進一步增大了誤差.因此,我們認(rèn)為本文所建更適用于中低緯地區(qū).

    圖10 VHF信號接收功率隨foEs/f變化(a)鏈路1;(b)鏈路2;(c)鏈路3Fig.10.Received power of VHF signal changed with foEs/f:(a)Link 1;(b)Link 2;(c)Link 3.

    4 結(jié)論

    本文基于實測數(shù)據(jù),結(jié)合電離層Es電子密度分布特點,建立包含散射、反射與散射共同作用以及反射的三段式動態(tài)電離層Es電波傳播模型,該模型中電離層Es對入射電波的作用隨foEs/f變化如下:

    1)當(dāng)foEs/f高于某一門限(完全可用反射作用機制解釋的門限值fr)時,VHF信號經(jīng)Es層的傳播符合反射作用機制;

    2)當(dāng)foEs/f低于某一門限(完全可用散射作用機制解釋的門限值fs)時,VHF信號經(jīng)Es層的傳播符合散射作用機制;

    3)當(dāng)foEs/f在fs與fr之間的過渡區(qū)域時,VHF信號經(jīng)Es層的傳播是反射和散射共同作用機制,且此區(qū)間隨著foEs/f增大,反射作用機制逐漸增強.

    對三條鏈路實驗數(shù)據(jù)統(tǒng)計可發(fā)現(xiàn),不同鏈路間fr和fs的值存在差異,尤其是鏈路1中fr,fs與鏈路2鏈路3中的fr,fs差異較大,其可能原因是鏈路1與鏈路2、鏈路3地理位置中心點距離較遠(yuǎn),電離層背景參數(shù)差異較大,如E層或F2層的散射強度等;同時由于采用的探測波頻率也不相同,也可能導(dǎo)致各鏈路模型參數(shù)存在差異.

    由于本文在模型建立過程中基于的實驗數(shù)據(jù)相對較少,模型參數(shù)確認(rèn)與優(yōu)化還需要后期通過大量實驗數(shù)據(jù)完成.但本文所建模型是依據(jù)中低緯度地區(qū)測量數(shù)據(jù),在后期模型優(yōu)化后可以用來預(yù)測VHF頻段信號經(jīng)Es傳播能量分布,具有較強的工程應(yīng)用價值.

    [1]Smith L G,Mechtly E A 1972 Radio Sci.7 367

    [2]Beynon W J G,Maude A D 1972 Planet.Space Sci.20 809

    [3]Kobayashi T 1964 Radio Res.11 181

    [4]Bramley E N 1972 J.Atmos.Terr.Phys.34 1495

    [5]Althouse E L,Davis J R 1972 Radio Sci.7 897

    [6]Miya K,Shimizu K 1978 Radio Sci.13 559

    [7]Kerblay T S,Makarenko S F 1980 Geomagn.Aeron.20 449

    [8]Sherstyukov O N,Akchurin A D 2009 Adv.Space Res.43 1835

    [9]Sherstyukov O N,Akchurin A D 2011 General Assembly and Scienti fi c SymposiumIstanbul,Aug.13–20,2011 p1

    [10]Sherstyukov O N,Akchurin A D,Sherstyukov R O 2015 Adv.Space Res.56 1169

    [11]Whitehead J D 1989 J.Atmosph.Solar-Terr.Phys.51 401

    [12]Tao K 1962 Ionospheric Sporadic E(New York:MacMillian Company)p235

    [13]Pan W Y 1981 Acta Phys.Sin.30 661(in Chinese)[潘威炎1981物理學(xué)報30 661]

    [14]Zeng Z,Sokolovskiy S 2010 Geophys.Res.Lett.37 1480

    [15]Huang C S,Li J 1994 Acta Phys.Sin.43 1476(in Chinese)[黃朝松,李鈞1994物理學(xué)報43 1476]

    [16]Zhang Q,Wu X J 2016 Acta Phys.Sin.65 038102(in Chinese)[張卿,武新軍2016物理學(xué)報65 038102]

    [17]Deng F,Zhao Z Y,Shi R,Zhang Y N 2009 Acta Phys.Sin.58 7382(in Chinese)[鄧峰,趙正予,石潤,張援農(nóng)2009物理學(xué)報58 7382]

    [18]Booker H G,Gordon W E 1950 Proceed.Ire.38 401

    [19]Norton K A 1956 Commun.Syst.Ire Trans.4 39

    [20]Malaga A 1986 Rome Air Development Center Air Force Systems Command Griffiss Air Force Base New York,Aug 1985 p13441

    [21]Booker H G,Gordon W E 1958 Proceed.Ire.46 298

    [22]Ovezgeldyev O G 1988 Geomagn.Aeron.28 1024

    (Received 4 January 2017;revised manuscript received 8 March 2017)

    PACS:94.20.Vv,94.30.TzDOI:10.7498/aps.66.119401

    F信號接收功率隨foEs/f的變化 Fig.4.

    power of VHF signal changed with foEs/f.

    Modeling of radio wave propagations under sporadic-E in fl uence at low and middle latitudes?

    Hao Shu-JiZhang Wen-Chao?Zhang Ya-BinYang Ju-TaoMa Guang-Lin
    (National Key Laboratory of Electromagnetic Environment,China Research Institute of Radiowave Propagation,Qingdao 266107,China)

    The sporadic-E(Es)layer is a thin layer of several kilometers existing at an altitude around 100 km and features extremely dense ionized irregularities,which can re fl ect or scatter high frequency(HF)and very high frequency(VHF)radio waves.The most popular theoretical explanation for mid-latitude Es formation is the wind shear theory.Measurements by rocket souding have shown that Es has high electron density and relatively sharp density gradient in the vertical direction.The one-hop propagation of VHF signal in Es can even reach as far as 2000 km.In this paper,we consider incident radio waves in fl uenced by Es via both re fl ecting and scattering processes at low and middle latitudes,the coefficients of which are related to and vary with the critical frequency of Es(foEs).Firstly,with a supposed parabolic density distribution and the autocorrelation function of the electron density given by Booker,HF and VHF radio wave propagations in Es are analyzed according to the re fl ection and scattering theory.Secondly,a numerical model for the combined re fl ecting and scattering processes is developed in the form of piecewise function,the contribution of which can be distinguished by the portion factor of re fl ection(kr).According to the model,there are two threshold ratios of the critical frequency to the wave frequencies frand fsrespectively.The incident radio waves are totally re fl ected by Es when foEs/f is higher than frand mostly scattered when foEs/f is lower than fs.A transition zone exists between two critical points,with the combined processes working together.Thirdly,HF/VHF radio wave propagations in low and middle latitudes of Es are are in the north-southern direction and east-western direction separately.The experiment link in the north-southern direction is from Kunming to Xi’an at distance of 1065 km,and the ionosonde used for Es observation is located at Chongqing.Two east-west links are arranged,one of which is from Dehong to Huaihua and the other is from Dehong to Chenzhou,with the ionosonde located at Kunming and the ground distance as far as 1240 km and 1590 km respectively.The measurement data are treated and parameters of the above mentioned model for wave propagation in Es are experimentally determined.Finally,our model is veri fi ed by comparing with ITU-R model.Our results are consistent with the results from the ITU-R model when the foEs/f is high(i.e.,the re fl ecting process plays a main role).When the scattering process dominates,the attenuation value of VHF signal is far less than that predicted by the ITU-R model,which is closer to actual measurements.It is concluded that our model is more preferable for HF and VHF radio wave propagations in Es at low and middle latitudes.

    sporadic-E,radio wave propagation,International Telecommunications Union model

    10.7498/aps.66.119401

    ?中國電科技術(shù)創(chuàng)新基金(批準(zhǔn)號:A171601C01)資助的課題.

    ?通信作者.E-mail:zhangwenchao0726@126.com

    ?2017中國物理學(xué)會Chinese Physical Society

    http://wulixb.iphy.ac.cn

    *Project supported by the Technology Innovate Fund of China Electronics Technology Group Corporation(Grant No.A171601C01).

    ?Corresponding author.E-mail:zhangwenchao0726@126.com

    猜你喜歡
    電子密度電波電離層
    永遠(yuǎn)的紅色電波
    中國記者(2021年11期)2021-12-11 09:05:04
    一種電離層TEC格點預(yù)測模型
    The Speed of Light
    Kalman濾波估算電離層延遲的一種優(yōu)化方法
    顧及地磁影響的GNSS電離層層析不等像素間距算法*
    不同GPS掩星電離層剖面產(chǎn)品相關(guān)性分析
    測繪通報(2019年11期)2019-12-03 01:47:34
    瞌睡電波
    幽默大師(2019年11期)2019-11-23 08:47:58
    等離子體電子密度分布信息提取方法研究
    一種適用于電離層電子密度重構(gòu)的AMART算法
    “電波衛(wèi)士”在行動
    99久久99久久久精品蜜桃| 一夜夜www| 亚洲美女黄片视频| 乱人伦中国视频| 日韩大尺度精品在线看网址 | 国产一区二区三区综合在线观看| 一级a爱片免费观看的视频| 亚洲精品中文字幕在线视频| 久久精品国产综合久久久| 多毛熟女@视频| 成人亚洲精品av一区二区 | 美女国产高潮福利片在线看| 亚洲午夜理论影院| 国产精品一区二区免费欧美| 欧美黄色淫秽网站| 精品一区二区三区四区五区乱码| a在线观看视频网站| 亚洲欧美一区二区三区黑人| 日韩欧美三级三区| 一级a爱片免费观看的视频| 国产深夜福利视频在线观看| 亚洲人成网站在线播放欧美日韩| 亚洲av片天天在线观看| 色尼玛亚洲综合影院| 欧美精品一区二区免费开放| 99在线视频只有这里精品首页| 午夜福利在线观看吧| 黄频高清免费视频| 交换朋友夫妻互换小说| 中文字幕人妻丝袜制服| 日韩三级视频一区二区三区| 岛国视频午夜一区免费看| 宅男免费午夜| 波多野结衣高清无吗| 黄色成人免费大全| 一区二区三区国产精品乱码| 久久国产精品人妻蜜桃| 国产成人欧美在线观看| 精品无人区乱码1区二区| 操出白浆在线播放| 热99re8久久精品国产| 一边摸一边抽搐一进一小说| 美女高潮到喷水免费观看| 久久亚洲精品不卡| av中文乱码字幕在线| 欧美一区二区精品小视频在线| 色老头精品视频在线观看| 久久精品亚洲精品国产色婷小说| 99精品在免费线老司机午夜| 人妻丰满熟妇av一区二区三区| 法律面前人人平等表现在哪些方面| 老司机福利观看| 精品欧美一区二区三区在线| 免费搜索国产男女视频| 午夜亚洲福利在线播放| 亚洲av熟女| 久久人妻熟女aⅴ| 日本精品一区二区三区蜜桃| 五月开心婷婷网| 国产成人精品久久二区二区91| 成人手机av| 国产欧美日韩一区二区精品| 丝袜人妻中文字幕| 久久久久国内视频| 大型av网站在线播放| 久久久精品欧美日韩精品| 叶爱在线成人免费视频播放| 中文字幕人妻丝袜制服| 香蕉丝袜av| 亚洲av日韩精品久久久久久密| 国产午夜精品久久久久久| 香蕉丝袜av| 一个人免费在线观看的高清视频| 亚洲av美国av| av视频免费观看在线观看| 日本欧美视频一区| 国产真人三级小视频在线观看| 国产成人系列免费观看| 久久精品国产清高在天天线| 日本黄色日本黄色录像| 久久久水蜜桃国产精品网| 免费高清视频大片| 日本精品一区二区三区蜜桃| 久久久国产精品麻豆| 国产精品自产拍在线观看55亚洲| 国产成人精品在线电影| 国产蜜桃级精品一区二区三区| 亚洲少妇的诱惑av| 三上悠亚av全集在线观看| 午夜激情av网站| 一本大道久久a久久精品| 91精品国产国语对白视频| 国产高清国产精品国产三级| 成年女人毛片免费观看观看9| 在线观看日韩欧美| 桃红色精品国产亚洲av| 夜夜爽天天搞| 久久中文字幕人妻熟女| 国产精品一区二区三区四区久久 | 精品久久久久久久毛片微露脸| 夜夜爽天天搞| 日本a在线网址| 女生性感内裤真人,穿戴方法视频| 黄色女人牲交| 国产一区在线观看成人免费| 午夜视频精品福利| 国产精品亚洲一级av第二区| 国产伦一二天堂av在线观看| 亚洲欧美日韩高清在线视频| 夫妻午夜视频| 校园春色视频在线观看| 亚洲成a人片在线一区二区| 久99久视频精品免费| 黄片大片在线免费观看| 1024视频免费在线观看| 国产精品自产拍在线观看55亚洲| 五月开心婷婷网| 天堂√8在线中文| 欧美老熟妇乱子伦牲交| 精品一品国产午夜福利视频| 99久久人妻综合| 久久久久亚洲av毛片大全| 午夜福利欧美成人| 99在线人妻在线中文字幕| 宅男免费午夜| 久久国产亚洲av麻豆专区| 久热爱精品视频在线9| 免费看a级黄色片| 黄色视频不卡| 一区在线观看完整版| 淫妇啪啪啪对白视频| 国产亚洲av高清不卡| 成人三级黄色视频| 亚洲精品一二三| 精品国产乱码久久久久久男人| 一本大道久久a久久精品| 50天的宝宝边吃奶边哭怎么回事| 国产片内射在线| 法律面前人人平等表现在哪些方面| 欧美一级毛片孕妇| 夜夜夜夜夜久久久久| 99热国产这里只有精品6| 久久精品亚洲av国产电影网| a级毛片在线看网站| 亚洲七黄色美女视频| 亚洲欧美一区二区三区黑人| 国产精品av久久久久免费| 国产一区二区激情短视频| 19禁男女啪啪无遮挡网站| 天堂动漫精品| 99在线人妻在线中文字幕| 久久香蕉国产精品| 99在线人妻在线中文字幕| 久久精品影院6| 免费高清视频大片| 亚洲av日韩精品久久久久久密| 女生性感内裤真人,穿戴方法视频| 在线天堂中文资源库| 19禁男女啪啪无遮挡网站| 亚洲av成人av| 欧美日韩中文字幕国产精品一区二区三区 | 两个人免费观看高清视频| 久热这里只有精品99| 午夜福利欧美成人| 日本 av在线| 80岁老熟妇乱子伦牲交| 丰满的人妻完整版| 国产蜜桃级精品一区二区三区| 日韩欧美三级三区| 美女高潮喷水抽搐中文字幕| 亚洲,欧美精品.| 一二三四社区在线视频社区8| 亚洲精品国产精品久久久不卡| 女性生殖器流出的白浆| 国产精品自产拍在线观看55亚洲| 国产有黄有色有爽视频| 欧美成人性av电影在线观看| 国产成人精品无人区| 男女下面插进去视频免费观看| 国产精品久久久久久人妻精品电影| cao死你这个sao货| 一级,二级,三级黄色视频| 久久天躁狠狠躁夜夜2o2o| 亚洲国产欧美一区二区综合| 午夜福利在线观看吧| 这个男人来自地球电影免费观看| 这个男人来自地球电影免费观看| 真人一进一出gif抽搐免费| 欧美一级毛片孕妇| 精品国产乱子伦一区二区三区| 久久 成人 亚洲| 曰老女人黄片| 国产深夜福利视频在线观看| av视频免费观看在线观看| 露出奶头的视频| av超薄肉色丝袜交足视频| 在线观看免费午夜福利视频| √禁漫天堂资源中文www| 好看av亚洲va欧美ⅴa在| 悠悠久久av| 在线免费观看的www视频| 99国产综合亚洲精品| 亚洲一区二区三区欧美精品| 99热国产这里只有精品6| 黑人猛操日本美女一级片| 国产成人精品久久二区二区免费| 亚洲精品一卡2卡三卡4卡5卡| 好男人电影高清在线观看| 在线视频色国产色| 午夜免费观看网址| 免费观看人在逋| 如日韩欧美国产精品一区二区三区| 日本vs欧美在线观看视频| 天堂影院成人在线观看| 99国产精品一区二区三区| 黄色 视频免费看| 少妇被粗大的猛进出69影院| 国产精品爽爽va在线观看网站 | 精品久久蜜臀av无| 9色porny在线观看| 亚洲av五月六月丁香网| 999久久久国产精品视频| 丰满的人妻完整版| 国产av精品麻豆| 欧美乱妇无乱码| av网站在线播放免费| 91老司机精品| 免费在线观看视频国产中文字幕亚洲| 国产亚洲av高清不卡| 精品久久久久久,| 91精品三级在线观看| 免费观看人在逋| 制服诱惑二区| 亚洲第一欧美日韩一区二区三区| 国产99久久九九免费精品| 又黄又粗又硬又大视频| 欧美亚洲日本最大视频资源| 亚洲片人在线观看| 免费观看人在逋| 日韩有码中文字幕| 亚洲情色 制服丝袜| 国产精品免费视频内射| av天堂在线播放| 欧美日本中文国产一区发布| 国产成人精品久久二区二区91| 国产国语露脸激情在线看| 久久久久精品国产欧美久久久| 丝袜人妻中文字幕| av福利片在线| 黄片播放在线免费| 欧美激情久久久久久爽电影 | 欧美日韩瑟瑟在线播放| 黄色成人免费大全| 嫩草影视91久久| 日韩 欧美 亚洲 中文字幕| 日本精品一区二区三区蜜桃| 69av精品久久久久久| 精品熟女少妇八av免费久了| 正在播放国产对白刺激| 丝袜美足系列| 免费在线观看视频国产中文字幕亚洲| 国产欧美日韩一区二区三区在线| 涩涩av久久男人的天堂| 欧美日本中文国产一区发布| 男女下面插进去视频免费观看| 久99久视频精品免费| 丰满人妻熟妇乱又伦精品不卡| 成人永久免费在线观看视频| 久久久久久人人人人人| 国产免费av片在线观看野外av| 久久久国产精品麻豆| 午夜免费激情av| 黄片播放在线免费| 亚洲欧美激情在线| 男女午夜视频在线观看| 日本黄色日本黄色录像| 十八禁人妻一区二区| 中文字幕色久视频| 在线国产一区二区在线| 一边摸一边做爽爽视频免费| 一级a爱片免费观看的视频| 亚洲一码二码三码区别大吗| 夫妻午夜视频| 一边摸一边做爽爽视频免费| 国产免费av片在线观看野外av| 久久99一区二区三区| 97超级碰碰碰精品色视频在线观看| xxxhd国产人妻xxx| 12—13女人毛片做爰片一| 亚洲精品美女久久av网站| 90打野战视频偷拍视频| 国产色视频综合| 91麻豆av在线| 在线国产一区二区在线| 精品国产超薄肉色丝袜足j| 日韩 欧美 亚洲 中文字幕| 国产三级黄色录像| 日韩人妻精品一区2区三区| 欧美成人性av电影在线观看| 久久午夜亚洲精品久久| 国产成人影院久久av| aaaaa片日本免费| 国产成人免费无遮挡视频| 久久欧美精品欧美久久欧美| 99精品欧美一区二区三区四区| 国产精华一区二区三区| 国产av一区在线观看免费| 变态另类成人亚洲欧美熟女 | 女同久久另类99精品国产91| 无遮挡黄片免费观看| 91精品三级在线观看| 黄频高清免费视频| 老司机午夜福利在线观看视频| 国产成+人综合+亚洲专区| 欧美在线黄色| 久久人妻熟女aⅴ| 精品一区二区三区四区五区乱码| 精品高清国产在线一区| av中文乱码字幕在线| 19禁男女啪啪无遮挡网站| 中文欧美无线码| 久久影院123| 免费日韩欧美在线观看| 女人高潮潮喷娇喘18禁视频| 成年人免费黄色播放视频| 一边摸一边抽搐一进一小说| 色尼玛亚洲综合影院| 国产精品乱码一区二三区的特点 | 少妇 在线观看| 久久国产精品男人的天堂亚洲| 亚洲av日韩精品久久久久久密| 青草久久国产| 热re99久久国产66热| 91成人精品电影| 俄罗斯特黄特色一大片| 亚洲第一欧美日韩一区二区三区| 国产乱人伦免费视频| 欧美老熟妇乱子伦牲交| 国产av又大| 国内久久婷婷六月综合欲色啪| 水蜜桃什么品种好| 欧美午夜高清在线| 黄片大片在线免费观看| 久久人妻av系列| 国产精品1区2区在线观看.| 99re在线观看精品视频| 国产精品av久久久久免费| 97超级碰碰碰精品色视频在线观看| 国产熟女xx| 国产又爽黄色视频| 在线看a的网站| 伦理电影免费视频| 老司机在亚洲福利影院| 亚洲第一av免费看| 久久中文字幕一级| 国产97色在线日韩免费| 韩国av一区二区三区四区| 日本五十路高清| 成人免费观看视频高清| 日本一区二区免费在线视频| 99精品在免费线老司机午夜| 亚洲一区二区三区色噜噜 | 久久久精品欧美日韩精品| 无人区码免费观看不卡| 久久国产亚洲av麻豆专区| 老司机靠b影院| 久久精品人人爽人人爽视色| 一区二区日韩欧美中文字幕| 国产精品爽爽va在线观看网站 | 丰满迷人的少妇在线观看| 黑人巨大精品欧美一区二区mp4| 亚洲,欧美精品.| 制服诱惑二区| av国产精品久久久久影院| 午夜视频精品福利| 国产精品九九99| 99久久国产精品久久久| 国产亚洲精品综合一区在线观看 | 午夜视频精品福利| 女生性感内裤真人,穿戴方法视频| 麻豆国产av国片精品| 国产激情欧美一区二区| 身体一侧抽搐| 中文字幕人妻丝袜一区二区| 波多野结衣一区麻豆| xxxhd国产人妻xxx| 成人国产一区最新在线观看| 国产一区二区三区综合在线观看| 国产精品美女特级片免费视频播放器 | 身体一侧抽搐| 十八禁网站免费在线| 女警被强在线播放| 亚洲av片天天在线观看| 一级a爱片免费观看的视频| 高清毛片免费观看视频网站 | 国产精品亚洲一级av第二区| 天堂√8在线中文| x7x7x7水蜜桃| 日韩国内少妇激情av| 久久精品aⅴ一区二区三区四区| 国产精品成人在线| 在线观看舔阴道视频| av视频免费观看在线观看| 国产欧美日韩一区二区精品| 国产熟女xx| 精品国产美女av久久久久小说| 两性夫妻黄色片| 亚洲专区字幕在线| 国产一卡二卡三卡精品| 欧美日韩福利视频一区二区| 成年版毛片免费区| 欧美午夜高清在线| 在线观看免费午夜福利视频| 免费日韩欧美在线观看| 99热只有精品国产| 黄色片一级片一级黄色片| 国产熟女xx| 免费av毛片视频| 免费在线观看亚洲国产| 国产精品亚洲一级av第二区| 精品久久久久久久毛片微露脸| 欧美一级毛片孕妇| 欧美日韩亚洲高清精品| 亚洲黑人精品在线| 一级毛片精品| 一进一出好大好爽视频| 黄色片一级片一级黄色片| 侵犯人妻中文字幕一二三四区| videosex国产| xxxhd国产人妻xxx| 国产欧美日韩一区二区三| 亚洲欧美激情综合另类| 国产精品98久久久久久宅男小说| 91麻豆精品激情在线观看国产 | 欧美精品亚洲一区二区| 人人妻人人爽人人添夜夜欢视频| 亚洲一卡2卡3卡4卡5卡精品中文| a级片在线免费高清观看视频| 亚洲中文av在线| 狂野欧美激情性xxxx| 少妇 在线观看| 黄色视频不卡| 国产单亲对白刺激| 亚洲中文字幕日韩| 90打野战视频偷拍视频| 久久久国产一区二区| 久久精品91蜜桃| 欧美日韩亚洲综合一区二区三区_| 熟女少妇亚洲综合色aaa.| 一区二区三区激情视频| 美女高潮到喷水免费观看| 自拍欧美九色日韩亚洲蝌蚪91| 夜夜看夜夜爽夜夜摸 | 一夜夜www| 久久人妻av系列| 男女午夜视频在线观看| 一边摸一边抽搐一进一小说| 久久 成人 亚洲| 男女下面进入的视频免费午夜 | 国产高清国产精品国产三级| 亚洲专区中文字幕在线| 久久久久久久精品吃奶| 免费观看精品视频网站| 亚洲国产毛片av蜜桃av| 大香蕉久久成人网| 国产亚洲精品一区二区www| 国产精品国产av在线观看| 女警被强在线播放| 9热在线视频观看99| 亚洲精品粉嫩美女一区| 人妻久久中文字幕网| 老汉色∧v一级毛片| 中文字幕人妻丝袜制服| 久久影院123| 涩涩av久久男人的天堂| www.www免费av| 在线观看舔阴道视频| 精品乱码久久久久久99久播| 亚洲人成伊人成综合网2020| 久久精品亚洲熟妇少妇任你| 国产精品一区二区三区四区久久 | 日日爽夜夜爽网站| 国产一区在线观看成人免费| 69精品国产乱码久久久| 久久久久国产精品人妻aⅴ院| 亚洲片人在线观看| 在线观看免费午夜福利视频| 欧美不卡视频在线免费观看 | 午夜福利,免费看| 成人国语在线视频| 亚洲精品中文字幕在线视频| 满18在线观看网站| 黄色怎么调成土黄色| 99re在线观看精品视频| 18禁裸乳无遮挡免费网站照片 | 欧美成人免费av一区二区三区| 老汉色∧v一级毛片| 69av精品久久久久久| 99国产精品99久久久久| 国产91精品成人一区二区三区| 女人精品久久久久毛片| 亚洲精品中文字幕在线视频| 亚洲一码二码三码区别大吗| 日日摸夜夜添夜夜添小说| 黄色女人牲交| 成人国语在线视频| ponron亚洲| 91国产中文字幕| 757午夜福利合集在线观看| 成年人免费黄色播放视频| 高清在线国产一区| 丁香欧美五月| 精品久久久久久久毛片微露脸| 国产片内射在线| 天天影视国产精品| 国产欧美日韩综合在线一区二区| 亚洲欧美日韩另类电影网站| 精品国产乱码久久久久久男人| 亚洲欧美日韩高清在线视频| 国产欧美日韩一区二区三| 欧美日韩亚洲国产一区二区在线观看| ponron亚洲| 高清毛片免费观看视频网站 | 少妇粗大呻吟视频| 中国美女看黄片| 一边摸一边抽搐一进一出视频| 91成年电影在线观看| 久久久久久久久免费视频了| 亚洲美女黄片视频| 久久精品亚洲熟妇少妇任你| 亚洲av日韩精品久久久久久密| 日韩人妻精品一区2区三区| 亚洲午夜理论影院| 国产精品av久久久久免费| 9191精品国产免费久久| 久久人妻熟女aⅴ| 日韩精品中文字幕看吧| 青草久久国产| 精品人妻1区二区| 怎么达到女性高潮| 欧美午夜高清在线| 看片在线看免费视频| 日韩欧美免费精品| 一边摸一边做爽爽视频免费| 美女 人体艺术 gogo| 天天躁狠狠躁夜夜躁狠狠躁| 国产aⅴ精品一区二区三区波| 99国产精品免费福利视频| 午夜91福利影院| 1024视频免费在线观看| 日韩欧美国产一区二区入口| 成人18禁在线播放| 91成年电影在线观看| 777久久人妻少妇嫩草av网站| 99香蕉大伊视频| 国产真人三级小视频在线观看| 免费在线观看日本一区| 免费女性裸体啪啪无遮挡网站| 国产亚洲av高清不卡| 99久久国产精品久久久| 精品一区二区三卡| 亚洲 欧美一区二区三区| 精品国产乱码久久久久久男人| 亚洲成人免费电影在线观看| 女人精品久久久久毛片| 精品乱码久久久久久99久播| 国产伦一二天堂av在线观看| 又大又爽又粗| 男女午夜视频在线观看| 丰满迷人的少妇在线观看| 国内毛片毛片毛片毛片毛片| 激情视频va一区二区三区| 成人国语在线视频| 国产精品一区二区三区四区久久 | 淫妇啪啪啪对白视频| 少妇裸体淫交视频免费看高清 | 满18在线观看网站| 悠悠久久av| 亚洲avbb在线观看| 色综合站精品国产| 少妇 在线观看| aaaaa片日本免费| 国产精品免费一区二区三区在线| 成人永久免费在线观看视频| 午夜精品在线福利| 美女午夜性视频免费| 日日爽夜夜爽网站| 嫁个100分男人电影在线观看| 一二三四在线观看免费中文在| 国内毛片毛片毛片毛片毛片| 日日摸夜夜添夜夜添小说| 成人亚洲精品一区在线观看| 麻豆一二三区av精品| 欧美日本中文国产一区发布| 亚洲情色 制服丝袜| 国产三级黄色录像| 国产免费男女视频| 久久久久久久久久久久大奶| 亚洲精品一卡2卡三卡4卡5卡| 国产成人影院久久av| 9热在线视频观看99| 久久久久久久精品吃奶| 亚洲精品一卡2卡三卡4卡5卡| 国产精华一区二区三区| 女人被躁到高潮嗷嗷叫费观| 热99re8久久精品国产| 91成人精品电影| 成熟少妇高潮喷水视频| 国产精品国产av在线观看| 国产日韩一区二区三区精品不卡| 老汉色av国产亚洲站长工具| 成人精品一区二区免费| 少妇裸体淫交视频免费看高清 | 亚洲专区国产一区二区| 欧美乱妇无乱码|