藍(lán)文威 陳 晨 張 金 張家琦 黎 峰 高軍峰*

1(中南民族大學(xué)生物醫(yī)學(xué)工程學(xué)院, 武漢 430074)

2(認(rèn)知科學(xué)國(guó)家民委重點(diǎn)實(shí)驗(yàn)室，武漢 430074)

引言

大腦功能狀態(tài)的研究一直是人類(lèi)探索自身奧秘的熱點(diǎn)，隨著近年腦科學(xué)、神經(jīng)科學(xué)及信息科學(xué)的快速發(fā)展，大腦功能狀態(tài)的研究日新月異，關(guān)于大腦情緒狀態(tài)的研究也逐漸深入。情緒作為人類(lèi)大腦的高級(jí)功能，綜合了個(gè)人的思想和行為，在一定程度上表現(xiàn)出對(duì)外界或者自身刺激的反應(yīng)，同時(shí)伴隨著一系列的生理以及心理上的變化。在情緒判別方面，之前人們對(duì)于情緒的理解皆基于自身對(duì)自身的內(nèi)在思想或者他人對(duì)自身的外在行為進(jìn)行主觀評(píng)估，實(shí)際上，通過(guò)生理信號(hào)對(duì)人類(lèi)大腦情緒狀態(tài)的客觀評(píng)估手段則更具有科學(xué)性，腦電信號(hào)(electroencephalogram, EEG)由于其高效率和低成本的優(yōu)勢(shì)，在近年來(lái)檢測(cè)情緒狀態(tài)的研究中得到了廣泛的應(yīng)用[1-2]。

國(guó)內(nèi)外已有大量關(guān)于大腦情緒狀態(tài)的研究。有報(bào)告指出，參與情緒的腦區(qū)主要集中于杏仁核、基底神經(jīng)節(jié)、前額皮層、扣帶回皮層和頂葉皮層[3]。另外，早在20世紀(jì)90年代Canli在使用正電子發(fā)射型計(jì)算機(jī)斷層顯像(positron emission computed tomography，PET)功能性成像技術(shù)對(duì)大腦情感的半球不對(duì)稱(chēng)研究時(shí)，已經(jīng)發(fā)現(xiàn)積極(或與進(jìn)近相關(guān))的情緒更傾向于偏向左半球，而消極(或與退縮相關(guān))的情緒則是更加強(qiáng)烈地偏向右半球[4]。Zibman等[5]進(jìn)行皮質(zhì)-皮質(zhì)配對(duì)聯(lián)合刺激(ccPAS)實(shí)驗(yàn)，結(jié)果發(fā)現(xiàn)，外側(cè)前額葉皮層的左右側(cè)的電生理和行為均顯示出半球不對(duì)稱(chēng)性的變化。Crawford等[6]在受試覺(jué)醒和催眠狀態(tài)下引發(fā)愉快和悲傷的情緒，并且對(duì)大腦前額區(qū)(F3F4)、中央?yún)^(qū)(C3C4)和頂區(qū)(P3P4)的EEG窄波頻帶進(jìn)行分析。研究發(fā)現(xiàn)，低頻alpha波(8.0～9.5 Hz)在頂區(qū)左右大腦半球存在差異。Lu等[7]使用左前(后)、右前(后)、左半球、右半球、前半球、后半球等8個(gè)不同大腦區(qū)域的腦電信號(hào)，基于動(dòng)態(tài)熵的模式學(xué)習(xí)來(lái)識(shí)別個(gè)體之間的積極和消極情緒。

