JP6515059B2 - 沈香の品質評価方法 - Google Patents

Description

本発明は、沈香の品質評価方法に関し、機器分析により沈香の品質を評価する。中でも、沈香の中でも最高品質である伽羅を特に見分けるための品質評価方法に関する。

沈香とは、アクイラリア（Ａｑｕｉｌａｒｉａ）属の材、特にその辺材の材質中に黒色の樹脂が沈着した香木であり、灰褐色〜黒褐色を呈し、所々に穴や溝を有するものもある。沈香は樹脂に富む部分は光沢のある黒点を有し、質は固く重い。また、僅かな香気があり、薫べると芳香を発する（非特許文献１）。

アクイラリア属は東南アジアに多く生息し、中でもベトナムやラオス、インドネシアに多く見られる。材質中に樹脂が沈着するためには落雷等の外傷や虫食い等、何らかの衝撃が必要とされるが、その詳細は明らかになっていない。
また沈香は、白檀や松等、樹齢を重ねると必ず心材に油が溜まるものとは異なり、必ずしも樹脂が沈着するものではないことから、非常に貴重な香木資源である。

沈香とひと口に言っても、その樹脂の種類や量は様々であり、その品質は様々である。沈香の中でも最高品質のものは伽羅（きゃら）と呼ばれ、非常に貴重である。
ことに近年は、沈香の限りある資源の貴重さから、財産的価値を見出す場合も多く、品質評価や真贋鑑定のニーズが強い。

一方、沈香は天然物であり、香木の種類、生育条件、樹齢や、樹脂沈着の要因、気象条件等、様々な要因が絡みあって沈香が得られることから、完全に同じ沈香は存在しないと言えるほどである。そのため、伽羅、沈香、人為的な沈香、偽物の沈香等と言った品質評価や真贋鑑定に際しての科学的分析手段による類型化は非常に困難であり、従来は熟練した鑑定者がその品質評価を担ってきた。

日本薬局方外生薬規格２０１５ 第４５頁

鑑定者による品質評価では、見た目（視覚）や、薫べた際の香り（嗅覚）、手触りや重み（触覚）等の五感を用いた官能評価が行われる。特に嗅覚を用いた評価においては、実際に沈香を加熱する必要があり、評価に用いた沈香は燃え尽きてしまう破壊試験となる。

一方で、沈香の濫獲等により資源そのものの数が激減し、沈香、特に伽羅は現在非常に入手が困難である。

本発明は上記実情に鑑みてなされたものであり、破壊試験、かつ熟練した技術が求められる鑑定者による官能評価を補助し、より確実な鑑定結果とするために、科学的分析技術として機器分析による新たな沈香の品質評価方法を提供することを目的とする。また、当該機器分析において、試料の破壊量をなくす（非破壊）又は極小にした新たな沈香の品質評価方法を提供することもさらなる目的とする。

本発明者らは、科学的分析技術における機器分析の一つである赤外分光分析を選択し、該分析により得られたスペクトルデータを比較、解析することによって、複数種の沈香又は偽物の沈香を類型化できることを見出し、本発明を完成するに至った。
すなわち本発明は、上記課題を解決するものであり、下記［１］〜［５］に関するものである。
［１］ 中赤外又は熱赤外の赤外分光分析を用いて、沈香の１８００〜９００ｃｍ^−１の領域における赤外吸収強度又は赤外透過強度を測定することを特徴とする沈香の品質評価方法。
［２］ 前記赤外分光分析が中赤外を用いたＡＴＲ法により行われる、前記［１］に記載の沈香の品質評価方法。
［３］ 中赤外又は熱赤外の赤外分光分析を用いて、さらに沈香の４０００〜１８００ｃｍ^−１の領域における赤外吸収強度又は赤外透過強度を測定し、
１６６０〜１５９０ｃｍ^−１の領域に二股に分かれた２つのピークが存在し、かつ
４０００〜９００ｃｍ^−１の全領域中において、前記二股に分かれた２つのピークのうちピーク強度が高い方のピークよりもピーク強度が高いピークが存在しない、又は、３つ以下存在する、前記［１］又は［２］に記載の沈香の品質評価方法。
［４］ ４０００〜９００ｃｍ^−１の全領域中に１０３０〜１０２８ｃｍ^−１の領域におけるピークよりもピーク強度が高いピークが４つ以上存在する、又は、１０３０〜１０２８ｃｍ^−１の領域にピークが存在しない、前記［３］に記載の沈香の品質評価方法。
［５］ ４０００〜９００ｃｍ^−１の全領域中、１６４０〜１６３０ｃｍ^−１の領域にピーク強度が最も大きなピークが存在する、前記［３］又は［４］に記載の沈香の品質評価方法。

