日本包鍵学会議ＶｂＬ２ＡｂＪａ９９ａ）

一般論文

フィルム包装ブロッコリーの呼吸モデル
佐藤博実＊石川豊＊平田学＊RespirationModeIforBroccoIiPackagedinPoIymericFilms
ＨｉｒｏｍｉＳＡＴＯｏ，YutakalSHⅨＡＷＡ.，TakashiHⅡﾐＡＴＡ．

Amethcdwasproposedtodevelop"spimtionmodelsoffineshproduce・

SamplesoffreshbmccoUweighmｇ６８－２２７ｇｗｅＩ℃packagedmdifferentpolymerlc
fihnｓ（oxygenpeImeability8920-7640cc／ｍ２・ｄａｙ・ａｔｍａｔ15℃）ａｎｄｓｔｏ妃ｄａｔ15℃

fbrlOdaysSurfacea1℃ａａｎｄｉｍｔｉａｌｖｏｉｄｖｏｌｕｍｅｏｆｆｉｌｍｐａＣｋagewe蝿alsovaIiedto
obtamdiffe1℃ntgasenvimnmentsduIingstorage､Respirationratesofthesampleswe1℃
calculatedfromcontinuouschangesinatmospheIiccompositionwithmpackages、Respimtionmodelswe1℃developedasfunctionsofstomgetiｍｅｓａｎｄｃｏｎｃｅｎｔｍｔｉｏｎｓｏｆ
Ｏ２ａｎｄＣＯ２ｕsingmultiPleI℃grCsslonanalysis・
MOdelsdevelopedsuccessfuUyp1℃dictedchangesmthegascompositionwithinpackagesofbmccoli．
Keywords：FICshprDduce,MApackagingOxygenconsumptionrateCarbondioxide
evolutjion垣te,Respimtionmodel，Bmccoli

青果物のＭＡ包装設計を適切に行うため、ブロッコリーを試料として呼吸モデルを開発した。
ガス透過性の異なる４種類のフィルムを用い、袋表面積と空容積、試料重量を変えてブロッコリーを包装し、１５℃で貯蔵した。各々のガス組成変化の実測値から酸素吸収速度および二酸化炭素生成速度を求め､対応する時間と酸素･二酸化炭素濃度のデータを抽出した。重回帰分析により、任意の包装条件における酸素吸収速度および二酸化炭素生成速度を、それぞれ酸素．二酸化炭素濃度および時間の関数として表すモデルを得た。並行貯蔵した別グループの試料による検証の結果、得られたモデルによって十分な精度で比較的長期間のガス濃度変化の予測が可能なことが判明した。
キーワード：ＭＡ包装、青果物、呼吸モデル、酸素吸収速度、二酸化炭素生成速度、ブロッコリー、農林水産省食品総合研究所（〒305茨城県つくば市観音台2-1-2）：NationalFoodResearchlnstitute,Tsukuba，Ibaraki305－２５－

フィルム包装ブロッコリーの辱』Zぞ天ル

1.緒言

１レートする手法が報告されている。いずれにおいても、鮮度予測を行うためには、呼吸

ガス透過性の高いプラスチック包装材料をを精度よく測定し、呼吸モデルを効率的に構用いて青果物を包装し、鮮度保持に適切な包装内ガス環境を得る技術、すなわちＭＡ包装

築する必要があるが、初期の研究では呼吸商を１と仮定したり2)、二酸化炭素濃度の影響はについては多くの報告がある')。

無いと仮定してs)構築したモデルが多い。

ＭＡ包装は、青果物の呼吸作用と包装材料

著者らは既に、ＭＡ包装したブロッコリーのガス透過性とをバランスよく調整し、鮮度を用い、包装材料のガス透過性、サイズ、ヘ保持に必要な低酸素・高二酸化炭素環境を得ッドスペースなどの初期条件から、袋の体積る手法であるから、対象青果物の特性に応じ変化、酸素吸収速度、二酸化炭素生成速度をた包装設計が必要であることはいうまでもな

簡便に精度よく測定する方法を報告した8)。

い。しかし、現状では経験則に基づいた実証

今回はこの手法を用い、初期条件のみでその試験的アプローチによる包装設計が大勢を占後の呼吸速度変化、ガス組成変化が予測できめており、多大の労力と時間が費やされている呼吸モデルの開発について検討したので報

る。

告する。

一方､ＭＡ包装における青果物の呼吸変化、

