Giới thiệu dự án

Bối cảnh thị trường và vấn đề thực tiễn

Thị trường chứng khoán Việt Nam (TTCK) trải qua hơn một thập kỷ hình thành và phát triển với tốc độ nhanh chóng nhưng cũng tiềm ẩn nhiều rủi ro mang tính chu kỳ sâu sắc. Sau giai đoạn bùng nổ đỉnh điểm 2006–2007 (khi VN-Index đạt đỉnh 921,07 điểm với tổng lượng vốn huy động kỷ lục lên tới 40.000 tỷ đồng), thị trường đã chứng kiến sự sụt giảm nghiêm trọng xuống dưới ngưỡng 600 điểm vào năm 2008 do tâm lý đầu tư "bầy đàn" và sự thiếu hụt các công cụ lượng hóa rủi ro độc lập.

Đến giai đoạn 2012–2014, cấu trúc thị trường bộc lộ nhiều điểm nghẽn kỹ thuật và cơ bản:

  • Chất lượng niêm yết phân hóa mạnh: 20 doanh nghiệp hàng đầu chiếm tới 74% tổng vốn hóa toàn thị trường, trong khi hơn 341 công ty có vốn điều lệ dưới 100 tỷ đồng (chiếm 50% số mã niêm yết) và gần 60% cổ phiếu giao dịch dưới mệnh giá.
  • Rủi ro vận hành & thanh khoản: Cơ chế khớp lệnh và thanh toán bù trừ áp dụng chu kỳ T+3 (thực tế là T+4) làm hạn chế năng lực quản trị rủi ro danh mục trong ngày, cùng với việc thiếu vắng các công cụ phái sinh phòng hộ giá.
  • Vị thế ngành Năng lượng điện: Theo Quy hoạch phát triển điện lực quốc gia giai đoạn 2011–2020 có xét đến năm 2030 (Tổng sơ đồ VII), nhu cầu điện thương phẩm dự kiến đạt 330–362 tỷ kWh vào năm 2020. Dù ngành điện sở hữu dòng tiền kinh doanh tương đối ổn định nhờ tính chất dịch vụ hạ tầng thiết yếu, cổ phiếu ngành điện trên sàn HOSE/HNX vẫn chịu tác động lớn bởi rủi ro biến động giá nguyên liệu (than, khí đối với nhiệt điện), yếu tố thủy văn (đối với thủy điện), đòn bẩy tài chính cao và chính sách tái cơ cấu thị trường phát điện cạnh tranh theo 3 giai đoạn của Chính phủ.

Vấn đề cốt lõi (Problem Statement)

Các nhà đầu tư và tổ chức tài chính tại Việt Nam đối mặt với 3 thách thức lớn:

  1. Hiện tượng phương sai thay đổi có điều kiện (Conditional Heteroskedasticity): Lợi suất chuỗi thời gian của cổ phiếu năng lượng không có phương sai bất biến (homoskedasticity), dẫn đến việc các mô hình hồi quy tuyến tính cổ điển OLS (Ordinary Least Squares) cho ước lượng sai lệch và đánh giá thấp rủi ro đuôi (tail risk).
  2. Sự tập trung của biến động (Volatility Clustering): Các biên độ biến động lớn thường có xu hướng kéo theo các biến động lớn tiếp theo, khiến các chỉ báo kỹ thuật truyền thống mất đi độ chính xác trong dự báo rủi ro biến động giá ngắn hạn.
  3. Mất cân đối giữa định giá cơ bản và rủi ro thị trường: Thiếu một khung phân tích tích hợp kết hợp giữa sức khỏe tài chính doanh nghiệp (khả năng thanh toán, hiệu quả vốn lưu động, biên lợi nhuận ROA/ROE) và các mô hình lượng lượng kinh tế lượng chuỗi thời gian để đo lường mức độ bền bỉ của biến động rủi ro (Volatility Persistence).