本発明に係る品質評価方法は、沈香の科学的センシングによる評価が可能であるという世界初の知見をもたらすものであり、破壊量が極小又は非破壊で沈香の品質や真贋を評価、判断することができる。さらには、鑑定者による鑑定技術を科学技術による評価手法で補助することにより、鑑定結果に客観性を与え、より確実なものとすることができる。本発明により、今後、沈香を始めとする香木の評価に共通の鑑定プラットフォーム基板技術としての発展が期待できる。

図１は、試料１である伽羅の中赤外を用いたＡＴＲ法によるスペクトルデータである。 図２は、試料２である伽羅の中赤外を用いたＡＴＲ法によるスペクトルデータである。 図３は、試料３である伽羅の中赤外を用いたＡＴＲ法によるスペクトルデータである。 図４は、試料４である沈香の中赤外を用いたＡＴＲ法によるスペクトルデータである。 図５は、試料５である沈香の偽物の中赤外を用いたＡＴＲ法によるスペクトルデータである。 図６は、試料１である伽羅の中赤外を用いたＡＴＲ法によるスペクトルデータの微分曲線を示すグラフである。 図７は、試料４である沈香の中赤外を用いたＡＴＲ法によるスペクトルデータの微分曲線を示すグラフである。 図８は、本発明に係る沈香の品質評価方法を示すフロー図である。 図９は、ＯＭＮＩＣ ｖｅｒ．６．２を用いて赤外吸収スペクトル（反射率）のグラフからピークを判定する手順のうち、閾値及び感度を調整し、ピークを検出するためのソフトウェア画面を示す図である。

以下、本発明を詳細に説明するが、本発明は以下の実施形態に限定されるものではなく、本発明の要旨を逸脱しない範囲において、任意に変形して実施することができる。
本明細書において数値範囲を示す「〜」とは、その前後に記載された数値を下限値及び上限値として含む意味で使用される。

本来、赤外分光分析は、分子の振動や回転の状態を変化させるのに必要なエネルギーを与え（赤外光を照射し）、透過又は反射した光を測定する分析方法である。得られた赤外吸収スペクトルは主として分子の固有振動数に基づくため、分子が異なれば必ずスペクトルも異なる。そのため、物質の構造解析（同定）や定性分析、定量分析等を行うことができる。
一方、沈香は木材に樹脂が沈着した香木であることから、赤外分光分析により得られるスペクトルは木材や樹脂の種類によって様々であり収束しないものと考えられていた。しかしながら本発明では、鋭意検討の結果、沈香の品質によって該スペクトルを類型化することが可能となった。

すなわち、本発明に係る沈香の品質評価方法は、中赤外又は熱赤外の赤外分光分析（ＩＲ）を用いて、沈香の１８００〜９００ｃｍ^−１の領域における赤外吸収強度又は赤外透過強度を測定することを特徴とする。
かかる方法により、試料となる沈香を非破壊で、又は極小の破壊量で沈香の品質を評価することができる。また、機器を用いた分析であることから、評価結果に客観性を与え、一般化することができる。なお、本明細書において、沈香の品質の評価とは、沈香の真贋の評価も含む意味である。

赤外分光分析（ＩＲ）における赤外光は中赤外又は熱赤外であればよいが、本明細書において、中赤外とは波長２．５〜２５μｍの波長域の光を意味し、熱赤外とは波長２５μｍ〜１ｍｍの波長域の光を意味する。中でも分子の振動の基本周波数を計測する観点から中赤外を用いたＩＲが好ましい。
また、赤外分光分析の中でも、フーリエ変換赤外分光法（ＦＴ−ＩＲ）が好ましい。

赤外分光分析において任意の測定法を用いることができ、透過法でも反射法でもよいが、試料を非破壊で測定できることからＡＴＲ（全反射）法がより好ましい。

