崔 仁，孟兆敏，譚巧林，邵長(zhǎng)高

（1.國(guó)土資源部海底礦產(chǎn)資源重點(diǎn)實(shí)驗(yàn)室 廣州海洋地質(zhì)調(diào)查局，廣東 廣州 510075；

2.上海工程技術(shù)大學(xué)，上海 201620；3.廣東省地質(zhì)測(cè)繪院，廣東 廣州 510800）

BSR測(cè)線導(dǎo)航數(shù)據(jù)入庫(kù)接口模型的設(shè)計(jì)與實(shí)現(xiàn)

崔 仁1，孟兆敏2，譚巧林3，邵長(zhǎng)高1

（1.國(guó)土資源部海底礦產(chǎn)資源重點(diǎn)實(shí)驗(yàn)室 廣州海洋地質(zhì)調(diào)查局，廣東 廣州 510075；

2.上海工程技術(shù)大學(xué)，上海 201620；3.廣東省地質(zhì)測(cè)繪院，廣東 廣州 510800）

BSR測(cè)線導(dǎo)航數(shù)據(jù)通常以單個(gè)文件形式進(jìn)行管理，數(shù)據(jù)的檢索查詢需手工進(jìn)行，管理復(fù)雜且易出錯(cuò)。采用數(shù)據(jù)庫(kù)形式對(duì)單個(gè)BSR測(cè)線導(dǎo)航數(shù)據(jù)進(jìn)行管理，設(shè)計(jì)了BSR測(cè)線導(dǎo)航數(shù)據(jù)入庫(kù)接口模型，并在VB環(huán)境下開發(fā)了該模型。模型包含BSR測(cè)線導(dǎo)航數(shù)據(jù)庫(kù)結(jié)構(gòu)、數(shù)據(jù)索引方式和導(dǎo)航數(shù)據(jù)入庫(kù)方法。模型對(duì)BSR測(cè)線導(dǎo)航數(shù)據(jù)進(jìn)行了統(tǒng)一管理，提高了數(shù)據(jù)管理的效率，改善了數(shù)據(jù)索引機(jī)制；同時(shí)整合的測(cè)線導(dǎo)航數(shù)據(jù)也為未來(lái)測(cè)線規(guī)劃提供了基礎(chǔ)，特別是在海洋地質(zhì)調(diào)查等數(shù)據(jù)量大的領(lǐng)域，可以較快的運(yùn)算速度提取所需數(shù)據(jù)。

海洋地質(zhì)調(diào)查；似海底反射波；導(dǎo)航數(shù)據(jù)；數(shù)據(jù)庫(kù)

天然氣水合物作為最有可能替代石油的下一代清潔能源，其研究越來(lái)越受重視[1-3]。2013年我國(guó)在南海鉆獲了高純度的天然氣水合物，為我國(guó)的商業(yè)化開采奠定了基礎(chǔ)。近年來(lái)，國(guó)家對(duì)天然氣水合物的調(diào)查產(chǎn)生了海量數(shù)據(jù)，如何有效、安全地管理這些數(shù)據(jù)資料成為亟待解決的問題。BSR指含水合物的地層在地震反射剖面上常常會(huì)出現(xiàn)一強(qiáng)振幅的連續(xù)反射波，大致與海底反射波平行，故又稱似海底反射波[1]。研究表明，BSR是目前天然氣水合物識(shí)別的重要標(biāo)志之一[4-7]。海洋地質(zhì)調(diào)查獲取的BSR測(cè)線導(dǎo)航數(shù)據(jù)通常以單個(gè)文件形式進(jìn)行管理，需手工進(jìn)行數(shù)據(jù)檢索與查詢，管理復(fù)雜且易出錯(cuò)；目前工作中尚無(wú)有效方式管理全部的BSR測(cè)線導(dǎo)航數(shù)據(jù)。本文提出了采用數(shù)據(jù)庫(kù)形式管理單個(gè)BSR測(cè)線導(dǎo)航數(shù)據(jù)的方法；設(shè)計(jì)了BSR測(cè)線導(dǎo)航數(shù)據(jù)的數(shù)據(jù)接口模型；并在VB環(huán)境下實(shí)現(xiàn)了BSR測(cè)線導(dǎo)航數(shù)據(jù)模型。模型為單條測(cè)線的索引、區(qū)域內(nèi)測(cè)線的索引以及任意屬性的索引均提供了基本性支持；同時(shí)整合的測(cè)線導(dǎo)航數(shù)據(jù)也為未來(lái)測(cè)線規(guī)劃提供了基礎(chǔ)，特別是在海洋地質(zhì)調(diào)查等數(shù)據(jù)量大的領(lǐng)域，可以較快的運(yùn)算速度提取所需數(shù)據(jù)[5]。