Mục tiêu dự án

  1. Khảo sát và đánh giá toàn diện bức tranh tài chính, khả năng thanh toán, cấu trúc nợ và hiệu quả sinh lời của nhóm doanh nghiệp đầu ngành năng lượng điện (tiêu biểu: VSH - CTCP Thủy điện Vĩnh Sơn - Sông Hinh và PPC - CTCP Nhiệt điện Phả Lại).
  2. Xây dựng và kiểm định chuỗi lợi suất cổ phiếu thông qua quy trình kinh tế lượng hiện đại: Kiểm định tính dừng nghiệm đơn vị (Augmented Dickey-Fuller - ADF), lọc nhiễu tự hồi quy tích hợp trung bình trượt ARIMA.
  3. Mô hình hóa phương sai thay đổi có điều kiện bằng họ mô hình ARCH/GARCH/IGARCH nhằm lượng hóa độ nhạy cảm và độ trễ của các cú sốc biến động giá trên thị trường chứng khoán Việt Nam.
  4. Cung cấp bộ chỉ số dự báo rủi ro hỗ trợ nhà đầu tư cá nhân và định chế tài chính tối ưu hóa danh mục đầu tư.

Phạm vi và Giới hạn nghiên cứu

  • Đối tượng: 2 mã cổ phiếu có quy mô vốn hóa lớn đại diện cho hai mảng trụ cột của ngành năng lượng điện: CTCP Thủy điện Vĩnh Sơn - Sông Hinh (Mã: VSH) và CTCP Nhiệt điện Phả Lại (Mã: PPC - vốn hóa 7.031,27 tỷ VNĐ).
  • Khung thời gian:
    • Báo cáo tài chính doanh nghiệp giai đoạn 2009–2012.
    • Dữ liệu chuỗi thời gian giá đóng cửa và lợi suất hàng ngày gồm 462 phiên giao dịch liên tục từ ngày 21/05/2012 đến 26/03/2014.

Phân tích và thiết kế giải pháp

Phân tích hiện trạng và So sánh giải pháp

Tiêu chí Phân tích cơ bản truyền thống (FA) Phân tích kỹ thuật (TA) & SMA Volatility Mô hình Kinh tế lượng ARCH-GARCH
Bản chất dữ liệu Báo cáo tài chính định kỳ (quý/năm) Giá & Khối lượng theo phiên Chuỗi thời gian Log-Return phân phối phi chuẩn
Khả năng dự báo biến động Kém, có độ trễ lớn theo chu kỳ công bố BCTC Trung bình, giả định phương sai cố định (Constant Variance) Rất cao, bắt trọn hiệu ứng Volatility Clustering
Xử lý cú sốc bất đối xứng Không thể lượng hóa Không có khả năng phân tách tin xấu/tin tốt Phân tách được qua mở rộng TGARCH/EGARCH
Đo lường Persistence ($\alpha + \beta$) Không hỗ trợ Không hỗ trợ Định lượng chính xác qua hệ số GARCH/IGARCH

Yêu cầu hệ thống theo mô hình MoSCoW

  • Must have: Kiểm định nghiệm đơn vị ADF cho chuỗi lợi suất; ước lượng phương trình trung bình ARIMA; kiểm định hiệu ứng ARCH-LM của phần dư; ước lượng tham số mô hình GARCH(1,1) và kiểm định dạng tích hợp IGARCH.
  • Should have: Phân tích đa tiêu chuẩn hệ số tài chính (Current Ratio, Quick Ratio, Vòng quay tồn kho, ROA, ROE, Nợ/VCSH).
  • Could have: Mở rộng kiểm định ảnh hưởng bất đối xứng đòn bẩy thông qua mô hình TGARCH và EGARCH.
  • Won't have (lần này): Tích hợp định giá hợp đồng quyền chọn phức tạp và giao dịch thuật toán tự động HFT (High-Frequency Trading).

Kiến trúc tổng thể hệ thống phân tích

flowchart TD
    A["Thu thập Dữ liệu Lịch sử (Giá đóng cửa P_t & Báo cáo Tài chính)"] --> B["Tiền xử lý: Tính Lợi suất Logarit r_t = ln(P_t / P_{t-1})"]
    B --> C["Kiểm định Tính dừng (Augmented Dickey-Fuller Test)"]
    C -- "Chuỗi dừng I(0)" --> D["Mô hình hóa Phương trình Trung bình ARIMA(p,d,q)"]
    D --> E["Trích xuất Phần dư u_t & Kiểm định Hiệu ứng ARCH (ARCH-LM Test)"]
    E -- "Tồn tại hiệu ứng ARCH" --> F["Ước lượng Mô hình GARCH(p,q) / IGARCH"]
    F --> G["Kiểm định Chẩn đoán (Ljung-Box Q-stat & Wald Test)"]
    G --> H["Tích hợp với Phân tích Chỉ số Tài chính Cơ bản"]
    H --> I["Dự báo Phương sai Điều kiện & Đưa ra Khuyến nghị Đầu tư"]

