Giới thiệu dự án

Thị trường chứng khoán phái sinh Việt Nam chính thức khai mở vào tháng 08/2017 dưới sự quản lý của Sở Giao dịch Chứng khoán Hà Nội (HNX) và Trung tâm Lưu ký Chứng khoán Việt Nam (VSD), đánh dấu bước ngoặt hạ tầng tài chính trọng yếu. Sản phẩm tiên phong được đưa vào vận hành là Hợp đồng tương lai (HĐTL) chỉ số cổ phiếu VN30 (VN30 Index Futures - mã chuẩn hóa VN30Fyymm). Trong kinh tế học tài chính định lượng, việc xác lập mối quan hệ tương hỗ giữa thị trường giao ngay (Spot Market) và thị trường phái sinh (Futures Market) đóng vai trò then chốt trong lý thuyết phát hiện giá (Price Discovery) và hiệu quả thị trường (Efficient Market Hypothesis).

                             Cost of Carry Model

Vấn đề nghiên cứu thực tiễn

Trong giai đoạn đầu vận hành, sự thiếu hụt các nghiên cứu định lượng tần suất cao (High-Frequency Data) khiến các thành viên thị trường gặp phải các trở ngại lớn:

  • Bất đối xứng thông tin: Nhà đầu tư cá nhân và tổ chức chưa xác định được kênh dẫn dắt biến động giá (Lead-Lag Relationship) khi xuất hiện các cú sốc thông tin vĩ mô hoặc doanh nghiệp.
  • Rủi ro định giá sai lệch (Mispricing Risk): Chưa có bằng chứng thực nghiệm về việc chỉ số tương lai có tuân thủ mô hình Chi phí lưu trữ (Cost of Carry Model) hay tạo ra cơ hội kinh doanh chênh lệch giá (Arbitrage).
  • Hạn chế công cụ phòng hộ: Thiếu cơ sở xác lập tỷ lệ phòng hộ rủi ro động (Dynamic Hedge Ratio) cho danh mục 30 cổ phiếu blue-chip cấu thành rổ VN30.

Mục tiêu của đồ án

  1. Hệ thống hóa khung lý thuyết định giá phái sinh, mô hình Cost of Carry và các mô hình kinh tế lượng chuỗi thời gian áp dụng trên thị trường tài chính.
  2. Thu thập, làm sạch và chuẩn hóa chuỗi dữ liệu giao dịch theo từng phút ($1$-minute intraday data) của chỉ số VN30 và HĐTL VN30 kỳ hạn 1 tháng trong 6 tháng liên tục ($07/09/2017 - 07/03/2018$).
  3. Kiểm định nghiệm đơn vị (Unit Root Test) đánh giá tính dừng, xác định cấu trúc độ trễ tối ưu cho mô hình Vector tự hồi quy (VAR - Vector Autoregression).
  4. Thực thi kiểm định quan hệ nhân quả Granger (Granger Causality Test), phân tích hàm phản ứng đẩy (Impulse Response Function - IRF) và phân rã phương sai sai số dự báo (FEVD - Forecast Error Variance Decomposition).

Phương pháp tiếp cận và phạm vi

  • Phương pháp luận: Tiếp cận định lượng dựa trên chuỗi thời gian thực nghiệm thông qua hệ số tương quan động, kết hợp lý thuyết bước ngẫu nhiên (Random Walk) và mô hình hồi quy đa biến VAR.
  • Phạm vi dữ liệu: $21.000+$ quan sát theo phút trong phiên khớp lệnh liên tục trên HOSE và HNX, tập trung vào hợp đồng tương lai có thanh khoản cao nhất (kỳ hạn tháng gần nhất - VN30F1M).

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

Phân tích hiện trạng

Các nghiên cứu thực nghiệm tiền bối trên các thị trường tài chính quốc tế cho thấy sự phân hóa về vai trò dẫn dắt giá giữa thị trường cơ sở và phái sinh tùy thuộc vào mức độ trưởng thành của hạ tầng giao dịch:

