基于同位模式的北京城市職住要素研究

陳思宇 孟斌 李夢天

[摘要]空間同位模式是指在一定空間范圍內(nèi)頻繁發(fā)生并且位置相鄰的事件集合，對于研究不同空間要素之間交互關(guān)系具有重要意義。城市中職住要素的研究是很典型的同位模式。基于Ripley′s K函數(shù)基礎(chǔ)上變形的雙變量K（r）函數(shù)和雙變量g（r）函數(shù)，對獲取到的北京市五環(huán)內(nèi)職住要素和交通要素進(jìn)行空間同位模式分析，最終得出：1）住宅區(qū)和寫字樓之間的同位模式較為復(fù)雜，低價位住宅區(qū)與寫字樓之間不存在同位模式;中等價位住宅區(qū)和高價

位住宅區(qū)與寫字樓存在同位模式，但尺度范圍不同;2）不同價位住

宅區(qū)與地鐵站之間都在0～2 km范圍上具有空間同位模式，并且隨著距離的增加，住宅區(qū)與地鐵站的空間同位模式強(qiáng)度在不斷減弱;3）寫字樓與地鐵站點同樣在中距離范圍內(nèi)具有空間同位模式，并且隨著距離的增加，空間同位模式的強(qiáng)度均在減小。

[關(guān)鍵詞]空間同位模式;Ripley′s K函數(shù);職住關(guān)系

[中圖分類號]TU 984[文獻(xiàn)標(biāo)志碼]A[文章編號]10050310（2021）01002009

Research on the Elements of Beijing′s Urban JobsHousing

Based on the Colocation Model

Chen Siyu， Meng Bin，Li Mengtian

（College of Applied Arts and Science， Beijing Union University， Beijing 100191， China）

Abstract： Spatial colocation mode refers to a collection of events that frequently occur within a certain spatial range and are located adjacent to each other， which is of great significance to the study of the interaction between different spatial elements. The research on the elements of the urban jobshousing is a typical colocation model. Based on the deformed bivariate K（r） function and bivariate g（r） stemming from Ripley′s K function， this paper analyzes the spatial colocation pattern of the acquired elements of jobshousing and traffic within the Fifth Ring Road of Beijing， and finally obtains： 1） The colocation pattern between residential buildings and office buildings is complicated. There is no colocation pattern between lowpriced residential areas and office buildings; mediumpriced residential areas and highpriced residential areas have colocation patterns relationship with office buildings， but the scale ranges are different. 2） There is a spatial colocation pattern between residential buildings of different prices and subway stations in the range of 0～2 km. As the distance increases， the intensity of the spatial colocation pattern between residential buildings and subway stations constantly weakens. 3） Office buildings and subway stations also have spatial colocation patterns in the middledistance range. As the distance increases， the intensity of the spatial colocation patterns decreases.

Keywords： Colocation pattern; Ripley′s Kfunction; Relationship of jobshousing

引言

空間同位模式是指頻繁地出現(xiàn)在同一區(qū)域的一組空間要素或事件[12]。例如，生態(tài)學(xué)中的尼羅鱷與埃及燕鸻、經(jīng)濟(jì)學(xué)中的供應(yīng)商與制造商的分布，均是典型的同位模式[3]?？臻g同位模式的挖掘是空間關(guān)聯(lián)性挖掘的一種，逐漸成為空間挖掘中的一項重要任務(wù)[45]，有助于揭示地理現(xiàn)象或事件的分布規(guī)律?？臻g同位模式挖掘的方法主要可分為空間統(tǒng)計方法和數(shù)據(jù)挖掘方法兩類[6]，其中，空間統(tǒng)計方法是通過關(guān)聯(lián)性指標(biāo)去描述不同類空間要素之間的關(guān)系，這些指標(biāo)包括交叉K函數(shù)、平均最近鄰距離、空間回歸模型、同位系數(shù)、Q 統(tǒng)計等[712]。