Implementation và kết quả

Mô hình Toán học và Kỹ thuật ước lượng

1. Phương trình Lợi suất và Phương trình Trung bình

Lợi suất logarit hàng ngày $r_t$ được tính theo công thức: $$r_t = \ln\left(\frac{P_t}{P_{t-1}}\right)$$

Phương trình trung bình dạng $\text{ARIMA}(p, d, q)$: $$r_t = \mu_t + u_t = \phi_0 + \sum_{i=1}^p \phi_i r_{t-i} + \sum_{j=1}^q \theta_j u_{t-j} + u_t$$ với $u_t = \sigma_t \epsilon_t$, trong đó $\epsilon_t \sim \text{i.i.d}(0, 1)$.

2. Cấu trúc Mô hình GARCH(1,1) và IGARCH

Phương sai có điều kiện $\sigma_t^2$ phụ thuộc vào tin tức biến động của kỳ trước (thành phần ARCH - $u_{t-1}^2$) và phương sai điều kiện dự báo kỳ trước (thành phần GARCH - $\sigma_{t-1}^2$): $$\sigma_t^2 = \alpha_0 + \alpha_1 u_{t-1}^2 + \beta_1 \sigma_{t-1}^2$$ Điều kiện ràng buộc: $\alpha_0 > 0$, $\alpha_1 \ge 0$, $\beta_1 \ge 0$, và $\alpha_1 + \beta_1 < 1$.

Khi $\alpha_1 + \beta_1 = 1$, mô hình trở thành mô hình tích hợp IGARCH(1,1): $$\sigma_t^2 = \alpha_0 + (1 - \beta_1) u_{t-1}^2 + \beta_1 \sigma_{t-1}^2$$ Trong mô hình IGARCH, tác động của một cú sốc biến động ở hiện tại sẽ tồn tại vĩnh viễn trong tương lai: $$\sigma_n^2(k) = (k - 1)\alpha_0 + \sigma_n^2(1)$$

import numpy as np
import pandas as pd
from arch import arch_model
from statsmodels.tsa.stattools import adfuller

def run_econometric_pipeline(price_series: pd.Series):
    """
    Quy trình tính toán kiểm định tính dừng và ước lượng GARCH(1,1)
    """
    # 1. Tính toán Log Returns
    log_returns = np.log(price_series / price_series.shift(1)).dropna()
    
    # 2. Kiểm định nghiệm đơn vị ADF
    adf_result = adfuller(log_returns)
    print(f"ADF Statistic: {adf_result[0]:.4f}, p-value: {adf_result[1]:.4e}")
    
    # 3. Ước lượng mô hình GARCH(1,1)
    garch_spec = arch_model(log_returns, mean='AR', lags=1, vol='Garch', p=1, q=1, dist='Normal')
    res = garch_spec.fit(disp='off')
    print(res.summary())
    
    # 4. Kiểm tra Volatility Persistence
    alpha1 = res.params['alpha[1]']
    beta1 = res.params['beta[1]']
    persistence = alpha1 + beta1
    print(f"Volatility Persistence (alpha + beta): {persistence:.4f}")
    
    return res

Kiểm định Thực nghiệm và Kết quả Ước lượng

1. Kiểm định Tính dừng chuỗi lợi suất

Kết quả kiểm định nghiệm đơn vị ADF trên mẫu 462 quan sát:

  • Cổ phiếu VSH: Thống kê ADF $t = -21,0000$ (Giá trị tới hạn $1%$ là $-2,616238$), $p\text{-value} < 0,01$. Giả thuyết $H_0$ (chuỗi có nghiệm đơn vị) bị bác bỏ hoàn toàn. Chuỗi lợi suất dừng $I(0)$.
  • Cổ phiếu PPC: Thống kê ADF $t = -19,0000$ (Giá trị tới hạn $1%$ là $-2,616238$), $p\text{-value} < 0,01$. Bác bỏ $H_0$, chuỗi lợi suất đạt tính dừng vững chắc ở mức ý nghĩa $1%$.

