شناسنامه علمی شماره

نویسندگان

1 کارشناس ارشد علوم اقتصادی دانشگاه شهید باهنر کرمان

2 استادیار بخش اقتصاد دانشگاه شهیدباهنر کرمان

3 دانشجوی دکتری مدیریت دانشگاه تهران

چکیده

این مقاله با رویکرد بهینه‌سازی دومرحله‌ای به برآورد تابع تقاضای انرژی در بخش صنعت و محاسبه کشش‌های جانشینی میان نهاده‌ها می‌پردازد. تخمین تقاضای انرژی در بخش‌های مختلف امکان تعیین سیاست‌های کلان انرژی را فراهم می‌کند. در این مطالعه در بخش اول با تخمین تابع تقاضای انرژی در بخش صنعت کشور به عنوان یکی از بخش‌های پرمصرف در زمینه انرژی، عوامل مؤثر بر این تقاضا شناسایی و در مرحله دوم با استفاه از تابع لاجیت و به روش رگرسیون معادلات به ظاهر نامرتبط به شناسایی اجزای تقاضای انرژی پرداخته و کشش‌های قیمتی و جانشینی میان این اجزاء محاسبه می‌شود.

کلیدواژه‌ها

عنوان مقاله [English]

Estimated Energy Demand Function and Substitution Elasticity between Inputs in Industry Section: Seemingly Unrelated Regression Estimation

نویسندگان [English]

  • Saeed Sadzada Mogadam 1
  • Zin al-Abedin Sadeghi 2
  • Ahmad Quds Elahi 3

1 Master of Economics, Shahid Bahonar University of Kerman

2 Assistant Professor of Economics, Shahid Bahonar University of Kerman

3 PhD student in Management, University of Tehran

چکیده [English]

This paper intends to estimate the energy demand function in industry section based on the two-stage optimization approach. Consequently, it defines the substitution elasticity between inputs. Energy demand estimation in various sections makes it feasible to define energy policies. At the first phase the key influential factors on the energy demand function estimation in industry section which is one of high consumption section, were defined and at the next phase the elements of energy demand were determined through logit function and seemingly unrelated regression estimation and price and substitution elasticity were calculated.

کلیدواژه‌ها [English]

  • Energy Demand
  • industry section
  • Seemingly Unrelated Regression Estimation
  • Substitution Elasticity
  • Price Elasticity

تخمین تابع تقاضای انرژی و کشش قیمتی و جانشینی نهاده‌ها در بخش صنعت: رگرسیون معادلات به ظاهر نامرتبط SUR

سعید صدرزاده مقدم*، دکتر زین‌العابدین صادقی** و احمد قدس الهی***

 

تاریخ دریافت: 21 اردیبهشت 1392         تاریخ پذیرش: 9 مهر 1392

 

این مقاله با رویکرد بهینه‌سازی دومرحله‌ای به برآورد تابع تقاضای انرژی در بخش صنعت و محاسبه کشش‌های جانشینی میان نهاده‌ها می‌پردازد. تخمین تقاضای انرژی در بخش‌های مختلف امکان تعیین سیاست‌های کلان انرژی را فراهم می‌کند. در این مطالعه در بخش اول با تخمین تابع تقاضای انرژی در بخش صنعت کشور به عنوان یکی از بخش‌های پرمصرف در زمینه انرژی، عوامل مؤثر بر این تقاضا شناسایی و در مرحله دوم با استفاه از تابع لاجیت و به روش رگرسیون معادلات به ظاهر نامرتبط به شناسایی اجزای تقاضای انرژی پرداخته و کشش‌های قیمتی و جانشینی میان این اجزاء محاسبه می‌شود.

 

واژه‌های کلیدی: تقاضای انرژی، بخش صنعت، رگرسیون به ظاهر نامرتبط، کشش جانشینی، کشش قیمتی.

طبقه‌بندی JEL: D24، C51، C31.

 

1. مقدمه

در دنیای کنونی حرکت به سمت توسعه اقتصادی بدون در نظر گرفتن انرژی غیرممکن به نظر می‌رسد. در کشورهای در حال توسعه، انرژی علاوه بر نقش خود به عنوان عامل اصلی تولید، به عنوان یک منبع درآمد ملی نیز محسوب می‌شود. لذا لزوم حرکت به سمت رشد اقتصادی در کشورهای در حال توسعه از یک سو و نقش اساسی انرژی در این راه از سوی دیگر، نشان‌دهنده اهمیت شناخت تقاضای انرژی در بخش‌های مختلف اقتصادی و به ویژه بخش صنعت می‌باشد؛ چرا که یکی از شاخص‌های توسعه اقتصادی کشورها افزایش سهم صنعت در تولید ناخالص ملی است. اصلی‌ترین مصرف‌کنندگان انرژی در جهان کشورهای پیشرفته هستند، بنابراین با رشد و توسعه اقتصادی و در نتیجه صنعتی شدن میزان تقاضای انرژی افزایش می‌یابد. در ایران مصرف انرژی در بخش صنعت از 7/11 میلیون بشکه معادل نفت خام در سال 1346 به 6/274 میلیون بشکه معادل نفت خام در سال 1389 رسیده است و بطور متوسط سالانه 22% مصرف انرژی را به خود اختصاص می‌دهد. مدل‌های تقاضای انرژی به عنوان ابزارهای تحلیل سیاستگذاران و برنامه‌ریزان جهت تصمیم‌گیری درخصوص آینده انرژی و پیش‌بینی مصرف انرژی طراحی می‌شوند. آگاهی از مصرف انرژی در بخش صنعت می‌تواند باعث بهبود روند تصمیم‌گیری درخصوص عوامل مؤثر در مصرف انرژی در این بخش و همچنین واکنش مناسب در مقابل تغییرات احتمالی این عوامل گردد.

در این مطالعه مصرف انرژی هم به صورت کل و هم به صورت جزء مورد بررسی قرار گرفته است. بدین منظور از یک روش دومرحله‌ای استفاده شده است. در مرحله اول مقدار کل مصرف انرژی و در مرحله دوم مقدار مصرف اجزای انرژی (فرآورده‌های نفتی، برق، ذغال سنگ، گاز طبیعی) مورد بررسی قرار گرفته است. در اکثر مطالعات مشابه از تابع هزینه ترانسلوگ[1] استفاده شده، با توجه به اینکه ویژگی‌های معمول تابع هزینه ترانسلوگ تنها در محدوده‌ای از قیمت‌های نسبی حفظ می‌شود و خارج از این محدوده خاصیت غیرمنفی بودن سهم‌های بازار و شبه مقعر بودن توابع هزینه تأمین نمی‌شود.