袋内ガス組成変化、鮮度変化などをシミュ

2.実験方法レーションにより予測して、適正包装設計を行う試みが検討されてきた。リンゴ2)、トマ

2.1包装フィルム

ト3)`)、モモ5)、ブロッコリー｡、バナナ7)等に供試した包装フィルムの表面積、厚さ、１５ついて、品質変化や袋内ガス組成変化をシミ℃におけるガス透過性をＴａｂＩｅｌに示した。

Ｔａｂｌｅ１．Conditionsofmodifiedatmosphe｢epackaging

No．ｍ

Nb０２ＣＯ２

（ｍ２）

（ｍｌ）

（９）

Fig.１

3155

210

920

27.8±1.0

1480

5400

ＬＤＰＥ

17.5±１．４

2150

7640

1365
1261

2173
1273

856
1008
８５１

｡ｍｅａｎ±Ｓ、．・・measmredat15℃

ＰＢＤ：Polybutadiene，ＯＰＰ：OrientedpolypropyleneoLDPE：Lowdensitypolyethylene．

－２６－

２

４５

30.1±0.2

ＬＤＰＥ

加皿瓢
釦１３
２３

ＯＰＰ

1260

●●●●●■●●●PBD２５．３±0.8１８００６６３０37200

８３６８８６０００

Conditionbroccoli

１１１１１１１１１

Weightof

volume０００００００００

Initialvoid

offilm澱躯麺諏報翠麺明閉（似、）Surfacearea

鰹却却却却却却却却

FnmThiClmeSs。

PermeabiUty.．
(ｍｌ／m2・ｄａｙ・at、）

３１

日本包鍵学会鯵VbL2jVO・Ia99a）

ガス透過性は、Gasperm-100（日本分光）

し、水圧による負の誤差は十分小さいものとで測定した。

して無視した。予備試験によりブロッコリー

2.2ブロッコリーの比重を0.96として、重量から体積を換算し、全容積からブロッコリーの体積を差し引

土浦公設卸売市場で購入した愛知県産のブいたものを袋内の空容積（ヘッドスペース）ロッコリー（みよ緑１号）を試料として用いとした。容積測定後、直ちに15℃、相対湿度

た。ブロッコリーは入手後、直ちに15℃、相90％の恒温恒温室で貯蔵した。空容積、ブロ

対湿度90％の恒温恒温室へ搬入し保管した。

ッコリーの重さ等の条件は各包装区ともに代表的なものをＴａｂＩｅｌに示した。

2.3包装と貯蔵一晩保管した約200～2509のブロッコリー

2.4袋内ガス濃度の測定を各包装フィルムで包装して熱シールで密封

約1ｍlのヘッドスペースガスを包装袋からした。包装袋の全容積は水浸法により測定

抜き取り、ガスクロマトグラフ（島津ＧＣ－

グロ●●■■■■|■D￣●￣￣●￣￣ロー￣￣Ⅱ■￣￣Ｃ■1いやＱ－－Ｃ－

０２４６８１０

雷釦幅扣５０.'．．､漣--..---…．…＝

０２４６Ｂ１０

０２４６８１０

Condition3

パ､．…－．－－．．－－．＿＿＿＿＿．＿＿＿．＿

５０５０５０
２２１１

蛮如帽扣５０

訳．■◎室何』一色の。■。。⑩⑩。

Ｃｏｎｄｉｔｉｏｎ２

印幅⑩頚如蚕如帽扣５０

露、帽、５０Ｃｏｎｄｌｔｉｏｎｌ

０２４６BＴＯ

CondiIion5に■●Ｃの■■￣●‐￣■￣￣￣●￣●｡￣●●ロロ－－

０２４６８１０

Ｓｔｏｒａｇｅｔｉｍｅ，day
FiglChangesingasconcentrationswithinfiImpackagesoffreshb｢occoIi-Oxygen．･･･……Ｃａ｢bondioxide－２７－

フィルム包装ブロッコリーの呼吸そ天ル

3AH）で酸素、二酸化炭素、窒素を経時的に重回帰分析に使ったグループとは別に並行し分析した。分析条件は前報8》の方法によって貯蔵したグループについて、ガス組成変化た。なお、各ガスのTCD検出器に対する相対の予測を行った。前報8)の理論により、包装

感度は、窒素100、酸素94.3、二酸化炭素

袋内における酸素量の変化速度dVo2／dtは、115.9として計算した，)。

塑二＝Ｋ●2.Ａ．(O21-Po2)－Bo2･Ｗ
ｄｔ

２．５計算

同様に二酸化炭素の変化速度dVco2／ｄｔ

計算は全て､PC-9801RA（日本電気）を用いてBASICプログラムにより行った。

は、

2.5.1袋の体積変化酸素吸収速度、二酸化dVco2

炭素生成速度の算出

dｔ

ガスクロマトグラフィーにより得られた袋＝－ＫＣＯ２．Ａ・Pco2＋ＢＣＣ2.Ｗ

内ガス組成分析データを用いて、前報帥の方

ここに、法により計算で求めた。

Ko2：包装袋の酸素透過度（cc／m2・ｄａｙ・ａｔ、）

2.5.2重回帰分析

Kco2:包装袋の二酸化炭素透過度

各袋のガス組成変化の実測値に当てはめた（cc／ｍ２・ｄａｙ・ａｔ、）

最小二乗曲線を0.1日毎に分解して種々の組

Ａ：包装袋の表面積（､?）み合わせの酸素、二酸化炭素濃度を抽出した

Po28包装袋内の酸素分圧（at、）(Fig.１)。そして酸素濃度（[02])、二酸化炭

Pco2:包装袋内の二酸化炭素分圧（at、）素濃度（[CO2])、時間(T）を独立変数とし、Bo2：単位雷蚤一日当りの酸素消費量

呼吸速度（酸素吸収速度、二酸化炭素生成速（cc／day.ｋg）

度）を従属変数とする重回帰分析（変数増加

Bco2:単位重量一日当りの二酸化炭素生成量法）を行った。酸素濃度、二酸化炭素濃度、（cc／day.ｋg）

時間、それぞれの二乗（[ｑ]２，［CO2]2,Ｔ２）、

Ｗ：ブロッコリー重量（kg）

平方根（[02]1忽)、立方根（[02]'/３，［CO2]1忽，Vo2：袋内酸素量（cc）

Ｔｌへ)、逆数（[02]-1)、指数（e｢〔･by)、三者間

Vco2:袋内二酸化炭素量（cc）の積（[02]・[CO2],Ｔ・[Ｏ２１Ｔ・[CO2])、除ｔ：貯蔵期間（day）

算した項（[02]／［CO2]）も独立変数として加えた。これは一次の項のみによる限界、すである。

なわち回帰方程式の非線形性や、－次の項同モデルで得られる酸素消費速度および二酸士の相互作用に対応できるようにするためで化炭素生成速度をBo2およびBco2に代入するある。

と、袋内の酸素量および二酸化炭素量の変化

2.5.3ガス組成変化の予測速度が求まる。これを用いて各ガス成分の袋

重回帰分析によって得られた酸素吸収速度および二酸化炭素生成速度のモデルを使い、内初発量からの増減を算出し、ガス組成変化の予測曲線を得た。

－２８－日本包鋳学会雄ＶＯｊＬ２ｊＶｂＬＺ｡,93）

従来の呼吸モデルの中にも、二酸化炭素濃度

3.結果と考察の影響を考慮していないものが報告されている2)s)。しかし、モモ、ブロッコリー、レタス

3.1ＭＡ包装系内のガス濃度変化
Ｔａｂｌｅｌの各包装条件下における袋内ガスのように、高二酸化炭素濃度で酸素吸収に顕濃度の変化をFig.１に示した。各コンディシ

著に影響を受けるものもあり、、包装条件に

ョン共に測定データに当てはめた最小二乗曲

よっては高二酸化炭素条件下に長期間さらさ

線で図示した。

れることになり、呼吸モデルの開発にあたっ

酸素濃度変化でみると、まず貯蔵開始後のても、二酸化炭素濃度の影響は無視できない

濃度低下の速さに差がみられた。コンディシ
ョン１と３で、フィルムの種類や表面積が共通と考えられる。したがって、青果物の呼吸モ
デル、すなわち酸素吸収速度および二酸化炭であるにもかかわらず、それぞれの最低酸素素生成速度を求める関数を構築する場合に濃度に達するまでの時間に明らかな差がみられたが、その理由は袋の表面積と空容積に対は、酸素濃度や時間の項のみでなく二酸化炭
素濃度の項も必要であると考えられる。

する試料重量の差であると考えられた'０)。ま得られたガス濃度変化は様々な酸素・二酸

た、各コンディションにおける最低酸素濃度

化炭素濃度の組み合わせを含んでおり、モデは、コンディション１を除くと包装フィルム

ル構築に必要な変化曲線として十分であるとの酸素透過度が低いほど小さな値となってい考えられた。

た。

一方、二酸化炭素濃度変化では、コンディ

3.2呼吸モデル構築に用いたガス組成の組み合わせション４において貯蔵２日後に20％を超える

急速な増加がみられ、その後も増加し続け呼吸モデルを構築するためには、酸素と二た。その他のコンディションでは1～２日後

酸化炭素の様々な組み合わせ下における呼吸に見かけ上、平衡に達した。貯蔵３日目以降

速度データを得なければならない。Ｔａｂｌｅｌの平衡二酸化炭素濃度は、フィルムの二酸化に示した条件で貯蔵した結果、酸素・二酸化炭素透過度の低いものほど高くなっていた

炭素濃度の組み合わせは、酸素：0.5～20％、が、酸素濃度と比べると各コンディションの二酸化炭素：２～20％の範囲をカバーするこ差は大きくなかった。

とができた（Fig.２)。図中で、右上側にくる

以上、包装条件を変えることによって、ような組み合わせ（酸素・二酸化炭素ともに様々な酸素・二酸化炭素濃度変化曲線が得ら

高濃度になるような組み合わせ）は、今回のれ、通常のＭＡ包装におけるガス濃度変化を実験系では得られなかった。この理由は、使カバーすることができたと考えられた。

用した各フィルムの二酸化炭素透過性が、いコンディション４にみられるような高二酸

ずれの場合も酸素透過性の数倍の値であり、化炭素状態が、青果物の呼吸生理に与える影

生成した二酸化炭素が袋外へ拡散・放出する響について、ニンジンやキャベツにおいては速さのほうが、酸素が袋内へ流入する速さよその影響を受けにくいという報告がある１１)。

りも大きく、高酸素・高二酸化炭素状態には－２９－

フィルム包装ブロッコリーの序iBvぞ完"

二酸化炭素環境は得られないと考えられる。

２０

したがって、通常のＭＡ包装を想定した場
合には、得られた組み合わせで必要十分であ

●

５

１

●

母集団とした。

●●

０１

3.3酸素吸収速度のモデル●

５

誤苣。｜］口』］Ｅ①⑨臣。○卸。。

ると判断し、これを呼吸モデル構築のための●．｡･･･9J：．｡、．

酸素吸収速度（Rojを従属変数とし、酸素●●●●●●●

五$．『.:.．．．.．． 弓８．

濃度（[O蜜])、二酸化炭素濃度（[CO2])、時間(T）およびそれらの組み合わせを独立変数と０

して、重回帰分析（変数増加法）を行った結

０５１０１５２０

ozconcentration,％

果をＴａｂＩｅ２に示した。第一ステップで最もFig.２oxygen-ca｢bondioxidecombinations
usedforprepa｢inglcspirationmodeIs．

有効な独立変数として選択されたのは［02］であった。独立変数が増えるにしたがって重相関係数（R）は増大したが、５変数以上になるとＲの増加が少なく、またモデルとしてはなりにくいためである。

ＹａｎｇとChinnanl2)は、ガスバリアー性の複雑になりすぎるため、ＡＩＣ(情報量規準）の高いフィルムでトマトを包装し、包装内の雰値からみて、５変数を使った回帰式が妥当で囲気を組成の異なる20通りの混合ガス（02：あると判断した。即ち、酸素濃度（[02])、酸

5～20％,CO2：Ｏ～20％）で置き換えて、呼素および二酸化炭素濃度の立方根（[OJI/３，

吸量を測定した。しかし、ＭＡ包装に使われ[CO2Ｗ)、二酸化炭素濃度の２乗（[CO2]2)、

る包装材料はＴａｂｌｅｌに示した物性に近いも酸素濃度と時間との積（Ｔ･[o2Dの５つであのが多く、、通常の条件では極端な高酸素･高る。酸素濃度の立方根（[OJI/s）および二酸

Table2MuItipIe「egIessionmodeIsfo「「eIationsbetweenO2consumption「ateandseveraIva｢iabIes

Step

Partial1℃g輻ssioncoefficient

ＡＩＣ

Ｒ

R２

0.7756[02］［CO2]'ね［02]'御［CO2]２Ｔ・[02］Const．
１２３４５

5.492

4.514

－３８．２２

1.933－３８．９９25.70-0.309－７７．１４

47.00

-1882

－９５．０３

60.82

28.76

493.150.8807101.57457.700.9386

0.881175.47

452.450.94590.89480.183-

111.33441.410.95690.91570241-0.563

129.02440.130.95930.9202

ＡＩＣ：Anmformationcriterion，Ｒ： Multiplecorrelationcoefficient，Ｒ２：Coefficientofdetermination．
T,Ｔ２,ＴＩ/３，［02]２，［CO21T・［CO2]， [Ｏｊ・[CO2]，［02］／［CO2］werealsodiscussedasmdependent
variableSbutnotselectedinmost suitable5variables．

－３０－日本包裟学会誌VbjL2jVUJd993）

化炭素濃度の立方根（[CO2]1忽）に対する偏呼吸モデルは、線形⑩，)'3)、非線形･)12)を含回帰係数をみてみると、前者は従属変数（酸めてこれまでにいくつか報告されている。こ素吸収速度）に対して正の寄与、後者は負のれらのうち時間を変数に組み入れないも寄与を示しており、低酸素・高二酸化炭素でのs)`)7)'3)は呼吸量の時間依存性が考慮されは呼吸が抑制されるという一般的な考え方'）ていないために数日間の変化についてのみ適との矛盾はない。また、この５変数を使った用できるモデルとなっており、長期貯蔵中に回帰式のＲ２(決定係数または寄与率）の値かおける呼吸変化を予測するのには適さないと

らこの式で９２％は説明できることを表して考えられる。また、酸素濃度のみを変数としおり、必要十分な精度であるといえる。

て酸素吸収量を表したモデルs)'３)では、比較