隨著計(jì)算機(jī)技術(shù)的蓬勃發(fā)展，人機(jī)交互技術(shù)被廣泛應(yīng)用于各大領(lǐng)域。虛擬現(xiàn)實(shí)技術(shù)(virtual reality, VR)是20世紀(jì)快速發(fā)展起來(lái)的一項(xiàng)新興的實(shí)用技術(shù)，通過(guò)計(jì)算機(jī)模擬虛擬環(huán)境，以自然的方式與虛擬三維環(huán)境中的物體或者虛擬人類(lèi)進(jìn)行交互，創(chuàng)造出身臨其境的感官體驗(yàn)。這種技術(shù)很適合作為認(rèn)知科學(xué)、心理學(xué)等領(lǐng)域的研究方法，衍生出更加新穎的實(shí)驗(yàn)方案。利用虛擬現(xiàn)實(shí)技術(shù)仿真技術(shù)，可以實(shí)現(xiàn)視覺(jué)與聽(tīng)覺(jué)上的雙重刺激，可以在大腦情緒狀態(tài)的研究中更加逼近現(xiàn)實(shí)場(chǎng)景的真實(shí)感受[8]。有效地解讀個(gè)體的情緒狀態(tài)能解決工作生活中的各種現(xiàn)實(shí)問(wèn)題，如判定人們工作生活中的精神亞健康狀態(tài)、了解客戶(hù)對(duì)產(chǎn)品的滿意度進(jìn)而調(diào)整營(yíng)銷(xiāo)策略、檢測(cè)車(chē)輛行駛過(guò)程中司機(jī)的精神狀況，在教育領(lǐng)域中可以通過(guò)解讀學(xué)生的學(xué)習(xí)狀態(tài)來(lái)調(diào)整學(xué)習(xí)內(nèi)容和方式等[9]。

本研究引入半球不對(duì)稱(chēng)性指數(shù)作為研究大腦情緒狀態(tài)的工具。實(shí)驗(yàn)使用VR中的3D影片刺激誘發(fā)積極、中性以及消極情緒，同時(shí)記錄腦電信號(hào)，從功率譜腦地形圖上觀察各腦區(qū)，在此基礎(chǔ)上重點(diǎn)分析差異較大的腦區(qū)，計(jì)算各頻段平均功率譜密度和半球不對(duì)稱(chēng)指數(shù)，從中探尋主要腦區(qū)的左右半球不同情緒下的差異和情緒生理活動(dòng)的顯著特征。

1 材料和方法

1.1 實(shí)驗(yàn)設(shè)計(jì)

由10/20系統(tǒng)標(biāo)準(zhǔn)64導(dǎo)腦電帽和德國(guó)BrainProducts公司的腦電圖儀一同采集腦電數(shù)據(jù)，使用Recorder數(shù)據(jù)采集軟件、BrainAmp信號(hào)放大器與隔離電源。VR使用的是美國(guó)Facebook公司輕量級(jí)VR頭顯一體機(jī)Oculus Go。實(shí)驗(yàn)中選擇雙側(cè)乳突作為參考電極，水平眼電電極放置于外眼角處，垂直眼電電極放置于左眼下約2 cm處，前額Afz電極接地，采樣頻率為500 Hz，各導(dǎo)聯(lián)阻抗均小于5 kΩ[10]。實(shí)驗(yàn)材料刺激采用的是YouTube上的3D影片以及網(wǎng)絡(luò)上一些VR團(tuán)隊(duì)所制作的影片，一共12個(gè)影片片段，如圖1所示，其中2個(gè)喜劇片和2個(gè)愛(ài)情片、4個(gè)風(fēng)景紀(jì)錄片、2個(gè)驚悚片和2個(gè)災(zāi)難片，分別對(duì)應(yīng)積極情緒、中性情緒、消極情緒，每個(gè)影片時(shí)長(zhǎng)4～6 min，12個(gè)片段依次播放，每觀看一個(gè)影片充分留出休息時(shí)間讓受試平復(fù)心情再繼續(xù)。圖2展示了受試在實(shí)驗(yàn)時(shí)的狀態(tài)。

圖1 部分刺激材料展示Fig.1 Display of some stimulating materials

圖2 受試實(shí)驗(yàn)中Fig.2 The subject in the experiment

1.2 受試者