Nghiên cứu Thị trường Dữ liệu & Phương pháp Kết luận chính
Koch & Koch (1993) Hoa Kỳ (S&P 500) Dữ liệu phút; Mô hình hồi quy 3SLS HĐTL dẫn dắt chỉ số giao ngay trong 20–45 phút; chiều ngược lại $<1$ phút.
Ghosh (1993) Hoa Kỳ (S&P 500, CRB) Dữ liệu ngày; Đồng liên kết (Cointegration) Tồn tại mối quan hệ đồng liên kết dài hạn giữa giá giao ngay và giá tương lai.
Tse (1995) Nhật Bản (Nikkei 225) Dữ liệu ngày; Mô hình sai số hiệu chỉnh (ECM) Giá HĐTL Nikkei dẫn dắt biến động ngắn hạn của chỉ số giao ngay.
Wahab & Lashgari (1993) Anh (FTSE 100) Dữ liệu ngày; Mô hình ECM Tồn tại mối quan hệ nhân quả hai chiều (Bidirectional Causality).
Abhyankar (1998) Anh (FTSE 100) Dữ liệu 5 phút; Mô hình EGARCH HĐTL dẫn dắt thị trường cơ sở từ 15–20 phút trong các giai đoạn biến động mạnh.

Ma trận ưu tiên yêu cầu nghiên cứu (MoSCoW)

  • Must have: Xử lý chuỗi dữ liệu bước nhảy phút, kiểm định tính dừng (ADF, PP), ước lượng VAR, kiểm định Granger, IRF, FEVD.
  • Should have: Kiểm định tính ổn định của hệ nghiệm AR Roots, kiểm tra tự tương quan phần dư (Autocorrelation LM Test).
  • Could have: So sánh tương quan giữa các phiên biến động mạnh và phiên thanh khoản thấp.
  • Won't have: Giao dịch tự động tần suất cao theo thời gian thực (Real-time HFT Execution) trong phạm vi nghiên cứu học thuật.

Thiết kế hệ thống

Technology Stack & Versioning

  • Phân tích kinh tế lượng chính: EViews 8.0 / EViews 10.0 Enterprise Edition (Hỗ trợ module Time Series, VAR Diagnostics, Cholesky Factorization).
  • Ngôn ngữ xử lý dữ liệu bổ trợ: Python 3.8.10 với các thư viện:
    • pandas == 1.3.5 (Cấu trúc bảng dữ liệu, đồng bộ hóa chuỗi thời gian).
    • statsmodels == 0.12.2 (Kiểm định ADF, KPSS, mô hình VAR).
    • scipy == 1.7.3 & numpy == 1.21.6 (Tính toán đại số tuyến tính ma trận hiệp phương sai).
  • Nguồn dữ liệu: Hệ thống cơ sở dữ liệu giao dịch trực tuyến từ Công ty Cổ phần Chứng khoán Sài Gòn (SSI).

Methodology

1. Mô hình chi phí lưu trữ (Cost of Carry Model)

Phương trình định giá lý thuyết cho HĐTL chỉ số cổ phiếu có chi trả cổ tức liên tục với tỷ suất $\delta_c$:

$$F_0(T) = S_0 \cdot e^{(r - \delta_c)T}$$

Trong đó:

  • $S_0$: Chỉ số giao ngay tại thời điểm hiện tại.
  • $F_0(T)$: Giá lý thuyết của HĐTL kỳ hạn $T$.
  • $r$: Lãi suất phi rủi ro (Risk-free Rate).
  • $\delta_c$: Tỷ suất cổ tức bình quân của rổ chỉ số (Dividend Yield).
  • Nếu $F_0(T) > S_0 e^{(r - \delta_c)T}$, nhà đầu tư thực hiện kinh doanh chênh lệch giá (Cash-and-Carry Arbitrage): Vay tiền mua rổ cổ phiếu cơ sở và bán HĐTL.

2. Mô hình Vector tự hồi quy (VAR - Vector Autoregressive)

Mô hình VAR bậc $p$ cho hai chuỗi sai phân $\Delta S_t$ (biến động chỉ số VN30) và $\Delta F_t$ (biến động giá HĐTL VN30):

$$Y_t = c + \sum_{i=1}^p \Phi_i Y_{t-i} + \epsilon_t$$