Fig.３（a）にモデルによる酸素吸収速度の的高二酸化炭素濃度における呼吸変化につい計算値と個々の実測値との関係を示した。吸ては検討されていない。本実験におけるモデ

収速度の大きい領域で多少ばらつきがみられ

ルは､酸素および二酸化炭素濃度､時間の3つた。これは、ＭＡ包装開始直後には袋内の酸のパラメータを含んでいる。このため、低酸素濃度が高いために吸収速度も大きく、急速素・高二酸化炭素濃度下における呼吸変化をな濃度変化とそれに伴う体積変化が起こるた

長期にわたり予測できると考えられる。

め、測定誤差が含まれたものと考えられる。
しかし、相関係数はｒ＝0.956と非常に高く、3.4二酸化炭素生成速度のモデル

計算値と実測値はほぼ一致し、重回帰分析に

Fig.４に酸素吸収速度(Rojと二酸化炭素

よるモデル構築が適切であることを示してい

生成速度（Rco2）との関係を示した。二酸化る。よって、これを酸素吸収速度のモデル、炭素生成速度を酸素吸収速度の一次関数とす`モデル１，とした。

る回帰式においてｒ＝0.869であり、両者の相

０
８

1５０

●●

(b）
０
６

Ｚ．

扣

ｍ●●

０
６

０
２

釦

:緯

〈』二・□畢芯⑨〉
ロ①肩一コ⑨一口。《①苗』こ◎室。－．シの両。。

卯

・・蜜

０
２

１

壱二・国》へ。。）つ⑭一口一コ⑨一口。
⑰肩』Ｅｏ－ａＥ．⑩臣ｏｏｎｏ(a）Ｙ=0.985Ｘ９ｒ=0.956●

Ｙ=0.987Ｘ，ｒ=0.904
００

０３０６０ｇｏ１２０１５q

O2consumptionrate,experimentaI（cc/kghr）

２０４０６０BO

CO2evoIutionratelexPerlmentaＩ（ccⅢ《９．ｍ）

Fig3CorreIogrambetweenexpe｢imentaIandcalcuIatedrespi｢ations．

－３１－

フィルム包装ブロッコリーの呼吸そ天ルて利用可能なことが示唆された。

1００

しかし、酸素・二酸化炭素濃度の組み合わ

Y=22.95+0.415Ｘ
０８ているとは限らない。異なる酸素・二酸化炭

素濃度の組み合わせが同一の酸素吸収速度を６０与える場合も考えられる。そこで酸素吸収速度のモデルと同様に、酸素および二酸化炭素

０４

濃度、時間、およびこれらを組み合わせた値

８０
２』二・ロエ｝。。（⑭］⑩』臣。一一ゴー。シ①側。。

せと酸素吸収速度は､必ずしも１対１に対応しＨけを独立変数、二酸化炭素生成速度を従属変数として重回帰分析を試みた（TabIe3)。変数増加法の第七ステップにおける寄与率はR2＝０

０３０６０９０１２０１５０0.844であり、酸素吸収速度のモデルと比べ

o2consumptionrate,ｃｃ/kghrてやや低い値であった。しかしFig.３（b）にFig、４ＣＯ｢｢elog｢ambetweenoxygenconsumption
rateandcarbondioxideevoIutio、rate．

みられるように、モデルによる計算値と実測値との相関係数はｒ＝0.904と十分高く、さ関はかなり高かった。また、一般に一定の酸らにモデルの独立変数の数やＡＩＣの値からみ素吸収速度以下においては呼吸商(RQ）が高て、二酸化炭素生成速度のモデルとしては十まり、二酸化炭素の生成が酸素の吸収を上回分な精度であると判断した。

る場合が考えられる。本回帰式の傾きは１よ二酸化炭素生成速度を酸素吸収速度の一次

り小さく、約40ｍl／h･k9以下の酸素吸収速関数として表したものを`モデル２，､別個の重度においてはＲＱが１より大きくなることを回帰分析によるものを`モデル３，とし、`モデ示し、二酸化炭素生成速度を表すモデルとし

ル１，と併せてＴａｂＩｅ４に示した。

TabIe3MultiplB「eg『essionmodelsfor「elationsbetwGenCO2evolution「ateandseve｢aIvariables
Step

PartialIEgressioncoefficient

Ｒ慾

ＡIＣＲ

[02]'/３［02]･[CO2］Ｔ･[CO2］ＴＴＩ/３［CO2]２Ｔ・[OjConst．

－８８．９１》画秘恋画趣》－８０．５５０１３０－１．４６４

０００００００0.063－

66.51

０００００００60.39極》》》》”躯

－３７．０５田、麹鍵釦銃、

７．５５

37.14

螺唖躯蝿奴認識-0.822

妬、印閖

-0.373

’一一一

52.30
55.10

0.802
0.142

１．５４
１５６１

37.69

０．９５
0.521
４２７６

29.70

－一一一一一

32.49

11.11

廻録翻麺翻懸36.39

００００００

１２３４５６７

22.44

AIC：Anmformationcritenon,Ｒ：MUltiplecomdationcoefficient，Ｒ２ ：Coefficientofdetermination[ｑ]，［ｑ]２，［02]-'，［02]!/2ｂｅ-[ob]（［CO2]，［CO2]'ねａｎｄ［02]／[CO2］ we妃alsodiscussedasindependent
vaZiables，butnotselectedmmostsuitable7variables．

－３２－

日本包装学会誌ＶｂＬ２ｊＶｂ・Iag93）

TabIe4Respi｢ationmodels
ＭＣｄｅｌ

Equation

ｌ＝O2consumptionrate

MRo2＝-1.882［02］－９５．０３［CO2]】/3＋60.82［α]'/3＋0.241［CO2]2-0.563Ｔ・[on＋129.O2