本発明においては、沈香の１８００〜９００ｃｍ^−１の領域における赤外吸収スペクトルから沈香の品質を評価できるものであるが、より詳細なスペクトルの解析を行うために、測定領域を４０００〜９００ｃｍ^−１とすることが好ましく、４０００〜６００ｃｍ^−１とすることがより好ましい。

以下、沈香の赤外分光分析により得られるスペクトルについて、品質評価方法を具体的に述べる。
なお、本明細書における赤外吸収スペクトルは、ＴｈｅｒｍｏＦｉｓｈｅｒ ＳＣＩＥＮＴＩＦＩＣ社製のＮｉｃｏｌｅｔ ８７００ＦＴ−ＩＲのＦＴ−ＩＲ装置とＳｅｎｓＩＲ社製のＤｕｒａＳｃｏｐｅのＡＴＲアクセサリを用い、ＡＴＲ法により下記の条件で測定した場合のものである。得られるピークトップの波数の誤差は約±１ｃｍ^−１である。

（試料調製）
試料となる沈香から、厚さ０．５ｍｍ程度、サイズ２〜３ｍｍ四方の破片を切り出した。
（測定条件）
・測定法：ＡＴＲ法
・光照射面積：直径１ｍｍ
・分解能：４ｃｍ^−１
・スキャン回数：６４回
なお、本明細書内で示す赤外吸収スペクトル（反射率）のグラフにおいて、反射率の値が急峻に下がっているところをピークと呼ぶ。ピークとして判定するか否かの詳細な判定方法については後述するが、反射率の微分と演算式とから判定される。スペクトルにおいてピークが見られるということは、その波数域で物質の吸収が見られることを表す。
また、本明細書においてはピーク強度の補正は行っていないが、ＡＴＲ法においては、ＡＴＲ結晶界面のもぐり込み深さが波長によって異なる。波数が低いほど（波長が長いほど）当該もぐり込み深さが大きくなるので、必要に応じて、ピーク強度が増大する傾向があることを勘案してスペクトルの解析を行う。

顕微−透過法で赤外分光分析を行う際の測定条件は上記と同様であり、上記とほぼ同様のスペクトルが得られる。

得られた赤外吸収スペクトル（反射率）のグラフの微分スペクトルと演算式から、ピークとするか否かの判定を行う。
まず、赤外吸収スペクトル（反射率）の微分スペクトルにおいて、長波数側から見た際に正から負の方向にゼロをまたぐものをピーク候補として検出する。なお、微分スペクトルの例として、伽羅の赤外吸収スペクトルを示す図１の微分スペクトルを図６に、沈香の赤外吸収スペクトルを示す図４の微分スペクトルを図７にそれぞれ示す。

ピークの判定はＯＭＮＩＣ ｖｅｒ．６．２を用いて、得られた赤外吸収スペクトル（反射率）のグラフのデータを読み込み、次いで、ピーク検出機能を選択し、図９に示したソフトウェア上の画面において閾値及び感度を調整し、閾値設定ラインと感度設定を行うことにより、上記で検出されたピーク候補からピークとして判定するものを自動検出する。なお、上記ソフトウェアを用いたピークの自動検出は、下記演算式を用いて判定されるものである。
Ａｂｓ Ｓｅｎｓ＝０．５＊（Ｍａｘ−ｍｉｎ）＊ｅｘｐ（−０．０１＊ｓｅｎｓｉｔｉｖｉｔｙ）
上記演算式中、「ｓｅｎｓｉｔｉｖｉｔｙ」は０〜１００の範囲で任意で設定できるものであるが、本明細書においては、試料間で差異が出ないように７０に固定する。
なお、得られたスペクトルによっては、スペクトル分析おいて、手動でピーク検出を行うなど、適宜調整を行ってもよい。

沈香の真贋は、４０００〜９００ｃｍ^−１の領域における赤外吸収スペクトルから判断することができる。具体的には、赤外吸収スペクトルにおいて下記（ａ）及び（ｂ）に該当すれば、本物の沈香であると判断することができる。
（ａ）１６６０〜１５９０ｃｍ^−１の領域に二股に分かれた２つのピークが存在する。
（ｂ）前記（ａ）における二股に分かれた２つのピークのうちピーク強度が高い方のピークよりもピーク強度が高いピークが存在しない、又は、３つ以下存在する。

