基于車載巡測的區(qū)域電磁環(huán)境水平與人口密度相關性分析

摘 要：對北京市2 個行政區(qū)開展了30 MHz～6 GHz 射頻車載巡測，獲取133 余萬個點位的電場強度數據記錄，通過地理位置信息歸類到2 個行政區(qū)各自的22 個街鄉(xiāng)鎮(zhèn)中并計算射頻總值均值，均值范圍分別為0. 43 ～ 1. 67V/ m 和0. 56～1. 51 V/ m;與各街鄉(xiāng)鎮(zhèn)的人口密度、基站密度的相關系數均高于0. 8，為強正相關。對北京市1 個商務區(qū)開展了一天24 小時30 MHz～6 GHz 的車載巡測，獲取24 組共14. 5 萬余個點位的射頻總值和基站總值數據，測值范圍分別為0. 70～1. 46 V/ m 和0. 67～1. 31 V/ m;對每小時截取的百度人流量熱力圖進行解析計算了小時人口密度值，范圍為6. 58～46. 54 人/ hm2 ，發(fā)現小時人口密度值與小時射頻總值和基站總值相關系數均高達0. 901。可見區(qū)域電磁環(huán)境水平與人口密度空間分布特點、區(qū)域電磁環(huán)境水平值24 小時動態(tài)波動與人流量聚集密度變化趨勢強相關。

關鍵詞：車載巡測;區(qū)域電磁環(huán)境;人口密度;基站密度;定量分析

中圖分類號：X837 文獻標志碼：A

近年來隨著車載巡測技術在電磁監(jiān)測領域的廣泛應用，國內北京、南通、成都及西班牙的卡塞雷斯、日本東京等城市均開展了基于車載巡測的局部區(qū)域電磁環(huán)境水平的監(jiān)測與評價，繪制了電磁地圖[1-6] ;有的對小范圍的居民區(qū)、商業(yè)區(qū)等功能區(qū)開展了24 小時不同時段電磁環(huán)境水平的巡測[4] ;電磁地圖展示的區(qū)域電磁環(huán)境水平在空間上呈不均勻分布，24 小時巡測電磁環(huán)境水平在時間上亦有波動。

國內外均有對區(qū)域電磁環(huán)境水平高低原因的研究分析，如段臨林等人[7] 在2014—2015 年對廈門市區(qū)開展射頻電場強度監(jiān)測，發(fā)現中心城區(qū)由于人口密集、商貿活動頻繁等而相應密集布設移動基站， 導致測值高于旅游景區(qū)等其他區(qū)域;Mazloum 等人[8] 在2016 年和2017 年分別對荷蘭的阿姆斯特丹和法國的安納西開展射頻電場強度監(jiān)測，認為前者由于人口密度大，相應的基站布設密度大等原因使得其電場強度測值高于后者;宋欣蔚等人[9] 于2019 年對北京市3 個典型商業(yè)區(qū)開展射頻電磁環(huán)境監(jiān)測發(fā)現，電磁環(huán)境水平高的區(qū)域移動基站安裝密集且有較大的人流量密度等。上述研究均認為區(qū)域電磁環(huán)境水平及時間的變化與人口密度、基站密度、人流量密度等相關。

本文利用車載巡測技術，一是對北京市兩個行政區(qū)（行政區(qū)A 和行政區(qū)B）共44 個街鄉(xiāng)鎮(zhèn)開展射頻車載巡測，對44 個街鄉(xiāng)鎮(zhèn)的人口密度、基站密度及射頻電磁水平值進行關聯性分析;二是選取北京市1 個商務區(qū)開展24 小時射頻電磁環(huán)境水平的車載巡測，每小時巡測中截取百度人流量熱力圖，計算每小時的人口密度值，分析每小時電磁環(huán)境水平與人口密度的關聯性。以上針對區(qū)域電磁環(huán)境水平空間上的變化規(guī)律及一天中不同時段變化規(guī)律的定量分析，對說清電磁環(huán)境空域時域變化原因、通過人口密度判斷區(qū)域電磁環(huán)境水平的空間分布特點、制定電磁環(huán)境的分區(qū)管理原則以及科學回應公眾疑問等都有重要意義。

1 研究基礎