實(shí)驗(yàn)選取了20名受試者，男受試10名，女受試10名，均來(lái)自本校在校大學(xué)生與研究生，年齡范圍在23～25歲，平均年齡23.9歲，身體健康，裸眼視力或矯正后視力正常，沒(méi)有色盲色弱缺陷，沒(méi)有心理疾病或精神疾病史，參加實(shí)驗(yàn)前睡眠良好，實(shí)驗(yàn)時(shí)精神狀態(tài)良好。在進(jìn)行腦電數(shù)據(jù)采集之前，讓受試了解實(shí)驗(yàn)無(wú)害性，告知實(shí)驗(yàn)流程及規(guī)則，并簽署實(shí)驗(yàn)知情同意書(shū)。

1.3 信號(hào)預(yù)處理

利用Matlab平臺(tái)中的EEGLAB工具箱進(jìn)行EEG預(yù)處理。順序?yàn)殡姌O通道定位，去除IO眼電電極。以雙側(cè)乳突(TP9、TP10)作為參考電極、進(jìn)行1～49 Hz帶通濾波去除工頻干擾，采樣頻率500 Hz保持不變，然后對(duì)其進(jìn)行獨(dú)立成分分析(independent component analysis, ICA)并使用AJUST插件輔助去除腦電信號(hào)中肌電、眼電等偽跡信號(hào)，得到干凈的腦電信號(hào)。

2 腦電分析方法

2.1 功率譜密度和半球不對(duì)稱(chēng)性指數(shù)

特征提取是對(duì)腦電信號(hào)進(jìn)行變換，找出具有代表性的特征。常用的腦電特征提取主要有3個(gè)方向：時(shí)域特征、頻域特征和時(shí)頻特征，本實(shí)驗(yàn)所采用的分析方法是功率譜估計(jì)分析。在對(duì)EEG信號(hào)進(jìn)行頻域分析時(shí)，功率譜估計(jì)是一種很常用的方法，可以將EEG信號(hào)轉(zhuǎn)換成功率隨頻率變化的譜圖，更加直觀的從功率譜圖中分析EEG信號(hào)不同波段的功率分布與變化情況[11]。本研究中采用的是Matlab中BrainStorm工具箱，使用Welch法對(duì)EEG信號(hào)中的各頻段進(jìn)行功率譜估計(jì)。

首先，分析大腦各腦區(qū)在積極情緒、中性情緒、消極情緒中的功率譜密度分布情況，64導(dǎo)腦電帽導(dǎo)聯(lián)分布如圖3所示。進(jìn)一步地，選取幾個(gè)腦區(qū)的相關(guān)導(dǎo)聯(lián)作為主要研究對(duì)象，這些腦區(qū)包括額區(qū)(F3-F4)、中央?yún)^(qū)(C3-C4)、頂區(qū)(P3-P4)以及枕區(qū)(O1-O2)，利用BrainStorm工具箱將功率譜數(shù)據(jù)導(dǎo)出，然后計(jì)算半球不對(duì)稱(chēng)性指數(shù)。

圖3 64導(dǎo)腦電帽電極分布Fig.3 Distribution of 64-lead EEG cap electrode

半球不對(duì)稱(chēng)性指數(shù)是計(jì)算同一頻率下的左右半球功率譜密度之間的相對(duì)差，差異表征大腦的半球不對(duì)稱(chēng)性，差異的大小表征偏大腦左半球或者右半球的不對(duì)稱(chēng)性程度[12]。本研究在alpha(8～13 Hz)、beta(14～30 Hz)、gamma(31～49 Hz)頻帶中分別計(jì)算不對(duì)稱(chēng)指數(shù)。首先，在每個(gè)EEG頻帶中，計(jì)算20名受試的某一腦區(qū)左側(cè)導(dǎo)聯(lián)平均功率譜密度(WLmn)和右側(cè)導(dǎo)聯(lián)平均功率譜密度(WRmn)。然后計(jì)算每個(gè)EEG頻段，通道和對(duì)象的半球間EEG不對(duì)稱(chēng)性指數(shù)為

(1)

式中，f表示某一頻段，m和n表示某一腦區(qū)左右半球?qū)?yīng)導(dǎo)聯(lián)。半球不對(duì)稱(chēng)性指數(shù)A越高則表示相對(duì)左半球越活躍，由此對(duì)不同腦區(qū)的兩個(gè)通道進(jìn)行不對(duì)稱(chēng)分析。