例えば伽羅の赤外吸収スペクトルを示す図１においては、４０００〜９００ｃｍ^−１の領域中、１６６０〜１５９０ｃｍ^−１の領域において、１６３３．６５ｃｍ^−１と１６０２．４５ｃｍ^−１に二股に分かれた２つのピークが存在し、１６３３．６５ｃｍ^−１のピークよりもピーク強度が高いピークは存在しない。
また、沈香の赤外吸収スペクトルを示す図４においては、４０００〜９００ｃｍ^−１の領域中、１６６０〜１５９０ｃｍ^−１の領域において、１６５６．０４ｃｍ^−１と１５９８．５６ｃｍ^−１に二股に分かれた２つのピークが存在し、１６５６．０４ｃｍ^−１のピークよりもピーク強度が高いピークは、３３５２．３９ｃｍ^−１、２９２３．２９ｃｍ^−１及び１０３０．３１ｃｍ^−１に３つ存在する。

一方、沈香の偽物の赤外吸収スペクトルを示す図５においては、４０００〜９００ｃｍ^−１の領域中、１６６０〜１５９０ｃｍ^−１の領域において、比較的ブロードではあるものの、１６３４．８０ｃｍ^−１と１５９７．３４ｃｍ^−１に二股に分かれた２つのピークが存在し、１６３４．８０ｃｍ^−１のピークよりもピーク強度が高いピークは、３３２９．２６ｃｍ^−１、２９５１．８４ｃｍ^−１、２９２１．９９ｃｍ^−１、２８５２．７３ｃｍ^−１及び１０２９．６８ｃｍ^−１に５つ存在する。

なお、本物の沈香であっても、沈香の種類によっては前記（ａ）及び／又は（ｂ）に該当しない場合がある。
また、偽物の沈香のうち、沈香となる木を用いて、そこに人工的に樹脂を注入し、さらに金属等を挿入して重みを持たせたものは、木そのものが沈香となる木であることから、赤外吸収スペクトルが本物の沈香のスペクトルと似る場合がある。その場合には、金属の有無を検出する等、他の方法により真贋を判定することが可能である。

前記（ａ）について、１６６０〜１５９０ｃｍ^−１の領域に二股に分かれた２つのピークがブロードとなっており、明確にピークと認識できない場合がある。この場合、沈香は偽物であることが多い。
また、１６６０〜１５９０ｃｍ^−１の領域における２つのピークの反射率（ピーク強度）の差が１％以上で、かつ、高波数側のピークの方が大きい場合に伽羅であることが多い傾向がある。一方、偽物の沈香の場合は、当該２つのピークの反射率（ピーク強度）の差が１％未満程度であることが多い。例えば図５に示す偽物の沈香の場合、１６３４．８０ｃｍ^−１における反射率は９１．５％であり、１５９７．３４ｃｍ^−１における反射率は９１．７％であることから、その差は約０．２％と、１％未満である。

前記（ａ）及び（ｂ）に該当する沈香に対し、赤外吸収スペクトルにおいて下記（ｃ）及び／又は（ｄ）により、その品質を評価することができる。なお、下記（ｃ）及び（ｄ）の少なくともいずれか一方に該当すれば伽羅である可能性が高い。なお参考として、図８に、特徴（ａ）〜（ｄ）を用いた本発明に係る沈香の評価方法を示すフロー図を示す。
（ｃ）４０００〜９００ｃｍ^−１の全領域中に１０３０〜１０２８ｃｍ^−１の領域におけるピークよりもピーク強度が高いピークが４つ以上存在する、又は、１０３０〜１０２８ｃｍ^−１の領域にピークが存在しない。
（ｄ）４０００〜９００ｃｍ^−１の全領域中、１６４０〜１６３０ｃｍ^−１の領域にピーク強度が最も大きなピークが存在する。

前記（ｃ）に関し、伽羅である図１、３の赤外吸収スペクトルでは、それぞれ１０２９．９４ｃｍ^−１、１０２９．８３ｃｍ^−１のピークよりもピーク強度が高いピークは４つ以上存在する。また、伽羅である図２の赤外吸収スペクトルでは、１０３０〜１０２８ｃｍ^−１の領域にピークが存在しない。
一方、沈香である図４の赤外吸収スペクトルでは、１０３０．３１ｃｍ^−１のピーク（誤差として１０３０〜１０２８ｃｍ^−１の領域のピークと判断）よりもピーク強度が高いピークは２９２３．２９ｃｍ^−１の１つのみである。