1 BSR測(cè)線導(dǎo)航數(shù)據(jù)結(jié)構(gòu)及數(shù)據(jù)索引

1.1 數(shù)據(jù)結(jié)構(gòu)

本文將BSR測(cè)線導(dǎo)航數(shù)據(jù)管理分為3個(gè)部分（見圖1）。①調(diào)查測(cè)線數(shù)據(jù)，負(fù)責(zé)管理海上野外勘查過(guò)程中所有的測(cè)線數(shù)據(jù)信息。調(diào)查測(cè)線數(shù)據(jù)表是空間圖形數(shù)據(jù)表，包含數(shù)據(jù)集編號(hào)、工區(qū)編號(hào)、測(cè)線編號(hào)、測(cè)線長(zhǎng)度、調(diào)查類別、起始日期、結(jié)束日期、起始經(jīng)度、起始緯度、終止經(jīng)度、終止緯度、起始點(diǎn)號(hào)、終止點(diǎn)號(hào)、調(diào)查比例、導(dǎo)航數(shù)據(jù)文件、數(shù)據(jù)文件路徑、坐標(biāo)系統(tǒng)等。②BSR測(cè)線信息，負(fù)責(zé)管理有BSR剖面圖像測(cè)線段的基本信息，是對(duì)調(diào)查測(cè)線信息的延伸，即調(diào)查側(cè)線上含有BSR信息的測(cè)線段信息。其內(nèi)容包含測(cè)線編號(hào)、BSR測(cè)線號(hào)、起點(diǎn)經(jīng)度、起點(diǎn)緯度、終點(diǎn)經(jīng)度、終點(diǎn)緯度、起點(diǎn)炮號(hào)、終點(diǎn)炮號(hào)、起點(diǎn)CDP號(hào)、終點(diǎn)CDP號(hào)、平均雙程時(shí)間、平均深度、BSR類型、波形組合、弱振幅、空白帶、備注等。③導(dǎo)航數(shù)據(jù)信息，即野外勘查測(cè)線和BSR測(cè)線段的導(dǎo)航信息，主要有測(cè)線編號(hào)、炮號(hào)、經(jīng)度、緯度4個(gè)方面信息。

圖1 BSR測(cè)線導(dǎo)航數(shù)據(jù)總體數(shù)據(jù)結(jié)構(gòu)

1.2 數(shù)據(jù)索引

調(diào)查測(cè)線表中主鍵為測(cè)線編號(hào)，外鍵為數(shù)據(jù)集編號(hào)，數(shù)據(jù)集編號(hào)負(fù)責(zé)與整個(gè)調(diào)查項(xiàng)目的索引；BSR測(cè)線表中主鍵為BSR測(cè)線號(hào)，外鍵為測(cè)線編號(hào)，測(cè)線編號(hào)負(fù)責(zé)與調(diào)查測(cè)線表的索引；地震測(cè)線導(dǎo)航數(shù)據(jù)表主鍵為炮號(hào)，外鍵為測(cè)線編號(hào)，測(cè)線編號(hào)負(fù)責(zé)與調(diào)查測(cè)線表的索引。BSR測(cè)線表與地震測(cè)線導(dǎo)航數(shù)據(jù)表之間索引較為復(fù)雜，首先通過(guò)調(diào)查測(cè)線表索引出調(diào)查測(cè)線編號(hào)，然后通過(guò)起始炮號(hào)、終止炮號(hào)字段索引出地震測(cè)線導(dǎo)航數(shù)據(jù)表中的數(shù)據(jù)。數(shù)據(jù)集編號(hào)采用廣州海洋地質(zhì)調(diào)查通用數(shù)據(jù)集編碼規(guī)則，按“項(xiàng)目代碼—年份—調(diào)查船—航次”編碼，項(xiàng)目代碼取自廣州海洋地質(zhì)調(diào)查局?jǐn)?shù)據(jù)庫(kù)的數(shù)據(jù)項(xiàng)目代碼表[3]。