2.2 統(tǒng)計(jì)學(xué)處理方法

對(duì)積極、中性、消極情緒下的半球不對(duì)稱(chēng)性指數(shù)運(yùn)用重復(fù)測(cè)量方差分析，采用Greenhouse-Geisser法校正，獲得3種情緒下半球不對(duì)稱(chēng)性指數(shù)上的差異，有

(2)

(3)

事后成對(duì)比較采用Bonferroni法校正，目的在于調(diào)整兩兩比較時(shí)的檢驗(yàn)水準(zhǔn)，有

(4)

式中，α表示之前的檢驗(yàn)水準(zhǔn)，m為比較次數(shù)。

3 結(jié)果

3.1 功率譜密度

對(duì)所采集的不同情緒的腦電信號(hào)進(jìn)行信號(hào)預(yù)處理，再通過(guò)EEGLAB將20名受試的3種情緒的功率譜密度根據(jù)不同情緒和不同波段分別疊加平均后，繪制alpha、beta和gamma頻段上的大腦地形圖，見(jiàn)圖4。將3種情緒的功率譜密度進(jìn)行對(duì)比，結(jié)果發(fā)現(xiàn)在各情緒下，積極情緒的額區(qū)的激活程度高于另兩種情緒，而在各頻段之間的對(duì)比以及不對(duì)稱(chēng)性上未觀察到有明顯差異。

圖4 不同情緒下各頻段功率譜密度(每列上中下分別為積極情緒、中性情緒和消極情緒)。(a)Alpha頻段；(b)Beta頻段；(c)Gamma頻段Fig.4 The power spectral density of each frequency band under different emotions (positive emotion, neutral emotion and negative emotion in each column). (a)Alpha band; (b) Beta band; (c) Gamma band

3.2 半球不對(duì)稱(chēng)性指數(shù)

圖5給出了在不同情緒下的各腦區(qū)通道之間的alpha、beta和gamma頻段的不對(duì)稱(chēng)性指數(shù)值。總體上看，在3種情緒下額區(qū)的各頻段不對(duì)稱(chēng)性指數(shù)均為正值，頂區(qū)以及枕區(qū)均為分?jǐn)?shù)較高的負(fù)值，表明3種情緒下額區(qū)的激活程度左側(cè)半球大于右側(cè)半球，頂區(qū)、枕區(qū)則相反。中央?yún)^(qū)的不對(duì)稱(chēng)性指數(shù)正負(fù)值均有。

圖5 不同情緒下各頻段不對(duì)稱(chēng)性指數(shù)。(a)積極情緒；(b)中性情緒；(c)消極情緒Fig.5 Asymmetry index of different frequency bands under different emotions. (a) Positive emotions; (b) Neutral emotions; (c) Negative emotions

進(jìn)一步集中觀察并對(duì)比積極情緒與消極情緒的不對(duì)稱(chēng)性指數(shù)，各頻段下的額區(qū)、頂區(qū)和枕區(qū)結(jié)果如下：

1)額區(qū)(F3-F4)的alpha、 beta和gamma頻段，積極情緒的指數(shù)值在0.1值上下，而消極情緒的指數(shù)值接近于0值；

2)頂區(qū)(P3-P4)的alpha、 beta和gamma頻段頻段在兩種情緒下均為負(fù)值(右半球更活躍)；

3)枕區(qū)(O1-O2)的3個(gè)頻段在兩種情緒下同樣均為負(fù)值(右半球更活躍)。

通過(guò)3種情緒下不對(duì)稱(chēng)指數(shù)的重復(fù)測(cè)量方差分析，統(tǒng)計(jì)結(jié)果非球形性時(shí)采用Greenhouse-Geisser法校正P值，k為3，n為20。結(jié)果發(fā)現(xiàn)，3種情緒下額區(qū)alpha頻段的半球不對(duì)稱(chēng)性指數(shù)的差異邊緣顯著，F(xiàn)(2,38)=2.879，P=0.069。事后成對(duì)比較采用Bonferroni法校正，α取0.05，m為3。發(fā)現(xiàn)額區(qū)alpha頻段的積極情緒與消極情緒兩組間有具顯著差異，P<0.05；積極與中性差異不顯著，P>0.05；消極與中性差異不顯著，P>0.05。而在3種情緒下額區(qū)的beta和gamma頻段，以及各頻段的中央?yún)^(qū)、頂區(qū)和枕區(qū)，均未發(fā)現(xiàn)顯著差異。