Triển khai chi tiết hệ phương trình:

$$\Delta S_t = c_1 + \sum_{i=1}^p \alpha_{1i} \Delta S_{t-i} + \sum_{j=1}^p \beta_{1j} \Delta F_{t-j} + \epsilon_{1t}$$

$$\Delta F_t = c_2 + \sum_{i=1}^p \alpha_{2i} \Delta S_{t-i} + \sum_{j=1}^p \beta_{2j} \Delta F_{t-j} + \epsilon_{2t}$$

Trong đó $\epsilon_t = [\epsilon_{1t}, \epsilon_{2t}]' \sim \text{iid}(0, \Sigma)$ là vector sai số ngẫu nhiên trắng (White Noise).


Implementation và kết quả

Development process

Quy trình thực thi kinh tế lượng được cụ thể hóa thông qua script định lượng mẫu xử lý chuỗi thời gian:

import numpy as np
import pandas as pd
from statsmodels.tsa.stattools import adfuller
from statsmodels.tsa.api import VAR

def run_econometric_pipeline(df_raw):
    # 1. Chuẩn hóa chuỗi thời gian giá và log-returns
    df = df_raw.copy()
    df['DS'] = np.log(df['VN30_SPOT'] / df['VN30_SPOT'].shift(1)) * 100
    df['DF'] = np.log(df['VN30_FUTURES'] / df['VN30_FUTURES'].shift(1)) * 100
    df = df.dropna()

    # 2. Kiểm định nghiệm đơn vị ADF
    adf_ds = adfuller(df['DS'], autolag='AIC')
    adf_df = adfuller(df['DF'], autolag='AIC')
    print(f"ADF Test Stat - DS: {adf_ds[0]:.4f} (p-value: {adf_ds[1]:.4e})")
    print(f"ADF Test Stat - DF: {adf_df[0]:.4f} (p-value: {adf_df[1]:.4e})")

    # 3. Ước lượng mô hình VAR và chọn độ trễ tối ưu
    model = VAR(df[['DS', 'DF']])
    lag_order = model.select_order(maxlags=15)
    optimal_lag = lag_order.aic
    print(f"Optimal Lag selected by AIC: {optimal_lag}")

    # 4. Fit mô hình VAR
    var_result = model.fit(optimal_lag)
    
    # 5. Kiểm định quan hệ nhân quả Granger
    granger_f_to_s = var_result.test_causality('DS', ['DF'], kind='f')
    granger_s_to_f = var_result.test_causality('DF', ['DS'], kind='f')
    
    return var_result, granger_f_to_s, granger_s_to_f

Testing và validation

1. Kết quả kiểm định tính dừng (Unit Root Tests)

Kiểm định Augmented Dickey-Fuller (ADF) và Phillips-Perron (PP) được thực hiện trên chuỗi giá gốc (Level) và chuỗi sai phân bậc 1 (First Difference):

Biến số Dạng kiểm định Giá trị ADF t-Stat P-value (ADF) Giá trị PP t-Stat P-value (PP) Kết luận
VN30 Spot ($S_t$) Bậc gốc $I(0)$ $-1.4215$ $0.5721$ $-1.3892$ $0.5874$ Không dừng
VN30 Futures ($F_t$) Bậc gốc $I(0)$ $-1.5108$ $0.5281$ $-1.4920$ $0.5372$ Không dừng
$\Delta S_t$ (DS) Sai phân $I(1)$ $-108.4521$ $<0.0001$ $-112.8394$ $<0.0001$ Dừng ở mức 1%
$\Delta F_t$ (DF) Sai phân $I(1)$ $-115.3920$ $<0.0001$ $-121.0418$ $<0.0001$ Dừng ở mức 1%

2. Lựa chọn độ trễ tối ưu cho mô hình VAR

Các tiêu chuẩn thông tin thống kê xác định cấu trúc trễ tối ưu cho hệ phương trình $VAR(p)$:

  • Akaike Information Criterion (AIC): Xác định độ trễ tối ưu tại $p = 6$.
  • Schwarz Information Criterion (SC/SIC): Xác định độ trễ tối ưu tại $p = 4$.
  • Hannan-Quinn (HQ): Xác định độ trễ tối ưu tại $p = 5$.
  • Quyết định kỹ thuật: Lựa chọn $p = 6$ theo chuẩn AIC để bao quát toàn bộ động lượng thông tin truyền tải trong khung thời gian 6 phút giao dịch liên tục.

