致密砂巖孔喉結(jié)構(gòu)分析與滲透率預(yù)測方法
——以川中地區(qū)侏羅系沙溪廟組為例

陳少云,楊勇強(qiáng),2,邱隆偉,2,王小娟,楊保良,葉熱杰普·哈布臘什木

1.中國石油大學(xué)(華東) 地球科學(xué)與技術(shù)學(xué)院,山東 青島 266580;2.深層油氣全國重點實驗室,山東 青島 266580;3.中國石油 西南油氣田公司 勘探開發(fā)研究院,成都 610041

隨著全球?qū)茉葱枨蟮牟粩嘣黾?非常規(guī)油氣已成為當(dāng)今世界油氣勘探開發(fā)的重要領(lǐng)域[1-7]。致密氣是一種分布廣泛、勘探潛力巨大的非常規(guī)油氣資源[8-11]。中國致密砂巖氣廣泛分布于四川、鄂爾多斯、塔里木、準(zhǔn)噶爾等多個盆地,有利勘探面積達(dá)32.46×104km2,地質(zhì)資源量為21.85×1012m3,是中國現(xiàn)階段十分重要的非常規(guī)天然氣勘探領(lǐng)域[11-13]。近年來,四川盆地中部地區(qū)中侏羅統(tǒng)沙溪廟組陸相致密砂巖取得了重大勘探進(jìn)展,累計提交探明儲量達(dá)千億立方米[14-16]。川中地區(qū)沙溪廟組屬于典型的“孔隙型”致密砂巖,孔隙類型多樣、孔喉結(jié)構(gòu)復(fù)雜,儲層質(zhì)量差異受孔喉結(jié)構(gòu)制約明顯,強(qiáng)烈的非均質(zhì)性嚴(yán)重制約著油氣開發(fā)工作的進(jìn)行[17-20]。

孔喉結(jié)構(gòu)包括孔喉的幾何形狀、大小、連通性以及分布特征等幾個方面[21]。表征儲層孔喉結(jié)構(gòu)的常用方法有:鑄體薄片、掃描電鏡、高壓壓汞、恒速壓汞、核磁共振、微米CT、氮氣吸附、二氧化碳吸附、圖像分析等[21-25]。近年來,分形理論表征孔喉結(jié)構(gòu)在非常規(guī)油氣研究中受到越來越多的關(guān)注[26-30]。前人利用分形維數(shù)(Df)在研究儲層孔喉非均質(zhì)強(qiáng)度、儲層分類、儲層評價等方面進(jìn)行了大量工作,但在利用分形維數(shù)表征滲透率方面研究較淺[27,31-32]。筆者從儲集空間類型、儲集物性研究出發(fā),結(jié)合高壓壓汞實驗和分形理論,重點分析了川中沙溪廟組致密砂巖孔喉結(jié)構(gòu)特征,探討了孔喉結(jié)構(gòu)、分形維數(shù)、儲層物性之間的關(guān)系,并進(jìn)一步分析孔喉尺寸對滲透率的貢獻(xiàn),建立了適用于研究區(qū)的滲透率預(yù)測模型。

1 地質(zhì)背景

四川盆地是在揚(yáng)子克拉通地臺基礎(chǔ)上發(fā)展起來的一個大型疊合盆地,現(xiàn)今已經(jīng)成為中國西部重要的含油、富氣的大型盆地之一[33-34]。四川盆地侏羅系自下而上發(fā)育自流井組、涼高山組、沙溪廟組、遂寧組和蓬萊鎮(zhèn)組5套地層[35-37]。研究區(qū)位于四川盆地中部,研究層位為中侏羅統(tǒng)沙溪廟組(圖1)。

圖1 四川盆地構(gòu)造分區(qū)與侏羅系巖性特征

沙溪廟組以中間一套分布廣泛的黑色“葉肢介頁巖”為分界線分為沙一段和沙二段[16]。沙一段發(fā)育2～4組灰色砂巖與灰色、紅色泥巖的韻律層,底部發(fā)育一套灰色厚層的砂巖與下伏涼高山組相區(qū)分;沙二段主要由灰色的中厚砂巖與紫紅色泥巖組成不等厚互層。地層沉積環(huán)境為淺水三角洲,三角洲前緣水下分流河道砂體分布廣泛,單期河道砂體分布面積超4 000 km2,河道水流方向主要呈北東—南西向[14]。