２＝CO2evO1utionrate

FIRco2＝0.415I6Ro2＋22.95

３＝COzevolutionrate

MRco2＝５５．１［02]'/3-0.366［ｑ]・[CO2］-0.822Ｔ・[02］-1.606Ｔ＋66.51Ｔｌ肉＋0.13［CO2]2-1.464Ｔ・[CO2］-88.91

炭素濃度変化が若干低く予測された。これ

3.5呼吸モデルの検証

モデル作成用の試料と並行して貯蔵した別は、重回帰分析に用いた独立変数の設定が、グループのブロッコリーを用いて、構築した

必ずしも完全ではなかったことによるものと呼吸モデルの検証を行った。Fig.５（a）にそ思われる。しかし、この試料における貯蔵後の中の１条件について、Ｔａｂｌｅ４のモデル1,期の酸素濃度の僅かな漸減傾向は捉えられて2を用い、ブロッコリー包装袋内の酸素･二酸

おり、モデルが適切かつ有効であると判断さ化炭素濃度変化の予測をし、実測値と比較しれた。これは、モデルが時間と二酸化炭素濃た結果を示した。予測曲線は実測値とよく一度の項を含むためであり、一定期間以上の貯致し、またモデル３を用いた場合も同様に実蔵においては両項が不可欠であると考えられ測値と予測値はよく一致した（Fig.５（b))。