2 數(shù)據(jù)接口設(shè)計(jì)和實(shí)現(xiàn)方法

2.1 數(shù)據(jù)入庫(kù)及獲取流程

按照BSR測(cè)線導(dǎo)航數(shù)據(jù)結(jié)構(gòu)，本文設(shè)計(jì)了BSR測(cè)線導(dǎo)航數(shù)據(jù)入庫(kù)及獲取流程圖，如圖2所示。

圖2 BSR測(cè)線導(dǎo)航數(shù)據(jù)入庫(kù)及獲取流程圖

在數(shù)據(jù)入庫(kù)方面，調(diào)查測(cè)線信息和BSR測(cè)線信息以基本信息形式進(jìn)行數(shù)據(jù)錄入或數(shù)據(jù)導(dǎo)入，要求信息完整，索引結(jié)構(gòu)字段主外鍵不能為空；地震測(cè)線導(dǎo)航數(shù)據(jù)以文件信息形式導(dǎo)入，先讀取導(dǎo)航數(shù)據(jù)的原始文件，再進(jìn)行測(cè)線編號(hào)的信息填充，填充完畢后入庫(kù)。

在數(shù)據(jù)獲取中，首先讀取BSR測(cè)線信息，通過(guò)測(cè)線編號(hào)在調(diào)查測(cè)線表中進(jìn)行信息索引查詢；再把BSR測(cè)線和調(diào)查測(cè)線兩個(gè)數(shù)據(jù)表信息進(jìn)行關(guān)聯(lián)；然后將關(guān)聯(lián)后的數(shù)據(jù)表與地震測(cè)線導(dǎo)航數(shù)據(jù)表數(shù)據(jù)進(jìn)行索引，尋找導(dǎo)航測(cè)線數(shù)據(jù)表中的測(cè)線信息和炮號(hào)信息；最終形成完整的BSR測(cè)線導(dǎo)航數(shù)據(jù)表。

2.2 數(shù)據(jù)原始文件讀取及入庫(kù)實(shí)現(xiàn)方法

根據(jù)BSR測(cè)線導(dǎo)航數(shù)據(jù)入庫(kù)及獲取流程，本文設(shè)計(jì)了一個(gè)在VB環(huán)境下的數(shù)據(jù)入庫(kù)方法。其難點(diǎn)主要為導(dǎo)航數(shù)據(jù)文件的導(dǎo)入和數(shù)據(jù)入庫(kù)（圖3）。

圖3 導(dǎo)航原始文件讀取

導(dǎo)航數(shù)據(jù)文件解析及導(dǎo)入具體實(shí)現(xiàn)代碼為：

Private Sub txtLNavFile_LostFocus() '自動(dòng)填充測(cè)線編號(hào)，起止點(diǎn)坐標(biāo)，把文本讀到內(nèi)存里

On Error GoTo ErrOccurs

Dim strTmp As String

strTmp = FileChk(txtLNavFile)

If strTmp <＞ "" Then

txtLineID.Text = strTmp '將文件名指定為默認(rèn)的測(cè)線編號(hào)

txtLFilePath.Text = txtLNavFile.Text '保存文本數(shù)據(jù)至內(nèi)存

Dim fso

Set fso = CreateObject("Scripting.FileSystemObject")

Dim tsFile As TextStream

Set tsFile = fso.OpenTextFile(txtLNavFile.Text, ForReading)

Dim curLine() As String, rLine As String, iNum As Integer, eLine As Integer

Dim curShot As String, curLon As Double, curLat As Double '當(dāng)前行之炮號(hào)，首尾兩點(diǎn)無(wú)條件保留，其余按200的間距抽稀

strLons = ""

strLats = ""

curShot = ""

iNum = -1

eLine = tsFile.Line

Do Until tsFile.AtEndOfStream

‘------原始文件讀取及空格信息處理-------------

rLine = Trim(tsFile.ReadLine)

rLine = Replace(rLine, " ", " ")

rLine = Replace(rLine, " ", " ")

rLine = Replace(rLine, " ", " ")

rLine = Replace(rLine, " ", " ")

rLine = Replace(rLine, " ", " ")

curLine = Split(rLine) '解析出的數(shù)據(jù)為(炮號(hào)，緯度，經(jīng)度)

iNum = iNum + 1

ReDim Preserve wShotID(iNum)

ReDim Preserve wLats(iNum)

ReDim Preserve wLons(iNum)