2023 年和2024 年利用車載巡測方式陸續(xù)開展了北京市2 個行政區(qū)共44 個街鄉(xiāng)鎮(zhèn)的區(qū)域電磁環(huán)境監(jiān)測及1 個商務區(qū)24 小時不同時段的電磁環(huán)境巡測。研究中采用的部分車載巡測系統(tǒng)、巡測車輛、監(jiān)測因子和巡測質量保證與文獻[4]一致。巡測、評估方法及電磁地圖繪制依據《射頻電磁輻射車載巡測技術規(guī)范》 （ DB11/ T 2017—2022） [10]和《區(qū)域電磁環(huán)境調查與評估方法（ 試行）》（HJ 1349—2024） [11] 。

1. 1 車載巡測系統(tǒng)

2 個行政區(qū)巡測采用的設備是DZER-100[4] ，配備一個30 MHz ～ 6 GHz 的射頻選頻探頭，2023年6 月經中國計量科學研究院校準符合電性能指標要求，有效期1 年。

商務區(qū)巡測設備是OS - 6S， 配備一個30 MHz～6 GHz 的射頻選頻探頭，電場強度測量范圍為1 mV/ m～300 V/ m;巡測采樣間隔為250 ms;巡測數據為電場強度（V/ m）。監(jiān)測設備于2024年2 月經中國計量科學研究院校準符合電性能指標要求，有效期1 年。

1. 2 巡測條件

1. 2. 1 2 個行政區(qū)

2023 年6 月—7 月和9 月—10 月工作日晝間9：00—18：00，基站等主要電磁輻射源流量高峰時段開展巡測。巡測沿區(qū)域內次干線以上的道路行駛，平均車速35 km / h;巡測中監(jiān)測數據和行駛軌跡在配套軟件系統(tǒng)中實時顯示并記錄，數據量分別為73 余萬條和60 萬條。

1. 2. 2 1 個環(huán)境功能區(qū)

2024 年9 月5 日—6 日對選定的商務區(qū)開展24 小時巡測（1 小時巡測1 次），數據間隔300 ms，并分別對原始數據進行小時數據統(tǒng)計。商務區(qū)內具備條件的所有道路行駛，平均車速25 km / h;每小時數據量為6 000 余條。

1. 3 巡測數據記錄

數據記錄包括監(jiān)測設備型號，探頭型號，監(jiān)測時間，每個監(jiān)測點位的經緯度（WGS84 地理坐標系）、30 MHz ～ 6 GHz 的綜合電場強度總值（以下簡稱射頻總值）、基站總值，4G、5G、移動、聯通、電信、廣播電視、調頻廣播和WIFI 等的選頻電場強度積分值等。

1. 4 行政區(qū)基本信息

街鄉(xiāng)鎮(zhèn)政府是我國基層行政機關，將行政區(qū)車載巡測的電場強度以街鄉(xiāng)鎮(zhèn)為單元進行分解，可為街鄉(xiāng)鎮(zhèn)提供監(jiān)測數據，有效支持屬地及網格精細化管理，推動城市—區(qū)—街鄉(xiāng)鎮(zhèn)級別的屬地監(jiān)管。

開展車載巡測的行政區(qū)A 和行政區(qū)B 常住人口數量分別約為180 萬和220 萬，面積分別約900km2 和1 300 km2 。行政區(qū)A 轄區(qū)內22 個街鄉(xiāng)鎮(zhèn)（其中11 個街道、10 個鎮(zhèn)、1 個鄉(xiāng)），行政區(qū)B 轄區(qū)內22 個街鄉(xiāng)鎮(zhèn)（其中8 個街道、14 個鎮(zhèn)）。兩個行政區(qū)均為北京市的近郊區(qū)，各街鄉(xiāng)人口密度范圍分別為0. 04 ～ 2. 55 萬人/ km2 和0. 006 ～ 3. 04 萬人/ km2 ;各街鄉(xiāng)基站密度范圍分別為1. 9 ～ 26. 8個/ km2 和0. 31 ～ 29. 4 個/ km2 [12] 。由于各街道面積相對較小、人口密集、商業(yè)娛樂活動多，而鄉(xiāng)、鎮(zhèn)相反，因此各街鄉(xiāng)鎮(zhèn)的人口密度相差較大。