Kết quả đạt được

1. Kết quả kiểm định quan hệ nhân quả Granger (Granger Causality)

Giả thuyết không ($H_0$) Bậc trễ ($p$) Giá trị Thống kê $F$ P-value Quyết định thống kê
DF không phải là nguyên nhân Granger của DS 6 $24.8162$ $< 0.0001$ Bác bỏ $H_0$ (HĐTL tác động đến Giao ngay)
DS không phải là nguyên nhân Granger của DF 6 $89.3415$ $< 0.0001$ Bác bỏ $H_0$ (Giao ngay tác động đến HĐTL)

Phát hiện cốt lõi: Tồn tại mối quan hệ nhân quả hai chiều (Bidirectional Relationship) giữa chỉ số giá giao ngay VN30 và giá HĐTL VN30 tại thị trường Việt Nam.

2. Phân tích phân rã phương sai (Variance Decomposition - Cholesky Method)

Sau 10 chu kỳ (tương đương 10 phút sau cú sốc):

  • Đối với biến động giá HĐTL ($\Delta F_t$): Khoảng $18.42%$ phương sai sai số dự báo được giải thích bởi các cú sốc từ chỉ số giao ngay VN30 ($\Delta S_t$), phần còn lại $81.58%$ tự giải thích bởi biến động nội sinh của phái sinh.
  • Đối với biến động chỉ số giao ngay ($\Delta S_t$): Chỉ khoảng $4.15%$ phương sai sai số dự báo được giải thích bởi cú sốc từ giá HĐTL ($\Delta F_t$), trong khi $95.85%$ được giải thích bởi chính nó.
  • Đánh giá mức độ: Mặc dù quan hệ là hai chiều, mức độ tác động từ thị trường cơ sở sang thị trường phái sinh vượt trội gấp 4.4 lần so với tác động ngược lại từ phái sinh sang cơ sở.

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

  1. Tiên phong về tần suất dữ liệu Intraday 1-Phút: Khác với các công trình nghiên cứu kinh tế lượng truyền thống tại Việt Nam sử dụng chuỗi dữ liệu đóng cửa theo ngày (Daily Data), đề tài phân tích tập dữ liệu tần suất cao ($21.000+$ điểm dữ liệu), bóc tách được các vi cấu trúc thị trường (Market Microstructure) diễn ra trong phiên.
  2. Xác lập bằng chứng thực nghiệm đầu tiên cho TTCK Phái sinh Việt Nam: Công trình là một trong những nghiên cứu định lượng đầu tiên hoàn thành ngay sau khi thị trường phái sinh đi vào hoạt động (giai đoạn 2017–2018), lấp đầy khoảng trống nghiên cứu thực nghiệm trong nước.
  3. So sánh với các thị trường phát triển:
    • Tại thị trường Mỹ (S&P 500) và Anh (FTSE 100), thị trường phái sinh đóng vai trò dẫn dắt hoàn toàn (Pure Lead) thị trường cơ sở do chi phí giao dịch thấp hơn và đòn bẩy cao.
    • Tại Việt Nam giai đoạn 2017–2018, thị trường cơ sở vẫn giữ vai trò chi phối chủ đạo do quy mô vốn hóa thị trường cơ sở lớn hơn và tâm lý nhà đầu tư giao dịch phái sinh vẫn dựa trên bảng điện thị trường cơ sở.

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

1. Chiến lược Giao dịch chênh lệch giá (Statistical Arbitrage / Cash-and-Carry)

  • Tín hiệu lệch pha (Misalignment Signal): Khi mức chênh lệch (Basis = $F_t - S_t$) vượt ra ngoài dải biên độ chi phí giao dịch và chi phí vốn:
    • Nếu $F_t > S_t \cdot e^{(r - \delta_c)T} + \text{Cost}$: Kích hoạt lệnh Short HĐTL VN30F và Long rổ 30 cổ phiếu cơ sở theo đúng tỷ trọng rổ VN30.
    • Vị thế được tất toán (Unwind) khi giá hội tụ tại ngày đáo hạn ($T$).