wShotID(iNum) = curLine(0)

wLats(iNum) = CLng(curLine(1) * 100) / 100

wLons(iNum) = CLng(curLine(2) * 100) / 100

If eLine Mod 200 = 1 Then '取點(diǎn)的行號(hào)為：1 201 401 601 ...，如果坐標(biāo)有誤，則跳過(guò)

curShot = curLine(0)

curLat = DegCon(CDbl(curLine(1)))

curLon = DegCon(CDbl(curLine(2)))

If curLat ＞= 0 And curLat <= 90 And curLon ＞= 0 And curLat <= 180 Then

strLats = strLats & curLat & " "

strLons = strLons & curLon & " "

End If

End If

eLine = tsFile.Line

DoEvents

Loop

txtLSPoint.Text = wShotID(0)

txtLSLatitude.Text = wLats(0)

txtLSLongitude.Text = wLons(0)

txtLEPoint.Text = wShotID(iNum)

txtLELatitude.Text = wLats(iNum)

txtLELongitude.Text = wLons(iNum)

If curShot <＞ curLine(0) Then

strLats = strLats & DegCon(CDbl(curLine(1)))

strLons = strLons & DegCon(CDbl(curLine(2)))

Else

strLats = Trim(strLats)

strLons = Trim(strLons)

End If

Set fso = Nothing

tsFile.Close

Set tsFile = Nothing

Else

txtLineID.Text = ""

txtLFilePath.Text = ""

ReDim wShotID(0)

ReDim wLats(0)

ReDim wLons(0)

ReDim curLine(0)

txtLSPoint.Text = ""

txtLSLatitude.Text = ""

txtLSLongitude.Text = ""

txtLEPoint.Text = ""

txtLELatitude.Text = ""

txtLELongitude.Text = ""

End If

strPreNavFile = txtLNavFile.Text

Exit Sub

ErrOccurs：

MsgBox "所用文件不是有效的導(dǎo)航數(shù)據(jù)文件！", , "原始文件讀取失?。?

End Sub

代碼通過(guò)解析導(dǎo)航數(shù)據(jù)文件的原始數(shù)據(jù)對(duì)數(shù)據(jù)信息進(jìn)行提取，再將提取出導(dǎo)航坐標(biāo)信息的屬性信息進(jìn)行數(shù)據(jù)入庫(kù)，空間信息進(jìn)行數(shù)據(jù)繪圖形成地圖點(diǎn)坐標(biāo)信息。

2.3 坐標(biāo)投影轉(zhuǎn)換

測(cè)線導(dǎo)航數(shù)據(jù)提取后要進(jìn)行WGS84投影下經(jīng)緯度坐標(biāo)與墨卡托投影的轉(zhuǎn)換。墨卡托投影正解公式為[8]：

式中，k0為原點(diǎn)處比例因子；θ1為墨卡托投影中第一標(biāo)準(zhǔn)緯線。

式中，E為墨卡托投影緯度坐標(biāo)；a為地球半徑，單位為m；FE為經(jīng)度偏移；λ為經(jīng)度；λ0為起點(diǎn)經(jīng)度。

N= FN+aK0ln{tan(π/4+θ/2)[(1-esinθ)/(1+esinθ)](e/2)}式中，N為墨卡托投影經(jīng)度坐標(biāo)；FN為緯度偏移；θ為測(cè)量經(jīng)度。

基于ArcEngine的坐標(biāo)轉(zhuǎn)換方法可詳見參考文獻(xiàn)[8]，墨卡托投影反解公式可詳見參考文獻(xiàn)[9]、[10]。

2.4 BSR測(cè)線生成

坐標(biāo)轉(zhuǎn)換后，再生成數(shù)據(jù)測(cè)線，其核心代碼為：

Case "BSR測(cè)線導(dǎo)航" '地震測(cè)線導(dǎo)航

If Trim(cmbELineID.Text) = "" Then

MsgBox "測(cè)線編號(hào)不能為空，請(qǐng)重新輸入！"

cmbELineID.SetFocus

ElseIf Trim(txtEShotID.Text) = "" Then

MsgBox "炮號(hào)不能為空,請(qǐng)重新輸入!"

txtEShotID.SetFocus

Else

strObject = "select OBJECTID from sde." & strCurSeaArea & "_調(diào)查項(xiàng)目基本信息 "