為進行數據歸類統(tǒng)計分析，將人口密度按大小劃分為三級，即一級： gt; 1 萬人/ km2 ;二級：0. 5萬～1 萬人/ km2 ;三級：lt;0. 5 萬人/ km2 。劃分后，行政區(qū)的22 個街鄉(xiāng)鎮(zhèn)按人口密度大小分別排序并編號，即A1 ～ A22 和B1 ～ B22 人口密度依次減少;按人口密度分級歸類后，行政區(qū)A 人口密度一級區(qū)包括A1 ～ A8， 二級： A9 ～ A12， 三級： A13～A22;行政區(qū)B 人口密度一級區(qū)：B1 ～ B6，二級：B7～ B10，三級：B11～ B22。

2 結果分析

2. 1 行政區(qū)電磁環(huán)境水平與人口密度關系分析

2. 1. 1 街鄉(xiāng)鎮(zhèn)區(qū)域電磁環(huán)境水平空間分布特點

對2 個行政區(qū)巡測的所有原始數據通過各點位經緯度信息與所屬街鄉(xiāng)鎮(zhèn)關聯后，分別對各街鄉(xiāng)鎮(zhèn)所有巡測點位的射頻總值進行平均后作為各街鄉(xiāng)鎮(zhèn)的射頻總值。經統(tǒng)計，行政區(qū)A 和行政區(qū)B 各街鄉(xiāng)鎮(zhèn)的射頻總值范圍分別為0. 43～1. 67 V/m 和0. 56～1. 51 V/ m，44 個街鄉(xiāng)鎮(zhèn)的射頻總值均明顯低于《電磁環(huán)境控制限值》 （ GB 8702—2014） [13] 中 30 MHz～6 GHz 頻段中最低電場強度限值12 V/ m。分別對兩個行政區(qū)各街鄉(xiāng)鎮(zhèn)人口密度、射頻總值繪制了電磁地圖，如圖1 和圖2所示。

由圖1 和圖2 可以看出，各街鄉(xiāng)鎮(zhèn)射頻總值的空間分布與人口密度分布特點基本一致，人口密集的街道區(qū)域電磁輻射水平較高，而人口稀少、商業(yè)活動不繁華的鎮(zhèn)、鄉(xiāng)等電磁環(huán)境水平較低。行政區(qū)A 人口主要集中在西北區(qū)域，主要是一級人口密度街道（圖1（a）），相應的電磁環(huán)境水平較高的街道也分布在西北側（圖1（b））;同樣，行政區(qū)B 人口密度一級區(qū)域主要為中心城區(qū)街道和行政區(qū)南側（圖2（a）），這些街道的電磁環(huán)境水平明顯高于其他街鄉(xiāng)鎮(zhèn)（圖2（b））。可見，從空間分布特點能夠凸顯出電磁環(huán)境水平與人口密度密切相關。

2. 1. 2 街鄉(xiāng)鎮(zhèn)區(qū)域電磁環(huán)境水平與人口密度關系

為進一步了解區(qū)域電磁環(huán)境水平與人口密度的關聯程度，對二者進行了量化分析。

圖3 為兩個行政區(qū)各街鄉(xiāng)鎮(zhèn)射頻總值與人口密度的關系對比圖，可以看出隨著街鄉(xiāng)鎮(zhèn)人口密度由高到低，射頻總值總體處于下降趨勢。行政區(qū)A 人口密度一級、二級和三級的街鄉(xiāng)鎮(zhèn)平均射頻總值分別為1. 44 V/ m、1. 06 V/ m 和0. 71 V/ m;行政區(qū)B 分別為1. 36 V/ m、1. 05 V/ m 和0. 72 V/m，與人口密度關聯性非常明顯。圖3 中部分街鄉(xiāng)鎮(zhèn)的射頻總值與人口密度的相關性較弱，如圖3（a）中街道A3、A10、A14 等的電場強度值較高，但人口密度較低，主要是由于基站架設較低等原因影響了測值的規(guī)律性。

徐輝等人[4] 在2023 年對北京市800 km2 區(qū)域開展的電磁環(huán)境車載巡測發(fā)現，基站電場強度在射頻總值中平均占比高達76. 7%，因此基站是環(huán)境中數量最多、貢獻最大的電磁輻射源。圖4 為各街鄉(xiāng)鎮(zhèn)的射頻總值與基站密度的對比，可以看出隨著街鄉(xiāng)鎮(zhèn)人口密度的減少，基站密度與射頻總值總體都呈下降趨勢。圖4（a）中，街道A17 基站密度高而射頻總值相對較低，結合圖3（a） A17人口密度也較低，主要原因是該街道主要為政府集中辦公區(qū)域，常住人口較少，但晝間工作人員較多，因此盡管人口密度低但基站架設相對較多;而街道A22 基站密度較低但射頻總值較高，主要是基站架設低或巡測點位處無樓房/ 樹木遮擋。