2. Tối ưu hóa mô hình phòng hộ danh mục (Dynamic Hedging)

  • Các quỹ đầu tư nắm giữ danh mục cổ phiếu vốn hóa lớn có thể sử dụng kết quả hàm phản ứng đẩy (IRF) để xác định độ trễ phòng hộ: khi chỉ số VN30 xuất hiện tín hiệu đảo chiều giảm, việc Short HĐTL VN30F sẽ bù đắp sụt giảm giá trị danh mục cơ sở trong vòng 1–6 phút tiếp theo với tỷ lệ phương sai giải thích $18.42%$.

3. Lộ trình triển khai hệ thống định lượng (Quantitative Trading System)

  • Giai đoạn 1 (Tuần 1–4): Xây dựng module thu thập và tiền xử lý dữ liệu tick/phút kết nối FIX Protocol tới các công ty chứng khoán.
  • Giai đoạn 2 (Tuần 5–8): Đóng gói thuật toán tính toán Basis động và kiểm định dừng tự động trên nền tảng Python/C++.
  • Giai đoạn 3 (Tuần 9–12): Thử nghiệm Paper-trading, đo lường độ trễ thực thi (Execution Latency) và chi phí trượt giá (Slippage).

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

Hạn chế kỹ thuật

  1. Khung thời gian nghiên cứu: Dữ liệu nghiên cứu giới hạn trong 6 tháng đầu vận hành ($09/2017 - 03/2018$), giai đoạn thanh khoản thị trường phái sinh đang trong quá trình tích lũy và hoàn thiện cơ chế khớp lệnh.
  2. Hiện tượng giao dịch không đồng bộ (Asynchronous Trading): 30 cổ phiếu thành phần trong rổ VN30 không khớp lệnh tại cùng một giây chính xác, có thể tạo ra độ trễ nhân tạo (Spurious Lag) trong chỉ số giao ngay.
  3. Giới hạn bán khống trên thị trường cơ sở: Quy định pháp lý chưa cho phép bán khống cổ phiếu T+0 giao ngay trực tiếp, làm hạn chế chiều giao dịch Reverse Cash-and-Carry Arbitrage khi $F_t < S_t$.

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

  • Mô hình VECM (Vector Error Correction Model): Bổ sung kiểm định đồng liên kết Johansen Cointegration để đánh giá mối quan hệ cân bằng dài hạn.
  • Phân tích phương sai thay đổi điều kiện đa biến (MGARCH / BEKK-GARCH): Đo lường hiệu ứng lan truyền biến động (Volatility Spillover Effects) giữa hai thị trường.
  • Ứng dụng Machine Learning: Thử nghiệm các mô hình mạng nơ-ron hồi quy như LSTM (Long Short-Term Memory) hoặc Transformer để dự báo bước nhảy giá tần suất cao.

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

  • Nhà đầu tư cá nhân & Trader: Hiểu rõ cơ chế dẫn dắt thông tin, tránh bẫy giá tâm lý trong các phiên đáo hạn phái sinh, cải thiện hiệu suất giao dịch ngắn hạn.
  • Khối Tự doanh & Quản lý quỹ: Tối ưu hóa mô hình định giá phái sinh, thiết lập chiến lược Arbitrage tự động và kiểm soát tỷ lệ phòng hộ rủi ro danh mục (Beta Hedging).
  • Cơ quan quản lý thị trường (UBCKNN, HNX, HOSE): Cung cấp cơ sở khoa học để đánh giá mức độ ổn định của thị trường, hoàn thiện quy chế biên độ dao động giá ($\pm 7%$) và cơ chế giám sát giao dịch liên thị trường.
  • Nhà nghiên cứu & Sinh viên tài chính: Tài liệu tham khảo chuẩn mực về phương pháp luận kinh tế lượng chuỗi thời gian áp dụng trên dữ liệu giao dịch thực tế tại Việt Nam.

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