2 實驗過程及研究方法

2.1 樣品來源和實驗方法

本次研究共收集川中地區(qū)沙溪廟組475塊致密砂巖樣品、2 421份物性測試數(shù)據(jù),樣品與數(shù)據(jù)均來源于中國石油西南油氣田勘探開發(fā)研究院。除收集到資料外,本次研究所有實驗測試均在中國石油大學(xué)(華東)深層油氣重點實驗室完成。

2.2 孔隙微觀形貌觀察

利用偏光顯微鏡對109塊藍(lán)色鑄體薄片的儲層空間類型與孔喉配置樣式進(jìn)行觀察;通過CoxemEM-30掃描電鏡對34塊砂巖樣品進(jìn)行孔隙微觀形貌觀察,并測量孔隙大小。

2.3 高壓壓汞實驗

參考石油工業(yè)標(biāo)準(zhǔn)《巖石毛管壓力曲線的測定:SY/T 5346—2005》,利用AutoPore Ⅳ 9500壓汞儀對33塊不同砂巖樣品進(jìn)行高壓壓汞實驗,取樣位置見圖1。在測試之前將樣品磨制成直徑2.5 cm、長度2.5 cm的兩端平整的圓柱體,并將樣品清洗干凈、烘干至恒重,測量樣品的幾何形態(tài)、孔隙度和滲透率。本次實驗最大壓力約為200 MPa,得到的主要實驗參數(shù)有排驅(qū)壓力(Pcd)、最大進(jìn)汞飽和度(Smax)、剩余飽和度(Sr)、最大孔喉半徑(Rp)、平均孔喉半徑(Ra)、孔喉分選系數(shù)(Sp)、滲透率貢獻(xiàn)度等。

2.4 計算分形維數(shù)

在高壓壓汞實驗分析的基礎(chǔ)上,采用Li Kewen Model[38]分形模型探究致密砂巖孔喉分形特征,具體實驗方法與計算過程見參考文獻(xiàn)[39]。分形維數(shù)的計算結(jié)果詳見表1。

表1 川中侏羅系沙溪廟組致密砂巖高壓壓汞實驗數(shù)據(jù)與分形維數(shù)

3 結(jié)果

3.1 儲層物性與儲集空間類型

研究區(qū)儲層以長石巖屑砂巖、巖屑長石砂巖和長石砂巖為主,含有少量的巖屑砂巖(圖2a)。砂巖整體致密,孔隙度主要分布在6%～14%之間(圖2b,c),平均孔隙度和滲透率分別為10.06%、1.999×10-3μm2,82.24%砂巖滲透率小于1×10-3μm2(圖2b,d),屬于典型的致密砂巖。

圖2 川中侏羅系沙溪廟組致密砂巖儲層巖石學(xué)與物性分布特征

孔隙類型主要包括粒間孔隙、粒內(nèi)孔隙、微裂隙及少量晶間孔。粒間孔包括剩余原生粒間孔(圖3a-c)和復(fù)合孔隙(圖3d),剩余原生孔隙邊緣平直,在鑄體薄片中呈現(xiàn)三角形或者多邊形,孔隙邊緣發(fā)育綠泥石薄膜(圖3a)。溶蝕孔主要是由長石(圖3e,f)、巖屑(圖3g)溶蝕形成的粒內(nèi)溶孔、沿濁沸石解理縫溶蝕形成的微裂縫(圖3h)及少量石英溶蝕形成的微小孔隙(圖3i)構(gòu)成。另外,黏土礦物晶間孔(圖3j,k)和云母層間孔(圖3l)也少量發(fā)育。喉道類型包括片狀喉道(圖3m)、縮頸喉道(圖3n)和孔隙縮小型喉道(圖3o),儲集空間類型多樣,孔喉大小從幾納米到幾十微米不等。

圖3 川中侏羅系沙溪廟組致密砂巖孔喉微觀特征

3.2 孔喉結(jié)構(gòu)參數(shù)特征

研究區(qū)樣品可分為4種類型,不同類型樣品的孔喉參數(shù)存在明顯差異(表1、圖4、圖5)。

圖4 川中侏羅系沙溪廟組致密砂巖不同類型樣品的毛細(xì)管壓力曲線、孔喉分布與孔喉含量特征

圖5 川中侏羅系沙溪廟組致密砂巖不同類型樣品的孔喉結(jié)構(gòu)參數(shù)與物性對比