た。

モデル３による予測では、貯蔵初期の二酸化

2５

2５

(a）

０

１

５

５

０

０２４６B10

StoragetIme,ｄａｙ

Ｆｉｇ．５

５

０

１

ＣＤヨゴコベ戸一L戸デー己一一５

１

訳（■。一一、』Ｅ①ｃＥｏｏの⑱。

５

０

影…－－…‐

０
２

CO■●ActuaI，…PredIcIed
１

訳（こ◎』』■』Ｅ①。Ｅ３⑪口③０
２０匁OActuaIo-PredIcted(b）

瞳SD-pD-o--丁。＝￣で~=￣

０２４６B10

Storagetime,ｄａｙ

P『edictedandactuaIchangesinoxygenandcarbondioxideconcentrations．
(a）：Ｐ｢ｅｄｉｃｔｅｄｂｙＭｏｄｅｌｌａｎｄ２ｉｎＴａｂＩｅ４(b）：ＰｒｅｄｉｃｔｅｄｂｙＭｏｄｅＩ１ａｎｄ３ｉｎＴａｂｌｅ４－３３－

フィルム包装ブロッコリーの呼吸ぞ天ル

4.結論

４）Ｃ,Ｃ,YangandMSChinnan,Trans・
ASAE・’３０(2)，５４８（1987）

新しく開発した呼吸測定法から得られる

５）Ｋ､Ｒ・ＤｅＵｙａｎｄＳ・SHRizvi,Ｊ､Food

データを基に、酸素濃度、二酸化炭素濃度、ProcessEngineeling,５，２３（1981）

時間を含む呼吸モデルを構築することができ６），.Ｓ,Lee,Ｐ.Ｅ､Haggar,Ｊ・ＬｅｅａｎｄＫた。前述したように、時間項を含まない呼吸

ＬＹａｍ,Ｊ､FoodSci.,５６(6)，１５８０(1991）

モデルでは、７日間以上の予測を行って実測

７）KHayakawa,Ｙ・ＳＨｅｎｉｇａｎｄＳ.Ｇ・値と比較した例は少ない。本モデルは時間のGilbert,lFoodSci.,４０，１８６（1975）項も含むため、９日間以上の比較的長期間に８）石川豊、佐藤博実、石谷孝佑、平田孝、日本

わたって精度よく予測を行うことができ、開発したモデルが適切かつ有用であることを示

包装学会誌、１(2)，１４３（1992）
９）Ｔ､Hirata,Ｔ､NishiyamaHSato,Ｙ、

すことができた。しかし、時間の項を含めた

Ishikawa，Ｔ・ＳｈｉｉｎａａｎｄＴ・Ishitani，L

モデルを構築するためには、実測する期間とPack・ScLTechnoLJapan,２(1)，１５

同じだけの実験期間が必要となり、迅速にモ（1993）

デルを構築することは困難である。促進試験

10）Ａ・CCameron,Ｗ・Ｂｏｙｌａｎ－ｐｅｔｔａｎｄＪ･法等を開発し、実験期間を短縮する必要があ

Lee,Ｊ､FoodSci.，５４(6)，１４１３（1989）

り、検討中である。

11）Y8Kubo,AInabaandRNakamura,Ｊ・
Ａｍｅｒ・ＳＯＣ､Hort､ＳＣＬ，１１５(6)，９７５<文献＞（1990）

１）AAKader,Ｄ､ＺａｇｏｒｙａｎｄＥ.Ｌ・

12）ＣＣＹａｎｇａｎｄＭＳＣｈｉｎｎａｎ,Trans、

Kelbel,ＣＲＣＣｒｉｔ・Rev・FoodScLNutr.，

２８(1)，１（1989）

ASAE.,３１(3),９２０（1988）

13）Ｐ・VeemjuandMKarel,Modem

２）Ｖ・ＪｕｒｉｎａｎｄＭ､Ｋａ1℃1,ＦｏｏｄＴもchno1．，

Packagmg,３９（12)，１６８（1966）

１７(6)，１０４（1963）

３）Ｙ､Ｓ・ＨｅｎｉｇａｎｄＳ.Ｇ,Gilbert,Ｊ・FoodSci.,40,1033（1975）(原稿受付1992年11月１３日）(審査受理1992年12月２４日）－３４－