2. Kết quả ước lượng mô hình biến động VSH

Sau khi lọc qua mô hình trung bình ARIMA(2,2,0) có tiêu chuẩn thông tin Akaike $\text{AIC} = -5,022239$ và Durbin-Watson $\text{stat} = 2,032855$, phần dư bình phương xuất hiện tự tương quan chuỗi đáng kể (kiểm định Q-stat cho $p < 0,001$), xác nhận sự tồn tại của hiệu ứng ARCH.

Ước lượng bằng phương pháp Cực đại Hợp lý (Maximum Likelihood - Thuật toán Marquardt):

  • Phương trình phương sai điều kiện: $$\sigma_t^2 = 0,00000217 + 0,0978 u_{t-1}^2 + 0,8689 \sigma_{t-1}^2$$
  • Hệ số Persistence: $\alpha_1 + \beta_1 = 0,0978 + 0,8689 = 0,9667 \approx 1$.
  • Kiểm định Wald: Kiểm định giả thuyết $H_0: \alpha_1 + \beta_1 = 1$ cho giá trị $F\text{-statistic} = 3,051280$ ($p > 0,05$). Do đó, không đủ cơ sở bác bỏ $H_0$, chứng minh chuỗi biến động của VSH tuân theo mô hình IGARCH(1,1).
  • Phương sai không điều kiện: $$\bar{\sigma}^2 = \frac{\alpha_0}{1 - \alpha_1 - \beta_1} = \frac{0,00000217}{1 - 0,0978 - 0,8689} = 0,00006516$$

Phân tích Sức khỏe Tài chính Doanh nghiệp (Case VSH)

Bên cạnh mô hình chuỗi thời gian, phân tích báo cáo tài chính giai đoạn 2010–2012 cho thấy bức tranh chi tiết về cấu trúc rủi ro nội tại của VSH:

Nhóm chỉ tiêu Năm 2010 Năm 2011 Năm 2012 Trung bình ngành Nhận xét rủi ro
Khả năng thanh toán hiện hành (lần) 0,92 0,98 1,04 0,98 Cải thiện vượt ngưỡng 1,0; tài sản ngắn hạn đủ bù đắp nợ ngắn hạn
Khả năng thanh toán nhanh (lần) 0,38 0,35 0,60 0,44 Trung bình 3 năm đạt 0,44 (< 0,5); vẫn tiềm ẩn rủi ro thanh khoản tức thời
Tỷ lệ Tổng Nợ / VCSH (lần) 1,43 1,75 1,38 1,60 Đòn bẩy tài chính cao, chi phí lãi vay chịu ảnh hưởng nặng từ lãi suất
Vòng quay hàng tồn kho (vòng) 2,92 5,44 5,51 4,42 Tốc độ luân chuyển tốt hơn trung bình ngành
Vòng quay khoản phải thu (vòng) 7,01 16,55 11,66 15,85 Thấp hơn ngành (15,85), chính sách tín dụng khách hàng còn chịu áp lực chậm trả
Hệ số sinh lời ROA (%) 1,56% 4,21% 6,92% - Hiệu quả khai thác tổng tài sản tăng trưởng vượt bậc
Hệ số sinh lời ROE (%) 4,32% 12,46% 20,07% - Vượt xa lãi suất tiết kiệm ngân hàng (10–13%), tạo giá trị vượt trội

Đổi mới và đóng góp

  1. Phương pháp tiếp cận tích hợp 2 chiều (Hybrid Analytical Framework): Kết hợp phương pháp phân tích cơ bản (Fundamental Analysis - đo lường các hệ số an toàn thanh khoản, cấu trúc đòn bẩy và sức sinh lời) với mô hình động lượng kinh tế lượng chuỗi thời gian (GARCH/IGARCH). Khắc phục triệt để điểm mù của phân tích kỹ thuật đơn thuần.
  2. Chứng minh hiện tượng "Trí nhớ dài hạn" của rủi ro (Long-memory Persistence): Phát hiện giá trị $\alpha_1 + \beta_1 = 0,9667$ của VSH kiểm định theo Wald test tuân theo mô hình IGARCH, chứng minh các cú sốc biến động giá trên thị trường chứng khoán Việt Nam có hiệu ứng lan tỏa kéo dài (persistence), không suy giảm nhanh chóng như trong các thị trường tài chính phát triển.
  3. Mô hình hóa chuyên biệt cho ngành Năng lượng điện: Đặt các chỉ số tài chính và chuỗi lợi suất trong bối cảnh tái cơ cấu ngành điện theo Tổng sơ đồ VII và cơ chế thị trường phát điện cạnh tranh giai đoạn 1 (2005–2014), cung cấp cơ sở định giá rủi ro sát thực tế vận hành của các nhà máy nhiệt điện/thủy điện niêm yết.