در این مقاله ضمن معرفی مدل‌های لاجیت به منظور تخمین از یک مدل لاجیت دومرحله‌ای استفاده شده است. این مقاله در 6 بخش تدوین شده است. در بخش دوم مبانی نظری مدل‌های تقاضای انرژی مورد بررسی قرار گرفته است. در بخش سوم مروری بر مطالعات داخلی و خارجی در حوزه توابع تقاضای انرژی مورد بررسی قرار گرفته است. در بخش چهارم مدل مورد استفاده در این تحقیق که یک مدل بهینه‌سازی دومرحله‌ای است معرفی و معادلات رگرسیون به ظاهر نامرتبط نشان داده شده است. برآورد مدل معرفی شده و نتایج آن در بخش پنجم آورده شده و سرانجام در بخش ششم کشش‌های جانشینی و قیمتی محاسبه شده و همچنین نتیجه‌گیری و پیشنهادات کاربردی ارائه گردیده است.

 

2. مبانی نظری

توابع تقاضای نهاده‌ها را می‌توان از دو روش مشتق‌گیری از تابع سود نسبت به قیمت نهاده‌ها و یا مشتق‌گیری از تابع هزینه نسبت به قیمت هر نهاده استخراج کرد که در روش اول تابع تقاضای مستقیم و در روش دوم توابع تقاضای غیرمستقیم (مشروط) برای نهاده‌ها بدست می‌آید. در اکثر مطالعات انجام شده به منظور استخراج توابع تقاضای نهاده‌ها از روش دوم استفاده شده است. در این روش ابتدا یک تابع تولید انتخاب و تابع هزینه همزاد آن مشخص می‌گردد، سپس با مشتق‌گیری از تابع هزینه نسبت به قیمت هر یک از نهاده‌ها، توابع تقاضا بدست می‌آید.

در بسیاری از مطالعات مربوط به برآورد تقاضای انرژی از توابع هزینه از نوع ترانسلوگ استفاده شده است. نشان داده‌اند که استفاده از مدل تابع هزینه ترانسلوگ با مشکلاتی همراه است که از اعتبار کار انجام شده می‌کاهد. در این توابع دو ویژگی ضروری غیرنزولی بودن و شبه مقعر بودن تابع هزینه به ازای تمام قیمت‌های عوامل ارضاء نمی‌شود. بنابراین نتایج حاصل از آنها نیز تنها در دامنه‌های خاصی معنادار می‌باشد و چنانچه توابع تقاضای نهاده‌ها از آنها استخراج شوند ممکن است ویژگی‌های مثبت بودن سهم هزینه و منفی نیمه معین بودن ماتریس کشش‌های جانشینی را دارا نباشند (کریستنسن و کیوز[2]،1980). به گونه‌ای که در توابع ترانسلوگ تقعر تنها زمانی به صورت مطلق و فراگیر برقرار است که تابع ترانسلوگ به حالت خاص تابع کاب- داگلاس[3] تقلیل یابد.[4] در توابع لاجیت محدودیت‌های توابع ترانسلوگ وجود ندارد و سهم هزینه‌ها در تمامی نقاط مثبت می‌باشد و انعطاف‌پذیری بسیاری وجود دارد.

مدل لاجیت خطی سهم هزینه[5]، یکی از مدل‌های با کشش متغیر با خواص فراگیر است. در تخمین بوسیله مدل لاجیت نیازی به تعیین تابع هزینه نمی‌باشد. با فرض بازدهی ثابت به مقیاس، دلیل تئوریک مبتنی بر نیاز به فرمول تابع هزینه وجود ندارد، چراکه تمامی اطلاعات لازم درخصوص ساختار تابع هزینه توسط  سهم‌های هزینه بدست می‌آید. همچنین به کمک مدل لاجیت ساختارهای هزینه غیرمتجانس را نیز می‌توان تخمین زد. به عبارت دیگر هرگاه فرض شود  تابع هزینه غیرمستقیم است در این صورت  قیمت‌های N عامل تولیدی و Y سطح تولید می‌باشد. اگر  سطوح نهاده‌های مورد استفاده باشند آنگاه رابطه (1) هزینه تولید محصول را با استفاده از نهاده‌های تولید نشان می‌دهد.

(1)                                                          

با استفاده از لم شفارد[6] تابع تقاضای مشروط برای نهاده iام متناسب با سطح مشخص تولید و حداقل نمودن هزینه به صورت زیر تعریف می‌شود.

(2)                                     

در صورتی که تابع تولید نسبت به نهاده‌ها همگن خطی[7] باشد می‌توان تابع هزینه را به صورت زیر بازنویسی کرد:

(3)                           

در این صورت معادلات تقاضای همزمان به صورت زیر می‌باشند:

(4)                                                   

با توجه به اینکه y یکی از پارامترهای معلوم در رابطه (3) می‌باشد، با توجه به همگن خطی بودن تابع تولید، خصوصیات تابع هزینه C به  منتقل می‌شود. همچنین معادلات سهم هزینه به صورت مستقل از متغیر تولید به صورت زیر تعریف می‌شوند:

(5)                                 

به جای اینکه معادلات فوق با تعیین برخی انواع توابع هزینه و مشتق‌گیری نسبت به قیمت‌ها (مانند روش ترانسلوگ) بدست آید، می‌توان برای قیمت‌های وزن داده شده به وسیله مشتق‌های جزیی در معادله (5)، شکلی به صورت تابع در نظر گرفت. به جای اینکه این شکل به صورت تابع خطی باشد از یک تابع نمایی استفاده می‌شود. معادلات سهم هزینه برای یک تابع نمایی به جای شکل خطی به صورت زیر تعریف می‌شوند (کانسیدین[8]، 1989):

(6)                            

که  تابعی از قیمت‌های n عامل تولیدی و  و  ضرایب هستند. کشش‌های سهم هزینه  به صورت زیر تعیین می‌شوند:

(7)                                     

به همین ترتیب کشش‌های قیمتی متقاطع عبارتند از:

(8)                                         

با توجه به شکل نمایی مدل لاجیت خطی، سطح نهاده‌های نمی‌تواند منفی باشد بنابراین برای هر دسته از قیمت‌های مثبت، سطوح داده‌های متناظر همگی مثبت می‌باشند. در نهایت کشش‌های جانشینی برحسب تابع هزینه به صورت زیر محاسبه می‌شوند:

(9)                                                   

سپس با فرض لم شفارد () و تعریف کشش متقاطع قیمتی خواهیم داشت:

(10)                                              

بنابراین در صورتی که  کشش متقاطع قیمتی و  سهم هزینه عامل jام باشد، کشش جانشینی جزیی به صورت زیر تعریف می‌شود.

(11)                                                  

تمامی کشش‌ها با توجه به روابط معرفی شده در مطالعات کانسیدین و مونت (1984 و 1989) می‌باشند.

 

3. پیشینه تحقیق

1-3. مطالعات داخلی