ObjectRs.Open strObject, Conn, adOpenKeyset, adLockOptimistic

While Not ObjectRs.EOF

m = m + 1

ObjectRs.MoveNext

Wend

If m <＞ 0 Then

ObjectRs.MoveFirst

End If

For n = 1 To m

If ObjectRs.Fields(0).Value = "" Then

Exit For

Else： ObjectRs.MoveNext

End If

Next

ObjectRs.Close

'查詢是否有重復(fù)記錄

sqlStr = "select 測(cè) 線 編 號(hào) , 炮 號(hào) from sde." & strCurSeaArea & "_BSR測(cè)線導(dǎo)航數(shù)據(jù) where trim(測(cè)線編號(hào))='" _& Trim(cmbELineID.Text) & "' and trim(炮號(hào))='" & Trim(txtEShotID.Text) & "'"

prjRS.Open sqlStr, Conn, adOpenKeyset, adLockOptimistic

If prjRS.BOF And prjRS.EOF Then

prjRS.Close

Set prjRS = Nothing

'允許插入一條新記錄

sqlStr = "insert into sde." & strCurSeaArea & "_BSR測(cè)線導(dǎo)航數(shù)據(jù) (OBJECTID,測(cè)線編號(hào), 炮號(hào), 經(jīng)度, 緯度) values ('" & n & "','"

sqlStr = sqlStr & Trim(cmbELineID) & "', '" & Trim(txtEShotID) & "', '"

sqlStr = sqlStr & Trim(txtELongitude) & "', '" & Trim(txtELatitude) & "')"

Conn.Execute sqlStr

openSDEWorkspace.ExecuteSql (sqlStr)

If Conn.State = adStateOpen Then Conn.Close

MsgBox "成功錄入一條記錄！", , "提示："

'********* (注：?jiǎn)蝹€(gè)導(dǎo)航點(diǎn)不生成圖形) *********

Else

prjRS.Close

Set prjRS = Nothing

If Conn.State = adStateOpen Then Conn.Close

MsgBox "數(shù)據(jù)庫(kù)中已有使用[" & Trim(cmbELineID.Text) & "]作為測(cè)線編號(hào), [" & Trim(txtEShotID.Text) & "]作為炮號(hào)的記錄存在！", , "重復(fù)錄入錯(cuò)誤："

End If

cmbELineID.SetFocus

End If

BSR測(cè)線導(dǎo)航數(shù)據(jù)生成前需進(jìn)行索引過(guò)程主鍵檢查，檢查通過(guò)后才能進(jìn)行打點(diǎn)成線。

3 模型應(yīng)用

實(shí)施天然氣水合物資源勘查試采工作、加快水合物的商業(yè)開發(fā)，將是解決我國(guó)能源供需矛盾、保證國(guó)家能源安全的有效途徑之一。為此，廣州海洋地質(zhì)調(diào)查局近年實(shí)施了大規(guī)模的水合物勘查工作，獲取了大量地球物理、地球化學(xué)、地質(zhì)學(xué)、環(huán)境科學(xué)等多個(gè)學(xué)科的數(shù)據(jù)。為了保證對(duì)勘查數(shù)據(jù)的有效安全管理，實(shí)現(xiàn)專業(yè)化角度的共享服務(wù)，廣州海洋地質(zhì)調(diào)查局建設(shè)了水合物數(shù)據(jù)庫(kù)管理系統(tǒng)（圖4）。目前系統(tǒng)已應(yīng)用于廣州海洋地質(zhì)調(diào)查局的實(shí)際工作中，持續(xù)為實(shí)際勘查工作服務(wù)。系統(tǒng)的數(shù)據(jù)結(jié)構(gòu)以及對(duì)BSR測(cè)線導(dǎo)航數(shù)據(jù)的管理均采用本文設(shè)計(jì)的BSR測(cè)線導(dǎo)航數(shù)據(jù)入庫(kù)接口模型。

