データの収集が行われない第1の期間と、前記第1の期間の後にデータの収集が行われる第2の期間とを有するパルスシーケンスを実行する磁気共鳴装置であって、前記パルスシーケンスは、以下の条件を満たす、磁気共鳴装置。1前記第1の期間は、隣接するRFパルスによって区切られたn個のサブ期間を含む2前記n個のサブ期間の各々は、1つ以上の勾配パルスを含む3前記n個のサブ期間に含まれている勾配パルスの極性は同じである41つのサブ期間に1つの勾配パルスが含まれている場合、前記1つの勾配パルスの面積はSであり、1つのサブ期間に複数の勾配パルスが含まれている場合、前記複数の勾配パルスの面積の総和はSである5前記n個のサブ期間のうちのp番目のサブ期間は、x個の勾配パルスGpii=1~xを含み、前記n個のサブ期間のうちのq番目のサブ期間は、y個の勾配パルスGqjj=1~yを含み、以下の関係式が成り立つここで、Spi:前記勾配パルスGpiの面積Δtpi:前記p番目のサブ期間の開始時点に印加される第1のRFパルスと、前記勾配パルスGpiとの間の時間Sqj:前記勾配パルスGqjの面積Δtqj:前記q番目のサブ期間の開始時点に印加される第2のRFパルスと、前記勾配パルスGqjとの間の時間前記x個の勾配パルスは、前記第1のRFパルスから第1の時間だけ離れた第1の勾配パルスを含んでおり、前記y個の勾配パルスは、前記第2のRFパルスから第2の時間だけ離れた第2の勾配パルスを含んでおり、前記第1の勾配パルスの面積と、前記第2の勾配パルスの面積は同じであり、前記第1の時間と前記第2の時間は異なっている、請求項1に記載の磁気共鳴装置。前記第1の勾配パルスの高さおよび幅は、前記第2の勾配パルスの高さおよび幅に等しい、請求項2に記載の磁気共鳴装置。前記第1の勾配パルスの高さおよび幅が、前記第2の勾配パルスの高さおよび幅とは異なっている、請求項2に記載の磁気共鳴装置。前記x個の勾配パルスは、前記第1のRFパルスから第1の時間だけ離れた第1の勾配パルスを含んでおり、前記y個の勾配パルスは、前記第2のRFパルスから前記第1の時間だけ離れた第2の勾配パルスを含んでおり、前記第1の勾配パルスの面積と、前記第2の勾配パルスの面積は異なっている、請求項1に記載の磁気共鳴装置。前記p番目のサブ期間には、1個の勾配パルスが含まれ、前記q番目のサブ期間には、複数個の勾配パルスが含まれている、請求項1~5のうちのいずれか一項に記載の磁気共鳴装置。前記p番目のサブ期間には、1個の勾配パルスが含まれ、前記q番目のサブ期間にも、1個の勾配パルスが含まれている、請求項1~5のうちのいずれか一項に記載の磁気共鳴装置。前記p番目のサブ期間には、前記第1のRFパルスと同時に印加される勾配パルスの一部が含まれている、請求項1~5のうちのいずれか一項に記載の磁気共鳴装置。前記q番目のサブ期間には、前記第2のRFパルスと同時に印加される勾配パルスの一部が含まれている、請求項8に記載の磁気共鳴装置。前記条件を満たす勾配パルスは、周波数エンコード方向、又はスライスエンコード方向に印加される、請求項1~10のうちのいずれか一項に記載の磁気共鳴装置。特許文献1の方法では、負の勾配パルスと正の勾配パルスとを用いる必要がある。したがって、3D高速スピンエコー3DFSEのスライスエンコード方向や周波数エンコード方向のように、正の勾配パルスしか使えない場合には適用することができず、血流信号を十分に抑制することができないことがある。したがって、正の勾配パルスしか使えない場合でも、血流信号を十分に抑制できることが望まれている。本発明の一態様は、データの収集が行われない第1の期間と、前記第1の期間の後にデータの収集が行われる第2の期間とを有するパルスシーケンスを実行する磁気共鳴装置であって、前記パルスシーケンスは、以下の条件を満たす、磁気共鳴装置である。1前記第1の期間は、隣接するRFパルスによって区切られたn個のサブ期間を含む2前記n個のサブ期間の各々は、1つ以上の勾配パルスを含む3前記n個のサブ期間に含まれている勾配パルスの極性は同じである41つのサブ期間に1つの勾配パルスが含まれている場合、前記1つの勾配パルスの面積はSであり、1つのサブ期間に複数の勾配パルスが含まれている場合、前記複数の勾配パルスの面積の総和はSである5前記n個のサブ期間のうちのp番目のサブ期間は、x個の勾配パルスGpii=1~xを含み、前記n個のサブ期間のうちのq番目のサブ期間は、y個の勾配パルスGqjj=1~yを含み、以下の関係式が成り立つここで、Spi:前記勾配パルスGpiの面積Δtpi:前記p番目のサブ期間の開始時点に印加される第1のRFパルスと、前記勾配パルスGpiとの間の時間Sqj:前記勾配パルスGqjの面積Δtqj:前記q番目のサブ期間の開始時点に印加される第2のRFパルスと、前記勾配パルスGqjとの間の時間上記の条件を満たすことにより、勾配パルスが同じ極性であっても、血流信号を十分に抑制できることができる。以下、発明を実施するための形態について説明するが、本発明は、以下の形態に限定されることはない。図1は、本発明の一形態の磁気共鳴装置の概略図である。磁気共鳴装置以下、MR装置と呼ぶ。MR:MagneticResonance100は、マグネット2、テーブル3、受信コイル4などを有している。マグネット2は、被検体12が収容されるボア21と、超伝導コイル22と、勾配コイル23と、送信コイル24とを有している。超伝導コイル22は静磁場を印加し、勾配コイル23は勾配パルスを印加し、送信コイル24はRFパルスを送信する。尚、超伝導コイル22の代わりに、永久磁石を用いてもよい。テーブル3は、クレードル3aを有している。クレードル3aは、ボア21内に移動できるように構成されている。クレードル3aによって、被検体12はボア21に搬送される。受信コイル4は、被検体12の頭部から頸部に渡って取り付けられている。受信コイル4は、被検体12からの磁気共鳴信号を受信する。MR装置100は、更に、シーケンサ5、送信器6、勾配磁場電源7、受信器8、中央処理装置9、操作部10、および表示部11を有している。シーケンサ5は、中央処理装置9の制御を受けて、パルスシーケンスの情報を送信器6および勾配磁場電源7に送る。送信器6は、シーケンサ5から送られた情報に基づいて、RFコイル24を駆動する駆動信号を出力する。勾配磁場電源7は、シーケンサ5から送られた情報に基づいて、勾配コイル23を駆動する駆動信号を出力する。受信器8は、受信コイル4で受信された磁気共鳴信号を信号処理し、信号処理により得たれたデータを中央処理装置9に出力する。中央処理装置9は、シーケンサ5および表示部11に必要な情報を伝送したり、受信器8から受け取ったデータに基づいて画像を再構成するなど、MR装置100の各種の動作を実現するように、MR装置100の各部の動作を制御する。中央処理装置9は、例えばコンピュータcomputerによって構成される。操作部10は、オペレータにより操作され、種々の情報を中央処理装置9に入力する。表示部11は種々の情報を表示する。MR装置100は、上記のように構成されている。次に、本形態において血流信号が抑制されたデータを取得するときのパルスシーケンスについて、図2~図11を参照しながら説明する。図2は、本形態におけるパルスシーケンスの一例の説明図である。パルスシーケンスは、3つの期間A~Cを有している。期間Aは、励起パルスP0とRFパルスP1とによって規定される期間である。期間Bは、データの収集が行われない期間であり、RFパルスP1とP4によって規定されている。更に、データの収集が行われない期間Bは、RFパルスP1~P4によって3つのサブ期間B1~B3に区分けされる。期間Cは、被検体のデータの収集が行われるデータ収集期間である。また、期間A~期間Cでは、勾配パルスG0~Gacqが印加される。勾配パルスG0~Gacqは、以下の条件1a~1dを満たすように設定されている。条件1a勾配パルスG0~Gacqは正のパルスである。条件1b勾配パルスG11、G21、G31の面積S11、S21、S31は以下の関係がある。S11=S21=S31=S条件1c勾配パルスG0、Gacqの面積S0、Sacqは以下の関係がある。S0=Sacq=S/2条件1d勾配パルスG11、G21、およびG31は、同じ高さHであり、同じ幅Wである。条件1e勾配パルスG0は、励起パルスP0から時間Δt0だけ離れている。勾配パルスG11は、RFパルスP1から時間Δt11だけ離れている。勾配パルスG21は、RFパルスP2から時間Δt21だけ離れている。勾配パルスG31は、RFパルスP3から時間Δt31だけ離れている。勾配パルスGacqは、RFパルスP4から時間Δtacqだけ離れている。条件1f勾配パルスG11の時間Δt11、勾配パルスG21の時間Δt21、および勾配パルスG31の時間Δt31は、以下の関係を満たしている。Δt11≠Δt21、Δt11≠Δt31、Δt21≠Δt311次に、励起パルスP0の送信時点t0からデータ収集時点tacqまでの間に静止組織が勾配パルスから受けるモーメントの累計と、励起パルスP0の送信時点t0からデータ収集時点tacqまでの間に血流が勾配パルスから受けるモーメントの累計について、順に説明する。1励起パルスP0の送信時点t0からデータ収集時点tacqまでの間に静止組織が勾配パルスから受けるモーメントの累計について先ず、期間Aの間に静止組織が勾配パルスG0から受けるモーメントMSAについて説明する。勾配パルスG0の面積はS0であるので、モーメントMSAは、以下の式で表される。MSA=S0rS2ここで、rS:静止組織の磁化の位置次に、データの収集が行われない期間Bの間の勾配パルスによって静止組織が受けるモーメントについて説明する。データの収集が行われない期間Bは、3つのサブ期間B1~B3を有しているので、各サブ期間の勾配パルスごとに、静止組織が受けるモーメントについて説明する。先ず、サブ期間B1の勾配パルスG11によって静止組織が受けるモーメントにMSB1について説明する。勾配パルスG11の面積はS11であり、静止組織の位置は時間変化しない。したがって、モーメントMSB1は、以下の式で表される。MSB1=S11rS3ここで、rS:静止組織の磁化の位置以下同様に、サブ期間B2の勾配パルスG21によって静止組織が受けるモーメントMSB2、およびサブ期間B3の勾配パルスG31によって静止組織が受けるモーメントMSB3は、以下の式で表される。MSB2=S21rS4MSB3=S31rS5最後に、期間Cacqの勾配パルスGacqによって静止組織が受けるモーメントMSacqについて説明する。勾配パルスGacqの面積はSacqであるので、モーメントMSacqは、以下の式で表される。MSacq=SacqrS6したがって、励起パルスP0の送信時点t0からデータ収集時点tacqまでの間に静止組織が勾配パルスから受けるモーメントの累計MStは、RFパルスP1~P4によってモーメントが反転することを考慮すると、以下の式で表される。MSt=-MSAMSB1-MSB2MSB3-MSacq=-S0rSS11rS-S21rSS31rS-SacqrS=0∵条件1bおよび1c7したがって、データ収集時点tacqにおいて静止組織の位相を揃えることができるので、静止組織を高信号にすることができる。2励起パルスP0の送信時点t0からデータ収集時点tacqまでの間に血流が勾配パルスから受けるモーメントの累計について先ず、期間Aの間に血流が勾配パルスG0から受けるモーメントMBAについて説明する。勾配パルスG0の面積はS0であるので、モーメントMBAは、以下の式で表される。MBA=S0rt018ここで、rt01:勾配パルスG0が印加された時点t01おける血流の磁化の位置血流が一定速度v0で等速直線運動をすると仮定すると、rt01は、以下の式で表される。rt01=r0v0t01=r0v0t0Δt09ここで、r0:血流の初期位置式9を式8に代入すると、期間Aの間に血流が勾配パルスG0から受けるモーメントMBAは、以下の式で表される。MBA=S0rt01=S0{r0v0t0Δt0}10次に、サブ期間B1の勾配パルスG11によって血流が受けるモーメントMBB1について説明する。勾配パルスG11の面積はS11であるので、モーメントMBAは、以下の式で表される。MBB1=S11rt1111ここで、rt11:勾配パルスG11が印加された時点t11おける血流の磁化の位置血流が一定速度v0で等速直線運動をすると仮定すると、rt11は、以下の式で表される。rt11=r0v0t1112式12を式11に代入すると、サブ期間B1の間に血流が勾配パルスG11から受けるモーメントMBB1は、以下の式で表される。MBB1=S11rt11=S11r0v0t1113以下同様に、サブ期間B2の勾配パルスG21によって血流が受けるモーメントMBB2、サブ期間B3の勾配パルスG31によって血流が受けるモーメントMBB3は、以下の式で表される。MBB2=S21r0v0t2114MBB3=S31r0v0t3115また、期間Cacqの勾配パルスGacqによって血流が受けるモーメントMBacqは、以下の式で表される。MBacq=Sacqr0v0t41=Sacq{r0v0tacq-Δtacq}16したがって、血流が受けるモーメントの累計MBtは、式10、13、14、15、16を用いて、以下の式で表される。MBt=-MBAMBB1-MBB2MBB3-MBacq=-S0{r0v0t0Δt0}S11r0v0t11-S21r0v0t21S31r0v0t31-Sacq{r0v0tacq-Δtacq}17S11=S21=S31=S条件1b、およびS0=Sacq=S/2条件1cであるので、式17は、以下の式に簡略化できる。MBt=Sv0-t0/2t11-t21t31-tacq/2Δt0/2-Δtacq/218ここで、Δt0=Δtacqとすると、式18は、以下の式で表される。MBt=Sv0-t0/2t11-t21t31-tacq/219MBt=0になると、等速直線運動をする磁化の信号はデータ収集時点tacqにおいて位相がそろい、血流の信号を抑制することができない。図3に、MBt=0に設定した場合のパルスシーケンスの例を示す。図3では、勾配パルスG11、G21、G31のΔt11、Δt21、およびΔt31が、Δt11=Δt21=Δt31=Δtになるように設定されている。この場合、MBt=0となるので、等速直線運動する血流の信号を十分に抑制することができない。そこで、血流の信号を十分に抑制するために、図2のパルスシーケンスでは、条件1fの式1を満たすように、時間Δt11、Δt21、およびΔt31が設定されている。このように時間Δt11、Δt21、およびΔt31を設定することによって、血流に関するモーメントの累計MBtが0ゼロにならないようにすることが可能となる。したがって、データ収集時点tacqにおいて血流の位相を分散させることができ、血流の信号を抑制することができる。また、図2のパルスシーケンスでは、Δt11、Δt21、およびΔt31は、式1の関係を満たすように設定されているが、式1の代わりに、以下の式1Aを満たすように設定されていてもよい。Δt11=Δt21≠Δt311A図4には、式1Aを満たすようにΔt11、Δt21、およびΔt31を設定した場合の一例が示されている。このように、Δt11=Δt21であっても、Δt21≠Δt31Δt11≠Δt31であれば、血流に関するモーメントの累計MBtが0ゼロにならないようにすることが可能となる。したがって、血流の信号を抑制することができる。尚、式1Aの代わりに、以下の式1Bを満たすように設定してもよいΔt11≠Δt21=Δt311B式1Bの場合、Δt21=Δt31であるが、Δt11≠Δt21Δt11≠Δt31であるので、血流に関するモーメントの累計MBtが0ゼロにならないようにすることが可能となる。また、上記のパルスシーケンスでは、データの収集が行われない期間Bの勾配パルスG11~G31は、同じ高さHに設定され、且つ同じ幅Wに設定されている条件1d。しかし、条件1dは満たしていなくてもよい。図5に、条件1dを満たしていないパルスシーケンスの一例を示す。図5では、勾配パルスG11~G31の高さは互いに異なっており、勾配パルスG11~G31の幅も互いに異なっている。しかし、勾配パルスG11~G31の面積は、図2と同じくSであるので、条件1fの式1が成り立つ。したがって、血流の信号を抑制することができる。尚、上記のパルスシーケンスでは、サブ期間B3には、勾配パルスが一つのみであるが、複数の勾配パルスが含まれていてもよい図6参照。図6は、サブ期間B3に複数の勾配パルスが含まれている一例を示す図である。図6では、図4に示すパルスシーケンスと比較すると、以下の点が異なっている。2aサブ期間B3には、2つの勾配パルスG31およびG32が含まれている。2b勾配パルスG31の面積をS31、勾配パルスG32の面積をS32とすると、S31とS32との間には、以下の関係がある。S=S31S32202c勾配パルスG31は、RFパルスP3から時間Δt31だけ離れている。勾配パルスG32は、RFパルスP3から時間Δt32だけ離れている。2d勾配パルスG21の面積S21および時間Δt21、勾配パルスG31の面積S31および時間Δt31、勾配パルスG32の面積S32および時間Δt32は、以下の式を満たしている。S21Δt21≠S31Δt31S32Δt3221図6に示すように、サブ期間B3に2つの勾配パルスG31およびG32が含まれている場合は、勾配パルスG31およびG32の面積の和をSにすることによって式20参照、静止組織に関するモーメントの累計MStを、MSt=0にすることができる。したがって、静止組織を高信号にすることができる。また、式21を満たすように、勾配パルスG31およびG32の面積および時間を設定することによって、血流に関するモーメントの累計MBtが0ゼロにならないようにすることが可能となる。図6では、サブ期間B2には、勾配パルスが一つのみ含まれているが、複数の勾配パルスが含まれていてもよい図7参照。図7は、サブ期間B2に複数の勾配パルスが含まれている一例を示す図である。図7では、図6に示すパルスシーケンスと比較すると、以下の点が異なっている。3aサブ期間B2には、2つの勾配パルスG21およびG22が含まれている。3b勾配パルスG21の面積をS21、勾配パルスG22の面積をS22とすると、S21とS22との間には、以下の関係がある。S21S22=S223c勾配パルスG21は、RFパルスP2から時間Δt21だけ離れている。勾配パルスG22は、RFパルスP2から時間Δt22だけ離れている。3d勾配パルスG21の面積S21および時間Δt21、勾配パルスG22の面積S22および時間Δt22、勾配パルスG31の面積S31および時間Δt31、勾配パルスG32の面積S32および時間Δt32は、以下の式を満たしている。S21Δt21S22Δt22≠S31Δt31S32Δt3223図7に示すように、両方のサブ期間B2およびB3に2つの勾配パルスが含まれている場合は、勾配パルスG21およびG22の面積の和もSにすることによって式22参照、静止組織に関するモーメントの累計MStを、MSt=0にすることができる。したがって、静止組織を高信号にすることができる。また、式23を満たすように、勾配パルスG21、G22、G31、およびG32の面積および時間を設定することによって、血流に関するモーメントの累計MBtが0ゼロにならないようにすることが可能となる。図7では、サブ期間B2およびB3に2つの勾配パルスが含まれている例について示されている。図7の説明を、サブ期間B2にx個の勾配パルスが含まれており、サブ期間B3にy個の勾配パルスが含まれていれる場合に一般化して考えると、以下のように説明できる。図8は、サブ期間B2にx個の勾配パルスが含まれており、サブ期間B3にy個の勾配パルスが含まれていれる場合の説明図である。図8では、図7に示すパルスシーケンスと比較すると、以下の点が異なっている。4aサブ期間B2には、x個の勾配パルスG2ii=1~xが含まれており、サブ期間B3には、y個の勾配パルスG3jj=1~yが含まれている。4b勾配パルスG2iの面積をS2ii=1~x、勾配パルスG3jの面積をS3jj=1~yとすると、S2iとS3jとの間には、以下の関係がある。4c勾配パルスG2iは、RFパルスP2から時間Δt2iだけ離れておりi=1~x、勾配パルスG3jは、RFパルスP3から時間Δt3jだけ離れているj=1~y。4d勾配パルスG2iの面積S2iおよび時間Δt2i、勾配パルスG3jの面積S3jおよび時間Δt3jは、以下の式を満たしている。図8に示すように、サブ期間B2にx個の勾配パルスが含まれており、サブ期間B3にy個の勾配パルスが含まれていれる場合は、勾配パルスG2iの面積の和をSとし、勾配パルスG3jの面積の和もSにすることによって式24参照、静止組織に関するモーメントの累計MStを、MSt=0にすることができる。したがって、静止組織を高信号にすることができる。また、式25を満たすように、勾配パルスG2i、G3jの面積および時間を設定することによって、血流に関するモーメントの累計MBtが0ゼロにならないようにすることが可能となる。尚、上記のパルスシーケンスでは、データの収集が行われない期間Bが3つのサブ期間B1~B3に区分けされた例について示されている。上記の説明を、データの収集が行われない期間Bがn個のサブ期間B1~Bnに区分けされた場合に一般化して考えると、以下のように説明できる。図9は、データの収集が行われない期間Bがn個のサブ期間B1~Bnに区分けされた場合の説明図である。図9のパルスシーケンスは、以下の条件5a~5cを満たすように設定されている。条件5a1つのサブ期間には、1つの正の勾配パルス、又は複数の正の勾配パルスが含まれている。1つのサブ期間に1つの正の勾配パルスのみが含まれている場合、この正の勾配パルスの面積はSである。また、1つのサブ期間に複数の正の勾配パルスが含まれている場合、複数の正の勾配パルスの面積の総和がSとなる。図9では、説明の便宜上、p番目のサブ期間Bpに含まれている勾配パルスと、q番目のサブ期間Bqに含まれている勾配パルスのみを示しており、その他のサブ期間に含まれている勾配パルスは図示省略されている。p番目のサブ期間Bpには、x個の勾配パルスGpii=1~xが含まれており、q番目のサブ期間Bqには、y個の勾配パルスGqjj=1~yが含まれている。勾配パルスGpiの面積をSpii=1~x、勾配パルスGqjの面積をSqjとするとj=1~y、SpiとSqjとの間には、以下の関係がある。条件5b勾配パルスGpiは、RFパルスPpから時間Δtpiだけ離れておりi=1~x、勾配パルスGqjは、RFパルスPqから時間Δtqjだけ離れているj=1~y。条件5c勾配パルスGpiの面積Spiおよび時間Δtpi、勾配パルスGqjの面積Sqjおよび時間Δtqjは、以下の式を満たしている。図9に示すように、サブ期間Bpにx個の勾配パルスが含まれており、サブ期間Bqにy個の勾配パルスが含まれている場合は、勾配パルスGpiの面積の和をSとし、勾配パルスGqjの面積の和もSにすることによって式26参照、静止組織に関するモーメントの累計MStを、MSt=0にすることができる。したがって、静止組織を高信号にすることができる。また、式27を満たすように、勾配パルスGpi、Gqjの面積および時間を設定することによって、血流に関するモーメントの累計MBtが0ゼロにならないようにすることが可能となる。尚、図9では、RFパルスが送信されていないときに勾配パルスが印加されているが、RFパルスが送信されているときに勾配パルスを印加してもよい図10参照。図10は、RFパルスが送信されているときに勾配パルスGp、Gp1、Gq、Gq1が印加される場合の説明図である。図10の場合は、勾配パルスGpのうち、p番目のサブ期間Bpに含まれている部分を、勾配パルスGp1と考え、勾配パルスGp1のうち、p番目のサブ期間Bpに含まれている部分を、勾配パルスGpxと考える。また、勾配パルスGqのうち、q番目のサブ期間Bqに含まれている部分を、勾配パルスGq1と考え、勾配パルスGq1のうち、q番目のサブ期間Bqに含まれている部分を、勾配パルスGqyと考える。そして、式26および27を満たすようにすればよい。以下に、図9、図10で示されるパルスシーケンスの具体的な例について、図11~図16を参照しながら説明する。図11~図15は、図9のパルスシーケンスの具体的な例であり、図16は、図10に示すパルスシーケンスの具体的な例である。以下に、図11~図16のパルスシーケンスについて、順に説明する。図11では、サブ期間の勾配パルスが以下の条件を満たすように設定されている。条件6aサブ期間B1には、2つの勾配パルスG11およびG12が含まれている。サブ期間B2には、2つの勾配パルスG21およびG22が含まれている。サブ期間B3には、2つの勾配パルスG31およびG32が含まれている。条件6bサブ期間B1~B3の各勾配パルスG11~G32の面積S11~S32は、以下の関係がある。S11=S12=S21=S22=S31=S32=S/2条件6c勾配パルスG11およびG12は、それぞれ、RFパルスP1から時間Δt11およびΔt12だけ離れている。勾配パルスG21およびG22は、それぞれ、RFパルスP2から時間Δt21およびΔt22だけ離れている。勾配パルスG31およびG32は、それぞれ、RFパルスP3から時間Δt31およびΔt32だけ離れている。条件6d各勾配パルスG11~G32の時間Δt11~Δt32は、以下の関係がある。Δt11≠Δt12≠Δt21≠Δt22≠Δt31≠Δt32図11では、各サブ期間の勾配パルスの面積の和は、S11S12=S21S22=S31S32=Sに設定されているので、静止組織に関するモーメントの累計MStを、MSt=0にすることができる。また、図11では、勾配パルスG11~G32の面積は同じ値S/2に設定されている条件6b。この場合は、各勾配パルスG11~G32の時間Δt11~Δt32を、互いに異なる値に設定することによって条件6d、図9の式27を満たすことができる。したがって、血流のモーメントの累計MBtが0ゼロにならないようにすることができる。次に、図12のパルスシーケンスについて説明する。図12では、サブ期間の勾配パルスが以下の条件を満たすように設定されている。条件7aサブ期間B1には、1つの勾配パルスG11が含まれている。サブ期間B2には、2つの勾配パルスG21およびG22が含まれている。サブ期間B3には、2つの勾配パルスG31およびG32が含まれている。条件7b勾配パルスG11~G32の面積S11~S32は、以下の関係がある。S11=S21S22=S31S32=S条件7c勾配パルスG11、G21、G31の面積S11、S21、S31は、以下の関係がある。S11≠S21、S11≠S31、S21≠S31また、勾配パルスG22、G32の面積S22、S32は、以下の関係がある。S22≠S32条件7d勾配パルスG11は、RFパルスP1から時間Δt11だけ離れている。勾配パルスG21およびG22は、それぞれ、RFパルスP2から時間Δt21およびΔt22だけ離れている。勾配パルスG31およびG32は、それぞれ、RFパルスP3から時間Δt31およびΔt32だけ離れている。条件7e各勾配パルスG11~G32の時間Δt11~Δt32は、以下の関係がある。Δt11=Δt21=Δt31=ΔtaΔt22=Δt32=ΔtbΔta≠Δtb図12では、各サブ期間の勾配パルスの面積は、S11=S21S22=S31S32=Sを満たしているので条件7b、静止組織に関するモーメントの累計MStを、MSt=0にすることができる。また、図12では、RFパルスから同じ時間ΔtaΔtbだけ離れている位置に、面積の異なる勾配パルスG11、G21、およびG31G22およびG32が設けられている条件7c、条件7e。このような条件を満たすことによって、図9の式27を満たすことができる。したがって、血流のモーメントの累計MBtが0ゼロにならないようにすることができる。尚、図11は、サブ期間B1~B3に2つの勾配パルスが含まれている例を示しており、図12は、サブ期間B2およびB3に2つの勾配パルスが含まれている例を示している。しかし、サブ期間B1~B3のうちの少なくとも1つのサブ期間に、3つ以上の勾配パルスが含まれるようにしてもよい。図13は、サブ期間B1~B3に3つの勾配パルスが含まれている例を示しており、図14は、サブ期間B2およびB3に3つの勾配パルスが含まれている例を示している。このように、1つのサブ期間に3つ以上の勾配パルスを含めてもよい。次に、図15のパルスシーケンスについて説明する。図15では、サブ期間の勾配パルスが以下の条件を満たすように設定されている。条件8aサブ期間B1には、2つの勾配パルスG11およびG12が含まれている。サブ期間B2には、2つの勾配パルスG21およびG22が含まれている。サブ期間B3には、2つの勾配パルスG31およびG32が含まれている。条件8bサブ期間B1~B3の各勾配パルスG11~G32の面積S11~S32は、以下の関係がある。S11S12=S21S22=S31S32=S条件8c勾配パルスG11、G21、G31の面積S11、S21、S31は、以下の関係がある。S11=S21=S31=Saまた、勾配パルスG12、G22、G32の面積S12、S22、S32は、以下の関係がある。S12=S22=S32=Sb条件8d勾配パルスG11およびG12は、それぞれ、RFパルスP1から時間Δt11およびΔt12だけ離れている。勾配パルスG21およびG22は、それぞれ、RFパルスP2から時間Δt21およびΔt22だけ離れている。勾配パルスG31およびG32は、それぞれ、RFパルスP3から時間Δt31およびΔt32だけ離れている。条件8e各勾配パルスG11~G32の時間Δt11~Δt32は、以下の関係がある。Δt11≠Δt21、Δt11≠Δt31、Δt21≠Δt31Δt12=Δt22=Δt32=Δt図15では、各サブ期間の勾配パルスの面積は、S11S12=S21S22=S31S32=Sを満たしているので条件8b、静止組織に関するモーメントの累計MStを、MSt=0にすることができる。また、図15では、勾配パルスG11、G21、G31は同じ面積Saに設定されている条件7c。この場合は、勾配パルスG11、G21、G31の時間Δt11、Δt21、およびΔt32を、互いに異なる値に設定することによって条件8e、図9の式27を満たすことができる。したがって、血流のモーメントの累計MBtが0ゼロにならないようにすることができる。また、勾配パルスG12、G22、G32を、RFパルスから同じ時間Δtだけずらすことによって条件8e、渦電流がRFパルスに与える影響を均一にすることができる。次に、図16に示すパルスシーケンスについて説明する。図16のパルスシーケンスは、RFパルスと同時に勾配パルスG1~G4が印加されている。図16の場合は、勾配パルスG1のうち、1番目のサブ期間B1に含まれている部分を、勾配パルスG11と考え、勾配パルスG2のうち、1番目のサブ期間B1に含まれている部分を、勾配パルスG13と考える。また、勾配パルスG2のうち、2番目のサブ期間B2に含まれている部分を、勾配パルスG21と考え、勾配パルスG3のうち、2番目のサブ期間B2に含まれている部分を、勾配パルスG24と考える。更に、勾配パルスG3のうち、3番目のサブ期間B3に含まれている部分を、勾配パルスG31と考え、勾配パルスG4のうち、3番目のサブ期間B3に含まれている部分を、勾配パルスG34と考える。そして、静止組織に関するモーメントの累計MStを、MSt=0にし、血流のモーメントの累計MBtが0ゼロにならないように、勾配パルスG11~G34の面積S11~S34および時間Δt11~Δt34を設定すればよい。また、上記のパルスシーケンスで使用されている勾配パルスは、正の勾配パルスである。しかし、正の勾配パルスの代わりに、負の勾配パルスを用いてもよい。尚、上記のパルスシーケンスは、1軸に印加される勾配パルスのみが示されている。次に、3軸に勾配パルスを印加する例について説明する。図17は、3軸に勾配パルスを印加するときのパルスシーケンスの一例の説明図である。図17には、3DFSEのパルスシーケンスの例が示されている。図17のパルスシーケンスでは、スライスエンコード方向および周波数エンコード方向に関しては、正の勾配パルスのみを用いて、静止組織に関するモーメントの累計MStをMSt=0にし、血流のモーメントの累計MBtが0ゼロにならないようにしている。したがって、スライスエンコード方向および周波数エンコード方向のように、負の勾配パルスと正の勾配パルスとの組み合わせが使えない場合でも、血流信号を十分に低減させることができる。