上述統(tǒng)計(jì)結(jié)果表明，研究只集中討論alpha頻段的3種情緒下的不對(duì)稱(chēng)指數(shù)即可。圖6顯示了20名受試在額區(qū)(F3-F4)的alpha頻段不對(duì)稱(chēng)性指數(shù)的分布情況。結(jié)果顯示3種情緒下的平均額區(qū)不對(duì)稱(chēng)性指數(shù)均為介于0～0.2的正值，說(shuō)明在3種情緒下的半球不對(duì)稱(chēng)性特征—左側(cè)額區(qū)比右側(cè)額區(qū)的活躍程度高。同時(shí)，額區(qū)出現(xiàn)積極(0.130±0.227)、中性(0.058±0.240)、消極(0.006±0.130)依次遞減的趨勢(shì)，箱型圖中的數(shù)值分布也顯示出這一點(diǎn)，表明從積極到中性、再到消極的情緒過(guò)程中，左側(cè)額區(qū)的活躍程度是逐漸降低的。

圖6 額區(qū)(F3-F4)alpha頻段在3種情緒下的不對(duì)稱(chēng)性指數(shù)對(duì)比(統(tǒng)計(jì)學(xué)顯著性表示為*P<0.05)Fig.6 Comparison of asymmetry index of frontal region (F3-F4) alpha band under three emotions(Statistically significant:*P<0.05)

4 討論

本研究聯(lián)合虛擬現(xiàn)實(shí)和腦電技術(shù)，首次將半球不對(duì)稱(chēng)性指數(shù)用于研究不同情緒下大腦功能狀態(tài)的不對(duì)稱(chēng)性。越來(lái)越多的神經(jīng)科學(xué)和認(rèn)知科學(xué)的研究顯示了虛擬現(xiàn)實(shí)和腦電技術(shù)結(jié)合的優(yōu)勢(shì)。臨床上，已經(jīng)有將其作為暴露療法治療心理疾病，主要包括自閉癥、創(chuàng)傷后應(yīng)激障礙、恐懼癥、社交焦慮癥、抑郁癥等疾病的治療。相關(guān)研究顯示該技術(shù)可以有效地減輕負(fù)面情感的癥狀，從而增強(qiáng)情緒調(diào)節(jié)能力[13]。因此，本研究在實(shí)驗(yàn)刺激材料上使用VR影片，更真實(shí)地模擬情感觸發(fā)的場(chǎng)景，然后記錄了20名受試在導(dǎo)致不同情緒的VR影片下的腦電信號(hào)，基于功率譜密度，計(jì)算不同腦區(qū)半球不對(duì)稱(chēng)性指數(shù)，統(tǒng)計(jì)分析后，集中對(duì)額區(qū)的alpha頻段進(jìn)行了深入分析。

從功率譜密度上來(lái)看，積極情緒與消極情緒下在額區(qū)存在差異，即積極情緒的額區(qū)能量高于消極情緒，但是在左右半球的不對(duì)稱(chēng)性上卻未見(jiàn)明顯差異。 進(jìn)一步地對(duì)半球不對(duì)稱(chēng)性指數(shù)的統(tǒng)計(jì)分析得出了額區(qū)alpha頻段激活程度為左側(cè)半球大于右側(cè)，從積極到消極情緒緩慢下降。體現(xiàn)了額區(qū)在表達(dá)情緒的過(guò)程中扮演著關(guān)鍵的角色。Bartolic等[14]在情緒狀態(tài)對(duì)額葉認(rèn)知任務(wù)影響的實(shí)驗(yàn)結(jié)論，也顯示出額葉對(duì)情緒表達(dá)和情緒調(diào)節(jié)的重要性。在3種情緒下，額區(qū)左半球的激活程度均高于右半球，可以認(rèn)為左半球在處理情緒方面占主導(dǎo)地位，在表達(dá)的感知中具有優(yōu)勢(shì)。