BSR測(cè)線導(dǎo)航數(shù)據(jù)是海洋地質(zhì)調(diào)查野外數(shù)據(jù)的重要內(nèi)容，針對(duì)目前處理BSR數(shù)據(jù)技術(shù)的缺乏，按照本文設(shè)計(jì)的BSR測(cè)線導(dǎo)航數(shù)據(jù)結(jié)構(gòu)和實(shí)現(xiàn)方法，在水合物數(shù)據(jù)庫(kù)管理系統(tǒng)中開發(fā)了一個(gè)數(shù)據(jù)錄入模塊。該模塊實(shí)現(xiàn)了BSR測(cè)線導(dǎo)航數(shù)據(jù)原始文件讀取、屬性數(shù)據(jù)錄入、空間數(shù)據(jù)成圖等功能，負(fù)責(zé)測(cè)線數(shù)據(jù)及BSR測(cè)線導(dǎo)航數(shù)據(jù)的管理及錄入。該模塊的實(shí)際應(yīng)用證明了本文設(shè)計(jì)模型的可靠性和可推廣性。

圖4 水合物數(shù)據(jù)庫(kù)管理系統(tǒng)界面

4 結(jié) 語(yǔ)

本文設(shè)計(jì)了BSR測(cè)線導(dǎo)航數(shù)據(jù)的數(shù)據(jù)庫(kù)結(jié)構(gòu)及索引機(jī)制，實(shí)現(xiàn)了對(duì)BSR測(cè)線導(dǎo)航數(shù)據(jù)的安全管理；還設(shè)計(jì)了BSR測(cè)線導(dǎo)航數(shù)據(jù)原始文件及數(shù)據(jù)庫(kù)之間的入庫(kù)接口，并在VB環(huán)境下實(shí)現(xiàn)了數(shù)據(jù)接口，提供了專業(yè)的代碼。目前BSR測(cè)線導(dǎo)航數(shù)據(jù)入庫(kù)接口模型已進(jìn)入實(shí)際應(yīng)用階段，應(yīng)用效果良好，具有解決實(shí)際工作的能力，值得借鑒。

[1] 金慶煥,張光學(xué),楊木壯,等.天然氣水合物資源概論[M].北京：科學(xué)出版社,2006：2-5

[2] 黃永樣,張光學(xué).我國(guó)海域天然氣水合物地質(zhì)：地球物理特征及前景[M].北京：地質(zhì)出版社,2009

[3] 陳宏文,王立飛,邵長(zhǎng)高,等.天然氣水合物屬性及空間數(shù)據(jù)庫(kù)結(jié)構(gòu)設(shè)計(jì)[J].海洋地質(zhì)前沿,2012,28(5)：59-63

[4] 沙志彬,楊木壯,梁金強(qiáng),等.BSR的反射波特征及其對(duì)天然氣水合物識(shí)別的應(yīng)用[J].南海地質(zhì)研究,2003(1)：55-61

[5] 徐華寧,楊勝雄,鄭曉東,等. 南中國(guó)海神狐海域天然氣水合物地震識(shí)別及分布特征[J].地球物理學(xué)報(bào),2010,53(7)：1 691-1 698

[6] 鄒大鵬,盧博,閻貧,等.南海北部海底沉積物在溫度變化下的三種聲速類型[J].地球物理學(xué)報(bào),2012,55(3)：1 017-1 024

[7] 宋海斌,吳時(shí)國(guó),江為為.南海東北部973剖面BSR及其熱流特征[J].地球物理學(xué)報(bào),2007,50(5)：1 508-1 517

[8] 邵長(zhǎng)高,譚建軍,荊麗梅,等.海洋小比例尺地圖精確測(cè)量及計(jì)算方法[J].地理與地理信息科學(xué),2009,25(2)：42-45

[9] International Association of Oil & Gas Producers. Coordinate Conversions and Transformation Including Formulas[M]. OGP Publicaiton,2013：35-39

[10] 邵長(zhǎng)高,崔仁.地震剖面解釋數(shù)據(jù)提取接口方法[J].測(cè)繪科學(xué),2015(專刊)：143-146

P208

B

1672-4623（2017）05-0025-04

10.3969/j.issn.1672-4623.2017.0050.8

崔仁，主要從事海洋地質(zhì)調(diào)查數(shù)據(jù)入庫(kù)、數(shù)據(jù)結(jié)構(gòu)制定、原始資料管理等方面研究。

2016-04-27。

項(xiàng)目來(lái)源：國(guó)家高技術(shù)研究發(fā)展計(jì)劃資助項(xiàng)目（2013AA092501、2013AA0925010401）；“127”國(guó)家專項(xiàng)資助項(xiàng)目（DD20160227、DD20160227-6、GZH201100312）；“729”國(guó)家專項(xiàng)資助項(xiàng)目（GZH201200512、GZH201200512-2）。