；優(yōu)選的，可W選取1個月或2個月內(nèi)污染物的后向軌跡 進(jìn)行計(jì)算，本實(shí)施例中取1個月內(nèi)污染物的后向軌跡執(zhí)行步驟S101至步驟S104,其中沒隔 1小時為污染物繪制一條后向軌跡，從而1個月內(nèi)共有720條污染物的后向軌跡，進(jìn)而可W 得到720條后向軌跡對應(yīng)的720個總輸送強(qiáng)度；
[0093] 步驟S106,取總輸送強(qiáng)度最大的前N個后向軌跡經(jīng)過的網(wǎng)格區(qū)域作為預(yù)測區(qū)域， 根據(jù)預(yù)測區(qū)域W及源清單確定污染源的地理位置；其中N為正整數(shù)，如取N= 10,即取總輸 送強(qiáng)度最大的10個后向軌跡所經(jīng)過的網(wǎng)格區(qū)域作為預(yù)測區(qū)域，然后結(jié)合源清單查找預(yù)測 區(qū)域內(nèi)污染物排放強(qiáng)度最大且對應(yīng)有污染源（如化工廠、排污廠等）的地理位置，即為污染 源的地理位置，如預(yù)測區(qū)域內(nèi)存在的某化工廠對應(yīng)的位置等。
[0094] 與現(xiàn)有技術(shù)相比，本發(fā)明基于后向軌跡的污染源追蹤方法通過后向軌跡計(jì)算污染 物在每一網(wǎng)格內(nèi)的輸送概率、距離權(quán)重W及時間權(quán)重，并通過源清單計(jì)算相應(yīng)網(wǎng)格內(nèi)污染 物的排放強(qiáng)度，最后計(jì)算得到每個網(wǎng)格內(nèi)污染物的輸送強(qiáng)度，進(jìn)而得到后向軌跡上污染物 的總輸送強(qiáng)度，然后通過計(jì)算多條后向軌跡的總輸送強(qiáng)度W及將總輸送強(qiáng)度最大的N個后 向軌跡經(jīng)過的網(wǎng)格區(qū)域作為預(yù)測區(qū)域，再將預(yù)測區(qū)域結(jié)合源清單確定污染源的地理位置， 實(shí)現(xiàn)了污染源的定位追蹤，且結(jié)果較為精確，為解決空氣污染問題提供了重要依據(jù)。
[0095] 需要說明的是，本實(shí)施例中對多條后向軌跡進(jìn)行計(jì)算來確定預(yù)測區(qū)域，相較于依 靠單條后向軌跡來確定污染源的預(yù)測區(qū)域，能夠有效解決一條后向軌跡線由于測量誤差、 內(nèi)插誤差、計(jì)算誤差等因素而導(dǎo)致測量結(jié)果存在較大誤差的問題。
[0096] 下面結(jié)合圖7說明求解污染物后向軌跡的方法，其包括W下步驟：
[0097]S201;從匪5模式生成的NetCDF文件獲取污染物在初始點(diǎn)P的初始速度V(P，t); [009引其中獲取初始速度V(P，t)的方式有多種，由于NetCDF文件中的數(shù)據(jù)網(wǎng)格大小為 3KM*3KM，因此在對初始點(diǎn)P進(jìn)行速度合成時，可W簡單地取初始點(diǎn)P周圍的四個網(wǎng)格點(diǎn)處 風(fēng)速的速度平均值作為初始速度，該結(jié)果非常接近污染物的實(shí)際運(yùn)行速度，同時，將風(fēng)速的 矢量平均值作為初始速度，運(yùn)算時間極小，使得所有程序在運(yùn)行時1秒內(nèi)運(yùn)算結(jié)束，大大縮 短了運(yùn)算時間。如圖8所示，步驟S201包括：
[0099]S2011;從NetCDF文件中讀取初始點(diǎn)P周圍的四個網(wǎng)格點(diǎn)處的風(fēng)速，如圖9所示， 初始點(diǎn)P周圍的四個網(wǎng)格點(diǎn)處的風(fēng)速分別記為V(P1)、V(P2)、V(P3)、V(P4);
[0100] 52012;計(jì)算風(fēng)速￥任1)、￥任2)、￥任3)、￥任4)的矢量平均值^作為初始速度￥任， t)，其中風(fēng)速包括U、V、W方向上的S個分量，U代表東、西方向風(fēng)速，V代表南、北方向風(fēng)速， W代表垂直地面風(fēng)速，從而計(jì)算得到的初始速度V(P，t)也包括S個方向上的分量，計(jì)算得 到的初始速度V(p，t)如圖9所示。
[0101] 在另一優(yōu)選實(shí)施例中，如圖10所示，為了使獲取的初始速度V(P，t)準(zhǔn)確度更高， 還可W通過反距離加權(quán)法計(jì)算初始速度V(p，t)。此時，步驟S201具體包括：
[0102]S2011';從NetCDF文件中讀取初始點(diǎn)P周圍四個網(wǎng)格點(diǎn)處的風(fēng)速；如圖11 所示，四個網(wǎng)格點(diǎn)對應(yīng)的坐標(biāo)分別為P1 (xl,yl,zl)、P2 (x2,y2,z2)、P3 (x3,y3,z3)、 P4 (x4,y4,z4)，初始點(diǎn)P處的坐標(biāo)為P(X，y，z);
[0103]S2012';分別計(jì)算四個網(wǎng)格點(diǎn)到初始點(diǎn)P的距離，如圖11所示，計(jì)算公式如下：
[0104]
【主權(quán)項(xiàng)】
1. 一種基于后向軌跡的污染源追蹤方法，其特征在于，包括： 51 :根據(jù)污染物的后向軌跡依次計(jì)算所述后向軌跡經(jīng)過的每一網(wǎng)格內(nèi)所述污染物的輸 送概率、距離權(quán)重以及時間權(quán)重； 52 :根據(jù)所述污染物的源清單計(jì)算所述網(wǎng)格內(nèi)所述污染物的排放強(qiáng)度； 53 :根據(jù)所述輸送概率、距離權(quán)重、時間權(quán)重以及排放強(qiáng)度計(jì)算所述網(wǎng)格內(nèi)所述污染物 的輸送強(qiáng)度，計(jì)算公式為： T_1 (i,j) =R_1 (i,j)*E(i,j)*ff_d(i,j)*ff_tl(i,j) 其中，T_l(i，j)為網(wǎng)格（i，j)內(nèi)所述污染物的輸送強(qiáng)度，R_l(i，j)為網(wǎng)格（i，j)內(nèi) 所述污染物的輸送概率，E(i，j)為網(wǎng)格（i，j)內(nèi)所述污染物的排放強(qiáng)度，W_d(i，j)為網(wǎng)格 (i，j)內(nèi)所述污染物的距離權(quán)重，W_tl(i，j)為網(wǎng)格（i，j)內(nèi)所述污染物的時間權(quán)重； S4:根據(jù)各所述網(wǎng)格內(nèi)的所述輸送強(qiáng)度計(jì)算所述后向軌跡上所述污染物的總輸送強(qiáng) 度； 55 :取所述污染物在不同時刻的所述后向軌跡重復(fù)步驟S1至S4以得到多條所述后向 軌跡對應(yīng)的多個所述總輸送強(qiáng)度；以及 56 :取所述總輸送強(qiáng)度最大的前N個所述后向軌跡經(jīng)過的網(wǎng)格區(qū)域作為預(yù)測區(qū)域，根 據(jù)所述預(yù)測區(qū)域以及所述源清單確定污染源的地理位置，其中N為正整數(shù)。
2. 如權(quán)利要求1所述的基于后向軌跡的污染源追蹤方法，其特征在于，步驟S5具體 為： 取1個月內(nèi)的720條所述后向軌跡重復(fù)步驟S1至S4以得到所述后向軌跡對應(yīng)的720 個所述總輸送強(qiáng)度，其中1個月內(nèi)每隔1小時為所述污染物繪制一條所述后向軌跡。
3. 如權(quán)利要求1所述的基于后向軌跡的污染源追蹤方法，其特征在于，根據(jù)污染物的 后向軌跡計(jì)算所述網(wǎng)格內(nèi)所述污染物的輸送概率具體為： R_(i,j) =T _1 (i,j)/n 其中T_l(i，j)為所述后向軌跡在所述網(wǎng)格（i，j)內(nèi)的停留時間，n為所述后向軌跡 的運(yùn)行總時間。
4. 如權(quán)利要求3所述的基于后向軌跡的污染源追蹤方法，其特征在于，所述停留時間 以所述后向軌跡在所述網(wǎng)格（i，j)內(nèi)的軌跡點(diǎn)個數(shù)表示，所述運(yùn)行總時間以所述后向軌跡 上的所有軌跡點(diǎn)個數(shù)表示。
5. 如權(quán)利要求1所述的基于后向軌跡的污染源追蹤方法，其特征在于，根據(jù)污染物的 后向軌跡計(jì)算所述網(wǎng)格內(nèi)所述污染物的距離權(quán)重具體為： ff_d(i,j) = 1/(d(i,j)/5+l) 其中d(i，j)為所述網(wǎng)格（i，j)與所述后向軌跡上所述污染物的目的地之間的距離。
6. 如權(quán)利要求1所述的基于后向軌跡的污染源追蹤方法，其特征在于，根據(jù)污染物的 后向軌跡計(jì)算所述網(wǎng)格內(nèi)所述污染物的時間權(quán)重具體為： ff_tl(i,j) =l/((t_l(i,j))/18+l) 其中t_l(i，j)為所述污染物由所述網(wǎng)格（i，j)移動至所述后向軌跡上所述污染物的 目的地所需的時間。
7. 如權(quán)利要求1所述的基于后向軌跡的污染源追蹤方法，其特征在于，根據(jù)所述污染 物的源清單計(jì)算所述網(wǎng)格內(nèi)所述污染物的排放強(qiáng)度具體為： 從所述源清單中讀取所述網(wǎng)格以及3個相鄰網(wǎng)格內(nèi)所述污染物在不同高度層處的排 放強(qiáng)度； 依次計(jì)算每一所述網(wǎng)格內(nèi)所述污染物的總排放強(qiáng)度，所述總排放強(qiáng)度為各個高度層相 應(yīng)的所述排放強(qiáng)度的總和； 將所述網(wǎng)格以及相鄰網(wǎng)格的所述總排放強(qiáng)度取平均值以作為所述網(wǎng)格內(nèi)所述污染物 的排放強(qiáng)度。
8. 如權(quán)利要求1至7任一項(xiàng)所述的基于后向軌跡的污染源追蹤方法，其特征在于，所述 后向軌跡的計(jì)算方法包括： (1) 從MM5模式生成的NetCDF文件中獲取所述污染物在初始點(diǎn)的初始速度； (2) 根據(jù)所述初始點(diǎn)的位置坐標(biāo)和所述初始速度計(jì)算預(yù)測點(diǎn)的位置坐標(biāo)，計(jì)算公式為： P，（t+At) =P(t)+V(P，t) *At; (3) 根據(jù)所述初始點(diǎn)和所述預(yù)測點(diǎn)的位置坐標(biāo)對所述預(yù)測點(diǎn)進(jìn)行修正，修正公式為： P(t+At) =P⑴ +0? 5[V(P，t)+V(P'，t+At) ] *At; (4) 將所述修正點(diǎn)作為所述初始點(diǎn)重復(fù)步驟（1)至步驟（3)以獲取所述所述污染物的 多個位置坐標(biāo)進(jìn)而得到所述所述污染物的所述后向軌跡； 其中，P'（t+At)為所述預(yù)測點(diǎn)的位置坐標(biāo)，P(t)為所述初始點(diǎn)的位置坐標(biāo)，V(P，t)為 所述初始速度，At為所述預(yù)測點(diǎn)與所述初始點(diǎn)的時間差，P(t+At)為對所述預(yù)測點(diǎn)進(jìn)行 修正后得到的修正點(diǎn)的位置坐標(biāo)，V(P'，t+A