基于复杂网络及其模体的京津冀及周边城市空气污染空间关联与城市协同治理.pdf

为了探究京津冀及周边31个城市空气污染的空间关联与季节演化情况，于2 0 15年1月1日至2 0 2 1年12 月31日选取该区域的空气质量指数（AQI)日均值作为样本，首先对31个城市AQI的时间变化特征和空间分布特征进行了分析；然后计算不同城市AQI的皮尔逊相关系数，结合引力模型，构建了该区域空气污染空间关联网络；最后对网络的整体特征与季节演化情况进行了分析。结果表明，该区域空气污染空间关联网络密度较高，关联比较紧密，度值与中心性较高的城市对空气污染有更高的传输贡献，石家庄、邢台及邯郸的中心性最高，对其他城市空气污染传输的控制能力最强；网络共分为3个凝集子群，石家庄、邢台及邯郸位于子群发生关联的中心位置；对4个季节空间关联网络图进行对比分析，春季网络密度更高，不同城市之间AQI的关联更为紧密；在模体A的四季关联网络中，石家庄、邢台和邯郸为主要传播城市；在模体B的四季关联网络中，邯郸和北京为枢纽城市，这些城市在空气污染的传播过程中起着重要作用。提出，加强石家庄、邢台、邯郸、北京这4个城市与其他城市空气污染协同治理力度，完善监督管理机制，促进环保产业发展；从行政、法律、经济等多方面人手，为解决京津冀大气污染问题提供充分保障；结合京津冀地区的季节气候条件，加强冬季煤改电、煤改气等管理力度，实行区域或集中采暖供热，减少燃料消耗和烟尘排放。关键词：复杂网络；模体；京津冀；空气污染；空间关联；城市协同治理中图分类号：X511;X823Spatial Correlation and Urban Collaborative Governance of Air Pollution in Beijing-Tianjin-Hebei and Nearby Cities Based on Complex Network and Its MotifZHANG Cuizhil-2,AN Haigang.1*(1.Hebei University of Geosciences,Shijiazhuang,Hebei 050031,China;2.Hebei Oriental University,Langfang,Hebei 065000,China)Abstract:In order to explore the spatial correlation and seasonal evolution of air pollution in Beijing-Tianjin-Hebei and itsnearby cities,the daily average AQI of Beijing-Tianjin-Hebei region and its nearby 31 cities from January 1,2015 to December 31,2021 was taken as a sample.Firstly,the regional air quality distribution characteristics are analyzed from the perspective of thetemporal and spatial distribution characteristics of the AQI index of 31 cities;Then,the Pearson correlation coefficient ofAQI in different cities is calculated,and the spatial correlation network of air pollution in this region is constructed based onthe gravity model;Finally,the overall characteristics and seasonal evolution of the network are analyzed.The results show thatthe spatial correlation network density of air pollution in this region is higher and the correlation is closer.The cities withhigher degree and centrality have higher transmission contribution to air pollution.Shijiazhuang,Xingtai and Handan havethe highest centrality and the strongest control ability to other cities on air pollution transmission.The network was divided intothree agglutination subgroups,Shijiazhuang,Xingtai and Handan were located in the center of the correlation of subgroups.A comparative analysis of the spatial correlation network maps of the four seasons shows that the network density is higher inspring,and the correlation between diferent cities is closer.In the four-season correlation network of model A,Shijazhuang,收稿日期：2 0 2 2-0 9-2 9；修订日期：2 0 2 2-11-0 3基金项目：河北省高等学校人文社会科学研究项目（SD191006）；河北省省级科技计划软科学研究专项（2 2 557 6 52 D）作者简介：张翠芝（19 9 5一），女，助教，硕士，主要研究方向为大数据分析。*通讯作者：安海岗E-mail：a n h a i g a n g 16 3.c o m一30 一文献标志码：A文章编号：16 7 4-6 7 32(2 0 2 3)0 4-0 0 30-0 8第15卷第4期张翠芝等基于复杂网络及其模体的京津冀及周边城市空气污染空间关联与城市协同治理2 0 2 3年7 月Xingtai and Handan are the main transmission cities.In the four-season correlation network of model B,Handan and Beijing arethe hub cities,which play an important role in the transmission process of air pollution.It is suggested that we strengthenthe coordination of air pollution control in Shijiazhuang,Xingtai,Handan and Beijing with other cities,improve thesupervision and management mechanism,promote the development of environmental protection industry,solve the Beijing-Tianjin-Hebeiregion air pollution problem from administrative,legal,economic and other aspects so as to provide long-termguarantee,strengthen the management of coal-to-electricity and coal-to-gas in winter by combining with the seasonal climateconditions in the Beijing-Tianjin-Hebei region,and implement regional or centralized heating to reduce fuel consumption and sootemissions.Key words:Complex network;Motif;Beijing-Tianjin-Hebei;Air pollution;Spatial correlation;Urban collaborative governance0引言京津冀地区人口较为密集，经济也比较发达，是目前我国北方地区最具经济活力的一个城市集群 。随着经济的快速发展，京津冀地区的空气污染问题也越来越严峻，引起了社会与政府的广泛关注。空气污染具有空间溢出性,某一地区的空气污染不仅与当地的能源消耗结构、经济结构、居民生活、交通结构等息息相关，而且还会受到附近区域的影响。近年来，国内外学者多运用空间计量经济学、多元统计、莫兰指数（MoransI）及灰色关联模型等 2-3 对我国城市集群的空气污染与空间关联问题进行研究。2 0 14年，国家提出了京津冀协同发展战略。自此以后,国内外学者对京津冀地区的生态环境与协同发展问题进行了大量的研究。一些学者对京津冀地区的能源结构与能源消耗问题进行了研究，并对其与空气污染的关联性进行了分析 4-5;部分学者对京津冀地区空气污染的空间关联进行了统计分析，对其时空演化进行了系统论述，发现该地区不同城市的空气污染具有较强的空间溢出效应 6-10 。在研究方法方面，国内学者也进行了一定的创新。基于复杂网络理论，马宇博等 对京津冀雾霾污染网络进行了统计分析。基于权重综合指数,乔建刚等 12 对京津冀地区不同污染物的重要性进行了排序。基于GIS空间技术和统计学方法，李慧杰等 13 对京津冀地区空气污染的影响因素进行了分析。“21世纪将是复杂性科学的世纪”，这是世界著名物理学家霍金对于复杂性科学的精彩论述。复杂网络是复杂性科学的核心内容，自从复杂网络的无标度特征和小世界性被提出以后 14-1，国内外掀起了研究复杂网络的热潮。当前，复杂网络的研究成果多集中于自然科学 16 、物理科学 17 、管理科学 18 、社会经济学 19-2 0 等领域。在生态环境领域，虽然有部分学者运用复杂网络理论对城市空气质量空间关联与演化进行了相应的研究，但是从模体 2 1视角出发对城市空气污染空间关联及其协同治理进行的探讨还鲜有涉及。为了探求京津冀及周边城市空气污染的空间关联、季节演化情况及模体视角下的城市环境协同治理，本研究于2 0 15年1月1日至2 0 2 1年12 月31日，选取京津冀及周边共31个城市的空气质量指数（AQI）日均值作为样本，通过计算不同城市AQI的皮尔逊相关系数，结合引力模型，构建京津冀及周边城市的空气污染关联网络，并对网络的整体特征与季节演化情况进行了分析；构建了春、夏、秋、冬4个季节的模体关联网络图，对模体A及模体B随季节演化情况进行了分析，并基于模体演化情况，对京津冀及周边城市的环境协同治理提出了相应的对策与建议。1数据来源与研究方法1.1数据来源1.1.1AQI 日均值2015年1月1日至2 0 2 1年12 月3 1日京津冀及周边31个城市的AQI日均值从中国空气质量在线检测平台（https:/wwwa q i s t u d y c n/historydata/）获取。31个城市包括北京、石家庄、天津、唐山、沧州、廊坊、邯郸、邢台、秦皇岛、保定、衡水、郑州、濮阳、鹤壁、安阳、焦作、开封、新乡、晋城、朔州、大同、晋中、忻州、济南、德州、长治、承德、太原、阳泉、张家口、聊城。1.1.2城市空间距离2个城市之间的空间直线距离利用Python软件通过经纬度计算获取。1.2研究方法1.2.1空气质量空间关联网络的建立为了使用复杂网络可视化软件Gephi构建一一31一第15卷第4期张翠芝等基于复杂网络及其模体的京津冀及周边城市空气污染空间关联与城市协同治理2 0 2 3年7 月个较为健壮的空气质量空间关联复杂网络，首先应用SPSS软件计算出31个城市中任意2 个城市的皮尔逊相关系数，构建了一个3131的相关矩阵。参考马宇博 2 2 研究中的使用方法，以所有相关系数的平均值作为阈值，如果2 个城市之间的相关系数大于阈值，则2 个城市之间存在1条边，如果2个城市之间的相关系数小于值,则边不存在。对于2 个城市之间的边的权重，除了考虑相关系数外,还应当考虑空间距离。权重由于受到量纲的影响，差异较大,需要进行归一化处理。1.2.2关联网络节点的度与中心性度是度量网络中节点重要性的最重要的指标。在无向网络中，度表示与这个节点直接相连的节点的数目。在构建的京津冀及周边城市空气质量空间关联网络中，如果1个节点度值越高，代表与这个城市节点空气质量存在空间关联的城市数量越多。节点的度的计算公式见式（1）。k=Zfermjjev式中：k;节点名称为i的度值；i和j任意节点，对应文中表示城市；V一一节点集合；m；节点i和j之间是否存在连边，若存在连边,则m；记为1,否则记为0。度量节点重要性的指标除了度以外，还有中介中心性与接近中心性。中介中心性 2 3 是度量网络中的某个节点的中心性的重要指标，代表该节点在多大程度上位于网络图中其他节点对的中间位置。如果1个节点处于多对节点对之间，那么该节点的度值一般比较低，但这个节点可能在网络连接中起到重要的“中介作用”。中介中心性的计算公式见式（2）。d,(i)BE;=Zd.r式中：BE,一一节点i的中介中心性；d.（i）一节点s到r的最短路径中经过i节点的数量；dsr一一节点s到r的最短路径数量。接近中心性 2 4 也是度量节点重要性的重要指标，一般用来度量网络中的某个节点到其他节点的最短路径的平均长度。即在复杂网络中,如果1个节点离其他节点距离越近，那么它的接近中心性就越高。接近中心性的计算公式见式（3）。n-1CL;=一3 2 一式中：CL；一一节点i的接近中心性；n节点数量；d一一节点i与i之间的最短距离。1.2.3季节关联网络及其模体演化分析为了清晰地展示京津冀及周边城市空气污染空间关联随季节演化的情况，构建了31个城市的空气质量空间关联四季网络图，并对网络图的度、平均度、网络密度及关键城市关联情况进行分析。为了揭示4个季节关联网络的构建模式，引入了“模体”的概念，使用fanmod软件对四季关联网络中反复出现的局部结构进行统计分析，,探究了季节关联网络模体的存在情况及其演化特征。模体是1种特殊的网络结构,通常由少数节点组成,是在复杂网络中出现多次的具有同样结构的网络子图。复杂网络分为有向图与无向图，在有向图中，如果模体由3个节点构成，共包含13种网络结构；如果模体由4个节点构成，共包含19 9 种网络结构。无向图中的模体结构则较为简单，由3个(1)节点构成的模体只包含2 种网络结构，见图1。图1无向图中3个节点的模体示意在复杂网络中，模体的频率是指该模体在复杂网络中出现的次数，计算公式见式（4）。n(V)f(V)=N式中：(V）一模体V在网络中出现的次数，即频率;n(V)一一由n个节点组成的模体V在网络中存在的次数;N1个网络中所有的由n个(2)节点组成的模体出现的总次数。2结果与分析2.131个城市空气质量分布特征2.1.131个城市AQI时间变化特征2015一2 0 2 1年京津冀及周边31个城市AQI季度均值时序变化见图2。由图2 可见，AQI总体上呈下降趋势，说明各个城市整体环境质量有所改善。从季节上看，冬季的AQI均值最大，随着气温上升，AQI均值逐渐呈现下降趋势，并在夏季或秋(3)季达到最低值。之后随着温度下降，AQI值又逐渐呈现上升趋势，并在冬季达到最高值。不同季节的(4)第15卷第4期张翠芝等基于复杂网络及其模体的京津冀及周边城市空气污染空间关联与城市协同治理2 0 2 3年7 月AQI均值对应的空气质量等级大部分集中在良（A Q I 值为50 10 0），说明2 0 15一2 0 2 1年该区域大气环境质量得到显著改善。呈现这种变化规律的原因可能是该区域冬季气温低，空气流通速度慢，不利于大气污染物扩散，而夏季气温高，降雨量大，秋季冷暖空气对流速度快，有利于污染物的沉降和扩散 2 160140105.989410080F60LRAR年份/季节图2 2 0 152 0 2 1年京津冀及周边31个城市AQI季度均值时序变化2.1.231个城市AQI空间分布特征2015一2 0 2 1年京津冀及周边31个城市AQI年均值变化见图3。由图3可见，空气质量整体呈现出北部区域较好，南部区域较差的现象。将31个城市按照AQI年均值分为3类：以承德、张家口、大同、秦皇岛为代表的北部地区，AQI年均值 10 0,对应空气质量等级为轻度污染，空气质量相对较差。这可能是因为北部地区重工业较少，污染物排放量较低，且北部地形较高,处于上风向，有利于大气污染物的扩散；而中部和南部地区重工业较多，污染物排放量较高，且中南部地形较低，处于下风向，不利于大气中污染物的扩散 2 51402021888980德家同京山 坊州津定州州图32 0 15一2 0 2 1年京津冀及周边31个城市AQI年均值变化2.2相关系数对城市AQI数据进行正态分布检验，发现总体符合正态分布特征，因此选用皮尔逊相关系数进行分析。选取31个城市的AQI日均值进行二元分析,得到城市之间的相关系数。2 0 15一2 0 2 1年京津冀及周边31个城市皮尔逊相关系数热力图见图4（横、纵坐标的城市排名按照纬度自北向南依次排列）。由图4可见，右下角区域颜色相对较深,表明这些城市之间的相关性较大。这一现象与各城市之间地域距离存在较大关联。以衡水、德州、邢台、邯郸等为代表的南部区域的相关性强于北部区域，说明空气污染比较严重的城市之间的关联性较强，溢出效应也较高，这一现象与区域间的空气流通和大气污染物的扩散密不可分。202020192018,201720162015廊朔大保阳衡德邢邮济聊长宏鹤濮普新售开郑原泉中泉州谷郸南城治阳壁阳城乡作封州城市承德张家口大同秦皇岛北京唐山廊坊朔州天津保定折州沧州石家庄太原阳泉晋中衡水德州邢台邯郸济南聊城长治安阳鹤壁濮阳晋城新乡焦作开封郑州德知洞岛?汕防州津淀洲洲驻源腺种冰州跆榔南城治湘幢湘城乡作街洲承家大皇北康廊天保怀沧家太阳晋衡德邢邯济斯长安鹤渊副新焦形郑图42 0 152 0 2 1年京津冀及周边31个城市皮尔逊相关系数热力图一3 3 一相关系数1.00.90.80.70.60.50.40.3第15卷第4期张翠芝等基于复杂网络及其模体的京津冀及周边城市空气污染空间关联与城市协同治理2 0 2 3年7 月2.3城市空气质量空间关联网络特征分析表12 0 15一2 0 2 1年京津冀及周边31个城市度值排序以31个城市为节点，城市之间的关联关系为地区邢台石家庄邯郸保定德州衡水聊城鹤壁长治安阳焦作度值25边，以城市相关系数与对应城市间的平均距离的比地区亲新乡沧州济南开封郑州濮阳廊坊晋中唐山天津晋城值为边的权重，以31个城市间的皮尔逊相关系数度值1616161515151414131313的均值为阈值（经计算均值为0.547 9）。由此，构地区阳泉忻州北京太原承德大同秦皇岛朔州张家口度值131212966建出了包含31个节点与2 31条边的空气污染空间关联网络。2.3.1空气污染空间关联网络特征与度分布2015一2 0 2 1年京津冀及周边31个城市空气污染空间关联社会网络示意见图5。此网络平均度值为14.9 0 3，网络密度为0.49 7，表明京津冀及周边城市空气污染空间关联网较为紧密。由图5可见，存在一些空气质量空间关联较为紧密的城市对,如邯郸与安阳、廊坊与北京、焦作与新乡、太原与晋中等，针对这些城市，解决空气污染问题时应尽量协同治理，制定合作共赢的方针政策。此外，在空间关联网络中,邢台与石家庄度值最高，表明这2 个城市与网络中多个城市空气质量存在空间关联，空气污染具有较强的空间溢出效应。新乡唐山北京晋城晋中保定张家只大同焦作邢台天津星德州朔州衡水濮阳沧州阳泉秦皇岛太原图52 0 152 0 2 1年京津冀及周边31个城市空气污染空间关联社会网络示意2015一2 0 2 1年京津冀及周边31个城市度值排序见表1。由表1可见,邢台、石家庄、邯郸、保定、衡水为度值最高的5座城市，这些城市度值均2 0,表明与其他城市空气质量关联较为紧密，在空气污染方面有较强的传播力和影响力。节点度值最低的城市是张家口，与其他城市相比，张家口全年空气质量较好，与其他城市的空间关联也较弱。一34一25242020201919 19 1918662.3.2中心性分析2015一2 0 2 1年京津冀及周边31个城市的中介中心性与接近中心性见图6（a)（b）。由图6 可见，中介中心性排名前5位的城市分别是邢台、石家庄、北京、邯郸及保定，这些城市的媒介作用较强,对其他城市有较大的溢出效应。相比于中介中心性,接近中心性在几何地理上更接近中心位置。接近中心性排名前3位的城市分别是邢台、石家庄及邯郸,这些城市距离其他城市的路径之和最短，一旦这些城市的空气受到严重污染，周边较大范围地区的空气质量也将在一定程度上受到负面影响。由上述研究可知,邢台、石家庄及邯郸这3个城市的中介中心性与接近中心性排名均靠前，说明它们处于网络中的核心位置，空气污染传播的能力最强。(a)中介中心性302550跆注京郸淀中州防冰州城城南津山源岛洞壁治阳德作泉州乡州封阳洲邢家北邯保晋忻廊衡德聊晋济天唐太皇大鹤长安承焦家阳沧新郑开濮朔浙州济南郑州4石城市(b)接近中心性0.90.80.70.60.50.40.30.20.10台庄郸淀阳冰州治阳城作坊中南乡州津山城州京州泉封原州洞岛德口邢家邯保濮衡德长安聊焦鹤廊晋济新忻天唐晋沧北郑阳开太朔大皇承家石秦张城市图6 2 0 152 0 2 1年京津冀及周边31个城市的中介中心性与接近中心性基于Gephi软件对网络图进行模块化处理，2015一2 0 2 1年京津冀及周边31个城市模块化与子群第15卷第4期张翠芝等基于复杂网络及其模体的京津冀及周边城市空气污染空间关联与城市协同治理2 0 2 3年7 月展示示意见图7。由图7 可见，京津冀及周边城市空为了探究季节对城市空气污染空间关联的影气污染空间关联网络共分为3个凝集子群，石家庄、响,对原始数据按照气温划分季节，春季为3、4、邢台及邯郸处于3个子群发生关联的中心位置，再次5月份，夏季为6、7、8 月份，秋季为9、10、11月份，验证了上述结论。冬季为12、1、2 月份。通过对4个季节的空间网络图的图密度和度值进行统计计算，可知春季图密度濮阳济南南德州开封聊城新乡衡水郑州安阳焦作修石家庄保定陈坊北余长治晋城图7 2 0 152 0 2 1年京津冀及周边31个城市模块化与子群展示示意2.4季节的空气污染空间关联网络与模体演化(a)模体A春季济南焦作邮郸聊城安阳郑州沧州晋中天津保定唐山长治衡水北京德州邢台新乡晋城承德阳泉秦皇岛大同(e)模体B春季郑州沧州焦作秦皇岛邮邦普中析州太原济南长治聊城廊坊唐山承德张家口鹤壁阳泉朔州麻坊君家连安阳石家庄普城邢台濮阳新乡大同天津保定朔州安阳北京薇水图8 2 0 15一2 0 2 1年京津冀及周边31个城市4个季节模体演化网络图由图8（a）一（d)可见，4个季节均以石家庄、邢台和邯郸为主要传播城市，其中春季模体关联密度明显高于其他3个季节，该网络的平均度值为13.2,平均加权度值为56.9 33。由图8（e）一沧州天津秦皇岛康山承德太源张家口大同朔州(b)模体A夏季焦作北京新乡大同晋城太原承德濮阳晋中邮郸州太原忻州衡水廊坊()模体B夏季新乡长治晋中濮阳济南唐山秦皇岛保定衡水德州沧州郑州太原天津城承德北京邢台焦作大同张家口忻州鹤为0.531，平均度值为15.9 35；夏季图密度为0.484,平均度值为14.516；秋季图密度为0.47 1，平均度值为14.12 9；冬季图密度为0.47 7，平均度值为14.32 3。可见4个季节的空气污染空间关联网络图差异不大，密度与平均度值均较高，不同城市之间的AQI关联较为紧密，春季网络密度最高，这可能与北方春季多大风天气有关。通过对4个季节的关联网络图的模体进行统计分析,发现每个季节中均存在2 个模体，模体A为3个节点两两关联的模体，模体B为3个节点中2个节点存在关联的模体。边代表2 个节点城市间存在关联关系，权重表示2 个节点间存在关联程度的大小。2 0 15一2 0 2 1年京津冀及周边31个城市4个季节模体演化网络图见图8（a）一（h）。由图8 可见，从整体来看，模体A网络图节点间的关联程度比模体B的更紧密。(c)模体A秋季(d)模体A冬季秦皇岛大同张家口新乡承德唐山邯郸阳泉唐山郑州沧州秦皇岛天津大同晋中北京太原新乡析州朔州阳泉邢台沧州唐山翔州热作淡阳安阳街水济南聊城保定康坊石家庄衡水济南天津北京德州晋中秦皇岛承德(g)模体B秋季(h)模体B冬季济南长治张家口郑州廊坊普城天津集作沧州新乡承德晋中阳菜石家店邮郸保定德州阳泉安阳普城天洋洗阳张家口北京濮阳秦皇岛大同太原朔州邢台聊城安阳析州衡水越生（h）可见，4个季节主要以邯郸和北京为枢纽城市，表明邯郸和北京在31个城市间空气污染的传播过程中起着重要作用，在治理中需要重点防控。春季模体的关联程度略高于其他3个季节，不同的网络一35一濮阳焦作晋城新乡焦作长治朔州保定衡水唐山沧州石家庄邢台聊城济南晋.中北京大同德州麻坊郑州桥州太原鹤壁秦皇岛承德参考文献第15卷第4期张翠芝等基于复杂网络及其模体的京津冀及周边城市空气污染空间关联与城市协同治理2 0 2 3年7 月中同一类模体的频率不同，模体频率越高，说明该的治理。地方政府须完善监督管理机制,发展环保模体在真实网络中出现的次数越多，该模体也越重产业，走可持续发展道路。要。4个季节模体B的出现频率为34%38%，其（2）依据当地不同的季节气候条件，制定不同中夏季和冬季较高，均 37%，说明夏、冬两季城市的空气污染防治政策。冬季采用煤改电、煤改气、之间的空气污染散播与关联程度高于春、秋两季。地热能等多种清洁取暖方式，实行区域或集中采暖供热，充分提高锅炉或天然气利用率，降低燃料消3结论与建议耗和烟尘排放量。3.1结论（3)石家庄、邢台、邯郸、北京在模体中处于中基于复杂网络，从京津冀及周边31个城市心位置，应与其他城市开展空气污染的协同治理。AQI的时间变化特征和空间分布特征的角度,构建从行政、法律、经济等多方面人手，为解决京津冀大了2 0 15一2 0 2 1年该区域的空气质量空间关联加权气污染问题提供充分保障。网络，分析了空气污染空间关联与空间溢出的基本状况；构建了春、夏、秋、冬4个季节的空气污染空间关联网络，探究了不同季节对不同城市空气污染的空间关联差异性，并对季节网络进行了模体演化分析。得出以下结论：（1)通过对京津冀及周边31个城市AQI的时间变化特征和空间分布特征进行分析可知,该区域环境质量总体趋于好转，且区域内北部地区空气质量较好，南部地区空气质量较差。（2）通过对该区域空气质量空间网络的度值、中介中心性、接近中心性及凝集子群的比较分析可知,邢台、石家庄、邯郸、保定及衡水在空气污染方面有较强的传播力和影响力，其中邢台、石家庄及邯郸的中心性相对较高,处于网络的中心位置。网络共包括3个凝集子群，石家庄与邯郸处于子群发生关联的核心位置。（3)通过对该区域4个季节空间关联网络图进行统计分析可知，4个季节的网络图密度与度值均较高,不同城市之间的AQI关联较为紧密。在不同季节,城市之间关联的紧密程度随着季节演变发生变化，春季网络图密度最高。（4)通过对该区域季节关联网络模体演化特征进行分析可知，在模体A的关联网络中，石家庄、邢台、邯郸为主要传播城市；在模体B的关联网络中，邯郸和北京为枢纽城市，邯郸和北京在31个城市间空气污染的传播过程中起着重要作用。3.2建议对京津冀及周边城市空气污染的治理，应在该区域制定协同治理对策，建议如下：（1)石家庄、邢台及邯郸处于3个凝集子群中心位置，污染传播能力最强，应加强对这3个城市一3 6 一1 范凤岩，王洪飞，樊礼军.京津冀地区空气污染的健康经济损失评估 J.生态经济，2 0 19,35（9）：157-16 3.2杨新平，徐登国.基于中小城市空气污染指标的时间序列模型构建 J.环境污染与防治,2 0 14,36（11）：111.3刘晓红，江可申，中国城市PM2.5的时空分异及影响因素分析一基于16 1个城市的实证研究 J.调研世界,2 0 18(1)：25-34.4刘晓红，江可申.我国雾霾污染影响因素的空间效应一基于PM2.5的实证分析J.科技管理研究，2 0 17,37（12）：247-252.5孟雪，狄乾斌，季建文.京津冀城市群生态绩效水平测度及影响因素 J.经济地理，2 0 2 0,40（1）：18 1-18 6,2 2 5.6买超平，连素琴，牛叔文，等.甘肃省能源消费结构与空气质量变化的相关性分析J.兰州大学学报（自然科学版），2014,50(4):508-512.7许燕婷，刘兴诏，王振波.基于AQI指数的中国城市空气质量时空分布特征 J.广西师范大学学报（自然科学版）,2 0 19，37(1):187-196.8程钰，刘婷婷，赵云璐，等.京津冀及周边地区“2+2 6”城市空气质量时空演变与经济社会驱动机理 J.经济地理，2 0 19，39(10):183-192.9彭功，赵媛媛，赵吉麟，等.京津冀大气污染传输通道区大气污染时空格局研究J.中国环境科学，2 0 19，39（2）：449-458.10刘华军，裴延峰.经济发展与中国城市雾霾污染一一基于空间关联网络情形下的考察J.城市与环境研究，2 0 18（3）：15-41.【11马宇博，高广阔.基于节点重要性评价的京津冀雾霾污染网络研究 J.环境科学学报，2 0 18,38（6）：2 2 8 7-2 2 9 6.12 乔建刚，李士宣，王英.京津冀空气质量状况研究 J.应用化工,2 0 19,48(11):2 7 9 1-2 7 9 3.13李慧杰，王秀兰，王计平，等.2 0 13一2 0 17 年间京津冀地区空气质量及影响因素分析 J.环境监测管理与技术,2 0 19,31(2):21-25.第15卷第4期张翠芝等基于复杂网络及其模体的京津冀及周边城市空气污染空间关联与城市协同治理2 0 2 3年7 月14 WATTS D J,STROGATZ S H.Collective dynamics of small-world networksJ.Nature,1998,393(6684):440-442.15 CAO Y J,WANG G Z,JIANG Q Y,et al.A neighbourhood evol-ving network modeJ.Physics Letters A,2006,349(6):462466.16陈刚，王信，肖伸平，等.复杂动态网络指数采样同步控制 J：中南大学学报（自然科学版），2 0 18，49（12）：2993-3001.17孔江涛，黄健，龚建兴，等.基于复杂网络动力学模型的无向加权网络节点重要性评估J.物理学报，2 0 18,6 7（9）：255-271.18】秦胜君.基于复杂网络的企业创新网络级联失效可靠性模型 J.科技管理研究，2 0 17,37（7）：19 9-2 0 4.19 樊斌锋，杨琼.社交复杂网络中基于无向图和聚类的社区检测关键技术研究（英文）J.机床与液压，2 0 19，47（12）：179-184.【2 0 游鸽，郭昊，刘向.复杂网络视角下的金融市场结构演化与风险传染 J.金融发展研究，2 0 2 0（1）：30-39.21】徐越，刘雪明.基于三元闭包模体的关键节点识别方法J.复杂系统与复杂性科学，2 0 2 2,11(2）：1-8.22马宇博.京津冀区域城市雾霾污染动态网络分析 D.上海：上海理工大学，2 0 19.23 FREEMAN L C.A set of measures of centrality based on be-tweennessJ.Sociometry,1977,40(1):35-41.24 ALEX B.Communication paterns in task-oriented groupsJ.Journal of the Acoustical Society of America,1950,22(6):725730.25梁银双，刘黎明，卢媛.基于函数型数据聚类的京津冀空气污染特征分析 J.调研世界，2 0 17（5）：43-48.征订启事欢迎订阅2 0 2 3年中国环境监测中国环境监测创刊于19 8 5年，是由生态环境部主管、中国环境监测总站主办的环保科技学术期刊，是全国中文核心期刊、中国科技核心期刊，被中国核心期刊（遵选）数据库、中国期刊全文数据库收录期刊，是中国学术期刊综合评价数据库统计源、中国科学引文数据库来源期刊，国际标准连续出版物编号为ISSN1002-6002、国内统一刊号为CN11-2861/X。期刊一贯坚持学术性、专业性与实用性相结合的办刊原则，宣传国家有关环境保护工作的方针、政策、法规，介绍国内外先进的环境管理理念和环境监测技术,交流环境监测科研成果。设有特约来稿、特别关注、监测管理、分析测试、调查评价、应急预警、自动监测、环境遥感、仪器设备等栏目。所载文章具有重要的指导作用，是从事环保、监测工作的管理干部和科技人员的必备工具书。期刊为大16 开，双月刊，国内外公开发行。全年共出版6 期，每双月月末出版，每期定价7 0 元，总价42 0 元。订阅方法：直接向编辑部汇款订阅单位名称：中国环境监测总站开户银行：中国工商银行北京市和平里支行银行账号：0 2 0 0 0 0 42 0 9 0 8 9 114334行号：10 2 10 0 0 0 0 42 3编辑部联系方式：地址：北京市朝阳区安外大羊坊8 号院乙中国环境监测编辑部邮编：10 0 0 12联系人：霍晓芹订阅电话：[phone]手机：15510 36 2 0 0 05邮箱：8 18 0 59 8 4 行微信：138 0 12 7 0 6 37中国环境监测编辑部一3 7 一