居住和就業(yè)是城市活動的兩大基本功能，城市居住與就業(yè)分布及關(guān)系的研究對城市的規(guī)劃和管理有極大的影響[13]。由于改革開放以來我國經(jīng)歷了快速的城市化，城市居民職住分離現(xiàn)象日趨嚴(yán)峻，職住過度分離對城市個體和整體都會產(chǎn)生嚴(yán)重影響[14]。城市職住空間是城市空間結(jié)構(gòu)與形態(tài)特征的重要方面，也是影響城市功能正常運行與居民生活質(zhì)量的重要內(nèi)容[15]。在職住關(guān)系研究中，以距離或通勤時間的研究成果較多，但對居住地和就業(yè)地及交通節(jié)點等職住要素的布局模式關(guān)注相對較少，張景秋等對北京寫字樓的空間集聚模式進(jìn)行了系統(tǒng)研究，得出各行業(yè)總體呈“大分散，小聚集”“大分散、大聚集”和“小分散、大聚集”3種分布模式[16];牛方曲等對住宅的空間布局進(jìn)行了分析，在綜合考慮

政策限制和市場規(guī)律模擬開發(fā)商決策行為后，構(gòu)建了FDM模型，為城市空間模擬分析提供了良好的解決方案[17]。研究居住地與工作地兩種要素在空間上是否同位，能夠有效揭示城市職住要素的分布規(guī)律，極大地優(yōu)化城市職住空間結(jié)構(gòu)，保障城市功能正常運行。利用空間模式研究中的同位模式，可以定量研究城市職住要素的同位模式，推進(jìn)城市職住要素的研究。

在特大型城市中，考慮到城市功能分區(qū)，居住和就業(yè)分區(qū)布局成為一種常態(tài)，一般通過交通節(jié)點將居住要素和就業(yè)要素進(jìn)行連接。本文基于空間同位模式，利用Ripley′s K函數(shù)對北京市五環(huán)內(nèi)的住宅區(qū)、寫字樓以及地鐵站點進(jìn)行分析，為促進(jìn)城市職住平衡、合理配置城市要素提供理論和實證支持。

1研究區(qū)概況與數(shù)據(jù)來源

1.1研究區(qū)概況

北京市五環(huán)內(nèi)區(qū)域是北京城市的中心區(qū)域，包括東城區(qū)、西城區(qū)的全部以及朝陽區(qū)、海淀區(qū)、石景山區(qū)、豐臺區(qū)的大部分區(qū)域。根據(jù)國家統(tǒng)計局在2015年發(fā)布的北京環(huán)路人口分布數(shù)據(jù)顯示，北京市五環(huán)內(nèi)有1 054萬常住人口，占全市人口的48.9%。北京市五環(huán)內(nèi)有CBD、金融街、中關(guān)村、三元橋、望京、總部基地、亞運村等就業(yè)集聚區(qū)，根據(jù)《2018年北京區(qū)域統(tǒng)計年鑒》，僅東城區(qū)、西城區(qū)、朝陽區(qū)和海淀區(qū)的生產(chǎn)總值就為17 746億元，占北京市當(dāng)年生產(chǎn)總值的63.3%。北京市五環(huán)內(nèi)常住人口眾多、地鐵網(wǎng)絡(luò)密集、經(jīng)濟(jì)發(fā)達(dá)、職住問題突出，是研究職住關(guān)系的重點區(qū)域，因此，本文將研究區(qū)選擇為北京市五環(huán)內(nèi)區(qū)域。

本論文以北京市五環(huán)內(nèi)的寫字樓、住宅區(qū)、地鐵站點為研究對象，其中，通過對北京市租房網(wǎng)站使用八爪魚軟件抓取2016年帶有經(jīng)緯度坐標(biāo)信息的寫字樓數(shù)據(jù)和住宅區(qū)數(shù)據(jù)，并結(jié)合2016年北京市五環(huán)內(nèi)地鐵站點數(shù)據(jù)進(jìn)行分析，共獲得北京市五環(huán)內(nèi)住宅區(qū)數(shù)據(jù)4 905條，寫字樓數(shù)據(jù)3 171條，地鐵站點數(shù)據(jù)186條（見圖1、圖2、圖3）。

由于對于住宅區(qū)中高低價位標(biāo)準(zhǔn)有不同劃分，本文采用均價正負(fù)一個標(biāo)準(zhǔn)差的方式將住宅區(qū)數(shù)據(jù)分為高價位住宅區(qū)、中等價位住宅區(qū)及低價位住宅區(qū)3類，其中，高價位住宅區(qū)價格≥59 614元/平方米，共544條數(shù)據(jù);中等價位住宅區(qū)價格為27 816元/平方米～59 613元/平方米，共3 938條數(shù)據(jù);低價位住宅區(qū)價格≤27 815元/平方米，共423條數(shù)據(jù)。