本研究的結(jié)果也提示額區(qū)可以作為情緒研究的重要腦區(qū)進(jìn)行研究。額區(qū)作為人類(lèi)大腦認(rèn)知功能最重要的腦區(qū)之一，其中背扣帶回、內(nèi)側(cè)前額葉皮質(zhì)參與情緒的產(chǎn)生[15]。大量的研究表明杏仁核是人體情緒產(chǎn)生的根源[16]。進(jìn)一步地，Davidson在其關(guān)于額葉不對(duì)稱(chēng)性研究的綜述中指出， alpha頻段的不對(duì)稱(chēng)性與背外側(cè)前額葉皮質(zhì)有關(guān)[17]。該作者認(rèn)為左側(cè)前額葉背外側(cè)的區(qū)域可能影響杏仁核的活動(dòng)，并推測(cè)該區(qū)域負(fù)責(zé)了部分的情緒產(chǎn)生，此外，還發(fā)現(xiàn)左前額葉活性較高的個(gè)體的應(yīng)激激素皮質(zhì)醇水平較低，有力地揭示了具有較高水平的左側(cè)前額葉活動(dòng)的個(gè)體確實(shí)具有更積極的外周生物學(xué)指標(biāo)。本研究的結(jié)果在一定程度上證實(shí)了該作者的上述推測(cè)。

本研究通過(guò)不對(duì)稱(chēng)性指數(shù)，僅發(fā)現(xiàn)額區(qū)在alpha頻段上不同情緒間的差異，不過(guò)文獻(xiàn)顯示情緒的產(chǎn)生和調(diào)節(jié)以及它的不對(duì)稱(chēng)性應(yīng)該不僅與前額葉相關(guān)，還可能與其它腦區(qū)存在重要的關(guān)系，本研究的結(jié)果不夠全面。另外，本研究?jī)H適用了頭皮腦電，未對(duì)皮層上的神經(jīng)元活動(dòng)信號(hào)進(jìn)行分析。今后可以通過(guò)溯源等方法進(jìn)行皮層腦電信號(hào)的分析，從而更精確地定位在不同情緒下激活的其他腦區(qū)。此外，還可以通過(guò)功能連接等方法進(jìn)一步了解左側(cè)前額葉背外側(cè)皮質(zhì)是如何調(diào)控其他腦區(qū)(如杏仁核、海馬體和丘腦)，從而導(dǎo)致在不同情緒下的大腦半球的神經(jīng)活動(dòng)的不對(duì)稱(chēng)性。

最后，本研究仍存在一定局限性，有些實(shí)驗(yàn)材料的選擇可能未能使受試達(dá)到理想的不同情緒，且個(gè)體間的生理差異也會(huì)影響實(shí)驗(yàn)結(jié)果的穩(wěn)定性等。需要擴(kuò)大試驗(yàn)人群數(shù)量及實(shí)驗(yàn)?zāi)挲g范圍，以驗(yàn)證本研究方法的普適性。在未來(lái)，需要在實(shí)驗(yàn)刺激上進(jìn)行更仔細(xì)的篩選，以期望可以更有效的模擬真實(shí)的腦電情緒狀態(tài)。

5 結(jié)論

本研究采用功率譜估計(jì)與半球不對(duì)稱(chēng)性指數(shù)分析了不同情緒下左右半球的整體和區(qū)域特性。在處理情緒的過(guò)程中，左側(cè)額區(qū)的主導(dǎo)程度在積極、中性、消極情緒下依次遞減，提示左側(cè)額區(qū)在人類(lèi)情緒狀態(tài)的表達(dá)中表現(xiàn)活躍，特別是在積極情緒中占據(jù)主導(dǎo)地位。本研究從功率譜的角度揭示了情緒與大腦半球不對(duì)稱(chēng)性的關(guān)系，對(duì)大腦情緒機(jī)制和情緒心理學(xué)的研究提供了一定的參考價(jià)值。