對行政區(qū)A 和B 各街鄉(xiāng)鎮(zhèn)的射頻總值與人口密度、基站密度分別進行了99% 置信區(qū)間皮爾遜相關性分析，結果顯示射頻總值與人口密度相關系數分別為0. 843 和0. 928，射頻總值與基站密度的相關系數分別為0. 834 和0. 906，均為強正相關，且人口密度相關性均略高于基站密度。由于基站數量、架設位置等基礎資料較難獲取，而街鄉(xiāng)鎮(zhèn)的人口密度為公開信息，因此，獲取一個行政區(qū)各街鄉(xiāng)鎮(zhèn)的人口密度，就可以了解該區(qū)域射頻電磁環(huán)境水平近似的空間分布特點。

2. 2 環(huán)境功能區(qū)電磁環(huán)境水平與人口密度的動態(tài)關系分析

2. 2. 1 電磁環(huán)境水平24 小時巡測

大尺度區(qū)域電磁環(huán)境水平在空間上具有不均勻性，與人口密度密切相關。聚焦小區(qū)域的電磁環(huán)境水平，一天24 小時中會隨時間變化而波動，徐輝等人[4] 2023 年對一個工業(yè)區(qū)進行24 小時電磁環(huán)境車載巡測后發(fā)現，電磁輻射水平時間分布特征主要與人員在一天中工作、生活的流動性規(guī)律相關。為進行電磁輻射水平與人員流動性相關的量化分析，選取一個一天中人員聚集變化較大的商務區(qū)于2024 年9 月開展了24 小時30 MHz～6 GHz 的射頻電磁輻射的車載巡測，每小時巡測路徑保持一致，巡測區(qū)域面積分別為1. 1 km2 ，數據量約6 000 條記錄/ 小時，24 小時共獲取145 000條原始數據，對原始監(jiān)測電場強度值按小時進行統(tǒng)計后獲得24 個小時均值記錄，每個記錄中包括射頻總值、基站總值等。圖5 為商務區(qū)1 個小時（16：00—17：00）的巡測電磁地圖。

對每小時巡測原始數據中的射頻總值、基站總值分別計算平均值作為小時射頻總值均值和基站總值均值，24 個小時均值范圍分別為0. 70 ～1. 46 V/ m 和0. 67 ～ 1. 31 V/ m，均明顯低于《電磁環(huán)境控制限值》 （GB 8702—2014） [13] 中30 MHz～6 GHz 頻段中最低電場強度限值12 V/ m。

2. 2. 2 基于百度熱力圖的人口密度解析

由于通信基站規(guī)劃的重點是覆蓋規(guī)劃和容量規(guī)劃，其中容量規(guī)劃與用戶數量即人口密度直接相關[1] ，且從2. 1. 2 節(jié)結論可知環(huán)境中電磁環(huán)境水平主要與人口密度相關，商務區(qū)24 小時射頻總值的波動也應與每小時人群聚集密度，即人口密度的變化相關。本文利用百度實時人流量熱力圖解析出了每小時的人口密度，并與小時電場強度值進行關聯性分析。

百度熱力圖是基于用戶訪問百度產品（如地圖、搜索、天氣、音樂等） 時所在的位置信息聚類，反映一個時間點的相對人口聚集情況[14-15] 。每小時巡測中獲取一次實時百度熱力圖，共24 張圖，圖6 為兩張例圖， 時間分別為19： 00 和凌晨00：00。然后在Matlab 中用Imread 將熱力圖讀取轉化為像素RGB 的函數，提取每個像素RGB 值，識別對應像素的人口密度后，計算所選區(qū)域內的人口密度。得到24 小時人口密度變化范圍為6. 58～46. 54 人/ hm2 ，人流活動最少和最多的時間點分別為凌晨4：00 和中午12：00。

2. 2. 3 電磁輻射水平與人口密度動態(tài)變化相關性分析

將24 小時射頻總值、基站總值與人口密度值進行對比分析，如圖7 所示。由圖7 可以看出：