本文采用R 3.5.3軟件對獲得的數(shù)據(jù)進(jìn)行分析和統(tǒng)計。本文中Ripley′s K（r）函數(shù)、g（r）函數(shù)、雙變量K（r）函數(shù)、雙變量g（r）函數(shù)皆采用R軟件中的ads軟件包進(jìn)行分析。

2研究方法

對實際的地理對象點集進(jìn)行集聚模式分析時，僅僅使用最鄰近距離會掩蓋結(jié)果中的其他模式[12]，為了解決這一問題，Ripley對事件之間的所有距離進(jìn)行研究，提出了K函數(shù)方法。與其他單一距離分析工具不同，Ripley′s K函數(shù)最大的優(yōu)勢在于可以研究不同尺度上的空間分布模式。本文使用基于Ripley′s K 函數(shù)變形的雙變量K（r）函數(shù)和雙變量g（r）函數(shù)進(jìn)行研究。

2.1雙變量K（r）函數(shù)分析法

Ripley′s K（r）用于分析不同空間尺度上的聚集程度是否一致，發(fā)現(xiàn)是否存在聚集及聚集的空間尺度。公式為

Kr=An2ni=1nj=1IrUijWiji≠j。（1）

式中，當(dāng)要素i與要素j之間的距離Uij≤r時，Ir為1，否則為0。Wij為權(quán)重，是以i為圓心

、Uij為半徑的圓周長在樣地面積A中的比例;A為面積，n為面積內(nèi)要素的總數(shù);r為空間尺度。

因為Ripley′s K（r）函數(shù)只適用于研究一類要素隨距離變化所呈現(xiàn)出的空間模式，所以Diggle在此基礎(chǔ)上進(jìn)行了改進(jìn)并提出雙變量K函數(shù)[18]，使改進(jìn)后的雙變量K函數(shù)能對兩類要素進(jìn)行分析。其公式為

K12r=An1n2n1i=1n2j=1IrUijWiji≠j。（2）

式中，n1和n2分別代表樣地面積內(nèi)兩類要素的總個數(shù);i和j分別代表第一類要素和第二類要素，而

A、Ir（Uij）、Wij的含義與公式（1）相同。

使用蒙特卡羅（Monte—Carlo）擬合檢驗，計算上下包跡線，即置信區(qū)間。當(dāng)實際的要素分布數(shù)據(jù)計算得到的K12（r）值高于置信區(qū)間時，則代表要素1與要素2之間在此尺度下具有空間同位模式。

2.2雙變量 g（r）函數(shù)分析法

Ripley′s K（r）函數(shù)及雙變量K函數(shù)在多尺度空間分布格局和空間關(guān)聯(lián)性的研究中有著廣泛的應(yīng)用，但其也存在著一些不足之處：Ripley′s K（r）函數(shù)是計算以指定距離為半徑的圓中的所用信息，但隨著距離的不斷增加，大尺度上的計算結(jié)果包含了小尺度的信息，這種累積效應(yīng)使大尺度上的效應(yīng)與小尺度上的效應(yīng)相混淆。

g（r）函數(shù)通過對Ripley′s K（r）函數(shù)的改進(jìn)，從而避免了以上問題。g（r）函數(shù)與Ripley′s K（r）函數(shù)最大的不同在于，Ripley′s K（r）函數(shù)是計算以r為半徑的圓內(nèi)的要素，而g（r）函數(shù)計算以r為半徑，dw為面積的圓環(huán)內(nèi)的要素（圖4和圖5）。g（r）函數(shù)能有效地消除大尺度格局在分析時容易受小尺度上信息影響的缺點，是與Ripley′s K（r）函數(shù)成對使用的重要函數(shù)。同理，雙變量g（r）函數(shù)公式為

g12r=

dK12（r）2πdr。（3）

與K（r）函數(shù)相同，采用蒙特卡羅（Monte—Carlo）模擬檢驗計算置信區(qū)間，以檢驗兩類要素之間是否具有顯著性。若實際的要素分布數(shù)據(jù)計算得到的g12（r）值高于置信區(qū)間，則代表要素1與要素2之間具有空間同位模式。

3結(jié)果與分析

3.1住宅區(qū)與寫字樓空間同位模式分析

