王飛余佳益劉顯龍蔡陽(yáng)健

1)(蘇州大學(xué)物理科學(xué)與技術(shù)學(xué)院,蘇州 215006)

2)(蘇州大學(xué),蘇州納米科技協(xié)同創(chuàng)新中心,蘇州 215006)

3)(山東師范大學(xué)物理與電子科學(xué)學(xué)院,光場(chǎng)調(diào)控及應(yīng)用中心,濟(jì)南 250358)

1 引 言

自由空間激光通信是以激光作為信息載體,在不同的平臺(tái)(地面平臺(tái)、空對(duì)地平臺(tái)等)之間開(kāi)展信息傳遞與交換的一種技術(shù).與傳統(tǒng)的無(wú)線電波通信相比,激光通信在很多方面具有一定的優(yōu)勢(shì):1)激光頻率高,具有超大帶寬,能極大地提高信息傳輸速率;2)激光方向性好,在一定程度上能提高信息傳輸?shù)陌踩?同時(shí)還具有較強(qiáng)的抗干擾能力;3)激光亮度高、發(fā)散角小,同時(shí),激光器體積小,易于集成,方便搭載于不同平臺(tái)上.除了自由空間光通信,以激光為載波的激光雷達(dá)、遠(yuǎn)程遙感也有獨(dú)特的優(yōu)勢(shì).然而,在這些應(yīng)用中,激光都是以自由空間(大氣)作為信息傳輸通道,大氣中的微粒和分子對(duì)激光的吸收和散射及湍流對(duì)激光造成的光斑漂移、光強(qiáng)閃爍等效應(yīng),嚴(yán)重影響通信系統(tǒng)的性能.因此,激光在大氣中的傳輸特性及大氣本身性質(zhì)的研究顯得尤為重要.

當(dāng)光波在大氣中傳輸時(shí),將受到大氣的影響發(fā)生吸收、散射及大氣折射率隨機(jī)變化導(dǎo)致的相位隨機(jī)起伏等現(xiàn)象.吸收和散射主要是由組成大氣的分子、大氣中的懸浮顆粒及氣溶膠等引起的,并使得光波在傳輸路徑上的能量發(fā)生衰減,吸收和散射的強(qiáng)度依賴于光波的波長(zhǎng).大氣不規(guī)則的隨機(jī)運(yùn)動(dòng),每一點(diǎn)上的壓強(qiáng)、速度、溫度等物理特性的隨機(jī)漲落通常稱為大氣湍流.它的存在導(dǎo)致大氣的折射率隨機(jī)起伏,湍流中傳輸?shù)墓獠ǖ恼穹拖辔灰搽S之發(fā)生無(wú)規(guī)則的起伏,是影響光波傳輸特性的主要因素.通常情況下,大氣湍流將引起光波發(fā)散角變大、光束展寬、光斑位置隨機(jī)抖動(dòng)、相干性退化以及到達(dá)角起伏等效應(yīng),嚴(yán)重影響自由空間光通信和激光雷達(dá)等系統(tǒng)的信息傳遞和交換性能.因此,國(guó)內(nèi)外研究人員開(kāi)展了大量激光束在大氣中的傳播行為及如何克服或抑制大氣湍流引起光波擾動(dòng)的研究工作.

從物理上來(lái)講,大氣湍流是一種非線性的隨機(jī)過(guò)程,其嚴(yán)格解可由著名的Navier-Stokes(N-S)方程導(dǎo)出,然而,從數(shù)學(xué)上嚴(yán)格求解N-S方程至今仍是一個(gè)難題.1941年,蘇聯(lián)學(xué)者Kolmogorov[1]發(fā)展了一種運(yùn)用統(tǒng)計(jì)學(xué)方法描述湍流的統(tǒng)計(jì)規(guī)律,雖然這種統(tǒng)計(jì)方法對(duì)實(shí)際復(fù)雜的湍流運(yùn)動(dòng)過(guò)程做了一些簡(jiǎn)化和近似,但是比起理論上嚴(yán)格求解N-S方程,這種方法顯得更為簡(jiǎn)易和方便.我們現(xiàn)在所知的光波在湍流大氣中的傳輸理論研究,幾乎都是基于Kolmogorov的湍流統(tǒng)計(jì)理論發(fā)展而來(lái)的.

另一方面,部分相干光理論經(jīng)過(guò)Zernike,Wolf,Gori及Friberg等幾代科學(xué)家的發(fā)展和完善,已經(jīng)形成了一套比較完善的理論框架.部分相干光束也已經(jīng)被廣泛地應(yīng)用在慣性約束核聚變、激光掃描、非線性光學(xué)等領(lǐng)域.高斯-謝爾模光束是部分相干光中最為熟知的一種模型(光強(qiáng)分布和相干度分布都滿足高斯分布),被廣泛應(yīng)用于各種理論和實(shí)際計(jì)算中.目前,研究部分相干光束傳輸?shù)闹饕椒ㄓ袕埩糠椒?、二階矩方法、高階矩方法、Wigner分布函數(shù)方法和數(shù)值計(jì)算等方法.

2000年以來(lái),受到自由空間光通信等應(yīng)用的推動(dòng),大氣湍流效應(yīng)抑制的研究成為一個(gè)重要的課題.目前已知的抑制技術(shù)和方法有:大孔徑接收技術(shù)、編碼調(diào)制技術(shù)、自適應(yīng)光學(xué)方法、部分相干光束作為信息載體傳輸方法等.作為一種重要的湍流效應(yīng)抑制方法,部分相干光束在湍流大氣中的傳輸特性的研究重新引起了國(guó)內(nèi)外研究人員的重視,已有大量關(guān)于部分相干光束在湍流大氣中傳輸特性的研究文獻(xiàn)發(fā)表,推動(dòng)了湍流理論和光傳輸理論的進(jìn)一步發(fā)展.特別是近幾年來(lái),光場(chǎng)調(diào)控逐漸成為光學(xué)領(lǐng)域的熱點(diǎn)研究課題之一,部分相干光束的相干結(jié)構(gòu)、相位和偏振調(diào)控理論及技術(shù)都得到了長(zhǎng)足的發(fā)展,多參量調(diào)控的部分相干光束與湍流大氣的相互作用也逐漸成為一個(gè)重要研究課題.

本文第2節(jié)回顧了部分相干光束在湍流大氣中傳輸研究的歷程和進(jìn)展;第3節(jié)簡(jiǎn)要介紹了部分相干光束在湍流中傳輸?shù)睦碚撃Ｐ图岸A統(tǒng)計(jì)特性,包括平均光強(qiáng)、相干度及光束發(fā)散角等二階矩的計(jì)算方法;第4節(jié)主要介紹部分相干光在湍流大氣中的相位屏數(shù)值模擬方法;最后簡(jiǎn)要展望了部分相干光束在湍流大氣中傳輸研究今后的發(fā)展方向.

2 部分相干光束在湍流大氣中傳輸研究回顧