فخرایی (1371) توابع تقاضای انواع حامل‌های انرژی برای بخش‌های مختلف اقتصادی کشور از جمله صنعت را ارائه نمود. در این مدل برای تعیین میزان فروش کل برق به بخش صنعتی، قیمت واقعی فروش برق به بخش صنعتی و تولید ناخالص داخلی در تصریح مدل آورده شده است.

کاظمی (1375) برای تخمین تابع تقاضای بخش صنعت از متغیرهای توضیحی قیمت واقعی برق فروخته شده به بخش صنعت، ارزش افزوده بخش صنعت، قیمت واقعی نفت کوره و نفت گاز و مقدار مصرف برق دوره قبل استفاده نمود که در نهایت تقاضای کل صنعت، تابعی از ارزش افزوده صنعت و مصرف دوره قبل بوده است. مدل ایشان نشان می‌دهد که تقاضا تابعی از قیمت کالاهای جانشین برق (نفت کوره و نفت گاز) نیست. در مدل بخش صنعت که توسط صادقی (1382) ارائه شده است، مقدار مصرف برق تابعی از متغیرهای قیمت واقعی برق، تعداد مشترکین برق در بخش صنعت، ارزش افزوده واقعی بخش صنعت و متغیر مجازی مربوط به دوران جنگ و صلح است. در مدل مذکور تمام ضرایب در سطح اطمینان 5 درصد معنی‌دار می‌باشند.

بهبهانی‌فر (1383) به بررسی عوامل مؤثر بر تقاضای انواع انرژی در بخش صنعت استان اصفهان و تخمین کشش­های جانشینی بین آنها پرداخته است.

قادری (1384) چهار مدل، تابع خطی، تابع نمایی، تابع درجه 2 و تابع به فرم S جهت تخمین تابع مصرف انرژی الکتریکی ارائه شده است. در این مطالعه از روند تابع مصرف انرژی الکتریکی در گذشته، جهت تحلیل و پیش­بینی مصرف برق برای آینده استفاده شده است.

در مقاله صمدی (1387) توابع تقاضای ایستا و پویا برق در ایران طی سالهای 1372-1346 در بخش‌های خانگی، تجاری و عمومی، صنعتی و کشاورزی به روش حداقل مربعات معمولی (OLS) برآورد شده و بدین وسیله کشش­های قیمتی و کشش­های درآمدی تقاضا در کوتاه­مدت و بلندمدت محاسبه شده ­است. ارقام محاسبه شده، نشان­دهنده کم­کشش بودن تقاضای مشترکین نسبت به تغییرات قیمت برق در کوتاه‌مدت و بلندمدت بوده که یکی از علت­های آن را پایین بودن قیمت واقعی برق عنوان کرده است.

دینایی (1378) به بررسی عوامل مؤثر بر بخش خانگی پرداخته است. نتایج حاکی از بی‌اثر بودن قیمت برق روی میزان مصرف در این بخش است بطوری که قیمت‌گذاری چندان اثری بر کاهش مصرف نداشته است.

 

2-3. مطالعات خارجی

برکسون[9] (1944) یک تابع لجستیک برای رابطه بین درصد حشرات مرده و لگاریتم مقدار داروی تجویز شده بکار برد. در این مسئله انتخابی گسسته، یک مدل لاجیت بکار رفته تا احتمالاتی را بیان کند که غیرمنفی و مجموعشان یک می‌باشد. این خصوصیات برای سهم هزینه نیز وجود دارد و استفاده از تابع لاجیت برای سهم هزینه‌ای مخارج امری معقول به نظر می‌رسد. همچنین با رفتار حداقل‌سازی هزینه نیز متناقض نیست.

تایل[10] (1969) مدل لاجیت خطی را برای تقاضای مصرف‌کننده بکار برد. مدل‌های خطی لاجیت که اغلب به مسائل انتخابی گسسته مربوط می‌شوند، مانند انتخاب مدل حمل و نقل مک فادن[11] (1974) ولی می‌توانند برای انواع مسائل کاربردی نیز بکار روند.

درخصوص تخمین تابع تقاضای انرژی به روش دومرحله‌ای مطالعات خارجی توسط ارنست برنت و دیوید وود[12] (1975)، روبرت هالورسن[13] (1977)، ملوین فوس[14] (1977)، روبرت پین دایک[15] (1979)، تی موتی کانسایدین و تی موتی مونت[16] (1984)، تی موتی کانسایدین[17] (1989)، سنفورد برگ و پراکاش لونگانی[18] (1990) انجام شده است. فرض لازم و کافی برای استفاده از تخمین دومرحله‌ای، خاصیت جدایی‌پذیری ضعیف[19] برای توابع تولید و هزینه می‌باشد که در تمامی مطالعات مذکور مورد توجه قرار گرفته است. به منظور  بدست آوردن توابع سهم هزینه[20] بایستی شکل تابع هزینه از نوع توابع انعطاف‌پذیر باشد. عمومی‌ترین شکل تابع هزینه استفاده شده در مطالعات گذشته تابع هزینه ترانسلوگ می‌باشد که توابع سهم هزینه براساس آن محاسبه شده است.

در مقاله ایشی و فرد[21] (2010) تقاضای برق خانگی و صنعتی فیلیپین مدل‌سازی شده است. در این مطالعه عنوان شده است که در بخش صنعت در یک دوره طولانی‌مدت، یک رابطه بلندمدتی بین مصرف برق صنعتی و GDP وجود دارد و به نظر می‌رسد که فاقد واکنش قابل توجه قیمتی در مقابل سیاست‌های توسعه‌ای دولت باشد. مهمت[22] (2009) و جمیل[23] (2010) نیز به بیان رابطه میان تقاضای برق و قیمت برق در ترکیه و پاکستان پرداخته اند.

همانطور که در بخش مرور بر مطالعات دیده شد، اکثر مطالعات انجام شده خارج از کشور در زمینه انرژی به ویژه مطالعات قدیمی‌تر و همچنین همه مطالعات داخلی که در این مقاله مورد بررسی قرار گرفته‌اند روش دومرحله‌ای را با استفاده از مدل ترانسلوگ انجام داده‌اند. در حالت عمومی نیز در بسیاری از مطالعات مربوط به انرژی از تابع ترانسلوگ استفاده می‌شود. چنانچه تابع هزینه یا مطلوبیت عوامل از نوع ترانسلوگ باشد ویژگی‌های قانونی توابع هزینه یا مطلوبیت فقط در محدوده‌ای از قیمت‌های نسبی عوامل حفظ می‌شود و در خارج از این محدوده خاصیت‌های غیرمنفی برای مقادیر سهم بازار و شبه مقعر بودن تابع هزینه ارضاء نمی‌شود. در این مقاله از تئوری رفتار مصرف‌کننده استفاده می‌شود که می‌تواند برای تئوری تولیدکننده نیز مورد استفاده قرار گیرد و از این نظر تفاوتی نخواهد داشت.

 