Ⅰ類砂巖的排驅(qū)壓力小于0.4 MPa,最大進(jìn)汞飽和度平均值為74.47%,最高值為94.22%,剩余汞飽和度在29.18%～57.09%之間;孔喉半徑差異較大,平均孔喉半徑為0.598～1.605 μm,最大孔喉半徑均值大于4 μm,孔喉分選系數(shù)位于3～4之間;該類砂巖平均孔隙度為13.61%,平均滲透率為4.624×10-3μm2,儲層物性好。

Ⅱ類砂巖的排驅(qū)壓力在0.4～1.0MPa之間,最大進(jìn)汞飽和度為45.24%～90.47%,平均剩余汞飽和度為46.31%;平均孔喉半徑為0.307 μm,最大孔喉半徑在1 μm左右,孔喉分選系數(shù)小于Ⅰ類砂巖;樣品平均孔、滲分別為9.72%、0.375×10-3μm2。

Ⅲ類樣品的排驅(qū)壓力在1.0～5.0 MPa之間,最大進(jìn)汞飽和度平均值為60.53%;平均孔喉半徑為0.081～0.161 μm,最大孔喉半徑小于0.6 μm,孔喉分選系數(shù)在2～3之間;該類樣品孔隙度較小,平均值僅為6.68%,平均滲透率為0.047×10-3μm2。

Ⅳ類樣品排驅(qū)壓力大于5 MPa,進(jìn)汞飽和度低至70%以下,剩余汞飽和度在28%～36%之間;平均孔喉半徑分布在0.015～0.046 μm之間,最大孔喉半徑小于0.2 μm,孔喉分選系數(shù)小于2;該類型砂巖平均孔隙度僅為5.03%,平均滲透率為0.196×10-3μm2。

3.3 孔喉半徑及分布特征

沙溪廟組致密砂巖平均孔喉半徑在0.015～1.605 μm之間,不同樣品孔喉分布具有明顯差異性。根據(jù)前人的孔喉分類標(biāo)準(zhǔn)[26,40],基于孔隙尺寸可劃分為大孔(R≥1 000 nm)、中孔(100 nm≤R<1 000 nm)、過渡孔(10 nm≤R<100 nm)和微孔(R<10 nm)。

Ⅰ類樣品孔喉分布具有單峰或者雙峰的特征,最大峰值在1 μm處,部分樣品在0.01～0.02 μm處存在次級峰(圖4a1)。大孔含量在20%～40%之間,大、中孔占比最大,微孔占比最小(圖4a2)。

Ⅱ類樣品在0.3～0.6 μm之間存在最大峰值,0.01μm處存在次級峰值(圖4b1)。Ⅱ類砂巖大孔喉含量減少,微孔喉比例增大,中孔和過渡孔是主要的孔隙類型(圖4b2)。

Ⅲ類砂巖缺失大孔,以中孔和過渡孔為主(圖4c2),孔喉半徑在0.1～0.2 μm范圍存在最大值(圖4c1),次級峰值出現(xiàn)在0.01 μm處,微孔喉含量較Ⅱ類樣品略有增加。過渡孔是Ⅳ類砂巖最主要的孔隙類型,孔喉半徑主要分布在0.01～0.1 μm之間(圖4d1,d2)。

3.4 分形維數(shù)

川中沙溪廟組致密砂巖孔喉分形曲線主要為兩段式(圖6a-d),部分為三段式(圖6e)和單段式(圖6f)。所有樣品lgSHg-lgPc圖分形曲線的擬合相關(guān)系數(shù)平方(R2)均大于0.9,確保不同段的分形維數(shù)的準(zhǔn)確性。

圖6 川中侏羅系沙溪廟組致密砂巖孔喉分形曲線

沙溪廟組致密砂巖孔喉分形維數(shù)主要在2.0～2.4之間,兩個樣品分形維數(shù)高,分別為2.741、2.785。Ⅰ類樣品分形維數(shù)在2.04～2.17之間,平均值為2.11,分形曲線有三段式和兩段式2種形式,以兩段式為主,拐點處半徑大約為1 μm(圖6a,e);Ⅱ類樣品分形維數(shù)在2.13～2.27之間,平均值為2.20,分形曲線為兩段式,拐點半徑略小于1 μm(圖6b);Ⅲ類樣品分型維數(shù)平均值為2.26,拐點處孔喉半徑在0.1～1 μm之間(圖6c)。Ⅳ類樣品分形曲線有兩段式和單段式2種形態(tài),分形維數(shù)較高,由Ⅰ類到Ⅳ類孔喉分形維數(shù)逐漸變大。