光波在湍流介質(zhì)中傳輸理論的研究可以追溯到激光器發(fā)明之前,在1960年和1961年,蘇聯(lián)學(xué)者Chernov[2]和Tatarskii[3]分別發(fā)表了討論光波在湍流介質(zhì)中的光束展寬、光強(qiáng)閃爍、相干性演化等特性的著作.從那時(shí)起,各國(guó)學(xué)者對(duì)光波在湍流大氣中的傳輸行為進(jìn)行了深入的研究,詳細(xì)討論了光波在湍流中的強(qiáng)度分布、光束漂移、光斑展寬等物理特性.這期間主要的理論研究方法包括:Rytov微擾法、差分方程法、梯度近似法、幾何光學(xué)近似方法等[4?6],研究對(duì)象主要是球面波、平面波和高斯光波.1971年,Lutomirski和Yura[7]引入了一種新的方法,就是現(xiàn)在熟知的廣義惠更斯-菲涅耳積分法,處理有限孔徑的平面波在湍流介質(zhì)中的傳輸問(wèn)題.1972年,Yura[8]運(yùn)用這種方法得到了高斯光束在湍流介質(zhì)中傳輸后的互相干函數(shù),在一定條件下可以退化到平面波和球面波情況.廣義惠更斯-菲涅耳積分法可以方便地處理不同模式的光束(包括部分相干光束)在湍流中的平均光強(qiáng)、光束展寬、相干性退化等特性,因此目前仍被研究人員廣泛采用并用來(lái)研究部分相干光束在湍流大氣中的傳輸行為.實(shí)際上,早在1967年,蘇聯(lián)學(xué)者Feizulin和Kravtsov[9]就獨(dú)立地發(fā)展和運(yùn)用廣義惠更斯-菲涅耳方法研究了高斯光束在湍流大氣中的光束展寬.

1972年,Kon和Tatarskiis[10]研究了部分相干光束(高斯-謝爾模光束)在湍流大氣傳輸過(guò)程中有效相干長(zhǎng)度及光束的初始相干長(zhǎng)度對(duì)光束平均光強(qiáng)的影響,這可能是最早的關(guān)于部分相干光束在湍流大氣中傳輸行為的研究文獻(xiàn).1978年,美國(guó)學(xué)者Leader[11]基于廣義惠更斯-菲涅耳原理的瑞利-索末菲積分公式得到了高斯-謝爾模光束在湍流大氣中傳輸?shù)幕ハ喔珊瘮?shù)和平均光強(qiáng)的解析表達(dá)式,并研究了相干特性和光強(qiáng)在湍流中的演化.隨后,Wang和Plonus[12]結(jié)合廣義惠更斯-菲涅耳原理和Rytov相位結(jié)構(gòu)二次近似方法得到了在弱湍流條件下高斯-謝爾模光束的互相干函數(shù)的普適表達(dá)式,在他們的研究中,還考慮了初始光源準(zhǔn)直特性.除了部分相干光束的光束展寬、平均光強(qiáng)、相干性等二階統(tǒng)計(jì)的建模,研究人員也開(kāi)始關(guān)注理論上更難處理的部分相干光束在湍流介質(zhì)中的光強(qiáng)起伏(閃爍)特性的研究.1981年,Leader[13]研究了在弱湍流條件下任意尺寸大小的準(zhǔn)直或聚焦的部分相干光束在傳輸過(guò)程中光強(qiáng)閃爍因子,研究中運(yùn)用了廣義惠更斯-菲涅耳積分和湍流中球面波四階關(guān)聯(lián)統(tǒng)計(jì)的方法.隨后,Wang等[14]運(yùn)用和Leader類似的方法進(jìn)一步考慮了在接收面為有限孔徑(孔徑平滑效應(yīng))的條件下部分相干光束經(jīng)過(guò)弱湍流大氣傳輸后的光強(qiáng)閃爍因子.蘇聯(lián)學(xué)者Banach等[15,16]在此期間也研究了部分相干光束在湍流大氣中的光強(qiáng)閃爍,其研究分成兩種類型,即探測(cè)器的響應(yīng)時(shí)間快于和慢于部分相干光束本身光強(qiáng)變化速度,發(fā)現(xiàn)兩種類型對(duì)實(shí)驗(yàn)探測(cè)到的光強(qiáng)閃爍有重要的影響.然而在此期間,學(xué)者們的研究重點(diǎn)放在理論研究上,實(shí)驗(yàn)研究成果較少.理論研究發(fā)現(xiàn)在滿足“慢探測(cè)器”的前提下,部分相干光束與完全相干光束相比,能有效地降低湍流引起的光強(qiáng)閃爍,具有重要的應(yīng)用價(jià)值.在1990—2000年,部分相干光束在湍流大氣中傳輸特性的研究似乎突然被研究人員忽視了,只有少數(shù)的研究成果發(fā)表[17,18].

從2002年開(kāi)始,受到自由空間光通信、激光雷達(dá)等應(yīng)用發(fā)展的影響,部分相干光束在湍流介質(zhì)中的傳輸行為又重新引起了研究人員的興趣.Gbur和Wolf[19]從理論上研究了部分相干光束在湍流介質(zhì)中光束平均光斑半徑的演化,并從這個(gè)角度揭示了部分相干光束受到湍流的影響比相同條件下的相干光束小.隨后,Dogariu和Amarande[20]利用熱板模擬大氣湍流效應(yīng),從實(shí)驗(yàn)上驗(yàn)證了部分相干光束的平均光斑半徑展寬受到湍流的影響比相干光束小,并隨著光束相干長(zhǎng)度的減小而進(jìn)一步減小.Shirai等[21]利用部分相干光束的相干模式分解,研究了該光束在湍流介質(zhì)中的光束展寬和平均光強(qiáng)分布,并從模式分解角度解釋了該光束受到湍流擾動(dòng)較小的物理原因.同期,美國(guó)陸軍實(shí)驗(yàn)室Ricklin和Davidson[22,23]著眼于部分相干光束在自由空間通信中應(yīng)用的可能性,對(duì)部分相干光束在弱湍流中的平均光強(qiáng)、相干性演化、光強(qiáng)閃爍等特性做了進(jìn)一步研究,并從理論上分析了誤碼率和光束相干性的關(guān)系,得到低相干性光源能大幅降低誤碼率的結(jié)論.Korotkova等[24]利用激光束經(jīng)過(guò)隨機(jī)擾動(dòng)屏模擬部分相干光束,得到了計(jì)算光強(qiáng)閃爍因子的解析形式,并分析了該光束在湍流大氣中的光強(qiáng)閃爍因子和誤碼率的大小,也得到了與文獻(xiàn)[23]類似的結(jié)論.此后,國(guó)內(nèi)外學(xué)者對(duì)部分相干光束及其他類型的光束在湍流介質(zhì)中傳輸?shù)难芯客度肓藰O大的熱情,每年有大量的研究文獻(xiàn)發(fā)表.