基于雙變量K、g函數(shù)計算結(jié)果得到，低價位、中等價位和高價位住宅區(qū)與寫字樓在不同尺度下的同位模式不同：低價位住宅區(qū)與寫字樓之間不存在同位模式;中等價位住宅區(qū)與寫字樓在0～2 km范圍內(nèi)有同位模式;高價位住宅區(qū)與寫字樓在4.8～6.3 km范圍內(nèi)有同位模式。

K、g函數(shù)圖，其中，縱坐標(biāo)是數(shù)據(jù)歸一化后的表達(dá);橫坐標(biāo)則是構(gòu)建的圓的邊界到住宅區(qū)的距離。

如圖6（a）、6（c）、6（e）所示，觀察不同價位住宅區(qū)與寫字樓的K（r）函數(shù)發(fā)現(xiàn)，雖然3種價位住宅區(qū)與寫字樓的實際值在一定距離內(nèi)與高置信值相距較近，并且實際值都有上升趨勢，但總體來看都沒有超過置信區(qū)間，因此都不具有同位模式;如圖6（b）、6（d）、6（f）所示，觀察不同價位住宅區(qū)與寫字樓的g（r）函數(shù)發(fā)現(xiàn)，中等價位住宅區(qū)與寫字樓在0～2 km范圍內(nèi)的實際值高于置信區(qū)間;高價位住宅區(qū)與寫字樓的實際值在4.8～6.3 km范圍內(nèi)高于高置信值，即中高價位住宅區(qū)與寫字樓在不同距離范圍內(nèi)具有空間同位模式。

3.2不同價位住宅區(qū)與地鐵站點空間同位模式分析

基于雙變量K、g函數(shù)計算結(jié)果得到，低價位、中等價位和高價位住宅區(qū)與地鐵站點在不同距離范圍下的同位模式不同：不同價位住宅區(qū)與地鐵站之間都在0～2 km范圍上具有空間同位模式，并且隨著距離的增加，不同價位住宅區(qū)與地鐵站的空間同位模式強(qiáng)度在不斷減少。

如圖7（a）、7（c）、7（e）所示，觀察不同價位住宅區(qū)與地鐵站的K（r）函數(shù)發(fā)現(xiàn)，低價位住宅區(qū)與地鐵站的實際值完全位于置信區(qū)間內(nèi)，因此不具有空間同位模式;中高

價位住宅區(qū)都無限趨近于高置信值，但并不能說明其具有空間同位模式。同樣，

如圖7（b）、7（d）、7（f）所示，觀察不同價位住宅區(qū)與地鐵站的g（r）函數(shù)發(fā)現(xiàn)，不同價位住宅區(qū)與地鐵站都會在小于2 km的某些距離范圍上的實際值高于高置信值，即在這些范圍內(nèi)具有同位模式。

3.3職住要素與地鐵站空間同位模式分析

基于雙變量K、g函數(shù)計算結(jié)果得到，住宅區(qū)與地鐵站和寫字樓與地鐵站都在中距離范圍內(nèi)具有空間同位模式，隨著距離的增加，空間同位模式的強(qiáng)度在減小。同時，寫字樓與地鐵站空間同位模式的強(qiáng)度要大于住宅區(qū)與地鐵站。

如圖8（a）、8（c）所示，觀察住宅區(qū)與地鐵站和寫字樓與地鐵站的K（r）函數(shù)發(fā)現(xiàn)，住宅區(qū)與地鐵站和寫字樓與地鐵站的實際值都只是接近于高置信值，但是并沒有超過，即住宅區(qū)與地鐵站和寫字樓與地鐵站均不具有同位模式。

如圖8（b）、8（d）所示，觀察住宅區(qū)與地鐵站和寫字樓與地鐵站的g（r）函數(shù)發(fā)現(xiàn)，住宅區(qū)與地鐵站和寫字樓與地鐵站的實際值都存在高于高置信值的現(xiàn)象，即具有空間同位模式。但是，對比住宅區(qū)與地鐵站和寫字樓與地鐵站的數(shù)據(jù)發(fā)現(xiàn)，兩者之間集聚的程度有較大的區(qū)別。寫字樓與地鐵站比住宅區(qū)與地鐵站在更短的距離上有空間同位模式，并且寫字樓與地鐵站存在同位模式的距離比住宅區(qū)與地鐵站遠(yuǎn)且持續(xù)，住宅區(qū)與地鐵站雖在0～4 km范圍內(nèi)具有空間同位模式，但也只是零星存在。