4 討論

4.1 分形維數(shù)與孔喉結(jié)構(gòu)的關(guān)系

孔喉結(jié)構(gòu)與半徑大小是影響砂巖儲層滲流能力的重要因素[20],孔喉結(jié)構(gòu)參數(shù)與分形維數(shù)的相關(guān)性研究是致密砂巖孔喉表征的基礎(chǔ)。研究區(qū)致密砂巖孔喉分形維數(shù)與孔喉半徑呈負(fù)相關(guān)關(guān)系,孔喉半徑越大,分形維數(shù)越小;最大孔喉半徑大于1 μm的樣品分形維數(shù)均小于2.2,而分形維數(shù)大于2.3的樣品最大孔喉半徑小于0.3 μm(圖7a,b)。分形維數(shù)與排驅(qū)壓力存在正相關(guān)關(guān)系,與最大進(jìn)汞飽和度存在負(fù)相關(guān)關(guān)系,孔喉分形維數(shù)越小的樣品排驅(qū)壓力越小、最大進(jìn)汞飽度越大,該類樣品孔隙連通性高、非均質(zhì)性弱、儲層物性好(圖7c,d)。另外,孔喉分形維數(shù)與顆粒粒度、顆粒分選系數(shù)存在著明顯的相關(guān)性,顆粒越大,分選越好的樣品,分形維數(shù)越小(圖7e,f),這可能與強(qiáng)水動力條件下原生孔保存較多、連通性良好、孔喉非均質(zhì)性弱有關(guān)。儲層物性與分形維數(shù)具有明顯的負(fù)相關(guān)關(guān)系(圖7g,h),分形維數(shù)越低、滲透率越高。

圖7 川中侏羅系沙溪廟組致密砂巖分形維數(shù)與儲層物性、孔喉結(jié)構(gòu)關(guān)系

孔喉相對分選系數(shù)為孔喉分選系數(shù)與孔喉半徑均值的比值,反映了孔喉大小分布的均勻程度。分形維數(shù)與孔喉相對分選系數(shù)存在明顯的正相關(guān)關(guān)系(圖7i),相關(guān)系數(shù)高達(dá)0.851 5,進(jìn)一步證實了分形維數(shù)表征孔喉結(jié)構(gòu)與非均質(zhì)性的可行性。

以0.1 μm為分界,將孔喉分為2部分,建立了不同半徑孔喉含量與分形維數(shù)的關(guān)系圖(圖8a)。該圖顯示分形維數(shù)與孔喉含量有明顯的線性關(guān)系(R2=0.943 2),分形維數(shù)會隨著半徑大于0.1 μm孔喉含量的減少而增大。分形維數(shù)小于2.15,半徑大于0.1μm的孔喉含量不小于70%;分形維數(shù)上升至2.15～2.25之間時,半徑大于0.1 μm孔喉含量下降至50%～70%之間;而分形維數(shù)大于2.25時,半徑大于0.1 μm孔喉含量下降至50%以下,半徑小于0.1 μm孔喉占主導(dǎo)作用(圖8a)?？缀眍l率與分形維數(shù)的強(qiáng)相關(guān)性,表明分形維數(shù)是指示孔喉組成的重要指標(biāo)。

圖8 川中侏羅系沙溪廟組致密砂巖孔喉分布頻率與滲透率貢獻(xiàn)度、分形維數(shù)的關(guān)系

4.2 孔喉半徑分布對滲透率的控制作用

孔喉大小及組成是控制致密砂巖儲層物性的關(guān)鍵因素,致密砂巖滲透率隨著孔喉半徑的增大而增大,滲透率貢獻(xiàn)度隨著大孔含量的增加而上升(圖8b,c)。

在一次完整的高壓壓汞實驗中,累積滲透率貢獻(xiàn)曲線都會隨著壓力的升高經(jīng)歷3個不同的階段(圖9)。第一階段發(fā)生在壓力到達(dá)排驅(qū)壓力之前,在此階段累積滲透率始終為0。壓力到達(dá)排驅(qū)壓力之后,滲透率累積量迅速增加至95%以上,與此同時,進(jìn)汞量也呈現(xiàn)為一種快速上升趨勢,但進(jìn)汞飽和度僅為25%～42%,遠(yuǎn)遠(yuǎn)未達(dá)到最大值,此階段稱為“快速增長階段”。隨著壓力的進(jìn)一步增大,累積汞飽和度曲線依舊呈現(xiàn)為一種快速上升過程,但累積滲透率貢獻(xiàn)曲線進(jìn)入緩慢增加階段,滲透率增加不明顯。