4. روش تحقیق

1-4. معرفی مدل

در این مطالعه از یک مدل بهینه‌سازی  دومرحله‌ای با فروض زیر استفاده شده است:

  1. تابع تولید یک بنگاه، بازده ثابت نسبت به مقیاس دارد.
  2. هر بنگاه اقتصادی به منظور مشخص نمودن سطح بهینه تولید و یا نهاده‌های مورد نیاز، تابع هدف خود (تابع سود و یا تابع هزینه) را بهینه‌سازی می‌نماید.
  3. در میان نهاده‌های تولید، نهاده انرژی نسبت به دیگر نهاده‌ها جدایی‌پذیر ضعیف می‌باشد بنابراین:
  • به دلیل استقلال میان نهاده‌های تولید مشکل همخطی در مدل کاهش می‌یابد.
  • امکان تخمین مقادیر نهاده‌های انرژی در روش دومرحله‌ای را فراهم می‌کند. به عبارتی در یک مرحله کل تقاضای انرژی برآورد شده و در مرحله بعدی اجزای تشکیل‌دهنده این تقاضا مورد بررسی قرار می‌گیرند.
  • به دلیل استفاده از روش دومرحله‌ای تعداد پارامترها نسبت به حالتی که کل تقاضای انرژی و اجزای تشکیل‌دهنده آن در کنار یکدیگر تخمین زده شوند، کاهش می‌یابد.

در مرحله اول با در نظر گرفتن عوامل مؤثر بر میزان تقاضای انرژی مانند فرآورده‌های نفتی (بنزین، نفت گاز، نفت سفید، نفت کوره)، برق، زغال سنگ و گاز طبیعی با استفاده از روش OLS تقاضای کل انرژی در بخش صنعت تخمین زده شده و برای این منظور از آمار و اطلاعات سری زمانی مربوط به سال­های 1367 تا 1389 بانک مرکزی جمهوری اسلامی ایران، مرکز آمار ایران، ترازنامه­های انرژی و آمار تفضیلی صنعت برق ایران برای متغیرهای مربوطه استفاده شده است. در مرحله دوم با استفاده از روش رگرسیون معادلات به ظاهر نامرتبط سهم بازار تک‌تک حامل‌های انرژی براساس قیمت‌های نسبی آنها مشخص شده است. با استفاده از این روش می‌توان کششهای قیمتی و متقاطع را برای حامل‌های انرژی با در نظر گرفتن معادلات سهم هزینه تخمین زد.

2-4. مرحله اول

در این مرحله تقاضای کل انرژی صنعت تابعی است از شاخص قیمت انرژی، ارزش‌فزوده واقعی بخش صنعت و میزان تقاضای انرژی و ارزش افزوده با یک وقفه زمانی در نظر گرفته شده است. میزان مصرف انرژی در بخش صنعت به عنوان شاخص تقاضای انرژی در نظر گرفته می‌شود. شاخص قیمت انرژی میانگین  وزنی قیمت‌های سوخت با در نظر گرفتن سهم کمی آنها در میزان مصرف انرژی می‌باشد. مدل مذکور در این مرحله یک مدل لگاریتمی و به صورت زیر می‌باشد:

(12)                 

: لگاریتم مصرف انرژی در زمان t

: لگاریتم ارزش افزوده واقعی بخش صنعت در زمان t

: لگاریتم شاخص قیمت انرژیدر زمانt

مدل نیازی به عرض از مبدأ ندارد.

 

3-4. مرحله دوم مدل