前記（ｄ）に関し、伽羅である図１〜３の赤外吸収スペクトルでは、それぞれ１６３３．６５ｃｍ^−１、１６３６．４２ｃｍ^−１、及び１６３４．９６ｃｍ^−１に最も大きなピークが存在する。一方、沈香である図４の赤外吸収スペクトルでは２９２３．２９ｃｍ^−１におけるピークが最も大きい。

また、上記（ｃ）、（ｄ）の他に、伽羅が示す傾向として、下記（ｅ）、（ｆ）を挙げることができる。
（ｅ）前記（ａ）における二股に分かれた２つのピークが１６４０〜１６３０ｃｍ^−１の領域と、１６０５〜１６００ｃｍ^−１の領域に１つずつ存在する。
（ｆ）１５００〜１０００ｃｍ^−１の領域において存在するピークが、シャープかつビジーに存在する。

上記（ｅ）について、沈香であって伽羅ではない図４の赤外吸収スペクトルでは、１６５６．０４ｃｍ^−１と１５９８．５６ｃｍ^−１にピークが存在しており、伽羅である図１〜３の赤外吸収スペクトルに比べて、二股に分かれた２つのピークが離れて存在している。

上記（ｆ）について、ピークがシャープであること、及びビジーに（密に）存在することの定義を定めることは難しいが、複数の試料を相対的に見た際に、赤外吸収スペクトルのグラフの微分スペクトルにおける縦軸の絶対値（振れ幅）が大きいほどシャープであると言え、特定波数領域中に検出されるピークの数が多いほどビジーであると言える。
例えば、ＡＴＲ法によるスペクトルデータとその微分曲線を見比べると、伽羅（図１〜３及び図６参照）に比べて沈香（図４及び図７参照）は、１５００〜１０００ｃｍ^−１の領域において検出されているピークがブロードになり、かつ密度も疎となっていることが分かる。

以下、実施例により本発明を具体的に説明する。ただし、本発明はこれらの実施例に限定されるものではない。
＜機器分析による品質評価：赤外分光分析−ＡＴＲ法＞
試料１〜５について、ＴｈｅｒｍｏＦｉｓｈｅｒ ＳＣＩＥＮＴＩＦＩＣ社製のＮｉｃｏｌｅｔ ８７００ＦＴ−ＩＲのＦＴ−ＩＲ装置とＳｅｎｓＩＲ社製のＤｕｒａＳｃｏｐｅのＡＴＲアクセサリを用いて、ＡＴＲ法により赤外吸収スペクトルを測定した。測定条件は下記に示したとおりである。

（試料調製）
試料である香木から、厚さ０．５ｍｍ程度、サイズ２〜３ｍｍ四方の破片を切り出した。なお、香木は色が均一ではなく濃淡があるが、色が濃い方が香木の特徴を有していると一般的にはされていることから、各試料のうち色の濃い部分を選択して試料の調製を行った。なお、香木は５種類用意し、試料１〜５とした。

（測定条件）
・測定法：ＡＴＲ法
・光照射面積：直径１ｍｍ
・測定範囲：４０００〜６５０ｃｍ^−１
・分解能：４ｃｍ^−１
・スキャン回数：６４回
ＴｈｅｒｍｏＦｉｓｈｅｒ ＳＣＩＥＮＴＩＦＩＣ社製のＦＴ−ＩＲ装置のソフトウェア（ＯＭＮＩＣ ｖｅｒ．６．２）を用い、得られたスペクトルデータから、先述した方法に基づいて、赤外吸収スペクトル（反射率）の微分スペクトルと演算式からピークを判定し検出した。なお、得られるピークトップの波数の誤差は±１ｃｍ^−１である。

得られた赤外吸収スペクトルから、下記（ａ）及び（ｂ）に該当するものを沈香又は伽羅とし、該当しないものを偽物の沈香とした。また、（ａ）及び（ｂ）に該当したもののうち下記（ｃ）及び（ｄ）の少なくともいずれか一方に該当したものを伽羅とした。以上により、伽羅、沈香、偽物の沈香の３種類に分類した。