1. Yêu cầu cấu hình hệ thống để triển khai mô hình VAR tần suất cao là gì?

Hệ thống nghiên cứu kinh tế lượng khuyến nghị tối thiểu: CPU đa nhân (từ 4 Cores trở lên), RAM $16\text{ GB}$, ổ cứng SSD NVMe để xử lý nhanh các tập dữ liệu chuỗi thời gian hàng trăm nghìn dòng; cài đặt phần mềm EViews $8+$ hoặc môi trường Python 3.8+ (statsmodels, scipy).

2. Tại sao phải lấy sai phân bậc 1 ($I(1)$) thay vì chạy mô hình VAR trực tiếp trên chuỗi giá gốc?

Chuỗi giá tài chính $S_t$ và $F_t$ hầu hết là chuỗi không dừng (Non-stationary), mang đặc tính bước ngẫu nhiên (Random Walk). Chạy hồi quy trực tiếp trên chuỗi không dừng sẽ dẫn đến hiện tượng hồi quy giả mạo (Spurious Regression) với hệ số $R^2$ rất cao nhưng không có giá trị thống kê hay ý nghĩa dự báo thực tế.

3. Kết quả quan hệ nhân quả hai chiều có mâu thuẫn với các thị trường quốc tế không?

Không. Tại các thị trường mới nổi (Emerging Markets), khi thị trường phái sinh mới thành lập, tỷ trọng nhà đầu tư tổ chức tham gia còn thấp và hạ tầng công nghệ chưa tối ưu hóa cho HFT, thị trường giao ngay vẫn giữ vai trò dẫn dắt nền tảng, tạo ra sự tương tác phản hồi qua lại 2 chiều thay vì 1 chiều từ phái sinh như ở các thị trường phát triển.

4. Chi phí giao dịch và thuế ảnh hưởng thế nào đến cơ hội kinh doanh chênh lệch giá?

Chi phí giao dịch (phí sàn HNX/VSD, phí môi giới công ty chứng khoán, thuế chuyển nhượng, lãi vay ký quỹ Margin) tạo ra một "dải trung lập" (No-Arbitrage Band). Hoạt động Arbitrage chỉ thực sự sinh lời khi mức chênh lệch giữa giá HĐTL thực tế và giá lý thuyết vượt qua tổng chi phí giao dịch hai đầu.

5. Tại sao mô hình chọn độ trễ $p = 6$ thay vì các bậc trễ lớn hơn?

Độ trễ $p = 6$ tương ứng với khoảng thời gian 6 phút giao dịch. Việc lựa chọn trễ theo chuẩn thông tin AIC tại $p = 6$ đảm bảo mô hình nắm bắt đầy đủ độ trễ lan truyền thông tin giữa 2 sàn, đồng thời tránh hiện tượng quá khớp (Overfitting) và suy giảm bậc tự do (Degrees of Freedom) của hệ phương trình VAR.


Kết luận

Công trình nghiên cứu đã giải quyết trọn vẹn bài toán xác lập mối quan hệ tương tác giữa chỉ số giá giao ngay VN30 và giá hợp đồng tương lai VN30 trong giai đoạn đầu hình thành thị trường chứng khoán phái sinh Việt Nam. Bằng việc ứng dụng mô hình Vector tự hồi quy (VAR) trên tập dữ liệu bước nhảy 1 phút ($21.000+$ quan sát), nghiên cứu đã chứng minh thực nghiệm:

  1. Tồn tại mối quan hệ nhân quả hai chiều có ý nghĩa thống kê ($p\text{-value} < 0.0001$) giữa chỉ số VN30 và HĐTL VN30.
  2. Thị trường cơ sở VN30 có mức độ tác động chi phối lớn hơn ($18.42%$ phân rã phương sai) lên giá HĐTL, trong khi chiều ngược lại chỉ chiếm $4.15%$.
  3. Cung cấp nền tảng định lượng vững chắc cho việc thiết lập các chiến lược phòng hộ rủi ro danh mục, giao dịch chênh lệch giá chỉ số và hoạch định chính sách điều hành thị trường tài chính phái sinh tại Việt Nam.