圖9 川中侏羅系沙溪廟組致密砂巖孔喉分布對滲透率的控制作用

致密砂巖滲透率主要受控于小部分半徑相對較大的孔喉[20],而相對小的孔喉對滲透率的貢獻(xiàn)度有限。川中沙溪廟組致密砂巖最小流動孔喉半徑平均值為0.059 μm,屬于過渡孔范圍;半徑大于0.1 μm的孔喉滲透率貢獻(xiàn)度達(dá)98%以上(圖8c、圖9、圖10)。半徑大于1μm大孔對滲透率有著絕對的控制作用,大孔含量越多,滲透率貢獻(xiàn)度越大(R2=0.972 5)(圖8b)。中孔滲透率貢獻(xiàn)度明顯高于過渡孔和微孔,過渡孔與微孔所占進(jìn)汞飽和度與滲透率貢獻(xiàn)度之間具有正相關(guān)關(guān)系,相關(guān)系數(shù)R2分別為0.59與0.79,但兩類孔隙滲透率貢獻(xiàn)度分別集中在0.01%～10%與0%～1%之間,暗示這類小孔隙對總體滲透率貢獻(xiàn)有限,但作用不可忽視。

圖10 川中侏羅系沙溪廟組致密砂巖孔喉相對百分含量與滲透率貢獻(xiàn)柱狀圖

不同類型砂巖樣品孔喉分布與滲透率貢獻(xiàn)特征差異性明顯。Ⅰ類樣品大孔以小于40%的含量貢獻(xiàn)了90%以上的滲透率;Ⅱ、Ⅲ類砂巖大孔的滲透率貢獻(xiàn)度降低,中孔成為最主要的滲透率貢獻(xiàn)者;Ⅳ類儲層中過渡孔含量、滲透率貢獻(xiàn)度均達(dá)到最高(圖10)。

4.3 基于孔喉半徑及分形維數(shù)的滲透率預(yù)測模型

利用巖心實測孔滲資料,建立滲透率與孔隙度的擬合關(guān)系(圖11a)。滲透率與孔隙度的線性、指數(shù)、對數(shù)擬合關(guān)系較差,相關(guān)系數(shù)小于0.5,而滲透率與孔隙度冪函數(shù)擬合R2雖然達(dá)到0.5,但依舊無法滿足儲層精細(xì)評價的要求[41]。前人研究中提出了多種基于孔隙度、壓汞參數(shù)以及其他參數(shù)的滲透率模型[42-44],但大多數(shù)模型并未考慮孔喉分形特征對滲透率的貢獻(xiàn)。

圖11 川中侏羅系沙溪廟組致密砂巖實測滲透率與孔隙度、預(yù)測滲透率交會圖

基于前面的討論,選取孔隙度、分形維數(shù)、不同孔喉半徑對滲透率進(jìn)行多元回歸擬合與相關(guān)性分析(表2),并選取相關(guān)性最好的公式作為滲透率預(yù)測模型。最終確定了基于最大孔喉半徑、分形維數(shù)、孔隙度的滲透率預(yù)測公式:

表2 川中侏羅系沙溪廟組致密砂巖滲透率與分形維數(shù)、孔隙度、孔喉半徑相關(guān)性

lgK=4.452Df+2.014lgRa+0.065Ф-10.588

(1)

式中:孔隙度(Φ)表示砂巖中孔隙含量,單位%;分形維數(shù)(Df)表征孔喉的非均質(zhì)性與組成特征;最大孔喉半徑(Ra)指示孔喉大小,單位μm。

通過公式(1)對研究區(qū)及鄰近研究區(qū)樣品進(jìn)行預(yù)測并驗證(圖11b),除兩個異常點之外,預(yù)測結(jié)果與實測結(jié)果整體相關(guān)性極高(R2大于0.9)。圖11b中有兩個預(yù)測滲透率明顯小于實測滲透率的紅色異常點,對兩個異常樣品薄片觀察發(fā)現(xiàn)發(fā)育明顯的微裂縫,造成滲透率偏高。由此表明,基于孔隙度、分形維數(shù)、最大孔喉半徑的回歸公式可以更準(zhǔn)確地計算致密砂巖“孔隙型”儲層的滲透率,但對微裂縫較為發(fā)育的砂巖預(yù)測效果不理想。