（1）小時人口密度與射頻總值、基站總值變化趨勢基本一致。人群聚集密度在一天中呈現低-高-低的規(guī)律性，與商務區(qū)工作時間吻合;而射頻總值和基站總值水平同樣也是低-高-低的變化趨勢。

（ 2）7：00 開始人流量增加，8：00 達到峰值，同樣電磁輻射水平也從7：00 開始上漲8：00 達到高值;中午12：00 工作人員外出就餐等因此人流量最高，而對應的電場強度為上午最高值;18：00 為下班時間，人流聚集量呈較高水平達到46. 54 人/hm2 ，此時電話聯絡、刷手機訪問基站量增多，因此電場強度達到一天中最大值1. 46 V/ m;之后夜間、凌晨由于基本沒有了辦公人員，人流量和電場強度都處于低水平狀態(tài)，凌晨3：00 和4：00 人口密度（6. 81 人/ hm2 和6. 58 人/ hm2 ） 和電場強度（0. 71 V/ m 和0. 73 V/ m）均為最低水平。

將24 小時人口密度與每個小時的射頻總值、基站總值進行99% 置信區(qū)間皮爾遜相關性分析，相關系數均為0. 901，呈強正相關。24 小時中區(qū)域電磁輻射水平與人群聚集密度的變化趨勢一致，區(qū)域射頻電場強度、基站電場強度等電磁輻射水平在時間上的高低變化規(guī)律與人口流動性變化規(guī)律密切相關。晝間強相關性的趨勢高于夜間，主要是商務區(qū)22：00 至凌晨2：00 人流量聚集程度明顯下降，人員活動少，雖然基站業(yè)務量下降，但基站的公共信道依然工作，而且二者之間并非線性關系，因此射頻總值下降幅度并沒有人口密度下降的幅度大，但下降的趨勢依然一致。

3 結論

本文開展了2 個行政區(qū)的區(qū)域電磁環(huán)境水平與人口密度相關性分析及1 個商務區(qū)24 小時電磁環(huán)境水平與基于百度人流量熱力圖的小時人口密度相關性分析，發(fā)現：

（1）行政區(qū)A 和B 各自22 個街鄉(xiāng)鎮(zhèn)，射頻總值范圍分別為0. 43 ～ 1. 67 V/ m 和0. 56 ～ 1. 51 V/m;商務區(qū)24 小時車載巡測的每小時射頻總值、基站總值范圍分別為0. 70 ～ 1. 46 V/ m 和0. 67～1. 31 V/ m，均明顯低于《電磁環(huán)境控制限值》（GB8702—2014） [13] 中30 MHz～ 6 GHz 頻段中最低電場強度限值12 V/ m。

（2）按人口密度大小將各街鄉(xiāng)鎮(zhèn)劃分為三級，行政區(qū)A、B 人口密度一級、二級和三級的街鄉(xiāng)鎮(zhèn)平均射頻總值分別為1. 44 V/ m、1. 06 V/ m、0. 71V/ m 和1. 36 V/ m、1. 05 V/ m、0. 72 V/ m，隨著人口密度級別的降低，街鄉(xiāng)鎮(zhèn)的區(qū)域電場強度值逐漸減小，相關系數分別為0. 843 和0. 928;電場強度值與基站密度相關系數分別為0. 834 和0. 906，均為強相關。

（3）商務區(qū)24 小時射頻總值、基站總值與人口密度相關系數均為0. 901，呈強正相關。

綜上，一個城市各街鄉(xiāng)鎮(zhèn)的區(qū)域電磁環(huán)境水平及空間分布特點與人口密度密切相關，因此可以通過一個區(qū)域的人口密度分布情況近似掌握電磁環(huán)境水平的空間分布，并為電磁污染防治立法中電磁環(huán)境分區(qū)管控提供依據;區(qū)域24 小時電磁環(huán)境水平的波動特點與人流量聚集密度的變化趨勢一致的規(guī)律，可以為日后通過百度人流量熱力圖實現區(qū)域電磁環(huán)境水平動態(tài)更新提供解決方案。

參考文獻：

[ 1 ] 武攀峰， 王國旗， 陸煒， 等. 基于車載監(jiān)測法的Ⅱ型大城市射頻公眾曝露探討[J]. 環(huán)境監(jiān)測管理與技術， 2017，29（5）： 60-63.