Ứng dụng thực tế và triển khai

Tình huống ứng dụng thực tế (Use Cases)

  • Đo lường Giá trị chịu rủi ro (Value at Risk - VaR): Sử dụng phương sai điều kiện dự báo $\sigma_{t+1}^2$ từ mô hình GARCH để tính toán ngưỡng tổn thất tối đa danh mục đầu tư ở mức tin cậy $95%$ hoặc $99%$: $$\text{VaR}{\alpha, t+1} = -(\mu{t+1} + z_{\alpha} \sigma_{t+1})$$
  • Thiết lập mức cắt lỗ động (Dynamic Stop-loss): Tự động điều chỉnh biên độ dừng lỗ dựa trên độ lệch chuẩn có điều kiện $\sigma_t$ thay vì sử dụng tỷ lệ phần trăm cố định.
  • Tối ưu hóa phân bổ tỷ trọng danh mục (Asset Allocation): Phân bổ tỷ trọng các mã ngành điện (như PPC, VSH) vào danh mục phòng thủ dựa trên ma trận phương sai - hiệp phương sai có điều kiện (DCC-GARCH).

Kế hoạch triển khai định lượng

flowchart LR
    M1["Giai đoạn 1: Data Ingestion (Tự động cập nhật giá OHLCV & BCTC)"] --> M2["Giai đoạn 2: Model Estimation (Fit ARIMA-GARCH hàng ngày)"]
    M2 --> M3["Giai đoạn 3: Risk Engine (Tính VaR, Stress Testing & Volatility)"]
    M3 --> M4["Giai đoạn 4: Dashboard & Alert (Khuyến nghị tỷ trọng & Cắt lỗ)"]

Hạn chế và hướng phát triển

Hạn chế kỹ thuật

  • Tính đối xứng của mô hình GARCH chuẩn: GARCH(1,1) giả định các cú sốc tích cực (tin tốt) và cú sốc tiêu cực (tin xấu) cùng độ lớn tạo ra mức biến động phương sai như nhau, trong khi thị trường thực tế thường ghi nhận hiệu ứng bất đối xứng đòn bẩy (leverage effect).
  • Mẫu dữ liệu thực nghiệm: Tập trung vào 462 quan sát hàng ngày của 2 doanh nghiệp đầu ngành (VSH, PPC); chưa bao phủ toàn bộ 20 mã doanh nghiệp trong toàn ngành năng lượng niêm yết.

Hướng phát triển mở rộng

  1. Ứng dụng các mô hình phi đối xứng bậc cao như EGARCH (Exponential GARCH) và TGARCH (Threshold GARCH) để đo lường tham số đòn bẩy $\gamma$.
  2. Áp dụng phân phối $t$-Student hoặc Generalized Error Distribution (GED) cho phần dư thay vì giả định phân phối chuẩn, giúp giải quyết triệt để vấn đề hiện tượng đuôi dày (fat tails) trong tài chính.
  3. Mở rộng sang mô hình đa biến DCC-GARCH (Dynamic Conditional Correlation) để phân tích mức độ lan truyền rủi ro (volatility spillover) giữa nhóm cổ phiếu năng lượng và chỉ số VN-Index.