5 結(jié)論

(1)川中沙溪廟組儲層具有巖性致密、孔喉結(jié)構(gòu)復(fù)雜、非均質(zhì)性強(qiáng)的特點。原生孔隙、次生溶蝕孔隙、微裂縫與晶間孔均有發(fā)育,孔隙半徑從幾納米到幾十微米均有分布。Ⅰ類砂巖排驅(qū)壓力小于0.4MPa,平均滲透率大于4×10-3μm2,半徑大于0.1 μm的大孔和中孔占絕對主導(dǎo)地位;Ⅱ類樣品排驅(qū)壓力在0.4～1.0 MPa之間,平均孔、滲分別為9.72%、0.375×10-3μm2,分形維數(shù)為2.20,半徑大于0.1 μm中孔含量上升;Ⅲ類砂巖滲透率僅為0.047×10-3μm2,以中孔和過渡孔為主;Ⅳ類砂巖排驅(qū)壓力大于5 MPa,孔隙度低,且因缺乏大孔導(dǎo)致滲透率較低。

(2)致密儲層最小流動孔喉半徑平均值為0.059 μm,處于過渡孔范圍內(nèi),半徑大于0.1 μm的孔喉滲透率貢獻(xiàn)度達(dá)98%以上。大孔貢獻(xiàn)了Ⅰ類砂巖90%以上的滲透率,中孔貢獻(xiàn)了Ⅱ、Ⅲ類砂巖絕大部分滲透率,Ⅳ類砂巖中過渡孔滲透率貢獻(xiàn)最大。大孔、中孔含量決定了滲透率的高低,過渡孔與微孔滲透率貢獻(xiàn)度有限,但其含量與滲透率貢獻(xiàn)關(guān)系明顯,不可忽視。

(3)分形維數(shù)是指示孔喉組成的重要指標(biāo)。沙溪廟組致密砂巖分形維數(shù)主要介于2.0～2.4之間,分形曲線呈現(xiàn)多段式的特征。分形維數(shù)與孔喉半徑、進(jìn)汞飽和度呈負(fù)相關(guān)關(guān)系,與排驅(qū)壓力呈正相關(guān)關(guān)系。半徑大于0.1 μm的孔喉含量不小于70%時,分形維數(shù)主要在2.00～2.15之間;分形維數(shù)在2.15～2.25之間時,半徑大于0.1 μm孔喉含量在50%～70%之間;而分形維數(shù)大于2.25時,半徑小于0.1 μm孔喉占主導(dǎo)作用,含量不小于50%。

(4)孔隙含量、孔喉組成、孔喉半徑是影響滲透率的關(guān)鍵因素,基于分形維數(shù)、孔隙度、最大孔喉半徑建立的滲透率預(yù)測模型在沙溪廟組“孔隙型”儲層中具有極高的適用性。

利益沖突聲明/Conflict of Interests

所有作者聲明不存在利益沖突。

All authors disclose no relevant conflict of interests.

作者貢獻(xiàn)/Authors’Contributions

陳少云、楊勇強(qiáng)、邱隆偉完成了研究思路制定、數(shù)據(jù)分析等方面的工作。王小娟提供了部分研究區(qū)的原始數(shù)據(jù),指出了研究區(qū)實際問題與研究內(nèi)容的聯(lián)系。楊保良、葉熱杰普·哈布臘什木對部分?jǐn)?shù)據(jù)進(jìn)行整理,并參與圖片繪制工作。陳少云、楊勇強(qiáng)全程參與論文寫作和修改。所有作者均閱讀并同意最終稿件的提交。

CHEN Shaoyun, YANG Yongqiang and QIU Longwei carried out the tasks of formulating research ideas and analyzing data. WANG Xiaojuan provided some of the raw data of the study area, pointed out the connection between the actual problems of the study area and the content of the study. YANG Baoliang and Erejep HABILAXIM organized part of the data and participated in the drawing of pictures. CHEN Shaoyun and YANG Yongqiang participated in writing and revising the paper. All the authors have read the last version of paper and consented for submission.