這一時(shí)期,研究人員關(guān)注的重點(diǎn)是光束的模式、偏振態(tài)、相位等與湍流介質(zhì)相互作用引起的一些特異的傳輸特性.各種模式的光束,如拉蓋爾-高斯光束[25]、雙曲余弦高斯光束[26]、厄米-余弦-高斯光束[27],橢圓高斯光束[28],徑向偏振光束[29]、貝塞爾-高斯光束[30]、艾里光束[31]、渦旋光束[32]等,在湍流大氣中的平均光強(qiáng)、光束展寬、光斑漂移等特性得到了詳細(xì)的研究.盡管初始的光波模式不盡相同,研究人員得到的結(jié)論可以歸納為:1)光束在湍流大氣中傳輸足夠遠(yuǎn)的距離時(shí),光強(qiáng)分布逐漸退化成高斯分布,與初始光束的光強(qiáng)分布無(wú)關(guān),且湍流強(qiáng)度越大,退化成高斯分布的傳輸距離越短;2)光束發(fā)散角和光束展寬隨湍流強(qiáng)度的增加而增大.與此同時(shí),研究人員還把目光轉(zhuǎn)向各種類型的光束陣列在湍流大氣中平均光強(qiáng)和發(fā)散角特性[33?39],研究發(fā)現(xiàn)與單束光波相比,陣列光束也能有效抑制湍流引起的光強(qiáng)閃爍[40,41].在部分相干光束湍流傳輸方面,不再是考慮單純的高斯-謝爾模光束的傳輸特性,與完全相干光束類似,研究人員把光束的偏振態(tài)及模式與相干性結(jié)合起來(lái)考慮在湍流大氣中的平均光強(qiáng)演化、光束發(fā)散角等二階統(tǒng)計(jì)特性.相干性和光波模式結(jié)合各類高階部分相干光束,包括帶扭曲相位的部分相干高斯-謝爾模光束[42,43]、部分相干高階高斯光束[44?51]、部分相干平頂及空心光束[52?55]等的平均光強(qiáng)和光束發(fā)散的統(tǒng)計(jì)特性得到了廣泛的研究.同時(shí),部分相干部分偏振光束(隨機(jī)電磁光束)在湍流大氣中的偏振特性演化、平均光強(qiáng)演化特性也得到了深入的研究[56?67].一個(gè)比較有意思的研究結(jié)果是隨機(jī)電磁光束在湍流大氣中傳輸時(shí),光束在遠(yuǎn)場(chǎng)的偏振度的值會(huì)恢復(fù)到光源面的偏振度的值[58],但這個(gè)結(jié)果至今還未得到實(shí)驗(yàn)證實(shí).