Đối tượng hưởng lợi

  • Nhà đầu tư cá nhân: Nắm bắt phương pháp quản trị rủi ro khoa học, loại bỏ tâm lý đầu tư bầy đàn; hiểu rõ mối quan hệ giữa chu kỳ kinh doanh của nhà máy điện và biến động thị giá cổ phiếu.
  • Chuyên viên phân tích định lượng (Quant / Equity Analyst): Tiếp cận quy trình kiểm định kinh tế lượng bài bản từ khâu lọc nhiễu ARIMA đến chẩn đoán tham số GARCH/IGARCH trên dữ liệu thực tế tại Việt Nam.
  • Doanh nghiệp ngành năng lượng: Có thêm góc nhìn về cách thức thị trường vốn đánh giá các rủi ro hoạt động, tồn kho và nợ vay để tối ưu hóa chiến lược quan hệ nhà đầu tư (IR) và cơ cấu vốn.
  • Cơ quan quản lý và Giới học thuật: Nguồn tài liệu tham khảo thực nghiệm phục vụ việc xây dựng các công cụ phái sinh phòng hộ rủi ro và hoàn thiện khung pháp lý giao dịch thị trường chứng khoán.

Câu hỏi thường gặp

1. Tại sao cần kiểm định tính dừng (ADF test) trước khi chạy mô hình ARCH/GARCH?
Chuỗi giá cổ phiếu thực tế hầu hết là chuỗi không dừng $I(1)$ (bước đi ngẫu nhiên). Việc sử dụng trực tiếp chuỗi không dừng sẽ dẫn đến hiện tượng "hồi quy giả mạo" (spurious regression). Kiểm định ADF giúp xác định chuỗi lợi suất logarit đạt tính dừng $I(0)$, đảm bảo các ước lượng thống kê có ý nghĩa và phương sai có thể hội tụ.

2. Ý nghĩa thực tiễn của hiện tượng IGARCH ($\alpha_1 + \beta_1 = 1$) đối với cổ phiếu VSH là gì?
Khi một cổ phiếu tuân theo mô hình IGARCH, các cú sốc biến động từ thị trường (như thay đổi chính sách giá điện, sự cố thiết bị, biến động thủy văn) sẽ có ảnh hưởng bền bỉ và kéo dài đến phương sai tương lai mà không bị triệt tiêu nhanh chóng. Nhà đầu tư cần dự trữ biên an toàn rủi ro rộng hơn khi giao dịch.

3. Khả năng thanh toán nhanh của VSH đạt 0,44 có phải là dấu hiệu phá sản?
Không. Dù hệ số thanh toán nhanh nhỏ hơn 0,5 thể hiện rủi ro thanh khoản trong ngắn hạn nếu phải trả nợ ngay lập tức, nhưng các doanh nghiệp hạ tầng điện như VSH có dòng tiền kinh doanh (operating cash flow) đều đặn từ hợp đồng mua bán điện (PPA) với EVN. Tỷ lệ ROE đạt 20,07% năm 2012 chứng minh khả năng sinh lời trên vốn chủ vẫn rất vững mạnh.

4. Mô hình GARCH có thể thay thế hoàn toàn việc phân tích báo cáo tài chính không?
Không. GARCH là công cụ định lượng biến động giá và rủi ro thị trường trong ngắn và trung hạn, trong khi phân tích cơ bản (FA) xác định giá trị nội tại và biên an toàn dài hạn của doanh nghiệp. Hai phương pháp phải được phối hợp bổ trợ cho nhau.

5. Cần chuẩn bị công cụ phần mềm gì để tái lập quy trình phân tích này?
Có thể sử dụng các phần mềm kinh tế lượng chuyên dụng như EViews (phiên bản 8.0 trở lên), Stata, R (sử dụng package rugarch / tseries), hoặc ngôn ngữ Python (thư viện arch, statsmodels, pandas).

Kết luận

Chuyên đề nghiên cứu đã xây dựng thành công khung phân tích rủi ro kết hợp giữa chỉ số tài chính cơ bản và mô hình kinh tế lượng chuỗi thời gian ARCH/GARCH/IGARCH đối với các cổ phiếu ngành năng lượng điện Việt Nam (trọng tâm là VSH và PPC). Nghiên cứu không chỉ giải quyết bài toán định lượng hiện tượng phân cụm biến động (volatility clustering) trên dữ liệu thực nghiệm 462 phiên giao dịch mà còn làm rõ tác động lan tỏa dài hạn của các cú sốc rủi ro trên thị trường chứng khoán Việt Nam. Kết quả này cung cấp nền tảng phương pháp luận vững chắc cho các nhà đầu tư và nhà quản trị danh mục trong việc tối ưu hóa hiệu quả sinh lời gắn liền với kiểm soát rủi ro hệ thống.