[ 2 ] 唐輝， 高鵬， 粟琨璞， 等. 車載式選頻測量（30 MHz ～ 6 GHz） 在城市電磁環(huán)境水平監(jiān)測的應用[J]. 輻射防護，2023， 43（4）： 311-317.

TANG Hui， GAO Peng， SU Kunpu， et al. Application of vehicle-mounted frequency-selection measurement in urbanelectromagnetic environment monitoring[J]. Radiation Protection， 2023， 43（4）： 311-317.

[ 3 ] 宋欣蔚， 韓孟琦， 陳依婕， 等. 北京城區(qū)道路電磁環(huán)境監(jiān)測及與人口熱力關聯分析[J]. 中國環(huán)境監(jiān)測， 2023， 39（6）： 204-209.

SONG Xinwei， HAN Mengqi， CHEN Yijie， et al. Monitoring of the electromagnetic environment at urban roads in Beijingand analysis of its correlation with population heat[J]. Environmental Monitoring in China， 2023， 39（6）： 204-209.

[ 4 ] 徐輝， 佟晶， 李飛， 等. 基于車載巡測的北京市局部區(qū)域電磁環(huán)境解析[J]. 輻射防護， 2024， 44（4）： 356-366.

XU Hui， TONG Jing， LI Fei， et al. Preliminary analysis on area ambient electromagnetic fields in part of Beijing byvehicle survey[J]. Radiation Protection， 2024， 44（4）： 356-366.

[ 5 ] Teruo O， Kaoru E， Kazuhiro T， et al. Large-Area monitoring of radiofrequency electromagnetic field exposure levels frommobile phone base stations and broadcast transmission towers by Car-Mounted measurements around Tokyo [ J ].Electronics， 2023， 12（8）： 1835.

[ 6 ] Paniagua-Snchez J M， García-Cobos F J， Rufo-Prez M， et al. Large-area mobile measurement of outdoor exposure toRadio frequencies[J]. Science of the Total Environment， 2023， 877： 162852.

[ 7 ] 段臨林， 楊傳俊， 唐超， 等. 廈門市電磁環(huán)境解析[J]. 中國環(huán)境監(jiān)測， 2018， 34（2）： 122-129.

DUAN Linlin， YANG Chuanjun， TANG Chao， et al. Preliminary studies on the electromagnetic environment of Xiamencity[J]. Environmental Monitoring in China， 2018， 34（2）： 122-129.

[ 8 ] Mazloum T， Aerts S， Joseph W， et al. RF-EMF exposure induced by mobile phones operating in LTE small cells in two different urban cities[J]. Annals of Telecommunications， 2019， 74（1）： 35-42.

[ 9 ] 宋欣蔚， 岳云濤， 朱心月， 等. 北京市中心商區(qū)射頻電磁環(huán)境測試與分析[J]. 中國環(huán)境監(jiān)測， 2021， 37（4）： 135-142.

SONG Xinwei， YUE Yuntao， ZHU Xinyue， et al. Measurement and analysis of the radio-frequency electromagneticenvironment in Beijing’s central business districts[J]. Environmental Monitoring in China， 2021， 37（4）： 135-142.

[10] 北京市核與輻射安全中心. 射頻電磁輻射車載巡測技術規(guī)范： DB11/ T 2017—2022[S]. 北京：北京市市場監(jiān)督管理局， 2022.

[11] 四川省輻射環(huán)境管理監(jiān)測中心站， 北京市核與輻射安全中心. 區(qū)域電磁環(huán)境調查與評估方法（試行）： HJ 1349—2024[S]. 北京：中國環(huán)境科學出版社， 2024.

[12] 紅黑人口庫[DB/ OL]. [2023-05-29]. http：/ / www. hongheiku. com/ sjrk/ 1088. html.

[13] 環(huán)境保護部輻射環(huán)境監(jiān)測技術中心. 電磁環(huán)境控制限值： GB 8702—2014[S]. 北京：中國環(huán)境科學出版社， 2014.

[14] 譚欣， 黃大全， 趙星爍， 等. 基于百度熱力圖的職住平衡度量研究[J]. 北京師范大學學報（自然科學版）， 2016，52（5）： 622-627.

[15] 郭翰， 郭永沛， 崔娜娜. 基于多元數據的北京市六環(huán)路內晝夜人口流動與人口聚集區(qū)研究[J]. 城市發(fā)展研究，2018， 25（12）： 107-121.