4結(jié)論與討論

4.1方法評述

通過對以上結(jié)果分析，可以發(fā)現(xiàn)K（r）函數(shù)與g（r）函數(shù)在進(jìn)行空間同位模式分析時結(jié)果差異明顯，這主要是由于方法本身的計算所導(dǎo)致的。K（r）函數(shù)與g（r）函數(shù)的最大不同之處在于，K（r）函數(shù)是用半徑不斷增加的圓去計算不同空間尺度下的格局，具有一定的累積效應(yīng)，這使得在不同距離內(nèi)的信息會相互干擾。g（r）函數(shù)則是計算給定距離周圍所構(gòu)成的圓環(huán)內(nèi)的密度比例，對其進(jìn)行空間同位模式分析，所以g（r）函數(shù)在3 km范圍內(nèi)反映格局的變化，而在10 km范圍內(nèi)波動較少?？傊?，使用K（r）函數(shù)去描述要素的空間格局全貌會較快捷，且較能反映其整體特征，但K（r）函數(shù)忽略了空間格局及其關(guān)聯(lián)性細(xì)節(jié)的變化特征。所以，若要研究要素空間格局與生態(tài)過程有聯(lián)系的“關(guān)鍵尺度”，即空間格局的局部變化特征時，g（r）函數(shù)分析得到的結(jié)果將更加接近實際情況[22]。在研究空間格局及空間同位模式時，結(jié)合K（r）函數(shù)與g（r）函數(shù)各自的長處，能夠更充分地探索空間格局。

4.2結(jié)論與建議

本文運用在Ripley′s K（r）函數(shù)的基礎(chǔ)上進(jìn)行改進(jìn)和提升的雙變量K（r）函數(shù)、雙變量g（r）函數(shù)，對北京市五環(huán)內(nèi)的職住要素和地鐵站點數(shù)據(jù)進(jìn)行空間同位模式分析，最終得出：

1） 職住要素空間同位模式分析的結(jié)果表明：低價位、中等價位和高價位住宅區(qū)與寫字樓在不同距離范圍內(nèi)的同位模式不同，低

價位住宅區(qū)與寫字樓不存在同位模式，中等價位住宅區(qū)與寫字樓在0～2 km范圍內(nèi)有同位模式，高價位住宅區(qū)與寫字樓在4.8～6.3 km范圍內(nèi)有同位模式，這個結(jié)果說明北京五環(huán)內(nèi)的寫字樓大多分布在中高價位住宅區(qū)附近，中高價位住宅區(qū)與寫字樓之間呈現(xiàn)集聚分布的現(xiàn)象，并且中等價位住宅區(qū)與寫字樓之間的距離比高價位住宅區(qū)更近。

2） 居住地與地鐵站點同位模式分析的結(jié)果表明：不同價位住宅區(qū)與地鐵站之間都在0～2 km范圍內(nèi)具有空間同位模式，并且隨著距離的增加，不同價位住宅區(qū)與地鐵站的空間同位模式強(qiáng)度在不斷減少，表明無論什么價位的居住地，地鐵站點都會是影響其選址的因素之一。

3） 職住及地鐵站點空間同位模式分析的結(jié)果表明： 雖然住宅區(qū)與地鐵站和寫字樓與地鐵站都在中距離范圍內(nèi)具有空間同位模式，并且

隨著距離的增加，空間同位模式的強(qiáng)度均在減小。但是，寫字樓與地鐵站點空間同位模式的強(qiáng)度要大于住宅區(qū)與地鐵站點，這說明北京五環(huán)內(nèi)的寫字樓分布比住宅區(qū)更加靠近地鐵站點，職業(yè)地的選擇比居住地更加傾向于地鐵站點。

綜上所述，北京市五環(huán)內(nèi)職住要素之間的同位模式以及職住要素與地鐵站點之間的同位模式研究，為全面理解職住關(guān)系提供了新的視角，未來北京市在城市規(guī)劃中應(yīng)在注重功能分區(qū)的同時，要有意識

地鼓勵在低價位住宅區(qū)附近配備寫字樓，通過職住要素布局的優(yōu)化配置，有效緩解職住不平衡的現(xiàn)象。

[參考文獻(xiàn)]

[1]LOTWICK H W，SILVERMAN B W. Methods for analyzing spatial processes of several types of points[J]. The Royal Statistical Society，1982（B44）：406-413.