به منظور برآورد توابع تقاضای حامل‌های انرژی (گاز طبیعی، فرآورده‌های نفتی، زغال سنگ و برق( از مطالعه تئوریک کانسیدین و مونت با استفاده از روش رگرسیون سیستم معادلات به ظاهر نامرتبط استفاده می‌شود. در این سیستم با وارد شدن متغیر زمان سیستم دینامیک مصارف و سهم حامل‌های انرژی قابل بررسی می‌باشد. جونز گلیفتون در راستای مطالعه کانسیدین و مونت سیستم معادلات ایستا و پویای هزینه را مورد توجه قرار داده‌اند. معادلات  [24]SUR با توجه به موارد ذکر شده به صورت زیر است:

(13)                          

(14)                          

(15) 

، ،  و  به ترتیب میانگین سهم هزینه سوخت i ام، سهم هزینه سوخت i ام، قیمت واقعی سوخت i ام و متغیر زمان می‌باشد. 1،2،3،4i=به ترتیب برای فرآورده‌های نفتی، برق، ذغال سنگ و گاز طبیعی می‌باشد.در هر یک از سال های مورد بررسی حامل‌های انرژی سهم‌های متفاوتی در تقاضای انرژی دارند. به عبارتی سهم هزینه هر یک از حامل‌های انرژی در دوره مورد بررسی (1389-1367) تغییر می‌کند. به منظور اعمال محدودیت، از میانگین سهم هزینه برای هر یک از متغیرها استفاده شده است.

 

5. تخمین مدل

1-5. آمارهای مورد استفاده

ابتدا لازم است در مورد آمارهای بکار رفته در تخمین مدلها توضیحاتی داده شود. در این مطالعه اطلاعات آماری سری زمانی دوره 1389-1367 شامل: تقاضای انرژی در بخش صنعت، ارزش افزوده بخش صنعت، سهم‌های هزینه حامل‌های انرژی، شاخص قیمت انرژی و قیمت‌های واقعی حامل‌های انرژی مورد استفاده قرار گرفته است. مجموع مصرف انرژی یکایک حامل‌های انرژی در صنعت (فرآورده‌های نفتی، گاز طبیعی، سوخت‌های جامد و برق) برحسب میلیون بشکه نفت خام به عنوان تقاضای کل انرژی در بخش صنعت در نظر گرفته شده است. شاخص قیمت انرژی جمع وزنی قیمت انواع سوخت‌ها براساس سهم کمی آنها در مصرف انرژی مورد محاسبه قرار گرفته است.

(16)                                           

 و  به ترتیب قیمت واقعی و سهم کمی حامل انرژی i ام در دوره t می‌باشد. قیمت فرآورده‌های نفتی برحسب لیتر به ریال، قیمت سوخت‌های جامد برحسب کیلوگرم به ریال، قیمت برق برحسب کیلووات ساعت به ریال و قیمت گاز طبیعی برحسب مترمکعب به ریال می‌باشد.

براساس آمارهای ترازنامه انرژی، مصرف انرژی برحسب میلیون بشکه نفت خام می‌باشد. بنابراین قیمت واقعی هر یک از حامل‌های انرژی براساس بشکه نفت خام مورد محاسبه قرار می گیرد. بدین منظور قیمت‌های جاری بر شاخص عمده‌فروشی (سال پایه 1376) تقسیم و طبق فرمول 16 شاخص قیمت انرژی محاسبه می‌شود. ارزش افزوده بخش صنعت به قیمت ثابت سال 1376 مورد استفاده قرار گرفته است. سهم‌های هزینه حامل‌های انرژی بوسیله نسبت حاصل ضرب قیمت واقعی و مقدار مصرف هر یک از سوخت‌ها بر هزینه کل انرژی مصرفی در صنعت بدست آمده است. هزینه کل انرژی مصرفی به وسیله مجموع حاصلضرب قیمت واقعی تک‌تک سوخت‌ها و مقدار کمی آنها محاسبه شده است.

 

2-5. نتایج تخمین

1-2-5. مرحله اول

در مرحله اول تقاضای کل انرژی به صورت تابعی از ارزش افزوده بخش صنعت، شاخص قیمت انرژی، ارزش افزوده و تقاضای انرژی با یک وقفه زمانی با استفاده از روش حداقل مربعات معمولی[25] تخمین زده می‌شود.

(17)       

در این مرحله با توجه به مقادیر آماره t (داخل پرانتز) تنها متغیر تقاضای انرژی با یک وقفه زمانی معنادار می‌باشد و دیگر متغیرها به لحاظ آماری معنادار نمی‌باشند. براساس تئوری تقاضا دو متغیر اصلی در تقاضا قیمت و درآمد می‌باشد. در مدل (17) این دو متغیر اصلی وارد شده‌اند. همچنین آزمون بروش- گادفری نشان می‌دهد که مدل با مشکل خودهمبستگی مواجه می‌باشد (مقدار آماره F در مدل برابر 36/0 می‌باشد).

به منظور دستیابی به یک مدل مناسب، متغیر عرض از مبدأ به علت عدم معناداری از مدل حذف می‌شود. نتیجه تخمین نشان می‌دهد که پس از حذف عرض از مبدأ تمامی متغیرها به جز متغیر ارزش افزوده با یک وقفه زمانی معنادار می‌شوند.

(18)       

مقدار آماره دوربین واتسون، کاهش یافته بنابراین مدل از خودهمبستگی برخوردار نیست. برای اطمینان بیشتر، مسئله خودهمبستگی به وسیله روش بروخ-گادفری نیز بررسی شد که احتمال پذیرش فرضیه عدم خودهمبستگی 51% است که در این صورت فرضیه پذیرفته می‌شود. متغیرهایی که به عنوان متغیر توضیحی بکار رفته‌اند همگی پایه تئوریک دارند.

در مرحله بعدی با توجه به عدم معناداری متغیر ارزش‌افزوده در بخش صنعت با یک وقفه، این متغیر از مدل حذف می‌شود. نتایج نشان می‌دهد که کلیه متغیرها معنادار می‌باشد. آماره دوربین واتسون نسبت به حالت قبلی افزایش داشته است. مقدار آماره F در روش بروخ- گادفری نشان می‌دهد که مدل با مشکل خودهمبستگی مواجه می‌باشد (مقدار آماره 47/0 است که بسیار به 5/0 نزدیک است).

(19)                          

می‌توان با وارد کردن یک متغیر (1)AR مشکل خودهمبستگی ایجاد شده در این قسمت را رفع نمود. نتایج نشان می‌دهد که آماره دوربین واتسون به مقدار حدی (مقدار 2) بسیار نزدیک شده و همچنین سطح معناداری هر سه متغیر در مدل افزایش یافته است. همچنین روش بروخ- گادفری نشان می‌دهد که مدل دچار مشکل خودهمبستگی نمی‌باشد (آماره F برابر است با 75/0) . همچنین در گزارش رگرسیون روش بروخ-گادفری متغیر جزء اخلال تأخیری در معادله معنادار نیست و این تأییدی بر مسئله عدم وجود خودهمبستگی باقیمانده‌ها است.

(20)              

بنابراین در نهایت مدل (20) هم به لحاظ معناداری متغیرها و کل رگرسیون و به دلیل مقدار آماره دوربین- واتسون بهتر از مدل‌های قبلی است.

 

2-2-5. مرحله دوم

هدف در مرحله دوم، تخمین تابع تقاضای هر یک از نهاده‌ها (حامل‌های انرژی) می‌باشد. در این مرحله از روش تخمین رگرسیون معادلات به ظاهر غیرمرتبط استفاده می‌شود. این روش مخصوص تخمین یک سری از معادلات است که ظاهراً ارتباطی با یکدیگر ندارند. هر یک از توابع تقاضای نهاده‌ها به صورت مستقل از یکدیگر اما همه آنها به وسیله لم شفارد و از یک تابع هزینه بدست می‌آیند. متوسط سهم هزینه حامل‌های انرژی برای فرآورده‌های نفتی، برق، زغال سنگ و گاز طبیعی به ترتیب 52/0، 00018/0، 13/0 و 34/0 می‌باشند. گزارش اولیه تخمین مدلهای مرحله دوم در جدول پیوست آمده است. همانطور که مشاهده می‌شود معیار برازش در معادلات مناسب است اما متغیرهای  و  از نظر آماری معنادار نمی‌باشند. به همین علت ضرایب اصلی و کشش‌های متناظر با این ضرایب نیز معنادار نمی‌باشند. ضرایب تخمین زده شده در این مرحله مقادیر  می‌باشند و به منظور محاسبه ضرایب اصلی () بایستی از روابط زیر استفاده نمود:

(21)                                                         

(22)                                                           

که d یک مقدار نامعلوم است. ضرایب اصلی محاسبه شده در جدول 1 آمده است.

 

 

جدول 1. ضرایب اصلی محاسبه شده از روی ضرایب برآورد شده

49108/0-

 

1253/0-*

 

53418/0-

 

3554/0

 

0001/0-*

 

00037/0-

 

0042/1

 

00018/0-

 

26064/0

 

7795/0

 

28056/0-

 

03279/0-*

 

7518/0

 

6538/0-

 

18947/0-*

 

032242/0-

 

* این ضرایب به لحاظ آماری معنادار نمی‌باشند

منبع: محاسبات محقق

 

6. نتیجه‌گیری

 با توجه به تخمین مرحله اول می‌توان نتیجه گرفت که چون مدل لگاریتمی است بنابراین ضرایب همان کششها هستند. کشش قیمتی کل انرژی 1/0-  است که نشان‌دهنده کم‌کشش بودن تقاضای انرژی است. یعنی افزایش قیمت، روی کاهش تقاضا تأثیر زیادی ندارد. ضرایب متغیر لگاریتم ارزش افزوده بخش صنعت ( و ) همان کشش‌های درآمدی یا تولیدی کوتاه‌مدت تقاضای بخش صنعت می‌باشند. کشش درآمدی اثر افزایش یک درصد در درآمد صنعت را روی تقاضای انرژی نشان می‌دهد. طبق نتایج بدست آمده، ارزش افزوده با یک وقفه تأخیر بر روی مصرف انرژی تأثیری نداشته و به لحاظ آماری بی‌معنا می‌باشد. به عبارتی تغییر در درآمد بخش صنعت در دوره قبل تأثیری در میزان مصرف انرژی در دوره حال نخواهد داشت. کشش اثر درآمدی دوره جاری برابر با 14/0 می‌باشد. بیشترین کشش در بخش اول مربوط به متغیر مصرف انرژی با یک وقفه زمانی است. به عبارتی افزایش میزان مصرف انرژی در دوره قبلی باعث افزایش میزان مصرف در دوره جاری خواهد شد.  با نگاهی به نتایج مرحله دوم می‌توان به این نکته رسید که عرض از مبدأ و متغیر زمان در تمام مدل‌ها بامعناست و این خود نشان می‌دهد تابع تقاضای نهاده‌ها در صنعت تابعی از متغیرهای برونزا است. همچنین یکی از کاربردهای چنین مطالعاتی محاسبه کشش‌های قیمتی و جانشینی است .محاسبه این کشش‌ها از روابط زیر میسر است.

(23)                                                    

(24)                                                

(25)                                                              

همانطور که مشاهده می‌شود محاسبه کشش‌های قیمتی و جانشینی نیازمند محاسبه کشش‌های سهم هزینه می‌باشد. از آنجا که مدل‌ها لگاریتمی هستند، ضرایب اصلی مدل () همان کشش‌های سهم می‌باشند. اکنون به جای مقادیر () از مقادیر محاسبه شده برای ضرایب اصلی در جدول 1 استفاده می‌کنیم.

با توجه به اینکه متغیرهای  و  از نظر آماری معنادار نمی‌باشند، بنابراین کشش‌های قیمتی سهم هزینه و جانشینی متناظر (کشش قیمتی متقاطع میان فرآورده‌های نفتی و ذغال سنگ و همچنین کشش قیمتی میان برق و گاز طبیعی) نیز از نظر آماری معنادار نیستند. نتایج مربوط به کشش قیمتی در جدول 2 آمده است.

 

جدول 2. کشش قیمتی سهم هزینه

کشش قیمتی سهم هزینه (خودی و متقاطع)

گاز طبیعی

ذغال سنگ

برق

فرآورده‌های نفتی

 

0196/0

0987/0*

6-10×48/4-

1282/0-

فرآورده‌های نفتی

1517/0*

1373/0-

0112/0-

0129/0-

برق

3096/0-

0949/0-

00019/0-

3949/0*

ذغال سنگ

0784/0

1184/0-

5-10×03/8*

03/0

گاز طبیعی

* این ضرایب به لحاظ آماری معنادار نمی‌باشند

منبع: محاسبات محقق

 

کشش‌های قیمتی برای هر کدام از نهاده‌های فرآورده‌های نفتی، برق و ذغال سنگ منفی و برای نهاده گاز طبیعی مثبت است. کمترین کشش متقاطع میان فرآورده‌های نفتی و برق می‌باشد. به عبارتی با تغییر در قیمت فرآورده‌های نفتی، مقدار مصرف برق افزایش چندانی نخواهد داشت. بیشترین کشش متقاطع میان ذغال سنگ و فرآورده‌های نفتی است یعنی با افزایش قیمت ذغال سنگ مقدار مصرف فرآورده‌های نفتی افزایش محسوسی خواهد داشت.

نتایج نشان می‌دهد که تمام نهاده‌ها در بخش صنعت نسبت به قیمت‌های خود کم‌کشش هستند. می‌توان نتیجه گرفت که یارانه‌های انرژی و سوبسیدهای ارائه شده باعث می‌شوند که انرژی با قیمت نسبتاً ارزانی بدست مصرف‌کننده برسد. بنابراین تغییر در قیمت حامل‌های انرژی تأثیر چندانی بر روی مصرف انرژی در بخش صنعت ندارد. پایین بودن کشش‌های قیمتی در بخش صنعت نشان‌دهنده حساسیت پایین تقاضای نهاده‌های انرژی به تغییرات قیمتی می‌باشد. بنایراین هرچند با تغییر الگوهای قیمتی امکان جانشینی میان نهاده‌ها وجود دارد اما کم‌کشش بودن نهاده‌ها نسبت به قیمت نشان می‌دهد که تقاضای انرژی صرفاً به قیمت‌ها بستگی نداشته و به شرایط تکنولوژی‌های مورد استفاده در بخش صنعت و عوامل محیطی دیگر به خصوص یارانه‌ها و سوبسیدهای دولتی در بخش صنعت نیز وابسته است.

نتایج جدول 2، نشان‌دهنده مثبت بودن کشش قیمتی نهاده گاز طبیعی می باشد. این موضع احتمالاً به علت روند قیمت‌گذاری گاز طبیعی در کشور می‌باشد. قیمت‌گذاری گاز طبیعی دیرتر از سایر نهاده‌های انرژی در کشور اجرایی گردید و گاز طبیعی نسبت به سایر حامل‌های انرژی هنوز دارای قیمت پایینی است. از طرفی نکته بسیار مهم درخصوص گاز طبیعی رابطه بسیار قوی میان مصرف گاز طبیعی و درجه حرارت هوا می‌باشد.[26] بنابراین مصرف‌کننده در فصول سرد توانایی کاهش مصرف به تناسب افزایش قیمت را ندارد. از سوی دیگر دولت در دهه اخیر به استفاده بیشتر از گاز طبیعی تلاش‌های بسیاری را در بخش‌های مختلف و به خصوص بخش صنعت انجام داده است. تمام این موارد باعث ناهمگونی روند قیمت واقعی گاز طبیعی گردیده به گونه‌ای که روند قیمت گاز طبیعی در طول سال‌های مورد بررسی صعودی  بوده است.

 

جدول 3. کشش جانشینی میان نهاده‌ها

کشش جانشینی

گاز طبیعی

ذغال سنگ

برق

فرآورده‌های نفتی

 

05779/0

7595/0*

0249/0-

2465/0-

فرآورده‌های نفتی

4463/0*

0568/1-

322/62-

 

برق

9108/0-

7302/0-

 

 

ذغال سنگ

2307/0

 

 

 

گاز طبیعی

 *این ضرایب به لحاظ آماری معنادار نمی‌باشند

منبع: محاسبات محقق

 

با توجه به کشش‌های جانشینی می‌توان گفت نهاده‌های فرآورده‌های نفتی و برق، برق و ذغال سنگ و همچنین ذغال سنگ و گاز طبیعی به صورت نهاده‌های مکمل رفتار می‌کنند. به عبارتی با افزایش قیمت یکی از این نهاده‌ها، میزان مصرف دیگر نهاده‌های در بخش صنعت کاهش می‌یابد. اما نهاده‌های فرآورده‌های نفتی و گاز طبیعی جانشین  یکدیگر هستند. همچنین با توجه به کشش‌های جانشینی، نهاده‌های برق و گاز طبیعی و فرآورده‌های نفتی و ذغال سنگ نیز جانشین هستند اما متغیرها به لحاظ آماری معنادار نمی‌باشند. با این حال با توجه به لزوم جایگزینی برق و گاز طبیعی و تلاش‌های صورت گرفته در این راه در سالهای اخیر بهتر است این متغیرها از مدل حذف نشوند، به خصوص در زمانی که از این نتایج به منظور پیش‌بینی استفاده شود. 

یکی از دلایل پایین بودن کشش‌های جانشینی میان نهاده‌ها در بخش صنعت سنتی بودن تکنولوژی مورد استفاده در این بخش است. به عبارتی به علت عدم وجود سرمایه و دانش فنی برای استفاده از روش‌های مدرن انعطاف‌پذیری فنی جهت جانشینی میان نهاده‌ها در این بخش بسیار پایین است. همچنین نهاده‌های برق و ذغال سنگ سهم هزینه بسیار اندکی (کمتر از یک درصد سهم هزینه) را در بخش صنعت به خود اختصاص داده‌اند و همین سهم هزینه اندک نیز می‌تواند باعث کاهش جانشینی میان نهاده شود.

 

منابع

الف- فارسی

بهبهانی‌فرد، پ. (1383)، «بررسی عوامل مؤثر بر تقاضای انواع انرژی در بخش صنعت استان اصفهان و تخمین کشش­های جانشینی بین آنها»، نهمین کنفرانس شبکه‌های تویع نیروی برق.

دینایی، م. (1378)، «تخمین تابع تقاضای برق خانگی در شهرستان کرمان»، سومین همایش ملی انرژی ایران.

صادقی، ن. (1382)، پیش‌بینی مصرف برق با استفاده از مدل‌های اقتصاد"، پایان‌نامه کارشناسی ارشد، دانشگاه تهران.

صمدی، س. (1387)، «تحلیل تقاضای برق در ایران با استفاده از مفهوم هم‌جمعی و مدل ARIM طی 1363تا 1387»،  مجلة دانش و توسعه.

فاکهی و همکاران (1386)، «مدل‌سازی تقاضای انرژی مفید در بخش‌های مختلف اقتصادی»، ششمین همایش ملی انرژی،22 و 23 خرداد ماه.

فخرایی، ح. (1371)، گزارش نهایی از پیش‌بینی تقاضای حامل‌های انرژی در بخش‌های مصرفی مختلف.

قادری، ف. و ش. بامداد (1384)، «مدل ریاضی جهت تخمین تابع مصرف انرژی الکتریکی با استفاده از روش سری زمانی»، مؤسسه پژوهش در مدیریت و برنامه‌ریزی انرژی، نشریه انرژی ایران، سال نهم، شماره 23.

کاظمی، ا. (1375)، تحلیل تقاضای انرژی الکتریکی در بخش خانگی و صنعتی، پایان‌نامه کارشناسی ارشد، دانشگاه تهران.

کشاورزحداد، غلامرضا و محمد میرباقری (1386)، «بررسی تابع تقاضای گاز طبیعی»، فصلنامه پژوهشهای اقتصادی ایران، سال نهم، شماره 32، ص 160-137.

لطفعلی­پور، م. و ا. لطفی (1383)، «بررسی برآورد عوامل مؤثر بر تقاضای برق خانگی در استان خراسان»، مجله دانش و توسعه (علمی ـ پژوهشی)، شماره 15.

معاونت امور برق و انرژی، مروری بر 23 سال آمار انرژی کشور، دفتر برنامه‌ریزی کلان برق و انرژی. اسفند1390

منظور، داوود (1388)، «روش‌شناسی الگوسازی تقاضای انرژی»، هفتمین همایش ملی انرژی، 1 و 2 دی ماه.

مهرآرا، م. و ع. عبدی (1384)، «برآورد توابع تقاضا برای نهاده های ساختمانی: مورد ایران»، فصلنامه پژوهشهایاقتصادی، سال پنجم شماره چهارم، زمستان.

 

ب- انگلیسی

Bilgili, Mehmet (2009), “Estimation of the Net Electricity Consumption of Turkey”, Cukurova University, Adana Vocational School of Higher Education, J.of Thermal Science and Technology, No. 29, Vol. 2, pp. 89-98.

Considine ,T. J. (1989), “Separability, Functional form and Regulatory Policy in Model of Interfuel Substitution”, Energy Economics.

Considine, T. J . and T. D. Mount (1984), “The Use of Logit Models for Dynamic Input Demand System”, Review of Economics and Statistics, No. 66, pp. 434-443.

Ishi, Kotaro and Joutz Fred (2010), “Modeling and Forecasting Electricity Demand in the Philippines”, The International Monetary Fund.

Jamil, Faisal and Ahmad, Eatzaz (2010), “The Relationship between Electricity Consumption, Electricity Prices and GDP in Pakistan”, Energy Policy.

 


پیوست

نتیجه رگرسیون  معادلات به ظاهر نامرتبط

 

System: SUR1

Estimation Method: Seemingly Unrelated Regression

Date: 06/12/13   Time: 10:36

Sample: 1367 1389

Included observations: 23

Total system (balanced) observations 69

Linear estimation after one-step weighting matrix

 

 

 

 

 

Prob.  

t-Statistic

Std. Error

Coefficient

 

 

 

 

 

 

0.0000

10.38554

4.187026

43.48452

C(1)

0.0000

-8.438601

0.121455

-1.024913

C(4)

0.5547

-0.594173

0.404623

-0.240416

C(5)

0.0000

-9.280017

0.101531

-0.942205

C(6)

0.0961

-1.692248

0.327193

-0.553692

C(7)

0.0008

-3.544397

0.539113

-1.910830

C(8)

0.0000

-11.00579

0.002926

-0.032205

C(10)

0.0000

5.578521

2.856594

15.93557

C(2)

0.0039

-3.013325

0.682589

-2.056861

C(9)

0.0000

-7.124047

0.001962

-0.013981

C(11)

0.0000

8.249538

12.97718

107.0558

C(3)

0.0000

-9.256504

0.008755

-0.081037

C(12)

 

 

 

 

 

   

1.80E-07

Determinant residual covariance

 

 

 

 

 

Equation: LW1_4=C(1)-(0.00018*C(4)+0.1346*C(5)+0.8634*C(6))

        *LP1_4+(C(4)-C(7))*0.00018*LP2_4+(C(5)-C(8))*0.1364*LP3_4

        +C(10)*T

Observations: 23

0.272928

 

    Mean dependent var

0.969990

R-squared

0.399463

 

    S.D. dependent var

0.958736

Adjusted R-squared

0.105352

 

    Sum squared resid

0.081145

S.E. of regression

     

0.906613

Durbin-Watson stat

 

 

 

 

 

Equation: LW2_4=C(2)-(0.5212*C(4)+0.1364*C(9)+0.34238*C(7))

        *LP2_4+(C(4)-C(6))*0.5212*LP1_4+(C(9)-C(8))*0.1364*LP3_4

        +C(11)*T

Observations: 23

-3.245368

 

    Mean dependent var

0.911862

R-squared

0.170322

 

    S.D. dependent var

0.878810

Adjusted R-squared

0.056250

 

    Sum squared resid

0.059293

S.E. of regression

     

1.088174

Durbin-Watson stat

 

 

 

 

 

Equation: LW3_4=C(3)-(0.5212*C(5)+0.00018*C(9)+0.4786*C(8))

        *LP3_4+(C(5)-C(6))*0.5212*LP1_4+(C(9)-C(7))*0.00018*LP2_4

        +C(12)*T

Observations: 23

-0.341133

 

    Mean dependent var

0.921910

R-squared

0.825710

 

    S.D. dependent var

0.892626

Adjusted R-squared

1.171318

 

    Sum squared resid

0.270569

S.E. of regression

     

0.790063

Durbin-Watson stat

 

 

 

 

 

 

 




* کارشناس ارشد علوم اقتصادی دانشگاه شهید باهنر کرمان                                                    Sadrzadeh.saeed@gmail.com

** استادیار بخش اقتصاد دانشگاه شهید باهنر کرمان                                                                          Abed.sadeghi@gmail.com

*** دانشجوی دکتری مدیریت دانشگاه تهران                                                                        Ahmad.ghodselahi@gmail.com

[1]. Translog Cost Function

[2]. Christensen and Caves (1980)

[3]. Cob-Douglas

[4]. منظور و همکاران (1388)

[5]. Linear Logit Cost Share Model

[6]. Shephard's Lemma

[8]. Considine T. J. (1989)

[9].  Berkson (1994)

[10].  Theil (1969)

[11].  Mc Fadden (1974)

[12].  Ernest r. Berndt and David D. Wood (1975)

[13].  Robert Halvorsen (1977)

[14].  Melvin A.Fuss (1977)

[15].  Robert S.Pindyck (1979)

[16].  Temothy J. Consdintimothy D. Mount (1984)

[17].  Temothy J. Consdintimothy (1989)

[18].  Sanford W. Bergprakash  Lovughani (1990)

[19].  Weakly Sepawrability

[20].  Cost Share

[21].  Ishi and Fred (2010)

[22]. Mehmet B. (2009)

[23].  Jamil F. (2010)

[24]. Simingly Unrelated Regression Estimation

[25]. Ordinary Least Square

[26]. کشاورز حداد (1386)

الف- فارسی
بهبهانی‌فرد، پ. (1383)، «بررسی عوامل مؤثر بر تقاضای انواع انرژی در بخش صنعت استان اصفهان و تخمین کشش­های جانشینی بین آنها»، نهمین کنفرانس شبکه‌های تویع نیروی برق.
دینایی، م. (1378)، «تخمین تابع تقاضای برق خانگی در شهرستان کرمان»، سومین همایش ملی انرژی ایران.
صادقی، ن. (1382)، پیش‌بینی مصرف برق با استفاده از مدل‌های اقتصاد"، پایان‌نامه کارشناسی ارشد، دانشگاه تهران.
صمدی، س. (1387)، «تحلیل تقاضای برق در ایران با استفاده از مفهوم هم‌جمعی و مدل ARIM طی 1363تا 1387»،  مجلة دانش و توسعه.
فاکهی و همکاران (1386)، «مدل‌سازی تقاضای انرژی مفید در بخش‌های مختلف اقتصادی»، ششمین همایش ملی انرژی،22 و 23 خرداد ماه.
فخرایی، ح. (1371)، گزارش نهایی از پیش‌بینی تقاضای حامل‌های انرژی در بخش‌های مصرفی مختلف.
قادری، ف. و ش. بامداد (1384)، «مدل ریاضی جهت تخمین تابع مصرف انرژی الکتریکی با استفاده از روش سری زمانی»، مؤسسه پژوهش در مدیریت و برنامه‌ریزی انرژی، نشریه انرژی ایران، سال نهم، شماره 23.
کاظمی، ا. (1375)، تحلیل تقاضای انرژی الکتریکی در بخش خانگی و صنعتی، پایان‌نامه کارشناسی ارشد، دانشگاه تهران.
کشاورزحداد، غلامرضا و محمد میرباقری (1386)، «بررسی تابع تقاضای گاز طبیعی»، فصلنامه پژوهشهای اقتصادی ایران، سال نهم، شماره 32، ص 160-137.
لطفعلی­پور، م. و ا. لطفی (1383)، «بررسی برآورد عوامل مؤثر بر تقاضای برق خانگی در استان خراسان»، مجله دانش و توسعه (علمی ـ پژوهشی)، شماره 15.
معاونت امور برق و انرژی، مروری بر 23 سال آمار انرژی کشور، دفتر برنامه‌ریزی کلان برق و انرژی. اسفند1390
منظور، داوود (1388)، «روش‌شناسی الگوسازی تقاضای انرژی»، هفتمین همایش ملی انرژی، 1 و 2 دی ماه.
مهرآرا، م. و ع. عبدی (1384)، «برآورد توابع تقاضا برای نهاده های ساختمانی: مورد ایران»، فصلنامه پژوهشهایاقتصادی، سال پنجم شماره چهارم، زمستان.
 
ب- انگلیسی
Bilgili, Mehmet (2009), “Estimation of the Net Electricity Consumption of Turkey”, Cukurova University, Adana Vocational School of Higher Education, J.of Thermal Science and Technology, No. 29, Vol. 2, pp. 89-98.
Considine ,T. J. (1989), “Separability, Functional form and Regulatory Policy in Model of Interfuel Substitution”, Energy Economics.
Considine, T. J . and T. D. Mount (1984), “The Use of Logit Models for Dynamic Input Demand System”, Review of Economics and Statistics, No. 66, pp. 434-443.
Ishi, Kotaro and Joutz Fred (2010), “Modeling and Forecasting Electricity Demand in the Philippines”, The International Monetary Fund.
Jamil, Faisal and Ahmad, Eatzaz (2010), “The Relationship between Electricity Consumption, Electricity Prices and GDP in Pakistan”, Energy Policy.