（ａ）１６６０〜１５９０ｃｍ^−１の領域に二股に分かれた２つのピークが存在する。
（ｂ）前記（ａ）における二股に分かれた２つのピークのうちピーク強度が高い方のピークよりもピーク強度が高いピークが存在しない、又は、３つ以下存在する。
（ｃ）４０００〜９００ｃｍ^−１の全領域中に１０３０〜１０２８ｃｍ^−１の領域におけるピークよりもピーク強度が高いピークが４つ以上存在する、又は、１０３０〜１０２８ｃｍ^−１の領域にピークが存在しない。
（ｄ）４０００〜９００ｃｍ^−１の全領域中、１６４０〜１６３０ｃｍ^−１の領域にピーク強度が最も大きなピークが存在する。

＜官能試験による品質評価＞
試料１〜５について、熟練した鑑定者による官能評価により沈香の真贋及び品質評価を行った。
官能評価は、見た目（視覚）や、薫べた際の香り（嗅覚）、手触りや重み（触覚）による総合判断を行った。嗅覚を用いた官能評価方法は以下のとおりである。
香炉に灰を入れ、火をつけた炭を前記灰の上にのせた。次いで前記炭が十分に熾っていることを確認したのち、炭の周りで熱くなった灰（約２００℃）の上に試料の極小片（約０．１ｇ）をのせて加熱した。試料から芳香する香りの種類や拡散性を確認した上で、鑑定者独自の評価基準に則り、伽羅、沈香、偽物の沈香の３種類に分類した。

試料１〜５の赤外分光分析結果の赤外吸収スペクトルを図１〜５に、試料１及び４の赤外吸収スペクトルの微分曲線を図６及び７にそれぞれ示した。また、表１に赤外吸収スペクトル解析の結果及びそれらによる品質評価結果と、官能試験による品質評価結果を示した。

その結果、試料１〜３は前記（ａ）〜（ｄ）のすべての特徴を有し、機器分析及び官能試験の結果、共に伽羅であると判定された。試料４は前記（ａ）及び（ｂ）の特徴を有し、（ｃ）及び（ｄ）のいずれの特徴も有さず、機器分析及び官能試験の結果、共に沈香であると判定された。試料５は前記（ｂ）の特徴を有さず、機器分析及び官能試験の結果、共に偽物の沈香であると判定された。

本発明に係る品質評価方法は、沈香の科学的センシングによる評価が可能であるという世界初の知見をもたらすものであり、破壊量が極小又は非破壊で沈香の品質や真贋を評価、判断することができる。そのため、鑑定者による鑑定技術と共に、当該科学技術による手法を用いることにより、鑑定に客観性を与え、より確実なものとすることができる。
さらには、沈香を始めとする香木の評価に共通の鑑定プラットフォーム基盤技術としての発展が期待でき、香木に関する科学的データや様々な知見を集積することで、鑑定ビジネスを視野に入れた新たなビジネスモデルの構築にもつながる可能性があることから、その技術的意義は極めて大きなものである。

Claims (4)

1. 中赤外又は熱赤外の赤外分光分析を用いて、沈香片の４０００〜９００ｃｍ^−１の領域における赤外吸収強度又は赤外透過強度を測定し、
   １６６０〜１５９０ｃｍ ^−１ の領域に二股に分かれた２つのピークが存在し、かつ４０００〜９００ｃｍ ^−１ の全領域中において、前記二股に分かれた２つのピークのうちピーク強度が高い方のピークよりもピーク強度が高いピークが存在しない、又は、３つ以下存在する場合に、本物の沈香であると判断する沈香の品質評価方法。
2. さらに、４０００〜９００ｃｍ^−１の全領域中に１０３０〜１０２８ｃｍ^−１の領域におけるピークよりもピーク強度が高いピークが４つ以上存在する、又は、１０３０〜１０２８ｃｍ^−１の領域にピークが存在しない場合に、伽羅であると判断する請求項１に記載の沈香の品質評価方法。
3. さらに、４０００〜９００ｃｍ^−１の全領域中、１６４０〜１６３０ｃｍ^−１の領域にピーク強度が最も大きなピークが存在する場合に、伽羅であると判断する請求項１又は２に記載の沈香の品質評価方法。
4. 前記赤外分光分析が中赤外を用いたＡＴＲ法により行われる、請求項１〜３のいずれか１項に記載の沈香の品質評価方法。