由于光強(qiáng)閃爍與信息傳輸?shù)恼`碼率直接相關(guān)(光強(qiáng)閃爍越小,誤碼率越低),因此光強(qiáng)閃爍因子的研究具有更實(shí)際的意義.從理論上來(lái)講,光束光強(qiáng)閃爍因子的計(jì)算比光束平均光強(qiáng)、發(fā)散角等二階統(tǒng)計(jì)要復(fù)雜得多,而且基于廣義惠更斯-菲涅耳原理和Rytov四階復(fù)相位統(tǒng)計(jì)理論得到的結(jié)果與從近軸波動(dòng)方程出發(fā)得到結(jié)果是不相符的.一般認(rèn)為從近軸波動(dòng)方程出發(fā)得到的結(jié)果更為接近實(shí)際結(jié)果,所以基于廣義惠更斯-菲涅耳原理得到的部分相干光強(qiáng)閃爍因子結(jié)果只能從一定程度上得到一些正確的結(jié)論.但由于這種方法能夠得到近似的解析解,還是被很多研究人員采用.基于此理論,人們研究了一些不同類型的部分相干光束,如部分相干平頂光束[68]、隨機(jī)電磁光束[69]、部分相干多高斯光束[70]、貝塞爾關(guān)聯(lián)光束[71]、帶扭曲相位部分相干光[72]等的光強(qiáng)閃爍因子.研究結(jié)果表明,這些部分相干高階光束的光強(qiáng)閃爍因子在相同條件下比高斯-謝爾模光束的低.但需要注意的是,部分相干高階光束的發(fā)散角在同等條件下比高斯-謝爾模光束的大,似乎是犧牲了接收端的能量達(dá)到降低閃爍因子的效果,這個(gè)結(jié)果在完全相干高階高斯光束中也適用.

近年來(lái),光場(chǎng)調(diào)控成為光學(xué)領(lǐng)域的熱點(diǎn)研究課題之一,部分相干光束的相干調(diào)控研究引起了國(guó)內(nèi)外眾多研究人員的興趣,一系列具有不同相干結(jié)構(gòu)的部分相干光束模型被理論建模[73].部分相干光束相干結(jié)構(gòu)與湍流介質(zhì)相互作用引發(fā)的一些新特性及新效應(yīng)的研究成為一個(gè)重要的課題.各類具有特殊相干結(jié)構(gòu)的部分相干光束在湍流大氣中傳輸?shù)钠骄鈴?qiáng)、質(zhì)量因子、光束擴(kuò)散等特性得到了深入的研究[74?87].由于受到湍流和相干結(jié)構(gòu)兩個(gè)因素的影響,這類部分相干光束在湍流介質(zhì)中的平均光強(qiáng)首先會(huì)演化成一些特殊的形狀,如平頂分布、空心分布等,然后逐漸退化為類高斯分布.值得指出的是,在同等條件下,特殊相干結(jié)構(gòu)的部分相干光束的發(fā)散角比傳統(tǒng)的高斯-謝爾模光束要大.這類特殊相干結(jié)構(gòu)的部分相干光束在湍流介質(zhì)中的光強(qiáng)閃爍因子更值得關(guān)注,尤其是非均勻關(guān)聯(lián)光束,它在自由空間中傳輸會(huì)出現(xiàn)類似于“自聚焦”的現(xiàn)象[88].Gu和Gbur[89]研究了弱湍流中非均勻關(guān)聯(lián)光束的光強(qiáng)閃爍因子,結(jié)果發(fā)現(xiàn)在合適的參數(shù)下,該光束比傳統(tǒng)的高斯光束不僅具有更低的閃爍因子,而且中心區(qū)域的光能量更大,這就解決了部分相干光束發(fā)散角過(guò)大導(dǎo)致光能量發(fā)散較快的問(wèn)題.我們課題組在理論上研究了多高斯關(guān)聯(lián)光束在弱湍流大氣中的光強(qiáng)閃爍因子,發(fā)現(xiàn)這種光束也能有效降低閃爍因子[90],該實(shí)驗(yàn)結(jié)果被美國(guó)邁阿密大學(xué)Korotkova及美國(guó)海軍學(xué)院Avramov-Zamurovic教授在戶外實(shí)驗(yàn)中驗(yàn)證[91,92].

總體而言,近幾年有關(guān)激光束經(jīng)過(guò)湍流大氣傳輸?shù)难芯课墨I(xiàn)主要關(guān)注不同初始參數(shù)(如偏振、相位和相干結(jié)構(gòu))的激光束與湍流介質(zhì)相互作用后的統(tǒng)計(jì)特性,研究的目的一方面是探討不同類型的光源與湍流相互作用的基本物理特性,更為重要的是期望調(diào)控光源的參數(shù)達(dá)到一些更優(yōu)的效果,比如更低的光強(qiáng)閃爍、更高的能量接收效率等.

3 部分相干光束經(jīng)過(guò)湍流大氣的傳輸模型

3.1 部分相干光經(jīng)過(guò)湍流大氣傳輸?shù)钠骄鈴?qiáng)和相干性統(tǒng)計(jì)特性

大氣層中空氣無(wú)規(guī)則的運(yùn)動(dòng)導(dǎo)致空氣的氣壓、溫度、密度等無(wú)規(guī)則起伏稱為大氣湍流.而大氣的折射率取決于上述參數(shù),因此大氣的折射率也隨空間和時(shí)間做無(wú)規(guī)則的變化,光波在湍流中的效應(yīng)主要表現(xiàn)為強(qiáng)度起伏(光強(qiáng)閃爍)、相位起伏(相位的空間、時(shí)間起伏)和光束漂移(方向起伏)等.一般而言,湍流屬于弱散射介質(zhì)范疇,因此光束的后向散射效應(yīng)及偏振退化(耦合)效應(yīng)可以忽略不計(jì);再者,我們還需要假設(shè)在光束傳輸方向上折射率起伏的關(guān)聯(lián)函數(shù)服從狄拉克函數(shù)分布.運(yùn)用上述假設(shè),從經(jīng)典電磁理論出發(fā),光波在湍流中的標(biāo)量亥姆霍茲方程表示為

式中,U(R)表示光束的某一個(gè)電場(chǎng)分量;R≡(x,y,z)是空間三維坐標(biāo),其中光束假定沿z方向傳播;k=2π/λ是光波的波數(shù),λ是光波波長(zhǎng);n(R)=n0+n1(R)表示湍流介質(zhì)的折射率分布,其中n0≈1,n1(R)表示折射率的起伏變化,?n1(R)?≈0,角括號(hào)表示對(duì)系綜(時(shí)間)的平均.根據(jù)湍流統(tǒng)計(jì)理論,折射率起伏在空間兩點(diǎn)R1和R2的協(xié)方差函數(shù)可以表示為

式中,R≡(r,z),r≡(x,y)表示與z軸垂直平面內(nèi)的笛卡爾坐標(biāo);Bn表示折射率起伏協(xié)方差(相關(guān))函數(shù),根據(jù)Kolmogorov湍流統(tǒng)計(jì)理論,在湍流慣性區(qū)域,協(xié)方差函數(shù)正比于兩點(diǎn)距離差的2/3次方.等式最右邊運(yùn)用了馬爾可夫近似(Markov approximation),即:在光波傳輸方向上,折射率起伏服從狄拉克函數(shù)分布,An(r1?r2)是二維的協(xié)方差函數(shù).

方程(1)最為人熟知的兩種處理方法是Born微擾法和Rytov擾動(dòng)方法.Born微擾法是把光波電場(chǎng)分量寫(xiě)成各級(jí)散射波的求和:

式中,U0(R)表示未被擾動(dòng)(散射)的電場(chǎng)分量,U1(R)和U2(R)分別表示一級(jí)和二級(jí)散射分量.與此相對(duì)應(yīng)的是,Rytov擾動(dòng)方法是考慮電場(chǎng)分量受到一個(gè)復(fù)相位的擾動(dòng),寫(xiě)成如下形式:

式中,ψ(R)表示由湍流引起的復(fù)相位擾動(dòng),也可寫(xiě)成如下級(jí)數(shù)形式,

ψ1(R)和ψ2(R)分別表示一級(jí)和二級(jí)復(fù)相位的擾動(dòng).對(duì)比(3)和(4)式,Born和Rytov方法的不同點(diǎn)是:Born方法把擾動(dòng)場(chǎng)表示成求和形式,而Rytov方法把擾動(dòng)場(chǎng)表示為“相乘”形式.以往的研究中,Born方法成功地用于各種弱散射物體散射場(chǎng)的計(jì)算,但是在湍流領(lǐng)域普遍使用的是Rytov方法研究光波經(jīng)過(guò)湍流介質(zhì)的統(tǒng)計(jì)特性.究其原因,可能是Born方法更適用于短距離的散射問(wèn)題,而Rytov方法更適用于光波長(zhǎng)距離傳輸?shù)奈_.早期的關(guān)于光波在湍流介質(zhì)中的傳輸研究大多是基于Kolmogorov湍流統(tǒng)計(jì)理論,借助于Rytov方法從亥姆霍茲(近軸傳輸方程)出發(fā)討論光波的統(tǒng)計(jì)性質(zhì),得到了一些非?；A(chǔ)性的結(jié)果[2,3],但是這種方法的數(shù)學(xué)要求較高,對(duì)于復(fù)雜高階光束以及部分相干光束,很難從數(shù)學(xué)上嚴(yán)格的求解.

進(jìn)入21世紀(jì)以來(lái),研究人員更傾向于用Lutomiriski和Yura[7]發(fā)展的廣義惠更斯-菲涅耳積分方法處理光波在湍流大氣中的統(tǒng)計(jì)特性.這種方法的優(yōu)點(diǎn)是可以方便地處理各類不同的光波模式,如高階高斯光束、平頂光束等復(fù)雜激光束在湍流介質(zhì)中的特性,而且在弱湍流和強(qiáng)湍流條件下,所得到的平均光強(qiáng)、光束發(fā)散角等性質(zhì)能較好地符合實(shí)驗(yàn)結(jié)果[93],甚至在做了二次方近似之后仍與實(shí)驗(yàn)結(jié)果較好地符合[94].

圖1是光波在湍流介質(zhì)中從入射面到出射面的示意圖.根據(jù)廣義惠更斯-菲涅耳原理,在湍流介質(zhì)中入射面和出射面上的光波電場(chǎng)滿足如下積分關(guān)系:

式中,L表示從入射面到出射面的傳輸距離;r和ρ是與傳輸軸垂直的入射面和出射面的位置矢量;U0(r,0)和U(ρ,L)分別表示入射面和出射面上光波的電場(chǎng)分量;ψ(r,ρ,L)表示的是從入射面到出射面由湍流介質(zhì)引起的Rytov復(fù)相位擾動(dòng),根據(jù)惠更斯-菲涅耳原理,在出射面上的光波是入射面上每一點(diǎn)發(fā)出的球面波的相干疊加,所以ψ(r,ρ,L)實(shí)際代表的是球面波情況下Rytov復(fù)相位擾動(dòng).在實(shí)際應(yīng)用中,人們更關(guān)注的是光波的平均光強(qiáng)、發(fā)散角等物理量,對(duì)(6)式取復(fù)共軛,然后乘(6)式,取系綜平均后得到

式中,是光波任意兩點(diǎn)之間的系綜平均,在部分相干光范疇內(nèi),通常稱為互相干函數(shù),包含了光波的相干信息,如果是完全相干光,不需要進(jìn)行系綜平均運(yùn)算;是出射面上的互相干函數(shù),包含光束的平均光強(qiáng)、相干性等信息;等式最后一項(xiàng)尖括號(hào)內(nèi)表示的是由湍流引起的復(fù)相位擾動(dòng)的二階統(tǒng)計(jì)平均;尖括號(hào)中的下標(biāo)“s”和“m”分別表示光源本身的振幅和相位的起伏和湍流引入的復(fù)相位二階統(tǒng)計(jì).根據(jù)Kolmogorov湍流統(tǒng)計(jì)理論和Rytov微擾(只考慮到二階微擾),(7)式最后一項(xiàng)可以表示為[95]

其中,ρd=ρ1? ρ2和rd=r1?r2分別表示出射和入射面上的位置矢量差;κ是空間頻率矢量κ≡(κx,κy)的模,κ與空間坐標(biāo)r組成傅里葉變換對(duì);J0是零階貝塞爾函數(shù);Φn(κ)表示湍流譜密度,是(2)式協(xié)方差函數(shù)的傅里葉變換.其中最為人熟知的湍流譜密度函數(shù)是Kolmogorov譜,它的表達(dá)式為

式中,是折射率結(jié)構(gòu)常數(shù),是衡量湍流強(qiáng)弱的一個(gè)參數(shù),通常情況下,當(dāng)時(shí),稱為“弱湍流”,當(dāng)時(shí),稱為“強(qiáng)湍流”這種對(duì)湍流強(qiáng)弱的劃分方法也不是絕對(duì)的,因?yàn)橥牧鞯膹?qiáng)弱也與在湍流中的傳輸距離密切相關(guān);L0和l0是湍流的外尺度和內(nèi)尺度常數(shù).(11)式適用的條件只是在湍流的慣性區(qū)域,即譜空間頻率值在內(nèi)尺度和外尺度數(shù)的倒數(shù)之間.如果假定外尺度為無(wú)窮大,湍流譜在κ非常小時(shí)發(fā)散(趨于無(wú)窮大),為了在計(jì)算上更好地應(yīng)用,Tataskii和von Karman對(duì)湍流譜進(jìn)行了數(shù)學(xué)上的擴(kuò)展,使之能包含全部的空間頻率,又不含有奇點(diǎn).現(xiàn)在常用的湍流為von Karman湍流譜,其表達(dá)式為

其中κm=5.92/l0,κ0=C0/L0,C0的值根據(jù)具體條件可以取不同的值,一般為2π—8π,其他湍流譜的模型可以參考文獻(xiàn)[96].在現(xiàn)實(shí)環(huán)境中,大氣湍流在不同的地理區(qū)域,同一天不同的時(shí)間段及不同的氣象條件下都是不同的,而且從物理上而言,它的統(tǒng)計(jì)特性并不是嚴(yán)格靜態(tài)的,統(tǒng)計(jì)特性隨時(shí)間的變化而變化,所以用固定的湍流譜模擬光波的統(tǒng)計(jì)性質(zhì)只能從一定程度上反映光波的變化.(12)式的湍流譜還可以擴(kuò)展到non-Kolmogorov情況.需要指出的是,在得到(9)和(10)式的前提是運(yùn)用了馬爾可夫近似條件及假定湍流是各向同性的.把(9)和(10)式代入(7)式,一般來(lái)說(shuō)是很難得到解析表達(dá)式的.為了得到近似的解析表達(dá)式,還需要引入一些近似條件,對(duì)(10)式做進(jìn)一步簡(jiǎn)化.一種常用的簡(jiǎn)化是認(rèn)為(10)式中的積分函數(shù)只在貝塞爾函數(shù)的變量κ|(1?ξ)ρd+ξrd|?1范圍內(nèi)有值,其他區(qū)域的值趨于0,在這個(gè)條件下,貝塞爾函數(shù)可以只取泰勒級(jí)數(shù)的前兩項(xiàng)作為近似值,即:J0(x)≈1?x2/4.把這個(gè)條件代入(10)式,積分后可以得到如下解析式:

運(yùn)用(13)式,大部分的部分相干光束模型在湍流介質(zhì)中的平均光強(qiáng)和相干性分布都能得到解析形式,可以比較方便地分析這類光束的平均光強(qiáng)及相干性的演化過(guò)程.值得注意的是,運(yùn)用貝塞爾函數(shù)級(jí)數(shù)展開(kāi)的前兩項(xiàng)作為近似是在強(qiáng)湍流條件下得到的[21],現(xiàn)有的部分文獻(xiàn)并不遵照這個(gè)條件使用級(jí)數(shù)展開(kāi)研究部分相干光束的二階統(tǒng)計(jì)特性,從一定程度上來(lái)說(shuō)會(huì)引入一定的誤差.另外一種在大氣湍流介質(zhì)中常用的二階復(fù)相位擾動(dòng)統(tǒng)計(jì)模型是基于球面波在湍流中的波結(jié)構(gòu)函數(shù)并結(jié)合二次近似條件[22]:

其中,D是從入射面到出射面的球面波結(jié)構(gòu)函數(shù);()是球面波在湍流大氣中傳輸?shù)南喔砷L(zhǎng)度.如果湍流譜為Kolmogorov譜的形式,從數(shù)學(xué)上得到的嚴(yán)格解是ρ0的5/3次方,而不是2次方的形式,但是為了計(jì)算的方便,大多數(shù)文獻(xiàn)采用了二次方近似條件.由二次方近似而引起的計(jì)算誤差一般都在可接受的范圍內(nèi).

圖1 光波從入射面經(jīng)過(guò)湍流大氣到達(dá)出射面的示意圖Fig.1.Schematic for the propagation of light beams through turbulent atmosphere from input plane to output plane.

3.2 部分相干光束在湍流大氣中的平均光斑寬度及發(fā)散角

光波在湍流介質(zhì)中傳輸?shù)墓馐骄獍邔挾燃鞍l(fā)散角的變化能提供光束能量發(fā)散快慢的信息,同時(shí)也能計(jì)算光束在傳輸過(guò)程中的質(zhì)量因子的演化.光束的平均光斑寬度和發(fā)散角根據(jù)定義可以從光束的平均光強(qiáng)信息中直接得到,更全面地計(jì)算部分相干光束在湍流大氣中階矩的統(tǒng)計(jì)信息,需要借助于部分相干光束在湍流介質(zhì)中的Wigner函數(shù).四川大學(xué)張彬課題組[97]得到了部分相干光束在湍流介質(zhì)中傳輸過(guò)程中光束二階矩矩陣,可以方便地得到各類光束的任意階矩信息.光束在傳輸過(guò)程中的平均光斑寬度定義為[21]

根據(jù)文獻(xiàn)[19],光束的平均半徑可以表示成如下普適的形式:

(17)式右端第一項(xiàng)表示光束的初始(光源面)的光斑寬度;第二項(xiàng)表示的是光束在自由空間中傳輸過(guò)程中由衍射引入的光斑發(fā)散,與傳輸距離的平方成正比,是自由空間中遠(yuǎn)場(chǎng)發(fā)散角;第三項(xiàng)表示的是由湍流介質(zhì)引入的光束的額外發(fā)散,與傳輸距離的三次方成正比,其中

(18)式表明第三項(xiàng)只與湍流特性有關(guān),而與光束的初始狀態(tài)無(wú)任何關(guān)聯(lián).這就是說(shuō),在傳輸距離足夠大或者湍流強(qiáng)度足夠強(qiáng)(F2足夠大),光束的平均光斑半徑只取決于第三項(xiàng),即湍流引起的光束擴(kuò)散,而且?guī)缀跛泄馐钠骄霃蕉际且恢碌?與光束的初始參數(shù)無(wú)關(guān).如果初始光斑是橢圓形狀,在湍流大氣中橢圓率會(huì)逐漸減小,最終變成近似圓形.

要得到部分相干光束任意階矩的信息,最方便的方法是借助于部分相干光束的Wigner分布函數(shù):

其中,ρs=(ρ1+ρ2)/2,Γ(ρs,ρd,z)是部分相干光束在出射面上的互相干函數(shù),θ是與部分相干光束平均發(fā)散角相關(guān)的一個(gè)量.從Wigner函數(shù)出發(fā),光束沿x軸和y軸的平均光斑半徑及發(fā)散角階矩定義為

P表示光束所攜帶的總能量.根據(jù)(20)式,部分相干光束的任意階矩都可以求得.利用Wigner函數(shù)分布具體計(jì)算部分相干空心光束、部分相干平頂光束及扭曲部分相干光束的二階矩的方法,可參考文獻(xiàn)[55,98,99].

3.3 部分相干光束在湍流大氣中傳輸?shù)墓鈴?qiáng)閃爍

光波在湍流介質(zhì)中的傳輸效應(yīng),除了光束發(fā)散角變大、相干性退化、光斑位置的隨機(jī)漂移等,最主要的是光強(qiáng)隨時(shí)間發(fā)生隨機(jī)變化(光強(qiáng)起伏).晚上我們觀察某些星星時(shí),會(huì)發(fā)現(xiàn)星星“眨眼睛”現(xiàn)象,這就是星星“發(fā)出”的光經(jīng)過(guò)地球大氣湍流引起的光強(qiáng)起伏.在自由空間光通信中,這個(gè)由湍流引起的隨機(jī)光強(qiáng)起伏是導(dǎo)致系統(tǒng)誤碼率上升的主要因素.通常,用閃爍因子定義光強(qiáng)起伏的強(qiáng)弱:

式中,是接收面上ρ點(diǎn)處的瞬時(shí)光強(qiáng);是平均光強(qiáng);尖括號(hào)表示對(duì)系綜的平均.從(22)式看出,閃爍因子是光強(qiáng)起伏的歸一化協(xié)方差.

對(duì)于部分相干光束來(lái)說(shuō),需要考慮三個(gè)特征時(shí)間的大小:1)部分相干光束自身光強(qiáng)擾動(dòng)的特征時(shí)間τs;2)大氣湍流介質(zhì)隨時(shí)間變化快慢的特征時(shí)間τa;3)探測(cè)器探測(cè)光強(qiáng)的積分時(shí)間τd.這三個(gè)特征時(shí)間的大小,將會(huì)大幅度影響探測(cè)到實(shí)際的光強(qiáng)閃爍因子的大小.如果τd?τs,探測(cè)器“感受”不到部分相干光束自身的光強(qiáng)擾動(dòng),如果τd?τs,那么探測(cè)器能感知部分相干光束自身引起的光強(qiáng)擾動(dòng).據(jù)此可以把探測(cè)器分成兩類,一類是“快探測(cè)器”,滿足條件τd?τs?τa,即探測(cè)器的探測(cè)時(shí)間要遠(yuǎn)小于部分相干光束自身光強(qiáng)擾動(dòng)時(shí)間和湍流介質(zhì)變化的時(shí)間,它可以同時(shí)探測(cè)到湍流引起的和光源自身擾動(dòng)導(dǎo)致的光強(qiáng)變化;另一類是“慢探測(cè)器”,滿足τs?τd?τa,即探測(cè)器的探測(cè)時(shí)間介于湍流介質(zhì)變化特征時(shí)間和部分相干光束自身光強(qiáng)擾動(dòng)時(shí)間,這里部分相干光束自身引起的光強(qiáng)變化時(shí)間要遠(yuǎn)遠(yuǎn)小于湍流介質(zhì)變化時(shí)間,所以探測(cè)器測(cè)量只是湍流介質(zhì)引起的光強(qiáng)起伏.現(xiàn)階段大部分研究部分相干光束在湍流介質(zhì)中的光強(qiáng)閃爍特性的文獻(xiàn)都是基于“慢探測(cè)器”原理而建立起來(lái)的理論框架進(jìn)行的.基于廣義惠更斯-菲涅耳原理,部分相干光束在湍流介質(zhì)傳輸過(guò)程中入射面和出射面上的四階統(tǒng)計(jì)特性的關(guān)系如下:

(23)式積分中尖括號(hào)下標(biāo)“s”和“m”分別表示對(duì)部分相干光束自身引起光強(qiáng)起伏的平均和湍流介質(zhì)復(fù)相位擾動(dòng)引起的統(tǒng)計(jì)平均.在慢探測(cè)器理論框架下,并假定部分相干光束的復(fù)振幅變化滿足高斯統(tǒng)計(jì),這樣我們可以把部分相干光束的四階統(tǒng)計(jì)簡(jiǎn)化為

其中Γs為入射面上光束的互相干函數(shù).湍流介質(zhì)引起的復(fù)相位四階統(tǒng)計(jì)可以簡(jiǎn)化為

E1和E2表達(dá)式見(jiàn)(9)和(10)式,E3表達(dá)式為

但是,(25)式過(guò)于復(fù)雜,利用(25)式很難得到部分相干光束在湍流介質(zhì)傳輸過(guò)程中的解析表達(dá)式.一種可行的方法是對(duì)(25)式簡(jiǎn)化,得到一些近似表達(dá)式.目前較多的是在弱湍流近似下,(25)式可以近似表示為

式中,1,2,3,4)是球面波的波結(jié)構(gòu)函數(shù);是振幅對(duì)數(shù)相位結(jié)構(gòu)函數(shù);是球面波對(duì)數(shù)振幅和相位的相干長(zhǎng)度;是對(duì)數(shù)振幅關(guān)聯(lián)函數(shù);是球面波對(duì)數(shù)振幅方差.基于(27)式,可以計(jì)算一些相對(duì)簡(jiǎn)單的部分相干光束模型的光強(qiáng)閃爍因子[41,72,100],但是對(duì)于一些復(fù)雜的部分相干光束,依靠(27)式還是很難得到解析式.

除了上述方法之外,另外一種計(jì)算部分相干光束在湍流中的閃爍因子的方法是在完全相干光束的出射面上放置一個(gè)滿足高斯統(tǒng)計(jì)的隨機(jī)相位屏,結(jié)合湍流介質(zhì)的譜特性,這種方法成功地得到了高斯-謝爾模光束和隨機(jī)電磁光束的光強(qiáng)閃爍因子[24,69].從原理上來(lái)說(shuō),這種方法可推廣到其他更為復(fù)雜的部分相干光束情況,但是這個(gè)推廣至今還未見(jiàn)文獻(xiàn)報(bào)道.2006年,Berman和Chumak[101]報(bào)道了另外一種光子流概率密度的方法計(jì)算部分相干光束的閃爍因子.

近幾年,本課題組開(kāi)展了一些部分相干光束經(jīng)過(guò)實(shí)驗(yàn)室模擬湍流后的光強(qiáng)閃爍因子的實(shí)驗(yàn)研究.相比于同等條件下的部分相干光束,附加一個(gè)渦旋相位(部分相干渦旋光束)后能有效地降低光束經(jīng)過(guò)湍流后的光強(qiáng)閃爍因子,并且光強(qiáng)閃爍因子還隨著渦旋拓?fù)浜蓴?shù)的增加而減小[102].除了相位的影響,我們還研究了部分相干徑向偏振光束經(jīng)過(guò)湍流介質(zhì)后的閃爍因子,與渦旋相位類似,疊加徑向偏振后,光強(qiáng)閃爍因子也進(jìn)一步減小[103].

4 部分相干光經(jīng)過(guò)湍流大氣傳輸?shù)南辔黄翑?shù)值模擬方法

4.1 隨機(jī)相位屏模擬完全相干光束在大氣湍流中的傳輸

在第3節(jié),我們主要回顧了基于廣義惠更斯-菲涅耳原理和Rytov微擾理論計(jì)算部分相干光束經(jīng)過(guò)湍流介質(zhì)的傳輸特性.由于Rytov方法是基于長(zhǎng)時(shí)間統(tǒng)計(jì)平均后的特性,所以只能計(jì)算光波在湍流介質(zhì)中的平均統(tǒng)計(jì)特性,即平均光強(qiáng)、光斑半徑等,而無(wú)法得知光束“瞬時(shí)”的強(qiáng)度特性.根據(jù)實(shí)驗(yàn)測(cè)量的數(shù)據(jù),高斯光束經(jīng)過(guò)湍流介質(zhì)傳輸后出射面上的光強(qiáng)是偏離高斯形狀的,會(huì)發(fā)生光斑變形、扭曲,在強(qiáng)湍流條件下,光斑會(huì)瞬時(shí)分裂成幾個(gè)或多個(gè)散斑.

多相位屏是目前運(yùn)用最廣泛的模擬光波經(jīng)過(guò)湍流大氣傳輸?shù)臄?shù)值模擬方法,該方法不僅可以得到光束的平均光強(qiáng)、相干性演化及光強(qiáng)閃爍因子等二階和四階統(tǒng)計(jì)特性,還能得到光束的瞬時(shí)光強(qiáng)分布.在這一節(jié)中,先介紹多相位屏方法的原理,然后再探討把多相位屏方法運(yùn)用到部分相干光束.相位屏方法的原理和方法在文獻(xiàn)[104]中已經(jīng)有詳細(xì)的描述和介紹,先簡(jiǎn)單回顧該方法的原理和流程.從(1)式所表達(dá)的湍流中的亥姆霍茲方程出發(fā),如果考慮的光束是沿z方向傳輸?shù)陌S光波,那么U(R)可以寫(xiě)成如下形式:

其中,R≡(x,y,z),V(R)是關(guān)于z的緩變函數(shù).把(28)式代入(1)式,運(yùn)用湍流折射率起伏n(R)=n0+n1(R)關(guān)系,經(jīng)過(guò)化簡(jiǎn)得到如下傍軸近似下光束傳輸方程

其中,(R)是湍流折射率起伏,其平均值為零.從(29)式可以看出湍流的擾動(dòng)來(lái)自于第3項(xiàng).如果只考慮自由空間光波的傳輸方程(第3項(xiàng)為零):

其中,r≡(x,y).對(duì)(30)式中的V(r,z)做關(guān)于r的傅里葉變換,得到

其中,κ≡(κx,κy)是r對(duì)應(yīng)的空間頻率;函數(shù)F是V的傅里葉變換.(31)式從z′面?zhèn)鬏數(shù)絲平面的一個(gè)解可以表示為

根據(jù)幾何光學(xué),在實(shí)空間相位隨機(jī)擾動(dòng)的關(guān)聯(lián)函數(shù)可以表示成

圖2 光波在多相位屏中傳輸?shù)氖疽鈭DFig.2.Schematic for light beams propagating in the multi-random phase screen.

把(2)式代入(33)式,利用狄拉克函數(shù)的積分特性得到:

An(r1?r2)是與傳輸軸垂直的平面內(nèi)湍流折射率擾動(dòng)的協(xié)方差(關(guān)聯(lián))函數(shù),與湍流譜密度的關(guān)系為

根據(jù)(33)—(35)式,可以得到湍流中相位關(guān)聯(lián)譜密度與折射率關(guān)聯(lián)譜密度(相位關(guān)聯(lián)函數(shù)的傅里葉變換)的關(guān)系:

根據(jù)(36)式,只要給定湍流的折射率譜密度函數(shù),就能通過(guò)數(shù)值方法得到相位屏上的隨機(jī)相位.具體計(jì)算方法如下:首先產(chǎn)生N×N行和列的滿足方差為1的高斯統(tǒng)計(jì)隨機(jī)數(shù),然后產(chǎn)生的隨機(jī)數(shù)矩陣與相位關(guān)聯(lián)譜密度的平方根相乘,其中?k=2π/(N?),?為離散的空間頻率的最小間隔(一個(gè)像素點(diǎn)大小).最后對(duì)相乘后的矩陣進(jìn)行二維離散傅里葉變換得到復(fù)數(shù)相位的矩陣θ1+iθ2.θ1和θ2任取一個(gè)都可以作為相位屏的隨機(jī)相位矩陣.這樣相位屏數(shù)值模擬方法的具體流程如下:

1)對(duì)發(fā)射面上光束的電場(chǎng)做傅里葉變化,通過(guò)自由空間傳輸?shù)竭_(dá)第1個(gè)相位屏前,

2)光束穿過(guò)第1個(gè)隨機(jī)相位屏,受到相位屏的相位擾動(dòng),

3)對(duì)電場(chǎng)做傅里葉變換,光束再通過(guò)自由空間傳輸?shù)竭_(dá)第2個(gè)相位屏前,

4)光束穿過(guò)第2個(gè)隨機(jī)相位屏,受到相位屏的相位擾動(dòng),

如此循環(huán),直到通過(guò)所有的湍流隨機(jī)相位屏.上述各式中,?z=z/m,z是光束的傳輸距離,m為湍流隨機(jī)相位屏的個(gè)數(shù);θi(ri)是第i張湍流隨機(jī)相位屏的擾動(dòng)相位.到達(dá)接收面的光波的電場(chǎng)分布認(rèn)為是通過(guò)湍流介質(zhì)后的瞬時(shí)電場(chǎng),要用數(shù)值模擬法得到光波的平均統(tǒng)計(jì)特性,需要生成n組(每一組包含m個(gè)相位屏)相位屏組合放置在傳輸路徑上,n組之間的電場(chǎng)認(rèn)為是不相干的,利用這n組電場(chǎng)就可以計(jì)算平均光強(qiáng)、相干性變化及光強(qiáng)閃爍因子等特性.

在數(shù)值模擬中,相位屏間隔的距離?z之內(nèi)屬于弱湍流的范疇,一般來(lái)說(shuō),是用平面波在?z距離之內(nèi)的Rytov方差(平面波的閃爍因子)0.1來(lái)簡(jiǎn)單的判定.并且,另一個(gè)需要滿足的條件是在傳輸距離?z內(nèi)Rytov方差是整個(gè)傳輸距離Rytov方差的10%以內(nèi).在合成相位屏幕過(guò)程中,由于Kolmogorov湍流譜:在κ→0時(shí)趨于無(wú)窮大,在零頻附近,湍流譜隨κ的變化而急劇變化,而數(shù)值模擬中一般是等間隔采樣的,這樣導(dǎo)致在零頻附近湍流譜的誤差增大,為了減小在零頻處的誤差,研究人員提出了一些采樣方法來(lái)修正零頻附近的誤差[104,105].

4.2 部分相干光束經(jīng)過(guò)湍流大氣傳輸?shù)南辔黄聊M方法

與完全相干光不同,部分相干光束一般是用交叉譜密度函數(shù)或者互相干函數(shù)來(lái)描述它的統(tǒng)計(jì)特性,而不是直接以光束的電場(chǎng)來(lái)描述,所以相位屏模擬方法無(wú)法直接應(yīng)用到部分相干光束的情況.一種可行方法是相干模式分解方法,即把部分相干光束看成是不同模式的完全相干光束的非相干疊加[106],

式中,W是部分相干光束的交叉譜密度;ω是光波的圓頻率;λn(ω)是非負(fù)的本征值,代表模式的權(quán)重;En(r,ω)表示其中一個(gè)模式的電場(chǎng).利用(41)式,把組成部分相干光束的每一個(gè)模式分量經(jīng)過(guò)同樣的n個(gè)相位屏組傳輸進(jìn)行數(shù)值模擬計(jì)算,然后非相干疊加,得到部分相干光束的統(tǒng)計(jì)特性.然而,對(duì)任意的部分相干光束進(jìn)行(40)式的模型分解并得到對(duì)應(yīng)的解析解還有一定的難度,目前只有少數(shù)幾種部分相干光束、如高斯-謝爾模光束,帶扭曲相位的高斯-謝爾模光束有解析形式的模式分解[107,108].

另外一種方法是在完全相干光束的光源面上設(shè)置一個(gè)滿足一定統(tǒng)計(jì)特性的隨機(jī)相位屏(隨機(jī)相位屏的制作方法與湍流隨機(jī)相位屏的方法一致,只是譜密度函數(shù)不同)來(lái)模擬部分相干光束[109],但是純實(shí)數(shù)的隨機(jī)相位屏只能用來(lái)模擬相干性分布是高斯函數(shù)的部分相干光束,研究人員利用這個(gè)類型的隨機(jī)相位屏開(kāi)展了高斯-謝爾模光束、部分相干平頂渦旋光束經(jīng)過(guò)湍流大氣傳輸?shù)南辔黄翑?shù)值模擬計(jì)算[110,111].純相位屏雖然可以在實(shí)驗(yàn)上通過(guò)空間光調(diào)制器方便地產(chǎn)生,但是由于純相位屏的相位關(guān)聯(lián)函數(shù)與所產(chǎn)生的部分相干光束的相干函數(shù)并不一致,導(dǎo)致在數(shù)值計(jì)算上只能模擬關(guān)聯(lián)函數(shù)為高斯函數(shù)的部分相干光束,存在一定的局限性.復(fù)數(shù)屏(光源的振幅和相位同時(shí)進(jìn)行調(diào)制)能很好地克服這一缺點(diǎn)[112],復(fù)數(shù)屏的關(guān)聯(lián)函數(shù)嚴(yán)格地等于部分相干光束的關(guān)聯(lián)函數(shù),因此可以用來(lái)模擬任意相干結(jié)構(gòu)的部分相干光束.假設(shè)一完全相干光束經(jīng)過(guò)復(fù)數(shù)屏調(diào)制后的電場(chǎng)為

其中,E0(r,0)是調(diào)制前的光場(chǎng)的復(fù)振幅;T(r)為具有一定空間關(guān)聯(lián)特性的隨機(jī)復(fù)透射函數(shù)(復(fù)數(shù)屏),假設(shè)T(r)滿足高斯統(tǒng)計(jì),其一階統(tǒng)計(jì)平均為零,二階統(tǒng)計(jì)只依賴于兩點(diǎn)間的距離變化,

?r=r1?r2,二階統(tǒng)計(jì)平均為對(duì)應(yīng)的光源的空間關(guān)聯(lián)(相干)函數(shù).

數(shù)值模擬產(chǎn)生復(fù)屏的方法與產(chǎn)生隨機(jī)相位屏的方法類似,首先我們需要獲知部分相干光場(chǎng)的空間關(guān)聯(lián)函數(shù)的功率譜密度函數(shù),對(duì)空間相干度函數(shù)進(jìn)行傅里葉變化,即

其中,f是空頻域矢量.在傅里葉域中,復(fù)屏又可以表示為功率譜密度函數(shù)的平方根和零均值、單位方差的復(fù)高斯隨機(jī)數(shù)的乘積,即

在湍流理論中,一般假設(shè)光源的振幅、相位變化時(shí)間要比探測(cè)器的積分時(shí)間小很多;探測(cè)器的積分時(shí)間要比湍流引起的相位變化時(shí)間小很多,即τs?τd?τa,也就是說(shuō)探測(cè)器僅僅只測(cè)量湍流引起的光束特征隨時(shí)間的變化.基于這個(gè)關(guān)系,部分相干光束在大氣湍流中的數(shù)值模擬邏輯過(guò)程表述如下:

1)湍流相位屏的合成,然后固定在傳輸路徑上(湍流相位屏的合成上文已提及);

2)通過(guò)K1張復(fù)屏合成部分相干光束,并使每一幀(一張復(fù)屏)電場(chǎng)單獨(dú)通過(guò)由多相位屏模擬的大氣湍流(上文提及)到達(dá)接收面;

3)在接收面接收K1幀經(jīng)過(guò)多相位屏的電場(chǎng),并對(duì)K1幀電場(chǎng)取平均,得到的平均值就是通過(guò)大氣湍流之后的一幀部分相干光束;

4)步驟1—3循環(huán)K2次,得到K2幀部分相干光束通過(guò)大氣湍流的結(jié)果,可以利用這K2幀結(jié)果分析計(jì)算部分相干光束通過(guò)大氣湍流之后的統(tǒng)計(jì)特性.

上述數(shù)值模擬方法也可以推廣用來(lái)處理部分相干矢量光束經(jīng)過(guò)湍流大氣傳輸.

5 總結(jié)與展望

本文回顧了部分相干光在湍流大氣中傳輸行為的研究進(jìn)展、基本傳輸模型理論及相位屏模擬的數(shù)值方法.與完全相干光束相比,部分相干光束在湍流大氣中傳輸能有效抵抗湍流帶來(lái)的光束展寬、光斑漂移、光強(qiáng)閃爍等負(fù)面影響.但是,也應(yīng)該注意到這種抵抗負(fù)面效應(yīng)的效果是犧牲了接收面上的能量接收率,因?yàn)橥葪l件下,部分相干光束的發(fā)散角比完全相干光束要大.惟一例外是非均勻關(guān)聯(lián)光束,在合適的參數(shù)下,它不僅具有較低的光強(qiáng)閃爍因子,而且有較高的能量接收效率,但是如何高效地產(chǎn)生這類非均勻關(guān)聯(lián)光束[113,114],目前還是一個(gè)難題.

近年來(lái),光場(chǎng)調(diào)控受到越來(lái)越多研究人員的關(guān)注,逐漸成為光學(xué)領(lǐng)域的熱點(diǎn)研究課題之一.光場(chǎng)的多維調(diào)控,包括振幅、相位、偏振、相干性等的方法和技術(shù)越來(lái)越成熟,這些特殊調(diào)制光場(chǎng)在湍流大氣中的傳輸行為及與湍流介質(zhì)的相互作用是值得重點(diǎn)研究的一個(gè)課題,特別是利用光束的某些特性變化反演大氣湍流的參數(shù).近幾年,關(guān)于部分相干光束在湍流中的研究文獻(xiàn)大多是關(guān)于理論計(jì)算和數(shù)值模擬研究光束的傳輸行為,相應(yīng)的實(shí)驗(yàn)研究文獻(xiàn)較少.開(kāi)展戶外部分相干光束在湍流大氣中的實(shí)驗(yàn)研究還需要進(jìn)一步加強(qiáng),這有待于部分相干光束的產(chǎn)生和調(diào)控技術(shù)的相對(duì)成熟.